- Open Access
Tissue dissociation for single-cell and single-nuclei RNA sequencing for low amounts of input material
Frontiers in Zoology volume 19, Article number: 27 (2022)
Recent technological advances opened the opportunity to simultaneously study gene expression for thousands of individual cells on a genome-wide scale. The experimental accessibility of such single-cell RNA sequencing (scRNAseq) approaches allowed gaining insights into the cell type composition of heterogeneous tissue samples of animal model systems and emerging models alike. A major prerequisite for a successful application of the method is the dissociation of complex tissues into individual cells, which often requires large amounts of input material and harsh mechanical, chemical and temperature conditions. However, the availability of tissue material may be limited for small animals, specific organs, certain developmental stages or if samples need to be acquired from collected specimens. Therefore, we evaluated different dissociation protocols to obtain single cells from small tissue samples of Drosophila melanogaster eye-antennal imaginal discs.
We show that a combination of mechanical and chemical dissociation resulted in sufficient high-quality cells. As an alternative, we tested protocols for the isolation of single nuclei, which turned out to be highly efficient for fresh and frozen tissue samples. Eventually, we performed scRNAseq and single-nuclei RNA sequencing (snRNAseq) to show that the best protocols for both methods successfully identified relevant cell types. At the same time, snRNAseq resulted in less artificial gene expression that is caused by rather harsh dissociation conditions needed to obtain single cells for scRNAseq. A direct comparison of scRNAseq and snRNAseq data revealed that both datasets share biologically relevant genes among the most variable genes, and we showed differences in the relative contribution of the two approaches to identified cell types.
We present two dissociation protocols that allow isolating single cells and single nuclei, respectively, from low input material. Both protocols resulted in extraction of high-quality RNA for subsequent scRNAseq or snRNAseq applications. If tissue availability is limited, we recommend the snRNAseq procedure of fresh or frozen tissue samples as it is perfectly suited to obtain thorough insights into cellular diversity of complex tissue.
Gene expression is a central molecular process that coordinates various aspects of organismal life, such as behavior  and development [2, 3]. Since differences in gene expression are often associated with variation in organismal phenotypes, comparative gene expression studies are powerful approaches to establish testable biological hypotheses . For instance, differences in the expression of the developmental transcription factor genes pitx1 and shavenbaby cause natural variation in armor plate formation in stickleback fish  and trichome formation in Drosophila , respectively. Similarly, natural variation in paternal care behavior in Peromyscus mice and density related stress behavior in zebrafish are tightly linked to differences in the expression of genes coding for the hormone vasopressin  and the neuropeptide Parathyroid hormone 2 (Pth2) , respectively. Advances in sequencing technologies have been facilitating extensive insights into the regulation of gene expression on a genome wide scale [9, 10]. A common observation of such studies is that gene expression strongly depends on the biological context. Spatial and temporal gene expression, for example, is tightly regulated throughout development resulting in tissue- and even cell type specific expression profiles [11,12,13,14]. In the light of this context-dependent gene regulation, it is becoming increasingly relevant to study gene expression on a cellular level.
Nowadays, multiple sequencing technologies are available allowing to quantitatively analyze the messenger RNA content of single cells . Single-cell RNA sequencing (scRNAseq) has been proven powerful to reveal the cell type composition of complex tissues or organs in model organisms, such as the fruit fly Drosophila melanogaster [16, 17], the nematode Caenorhabditis elegans  and mouse . Also, biological processes, such as development of the optic lobe of the fly brain , cell–cell communication in tumors  and immunity [22, 23] have been successfully studied. Since the analysis of scRNAseq data does not require prior knowledge of the tissue of interest, this method is exceptionally well-suited to study the cell type composition of emerging model organisms, such as sponges , the cnidaria Nematostella vectensis , Hydra vulgaris  and Clytia hemisphaerica , the annelid Platynereis dumerilii  and the planarian Schmidtea mediterranea [29, 30], the ant Harpegnathos saltator  and multiple vertebrates [32, 33]. Comparative studies have been performed to reveal divergent and conserved aspects of the motor cortex in human, marmoset, and mouse  and during early embryonic development in pigs, humans and cynomolgus monkeys .
scRNAseq protocols are composed of the following key steps [36, 37]: The tissue of interest is dissociated, and individual cells are captured either in microwell plates  or in micro-droplets . Individual captured cells are lysed in the microwell or droplet and the released polyadenylated RNA (mRNA) is captured using poly-T oligos. The mRNA is reverse transcribed into complementary DNA (cDNA) and cell and molecule specific barcodes are added. Subsequently, sequencing libraries are generated by fragmentation, Illumina sequencing adapter ligation and amplification. The amplified libraries are eventually sequenced using next generation sequencing technologies (e.g. Illumina).
While current scRNAseq technologies allow sequencing up to 10,000 cells in one run , many more cells are needed as input material. For instance, mechanical stress during dissociation of complex tissue leads to increased cell death . Additionally, harsh dissociation conditions using enzymes, such as Trypsin, contribute to cell damage [42, 43], altered gene expression [43, 44] and RNA degradation . Due to the high cell loss during dissociation current scRNAseq methods are limited if small tissue samples are analyzed because tissue from multiple animals must be collected to obtain sufficient starting material.
Larval imaginal discs of the fruit fly Drosophila melanogaster are such tiny tissues. These flat epithelial sac-like tissues are specified as about 20 embryonic cells and they grow extensively during larval development to up to 44,000–60,000 cells [46,47,48,49]. During pupae stages, imaginal discs evert and give rise to external adult organs, such as wings, walking legs, genitals and compound eyes . Imaginal discs are excellent models to study fundamental developmental and cellular processes, such as cell proliferation, tissue patterning and morphogenesis [51, 52]. Due to its highly heterogeneous cell type composition, the eye-antennal disc that gives rise to the compound eye, the dorsal ocelli, the antennae, and most of the head capsule [53, 54] (Fig. 1) is especially interesting for scRNAseq applications. Moreover, recent comparative work on the evolution of compound eye size and head morphology in Drosophila species revealed pervasive variation in these adult traits [55,56,57,58,59,60,61,62]. Accordingly, inter- and intraspecific comparisons of eye-antennal disc development have been successful in revealing underlying developmental and molecular mechanisms [59, 60, 63,64,65]. While gene expression in late eye-antennal disc have been studied at single cell resolution [66, 67], earlier stages are less accessible due to low cell numbers. Therefore, we evaluated different dissociation, tissue preservation and sequencing methods to establish an efficient protocol for single-cell transcriptomics in eye-antennal discs.
We show that a combination of mechanical and chemical dissociation works best to obtain sufficient and representative cells for single-cell RNA sequencing (scRNAseq). However, we observed artificial expression of stress related genes, which was most likely due to rather harsh dissociation and cell-sorting conditions. As an alternative, we tested different protocols to isolate single nuclei from fresh and frozen tissue and we show that single-nuclei RNA sequencing (snRNAseq) successfully allowed identifying key cell types without the drawback of stress-response. We find differences in the relative contribution of scRNAseq and snRNAseq to common cell types and we discuss the advantages and disadvantages of both methods. Our work provides an excellent overview of different single cell sequencing approaches when accessibility to tissue samples is limited.
Results and discussion
Tissue dissociation for scRNAseq with low amount of input material
For RNA sequencing of single cells (scRNAseq), heterogenous tissue samples need to be dissociated into live and intact cells. Since about 10,000 cells can be analyzed using the 10 × Genomics Chromium System and about 50% of input cells are lost throughout the preparatory steps, we first tested different dissociation protocols to obtain about 20,000 cells from entire larval organs or about 30 late third instar eye-antennal imaginal discs.
The success of different tissue dissociation protocols was evaluated by estimating the ratio of dead and live cells, as well as the final number of live cells. A dead cell staining with Trypan blue is well-established in homogeneous cell suspensions obtained from cell culture [68, 69]. However, we experienced unreliable dead/live cell ratios with our complex cell suspensions, which was most likely due to Trypan blue positive debris. Therefore, we applied a live-dead assay based on propidium iodide (PI) and Calcein green/violet to identify dead and live cells, respectively. This method allows enrichment of live cells via fluorescence activated cell sorting (FACS), which efficiently also removed debris (Additional file 1: Fig. S1A). Note that the combination of PI and Calcein violet resulted in the most efficient separation of live and dead cells due to a lower spectral overlap of both dyes during FACS. Sorted cells were examined by fluorescent microscopy to confirm that they were mostly Calcein positive and PI negative.
First, we tested purely enzymatic or mechanical dissociation protocols, respectively. Incubation of eye-antennal discs in 10 × TrypLE and 2.5 mg/ml Collagenase even for 2 h did not result in single cell solutions based on visual assessment. Imaginal discs ground with a Dounce homogenizer showed a high proportion of debris and what appeared to be single-nuclei suspensions. Additionally, different attempts resulted in inconsistent dissociation because the low amount of input tissue was barely visible and due to the manual component, it was difficult to balance complete dissociation with the destruction of cells. Based on these observations we reasoned that efficient tissue dissociation required a combination of enzymatic dissociation with gentle mechanical force.
The basic protocol was based on treatment of the tissue with TrypLE and Collagenase on a shaker at 300 rpm with pipet strokes (1000 µl pipet tips) during and after the incubation. We varied the following parameters (see Additional file 1: Table S1): enzyme concentration (1 × and 10 × TrypLE; 2.5 mg/ml and 10 mg/ml Collagenase), incubation time (10–60 min), incubation temperature (37 °C and 30 °C), number of pipet strokes (5 strokes during the incubation and 17–20 strokes after the incubation) and filtration of the cell suspension (no filter, 20 µm and 35 µm filters). 1 × TrypLE was insufficient to achieve complete dissociation in a timely manner and the addition of 10 mg/ml Collagenase resulted in an increased yield, as well as fewer cell aggregates (visual assessment). Incubation for up to 60 min at 30 °C resulted in comparable or slightly more live cells compared to a digestion at 37 °C. Filtration with a filter of 35 µm mesh size did not drastically reduce the proportion of live cells but decreased the amount of debris. The number of pipet strokes after incubation had the highest impact on cell survival with significantly reduced cell survival after more than 17 strokes. We obtained the best results with 16,208 live cells (58% survival rate) from 28 eye-antennal discs after 60 min incubation at 30 °C in 10 × TrypLE and 10 mg/ml Collagenase and 5 pipet strokes during and 17 pipet strokes after the incubation (Additional file 1: Fig. S1B). RNA extracted from this sample was of high quality (Additional file 1: Fig. S2) and suitable for 10X Genomics scRNAseq.
In summary, for low amount of input material, such as < 50 late L3 eye-antennal discs we propose a protocol that combines enzymatic dissociation in conjunction with slight mechanical disruption.
scRNAseq reveals relevant cells and a major impact of heat shock and ribosomal genes
We subjected cells obtained after FACS to a 10X Genomics Chromium run to test if the established dissociation protocol resulted in representative cell types expected in the eye-antennal disc. After droplet-based isolation of RNA from individual cells and subsequent Illumina sequencing, we obtained almost 200 million reads from about 14,500 cells with 13,303 reads and 537 genes per cell (Table 1). 12,000 cells showed less than 10% mitochondrial gene expression (Additional file 1: Fig. S3A) confirming that we mostly isolated live cells. Many reads of the scRNAseq dataset mapped to genes coding for heat shock proteins (Fig. 2A) and among the top ten genes with most variable expression across cells, we found two heat shock related genes (Hsp23 and lncRNA:Hsromega, Fig. 2B). To evaluate whether our rather high dissociation temperature of 30 °C may impose stress on the cells, we re-evaluated a previously published scRNAseq dataset for eye-antennal discs that is based on dissociation at room temperature  (personal communication). As the list of genes expressed in different cell types contains 19 heat shock related genes of which three (Hsp27, Hsp26 and Hsp83) show moderate to high expression in almost all cell types [66 Supplementary Data 2], we suggest that other steps of the dissociation protocol may be responsible for the observed stress response. The distribution of reads also showed a high expression of cytoplasmic genes, such as eEF1alpha1 and eukaryotic elongation factors (Fig. 2A). Additionally, a lot of genes coding for ribosomal proteins were expressed in our dataset (Fig. 2A). The high content of ribosomal genes is expected for scRNAseq because cytoplasmic mRNA is extracted and ribosomal mRNAs are known to be very stable [70, 71]. However, they are often considered uninformative.
We performed an unbiased cluster analysis based on variable gene expression and performed dimension reduction using Uniform Manifold Approximation and Projection (UMAP) . We obtained 22 cell clusters which we annotated based on marker genes used in previous scRNAseq analyses [66, 67] and based on prior knowledge from the literature (see Materials and Methods for details, Additional files 2, 3, and 4: Tables S2 and S3 and Fig. S4). Gene ontology (GO) term enrichment analyses for marker genes defining each of these 22 clusters (Additional file 4: Table S4) suggest that cells in each cluster expressed genes involved in relevant biological processes (Additional file 5: Table S5).
For the sake of comparability (see below), we manually combined similar cell types (see Additional file 3: Table S3) to obtain 11 clusters (Fig. 3A) representing most cell types that have been previously described in scRNAseq data for eye-antennal discs [66, 67] (Fig. 1). To validate our automatic cluster annotation, we tested whether relevant genes were among the top four genes that define a certain cluster. For instance, cut (ct) [73, 74] and homothorax (hth)  have been shown to be involved in antennal development and they are among the top four genes expressed in antennal and eye-antennal border (“EAB”) clusters (Fig. 3B). In line with previously reported roles in photoreceptor differentiation [75, 76], we found amyloid protein precursor-like (Appl) and scratch (scrt) among the top four genes in photoreceptors (“PhR”) (Fig. 3B). And in the morphogenetic furrow (“MFurrow”) cluster, we observed a member of the enhancer of split gene complex (E(spl)m4-BFM) (Fig. 3B), which are broadly expressed in this tissue [77, 78]. Besides these relevant biological findings, the potential stress response of the cells was also evident in our cluster analysis because three heat shock genes (Hsp23, Hsp26 and Hsp68) were among the top four cluster defining genes in the “MFurrow” cluster and those three genes were expressed in most cells of all clusters (Fig. 3B). Note that the “Other” cluster was composed of diverse cell types, such as antennal, dorsal and second mitotic wave (“SMW”) cells (see Fig. 4E, left panel) and accordingly, this cluster was not clearly defined by a set of marker genes (Fig. 3B).
In summary, our tissue dissociation protocol successfully resulted in a cell suspension containing major cell types of the eye-antennal disc but may pose stress on the cells which is detectible through high expression of heat shock genes. Additionally, the high level of ribosomal genes may introduce a bias during further analysis of such data.
Cryo-preservation of imaginal discs for efficient isolation of single nuclei
Our scRNAseq data suggested that the applied dissociation conditions were still stressful for the cells. Additionally, the protocol relies on the processing of fresh tissue samples hampering the analysis of even smaller tissue samples. For instance, eye-antennal discs at the late L3 larval stage contain 44,000–60,000 cells, while discs at the transition from the L2 to the L3 stage are only composed of about 5000 cells . Therefore, about 12 times more discs are needed to obtain sufficiently high cell numbers for scRNAseq applying our single cell dissociation protocol. As tissue growth is an integral part of developmental processes, more efficient protocols are needed to harness the full potential of single cell sequencing methods for developmental biology. To this end, we tested two main approaches: First, we evaluated the use of single nuclei for RNA sequencing (i.e. snRNAseq) as snRNAseq has been shown to result in comparable data, especially for tissue samples that are difficult to dissociate into single cells [80,81,82,83]. Second, we tested the effect of cryo-preservation on the subsequent isolation of single nuclei and RNA integrity as this step allows collecting small tissue samples over time.
For the isolation of single nuclei, we tested two main protocols: One protocol suggested by 10X Genomics is based on NP40 as detergent and a small number of centrifugation and pipetting steps . The other protocol had been established for human heart tissue and is based on using Triton X-100 as a detergent and a variety of RNAse inhibitors to preserve RNA in single nuclei . When 30–50 freshly dissected eye-antennal discs at late L3 stage were used for nuclei isolation, both protocols resulted in more than 20,000 nuclei and extracted RNA was of high quality suitable for snRNAseq (Additional file 1: Fig. S5).
We next dissected imaginal discs, snap-froze them in liquid nitrogen and stored them at − 80 °C for at least one day, or up to four weeks. All applied protocols allowed us to isolate more than 20,000 intact nuclei from about 30 cryo-preserved eye-antennal discs. RNA extracted from nuclei isolated with the 10X Genomics protocol resulted in low RNA quality, suggesting a high level of RNA degradation (Additional file 1: Fig. S6, lanes 7 and 8). The addition of Citric acid to the dissociation buffer has been shown to preserve RNA integrity in human pancreatic cells . However, the use of Citric acid in the 10X Genomics protocol did only marginally improve the quality of RNA extracted from cryo-preseved samples (Additional file 1: Fig. S6, lanes 10 and 11). In contrast, we observed almost no RNA degradation and high RNA quality when we used the protocol that employs RNAse inhibitors (Additional file 1: Fig. S6, lanes 1, 2, 4 and 5). RNA integrity was preserved even when eye-antennal discs were thawed for 3.5 h and frozen again prior to nuclei isolation and RNA extraction (Additional file 1: Fig. S5), showing that the use of RNAse inhibitors are highly efficient to prevent RNA degradation when processing of cryo-preserved tissue samples. Based on the high yield and the high RNA quality, we conclude that the combination of cryo-preservation and nuclei isolation employing RNAse inhibitors is highly efficient to process low input material for snRNAseq.
snRNAseq identifies eye-antennal disc gene expression and reduces technical biases
To test if the nuclei obtained after cryo-preservation are suitable for snRNAseq and represent main cell types of the eye-antennal disc, we subjected the nuclei to a 10X Genomics run to obtain 120 million reads from about 9000 cells with 13,334 reads and 812 genes per cell (Table 1). For the analysis of the data, we applied the same pipeline and settings as for the scRNAseq dataset, with the exception that intronic reads were included because pre-mRNA is expected in nuclei . Among the 10 genes with most variable expression in the dataset (Fig. 2C), we found some with known functions and expression in late L3 eye-antennal discs. For instance, the homophilic cell adhesion molecule Klingon (Klg) is strongly expressed during R7 photoreceptor development  and rotund (rn) that codes for a Kruppel zinc-finger transcription factor is expressed in large parts of the antennal field . Importantly, we did not observe genes associated with heat shock response among the top 10 variable genes (Fig. 2C) and no bias of reads originating from heat shock genes was observed (Fig. 2A), suggesting that they will not impact subsequent clustering analyses as observed for the scRNAseq data. As expected for snRNAseq data [80, 82, 89, 90], we found only two cells with more than 10% mitochondrial gene expression (Additional file 1: Fig. S3B) and 28% fewer reads originating from ribosomal genes (Fig. 2A) compared to scRNAseq, suggesting lower impact of such uninformative reads on the entire dataset.
An unbiased clustering of the snRNAseq data using the same resolution as for scRNAseq (see above) resulted in 24 unique clusters, which we annotated using the same pipeline as described for the scRNAseq data (Additional file 1: Fig. S7). GO term enrichment analyses for marker genes defining each of these 24 cell clusters (Additional file 6: Table S6) revealed biological processes relevant for eye-antennal disc cells (Additional file 7: Table S7).
To improve comparability of the snRNAseq and scRNAseq data (see below), we combined similar clusters (see above & Additional file 3: Table S3), and we obtained the same clusters as for scRNAseq data (Fig. 3B). Additionally, we found clearly assigned clusters for second mitotic wave (“SMW”) and preproneural (“PPN”) cells, which were missing from the scRNAseq data (compare Fig. 3A to C). A closer examination of the top four genes that define a certain cluster revealed for instance Sp1 , disco-related (disco-r) , Distalless (Dll)  and spineless (ss)  in the antennal cell cluster (Fig. 3D) and indeed all four genes have been implicated in antennal development. In line with their role in retinal differentiation and their expression in retinal progenitors we found optix/six3  and twin of eyeless (toy)  among the top four genes in the preproneural (“PPN”) cell cluster. glass (gl) [97,98,99], scratch (scrt)  and anterior open (aop) [100,101,102] have been implicated in late retinal differentiation processes and accordingly, we observed highly specific expression in cells posterior of the morphogenetic furrow (i.e. second mitotic wave, “SMW”; photoreceptors, “PhR” and interommatidial cells, “Interom”). Note that the “Other” cluster was largely defined by ribosomal gene expression (Fig. 3D). As ribosomal mRNAs are more stable in the cytoplasm, we hypothesize that all nuclei that still had remnants of cytoplasm attached show higher levels of ribosomal mRNAs compared to pure nuclei and thus lead to the distinct expression profile and clustering result.
Overall, we conclude that the snRNAseq dataset captured gene expression profiles of eye-antennal disc cells and at the same time, technical artifacts, such as heat shock gene expression and an excess of reads from ribosomal genes were diminished.
scRNAseq and snRNAseq show largely overlapping gene expression profiles
Besides the reduced expression of mitochondrial, ribosomal and heat shock genes in the snRNAseq dataset, we wanted to compare the snRNAseq and scRNAseq datasets more thoroughly. The comparison of the top four genes defining comparable cell clusters in the scRNAseq and snRNAseq datasets revealed five shared genes (arrows in Fig. 3B, D), suggesting that cell-type specific expression is variable in our two datasets. To characterize the overlap of gene expression profiles further, we first asked to what extent the most variable genes overlapped in both datasets. Most of the top ten variable genes in the scRNAseq dataset were also among the most variable genes in the snRNAseq dataset and vice versa (Fig. 2B, C; Additional file 8: Table S8). Among the top 3000 most variable genes, 1480 were present in both datasets, while 1520 genes were unique for the scRNAseq or the snRNAseq dataset, respectively (Fig. 2D; Additional file 9: Table S9). The genes unique for scRNAseq or snRNAseq predominantly represented general cellular and metabolic processes (Additional files 1 and 9: Fig. S8A, B, Table S9), while the shared genes were highly enriched for developmental processes relevant for larval eye-antennal discs (Additional files 1 and 9: Fig. S8C, Table S9). This finding suggests that cell type specific genes may largely be shared between datasets. To confirm this, we compared two types of gene sets.
First, a direct comparison of the manually curated marker genes used for automatic cluster annotation (Additional file 2: Table S2) showed similar expression in comparable clusters between both datasets (Fig. 4A). Exceptions were the antennal marker genes Dll, ss, salr, salm, zfh2, ct and disco-r, which were highly expressed in the antennal (“Antenna”) and eye-antennal boarder (“EAB”) clusters in the snRNAseq and scRNAseq dataset, respectively (Fig. 4A). This observation suggests that these two cell types were not clearly distinguishable in the scRNAseq data.
Second, we asked for each of the comparable clusters (i.e. clusters defined in Fig. 3A and C) how many of the differentially expressed genes (FDR 0.05 and log2-fold change > 0.25) were shared between datasets. Depending on the cell type, we found 6% (e.g. “Antenna” and “Ventral PE”) to 26% (e.g. “PhR” and “Hemocytes”) overlap (Additional files 1 and 10: Fig. S9, Table S10). Among the genes that were shared between datasets, we observed factors with well-established functions in the respective cell types, such as glass (gl) [97,98,99], scratch (scrt)  and amyloid protein precursor-like (Appl)  in photoreceptors, cut (ct) [73, 74] and homothorax (hth)  in antennal and eye-antennal border cells and members of the enhancer of split gene complex [77, 78] in the morphogenetic furrow. Our observation is well in line with a comparison of scRNAseq and snRNAseq data in heterogeneous human liver tissue that also showed cell type specific correlation coefficients in log2-fold changes among the two datasets (coefficients ranging from 0.06 to 0.66) . Therefore, despite fundamental differences in the two applied sequencing approaches, we found that global and cluster-specific gene expression was considerably comparable between our scRNAseq and snRNAseq datasets.
scRNAseq and snRNAseq show differences cell type composition
To directly compare the snRNAseq and the scRNAseq datasets, we combined both datasets by clustering cells and nuclei based on expression profiles and performing dimension reduction. The obtained 21 clusters were annotated based on our manually curated marker gene list (Additional files 1, 2 and 3: Fig. S10, Table S2 and Table S3) and similar clusters were combined as described above for the individual datasets (Fig. 4B; Additional file 3: Table S3). Based on our cluster (i.e. cell type) annotation, we found that the combined dataset contained all major cell types that have been previously described in scRNAseq data for eye-antennal discs [66, 67] (Fig. 4B, Fig. 1).
The UMAP of the integrated data with cells color-coded by the dataset suggested that all major cell types were represented in both datasets and that the different datasets contributed unequally to some cell types (Fig. 4C). Therefore, we quantified the number of cells of each cluster originating from scRNAseq and snRNAseq data, respectively (Fig. 4D). This analysis confirmed the biases in the data composition as 8 of the 13 clusters were predominantly defined by snRNAseq data and 5 clusters were defined by scRNAseq data (Fig. 4D). Clusters with a strong scRNAseq-bias were the ones with ambiguous cell type assignment (“Other”), second mitotic wave (“SMW”) and photoreceptors (“PhR”) (Fig. 4D). The significant contribution of the scRNAseq dataset to second mitotic wave (“SMW”) cells was particularly unexpected as we were unable to unequivocally identify this cell type in the individual scRNAseq dataset (Fig. 3A). To test, which other cell type in the scRNAseq data may contribute to the “SMW” cluster in the integrated dataset, we used the cell identifiers to compare the cluster assignment of each cell in the individual scRNAseq dataset to the integrated dataset (Fig. 4E, left panel). This analysis revealed that cells defined as second mitotic wave (“SMW”) cells in the integrated dataset were part of the large “Other” cluster in the scRNAseq data. Moreover, this analysis showed that the “Other” cluster also contained preproneural (“PPN”) cells (Fig. 4E, left panel), which we could not find in the individual scRNAseq dataset (Fig. 3A). This observation suggests that the scRNAseq data did not contain sufficient information to clearly define these two cell types.
Since we were unable to identify a clear “Dorsal” cluster in our snRNAseq data, we used the specific nuclei identifiers to track individual nuclei between the snRNAseq data and the integrated data (Fig. 4E, right panel). The alluvial plot showed that “Dorsal|Ocelli”, ventral peripodial epithelium (“Ventral PE”), preproneural (“PPN”) and almost all “Other” nuclei contributed to the “Dorsal” cluster in the integrated dataset. This heterogeneous contribution is expected as we combined multiple clusters containing different cell types with dorsal identity for the comparative analysis (Additional file 3: Table S3). For instance, the peripodial epithelium covers the entire eye-antennal disc  and preproneural cells are present along the entire dorsal–ventral axis . Accordingly, both cell types may contribute to dorsal cells.
Generally, the alluvial plots showed that some cell types could be clearly defined in the two different datasets. For instance, photoreceptors (“PhR”) and ventral peripodial epithelium (“Ventral PE”) cells of the scRNAseq dataset, eye-antennal border (“EAB”) and antennal cells of the snRNAseq dataset and interommatidial cells (“Interom”) of both datasets were almost fully recovered as such in the integrated data. For other clusters the cell type assignment in the individual datasets was less clear, which may be caused by the observed differences in cluster-specific gene expression (see above).
In summary, our direct comparison of scRNAseq and snRNAseq data showed that both methods identified major cell types present in the eye-antennal disc, but we observed differences in the relative contribution of the datasets to different cell types. While previous comparisons of scRNAseq and snRNAseq also showed such biases [80, 103, 106, 107], the observed discrepancy between our two datasets could partially be explained by experimental difference. Specifically, larvae for scRNAseq were staged based on morphology and behavior (i.e. late wandering larvae at the transition to white prepupae) and discs for snRNAseq were dissected from larvae staged by developmental time (120 h after egg laying). Therefore, we assume that the larvae used for scRNAseq had been slightly older than those obtained for snRNAseq.
Assessing genome wide gene expression for individual cells has proven powerful to describe the heterogeneity of complex tissues, identify novel cell types and to study biological processes, such as immunity and cell–cell interactions at unprecedented detail. Despite the technological advances, single-cell RNA sequencing (scRNAseq) methods still require many cells as starting material. Therefore, we evaluated different dissociation protocols and compared scRNAseq to single-nuclei RNA sequencing (snRNAseq) with special emphasis on low-input material. Based on data obtained for eye-antennal imaginal discs of Drosophila melanogaster, we found snRNAseq superior to scRNAseq for the following reasons: (1) The isolation of nuclei requires fewer experimental steps compared to tissue dissociation into live cells, increasing reproducibility across experiments. This feature is especially relevant if gene expression comparisons are needed on the level of individual cells, for example to assess the effect of experimental manipulations, to study different developmental stages or to compare species/populations. (2) We observed significantly reduced stress-related expression responses and in line with a previous report in human liver data  we found reduced ribosomal and mitochondrial gene expression in snRNAseq data, suggesting that more informative reads contribute to biological insights. (3) We showed highly efficient nuclei isolation and high-quality RNA extraction from frozen tissue [see also e.g. 108]. It is a major advantage to have the opportunity to collect tissue over time and process samples simultaneously, especially for low-input material. (4) In line with previous reports [80, 82, 83, 89, 90, 107, 109, 110], our snRNAseq dataset contained sufficient expression information to unravel almost all major cell types expected in eye-antennal imaginal discs. We only lacked a clearly defined cluster of dorsal cells in our snRNAseq data, while we were unable to unequivocally identify important second mitotic wave (“SMW”) and preproneural (“PPN”) cells in our scRNAseq. (5) While scRNAseq has been shown to result in biased cell composition, due to different cell sizes, shapes, and survival rate upon dissociation , the more streamlined nuclei isolation procedure should ensure a more representative assessment for snRNAseq, especially for complex organs, such as nervous tissue . For instance, we found indications that snRNAseq may be more efficient in capturing rather complicated cell types, such as the large polyploid cells of the peripodial epithelium (“Ventral PE” in Fig. 4D).
It is important to consider major differences in the analysis and interpretation of scRNAseq and snRNAseq data. For instance, snRNAseq data contains intronic reads originating from immature nuclear RNA . Accordingly, well-annotated genome resources are advantageous and analyses pipelines need to be adjusted to also include reads mapped to introns in subsequent read quantification. snRNAseq data captures only rather transient nuclear RNA, while scRNAseq also includes cytoplasmic mature mRNA. Hence, gene regulation events acting on the level of nuclear export [112, 113], splicing  or mRNA maturation [115, 116] may contribute to differences in expression information derived from nuclei and cells, respectively. If cytoplasmic RNA molecules are of special interest and thus single cells need to be isolated, we strongly suggest a dissociation protocol combining chemical and mechanical treatment of tissue samples in conjunction with FACS-aided life cell selection based on fluorescent life-dead cell staining. Our direct comparison of scRNAseq and snRNAseq data showed that the different datasets contributed differently to the obtained cell types. While some of these differences can be explained by slightly different larval staging procedures in our experiments, it remains to be established, which exact cellular or molecular features influence the more efficient recovery of certain cell types in sc/snRNAseq. In the light of our findings, it will be important to stick to either scRNAseq or snRNAseq if comparative questions are tackled.
In summary, based on a thorough evaluation of different dissociation and sequencing protocols we suggest a highly efficient snRNAseq procedure to obtain high-quality expression data for individual nuclei. Our procedure is specifically tested for low-input material and will therefore be perfectly suited for future studies with limited access to tissue samples.
Fly stock keeping and tissue dissection
Flies of the Oregon R strain of Drosophila melanogaster were kept in fly food vials (400 g of malt extract, 400 g of corn flour, 50 g of soy flour, 110 g of sugar beet syrup, 51 g of agar, 90 g of yeast extract, 31.5 ml of propionic acid, and 7.5 g of Nipagin dissolved in 40 ml of Ethanol, water up to 5 L) in incubators at 25 °C prior to the experiment for at least one generation. The incubators maintain 12 h light and dark cycles. For single-cell RNA sequencing (scRNAseq), eye-antennal discs were dissected from late wandering L3 larvae. For single-nuclei RNA sequencing (snRNAseq), eye-antennal discs were dissected from late L3 larvae at 120 h after egg laying. To control for larval density, eggs were deposited on yeast-coated apple agar plates for 1–2 h and incubated for 24 h before 40 first instar larvae were transferred to food vials for further incubation.
For each dissociation experiment, at least 15 larvae were dissected for no more than 1 h in 400 µl ice-cold 1 × PBS. When larval tissue was used for tests, the larvae were everted, and a mix of inner organs (e.g. imaginal discs, gut, brain etc.) was isolated. When eye-antennal discs were used, ~ 30 eye-antennal discs were dissected (generally from about 15 larvae). All organs were transferred into a microcentrifuge tube containing Storage Buffer (4% BSA, 0.2U Protector RNAse inhibitor (Merck; 3,335,399,001) in PBS). If the sample was to be frozen for later nuclei extraction, the tube was submerged in liquid nitrogen for 2 min and stored at − 80 °C until further processing.
Recommended dissociation protocol to obtain single-cell suspensions for scRNAseq
In the following, the dissociation protocol is described that was used to obtain the live cells used for scRNAseq. Additional file 1: Table S1 contains detailed information about the protocol steps that had been varied and tested to achieve efficient dissociation.
About 30 eye-antennal discs were dissociated in 10 × TrypLE (Thermo Fisher Scientific; A1217701) containing 2.5 mg/ml Collagenase (Invitrogen; 17,100,017) for 30 min on the shaker at 30 °C and 300 rpm. Every 15 min, or, if the digestion time was only 15 min, once after 7.5 min, the discs were pipetted up-and-down 5 times using a 1000 µl pipet tip to dissociate cell clumps efficiently. The reaction was stopped using Schneider’s supplemented medium (SSM, 0.02 mg/ml Insulin in Schneiders Medium (Thermofisher/Gibco; 21,720-024)). The suspension was gently pipetted up and down ~ 17 × with a 1000 µl pipet tip and passed through a 35 µm cell strainer (Corning; 352,235). The suspension was centrifuged for 5 min at 1000 rcf. For low amounts of tissue, the pellet might be small and barely visible on the side wall of the tube. Therefore, it is advantageous to use a swing bucket centrifuge to ensure that the pellet accumulates in the center at the bottom of the tube. The supernatant was removed, and the pellet was resuspended in 1 × PBS. The suspension was centrifuged again (see above), the supernatant was removed, and the pellet was resuspended in 0.04% BSA (Invitrogen; 17,100,017) and 0.2 U/µl Protector RNAse Inhibitor (Sigma-Aldrich; 3,335,402,001) in 1 × PBS.
The cells were stained depending on the application: For testing non-fluorescent live-dead assays, 10 µl of a cell suspension were mixed with 10 µl Trypan Blue (Invitrogen; 15,250,061). 10 µl of this solution were transferred onto a counting chamber and cells were counted using a Zeiss Telaval 31. For fluorescence-based assays, Calcein-AM (green: Sigma-Aldrich; 56,496 or violet: Sigma-Aldrich; ThermoFisher Scientific; C34858) was used to stain live cells at a final concentration of 0.5 µg/ml. The suspension was incubated for 30 min to 1 h in the dark at room temperature on a shaker. Either DAPI or Propidium Iodide were used to stain nuclei at a final concentration of 1 μg/ml each and incubated for 10–30 min. The cell suspension was then immediately processed by Fluorescence Activated Cell Sorting (FACS) at the Universitätsmedizin Göttingen or at the Center for Molecular and Cellular Bioengineering Dresden using a Becton Dickinson BD FACSAria™ II Cell Sorter or BD FACSAria™ III Cell Sorter. In consecutive gating steps, living cells were selected out from debris, damaged cells, and doublets. Events which were positive for Calcein, as well as negative for Propidium Iodide were interpreted as live, undamaged cells. After FACS, cells were visually inspected under the microscope, counted and the volume of the suspension was adjusted with PBS and 0.04% BSA to achieve a concentration of ~ 1,000 cells per µl to match the optimal requirements for 10X Genomics scRNAseq.
Recommended dissociation protocol to obtain single nuclei for snRNAseq
Frozen tissue was thawed at 4 °C and kept on ice for the following steps unless specified otherwise. The tissue was transferred into a precooled Dounce homogenizer (2 ml) and 500 µl of Homogenization Buffer (HB) (0.4 U/µl RiboLock RNase Inhibitor (ThermoFisher Scientific; EO0381), 0.2 U/µl SUPERase In™ RNase Inhibitor (ThermoFisher Scientific; AM2694), 0.10% (v/v) Triton X-100 in NIM2; Nuclei isolation buffer 2 (NIM2): 1 µM DTT, 1 × Protease Inhibitor (Promega; G6521) in NIM1; Nuclei isolation buffer 1 (NIM1): 250 mM Sucrose, 25 mM KCl, 5 mM MgCl2, 10 mM Tris HCl, ph 8 in nuclease free water) was added. The tissue was homogenized with 8 strokes of the tight pestle and kept on ice whenever possible. If the homogenization seemed incomplete after visual inspection, 1 stroke was added at a time up to a maximum of 11 strokes. The homogenized tissue was filtered through a 30 µm MACS SmartStrainer (Miltenyi; 130-098-458) to exclude larger debris. The homogenizer was furthermore washed with 2 × 500 µl of HB to transfer as much of the tissue as possible to the cell strainer. The nuclei suspension was centrifuged at 500 g for 5 min at 4 °C in a swing bucket centrifuge to obtain a nuclei pellet. The supernatant was removed, and the pellet was resuspended in 500 µl Storage Buffer.
For subsequent FACS either 5 µl of a 100 μg/ml DAPI solution (Carl Roth; 6335.1) or 1 drop of NucBlue™ (Hoechst 33,342; Invitrogen Live ReadyProbes™; R37605) was added and the nuclei were incubated for 10–20 min for the staining to occur. During the exposure time of the staining, the sample was immediately transferred to FACS (Becton Dickinson (BD™) FACS Aria III Flow Cytometry Cell Sorter) to collect intact nuclei into a 1.5 ml microcentrifuge tube pre-coated with 1% BSA containing 0.04% BSA in 5 µl PBS. The concentration should be ~ 1,000 nuclei per µl to match the optimal requirements for 10X Genomics snRNAseq. The gates were set to select for DAPI positive nuclei. Particles smaller than 1 µm were excluded to remove small debris and damaged nuclei. Doublets and irregular shaped debris were also filtered out through gating as much as possible. Nozzle size was 100 µm. FACS was performed at the Universitätsmedizin Göttingen or at the Dresden Concept Genome center using a Becton Dickinson BD FACSAria™ II Cell Sorter or BD FACSAria™ III Cell Sorter. The settings were adjusted using unstained and stained samples.
Library preparation and 10 × genomics sequencing
scRNAseq and snRNAseq were performed at the Dresden Concept Genome Center on a 10 × Genomics Chromium sequencing system. The viability of the sorted cells or quality of nuclei were visually inspected under a light microscope (with 200 × magnification) from a small aliquot of cells or nuclei stained with Trypan blue.
Up to 20,000 cells/nuclei were carefully mixed with reverse transcription mix using the Chromium Single Cell 3’ Library & Gel beads chemistry v3 (10X Genomics, PN 1,000,075) and loaded into a Chromium Single Cell B Chip (10X Genomics, PN 1,000,073) on the 10X Genomics Chromium system .
Following the guidelines of the 10X Genomics user manual, the droplets were directly subjected to reverse transcription, the emulsion was broken, and cDNA was purified using Dynabeads MyOne Silane (10X Genomics). After cDNA amplification (11 cycles for cells, 12 cycles for nuclei), the sample was purified and underwent a quality control check on the Fragment Analyzer.
Preparation of single-cell or -nuclei RNA-seq libraries (fragmentation, dA-Tailing, adapter ligation and an indexing PCR step with 12 cycles (cells) or 15 cycles (nuclei)) followed the manufacture’s recommendations. After quantification, the libraries were sequenced on an Illumina NextSeq 500 using a high-output flowcell in PE mode (R1: 28 cycles; I1: 8 cycles; R2: 56 cycles) or on the Illumina Novaseq 6000 system with a S2 flowcell in PE mode (R1: 28 cycles; I1: 8 cycles; R2: 94 cycles). An average of 13,000 fragments per cell were sequenced.
The obtained sequencing data from scRNAseq/snRNAseq were mapped to a genome of the D. melanogaster strain Oregon-R (OreR) (FBsn0000276) and reads mapped to individual genes were counted using 10 × Genomics Cellranger 5 using default settings for mapping single-cell data. For mapping single-nuclei data, the option –include-introns was added. The OreR genome was annotated by transferring the annotation of D. melanogaster genome r6.37 to a previously sequenced genome of Oregon-R [64, 118] using Liftoff .
Further data analyses were performed using R version 4.1.1 (2021-08-10). Specifically, the package Seurat  was used for single-cell specific applications. This includes quality control steps such as calculating the percentage of mitochondrial, ribosomal and heat shock related genes and removing doublets and cells or nuclei of poor quality. Cells of poor quality were defined as expressing more than 3000 or less than 300 genes. Nuclei of poor quality were defined as the top 1% of nuclei expressing the highest number of genes or less than 300 genes. As a threshold of 5–10% of mitochondrial reads per cell resulted in a single cluster defined by mitochondrial genes (in scRNAseq, none in snRNAseq), we decided to exclude cells and nuclei with more than 10% of mitochondrial reads to include as much data as possible for further analyses. Genes were kept if they were expressed in at least 5 cells (for scRNAseq) or 3 nuclei (for snRNAseq). Normalization was performed using the SCTransform method . Unbiased clustering was performed in Seurat  and marker genes enriched in each cell cluster were identified by differential expression analyses (i.e. genes expressed in a cluster vs. all other clusters) followed by a cutoff of log2-fold change > 0.25 and an adjusted p-value < 0.05 (list of cell cluster markers for scRNAseq: Additional file 4: Table S4; list of cell cluster markers for snRANseq: Additional file 6: Table S6). Marker genes for each cell cluster were used to test for enrichment of gene ontology (GO) terms (i.e. Biological Process) using the R package gprofiler2 [123, 124] (GO enrichment for scRNAseq: Additional file 5: Table S5; GO enrichment for snRNAseq: Additional file 7: Table S7). The top four genes defining each cell cluster were chosen by the lowest adjusted p-value.
To compare the most variable genes between scRNAseq and snRNAseq, the top 3000 variable genes for both datasets were obtained based on the differential expression analysis (see above). Both lists were compared to identify those genes that were unique for each dataset and those that were shared. A gene ontology enrichment analysis was performed for the two lists of unique genes and the list of shared genes, respectively, using Gene Ontology. Gene ontology enrichment plots were created using ShinyGO (version 0.76) . Radar plots were generated using Excel.
For the analysis of integrated data from both single cells and single nuclei, we applied the standard workflow in Seurat . We used 3000 integration features and used FindIntegrationAnchors and IntegrateData adapted for datasets normalized using SCtransform. Barplots of dataset-specific celltype proportions were created using the R package dittoSeq . The percentages of cells or nuclei in these plots were corrected for absolute number by multiplying the percentages in cells by 0.74, the ratio of filtered nuclei to filtered cells. UMAP plots were created using Seurat. Each dot represents a single cell or nucleus. They are positioned based on their relative transcriptional similarity to each other. Clusters were identified using a nearest neighbor clustering algorithm and the resolution for the depicted number of clusters was chosen based on visual inspection using the R package clustree . Clusters were annotated by performing a differential gene expression analysis and scoring the overexpressed differentially expressed genes using a matrix of marker genes (Additional file 2: Table S2). To obtain comparable clusters across datasets, similar clusters were manually merged (Additional file 3: Table S3). Variable genes were obtained from the Seurat object and plotted using ggvenn (version 0.1.9.) . Alluvial plots were created using the R package ggalluvial (version 0.12.3) .
The ggplot2 package was used to create plots if not otherwise stated. Note that all scripts and the entire analysis pipeline are available online (https://doi.org/10.25625/YHG4ET).
Availability of data and materials
All Additional files tables and figures are part of this submission. All scripts and the entire analysis pipeline are available online (https://doi.org/10.25625/YHG4ET).
Whitfield CW, Cziko A-M, Robinson GE. Gene expression profiles in the brain predict behavior in individual honey bees. Science. 2003;302:296–9. https://doi.org/10.1126/science.1086807.
Bakken TE, Miller JA, Ding S-L, Sunkin SM, Smith KA, Ng L, et al. A comprehensive transcriptional map of primate brain development. Nature. 2016;535:367–75. https://doi.org/10.1038/nature18637.
Cardoso-Moreira M, Halbert J, Valloton D, Velten B, Chen C, Shao Y, et al. Gene expression across mammalian organ development. Nature. 2019;571:505–9. https://doi.org/10.1038/s41586-019-1338-5.
Buchberger E, Reis M, Lu T-H, Posnien N. Cloudy with a chance of insights: context dependent gene regulation and implications for evolutionary studies. Genes (Basel). 2019. https://doi.org/10.3390/genes10070492.
Shapiro MD, Marks ME, Peichel CL, Blackman BK, Nereng KS, Jónsson B, et al. Genetic and developmental basis of evolutionary pelvic reduction in threespine sticklebacks. Nature. 2004;428:717–23. https://doi.org/10.1038/nature02415.
McGregor AP, Orgogozo V, Delon I, Zanet J, Srinivasan DG, Payre F, Stern DL. Morphological evolution through multiple cis-regulatory mutations at a single gene. Nature. 2007;448:587–90. https://doi.org/10.1038/nature05988.
Bendesky A, Kwon Y-M, Lassance J-M, Lewarch CL, Yao S, Peterson BK, et al. The genetic basis of parental care evolution in monogamous mice. Nature. 2017;544:434–9. https://doi.org/10.1038/nature22074.
Anneser L, Alcantara IC, Gemmer A, Mirkes K, Ryu S, Schuman EM. The neuropeptide Pth2 dynamically senses others via mechanosensation. Nature. 2020. https://doi.org/10.1038/s41586-020-2988-z.
Wang Z, Gerstein M, Snyder M. RNA-seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63. https://doi.org/10.1038/nrg2484.
Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Meth. 2013;10:1213–8. https://doi.org/10.1038/nmeth.2688.
Cusanovich DA, Reddington JP, Garfield DA, Daza RM, Aghamirzaie D, Marco-Ferreres R, et al. The cis-regulatory dynamics of embryonic development at single-cell resolution. Nature. 2018;555:538–42. https://doi.org/10.1038/nature25981.
Kvon EZ, Kazmar T, Stampfel G, Yáñez-Cuna JO, Pagani M, Schernhuber K, et al. Genome-scale functional characterization of Drosophila developmental enhancers in vivo. Nature. 2014;512:91–5. https://doi.org/10.1038/nature13395.
Kim-Hellmuth S, Aguet F, Oliva M, Muñoz-Aguirre M, Kasela S, Wucher V, et al. Cell type-specific genetic regulation of gene expression across human tissues. Science. 2020. https://doi.org/10.1126/science.aaz8528.
The GTEx Consortium. The GTEx consortium atlas of genetic regulatory effects across human tissues. Science. 2020;369:1318–30. https://doi.org/10.1126/science.aaz1776.
Tang F, Barbacioru C, Wang Y, Nordman E, Lee C, Xu N, et al. mRNA-seq whole-transcriptome analysis of a single cell. Nat Meth. 2009;6:377–82. https://doi.org/10.1038/nmeth.1315.
Li H, Janssens J, de Waegeneer M, Kolluru SS, Davie K, Gardeux V, et al. Fly cell atlas: a single-nucleus transcriptomic atlas of the adult fruit fly. Science. 2022;375:eabk2432. https://doi.org/10.1126/science.abk2432.
Davie K, Janssens J, Koldere D, de Waegeneer M, Pech U, Kreft Ł, et al. A single-cell transcriptome atlas of the aging drosophila brain. Cell. 2018;174:982-998.e20. https://doi.org/10.1016/j.cell.2018.05.057.
Cao J, Packer JS, Ramani V, Cusanovich DA, Huynh C, Daza R, et al. Comprehensive single-cell transcriptional profiling of a multicellular organism. Science. 2017;357:661–7. https://doi.org/10.1126/science.aam8940.
The Tabula Muris Consortium. A single-cell transcriptomic atlas characterizes ageing tissues in the mouse. Nature. 2020;583:590–5. https://doi.org/10.1038/s41586-020-2496-1.
Özel MN, Simon F, Jafari S, Holguera I, Chen Y-C, Benhra N, et al. Neuronal diversity and convergence in a visual system developmental atlas. Nature. 2021;589:88–95. https://doi.org/10.1038/s41586-020-2879-3.
Kumar MP, Du J, Lagoudas G, Jiao Y, Sawyer A, Drummond DC, et al. Analysis of single-cell RNA-seq identifies cell-cell communication associated with tumor characteristics. Cell Rep. 2018;25:1458-1468.e4. https://doi.org/10.1016/j.celrep.2018.10.047.
Stubbington MJT, Rozenblatt-Rosen O, Regev A, Teichmann SA. Single-cell transcriptomics to explore the immune system in health and disease. Science. 2017;358:58–63. https://doi.org/10.1126/science.aan6828.
Shalek AK, Satija R, Adiconis X, Gertner RS, Gaublomme JT, Raychowdhury R, et al. Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells. Nature. 2013;498:236–40. https://doi.org/10.1038/nature12172.
Musser JM, Schippers KJ, Nickel M, Mizzon G, Kohn AB, Pape C, et al. Profiling cellular diversity in sponges informs animal cell type and nervous system evolution. Science. 2021;374:717–23. https://doi.org/10.1126/science.abj2949.
Sebé-Pedrós A, Saudemont B, Chomsky E, Plessier F, Mailhé M-P, Renno J, et al. Cnidarian cell type diversity and regulation revealed by whole-organism single-cell RNA-seq. Cell. 2018;173:1520-1534.e20. https://doi.org/10.1016/j.cell.2018.05.019.
Klimovich A, Giacomello S, Björklund Å, Faure L, Kaucka M, Giez C, et al. Prototypical pacemaker neurons interact with the resident microbiota. Proc Natl Acad Sci U S A. 2020;117:17854–63. https://doi.org/10.1073/pnas.1920469117.
Chari T, Weissbourd B, Gehring J, Ferraioli A, Leclère L, Herl M, et al. Whole-animal multiplexed single-cell RNA-seq reveals transcriptional shifts across Clytia medusa cell types. Sci Adv. 2021;7:eabh1683. https://doi.org/10.1126/sciadv.abh1683.
Vergara HM, Bertucci PY, Hantz P, Tosches MA, Achim K, Vopalensky P, Arendt D. Whole-organism cellular gene-expression atlas reveals conserved cell types in the ventral nerve cord of Platynereis dumerilii. Proc Natl Acad Sci U S A. 2017;114:5878–85. https://doi.org/10.1073/pnas.1610602114.
Plass M, Solana J, Wolf FA, Ayoub S, Misios A, Glažar P, et al. Cell type atlas and lineage tree of a whole complex animal by single-cell transcriptomics. Science. 2018. https://doi.org/10.1126/science.aaq1723.
Fincher CT, Wurtzel O, de Hoog T, Kravarik KM, Reddien PW. Cell type transcriptome atlas for the planarian Schmidtea mediterranea. Science. 2018. https://doi.org/10.1126/science.aaq1736.
Sheng L, Shields EJ, Gospocic J, Glastad KM, Ratchasanmuang P, Berger SL, et al. Social reprogramming in ants induces longevity-associated glia remodeling. Sci Adv. 2020;6:eaba9869. https://doi.org/10.1126/sciadv.aba9869.
Tosches MA, Yamawaki TM, Naumann RK, Jacobi AA, Tushev G, Laurent G. Evolution of pallium, hippocampus, and cortical cell types revealed by single-cell transcriptomics in reptiles. Science. 2018;360:881–8. https://doi.org/10.1126/science.aar4237.
Chen D, Sun J, Zhu J, Ding X, Lan T, Wang X, et al. Single cell atlas for 11 non-model mammals, reptiles and birds. Nat Commun. 2021;12:7083. https://doi.org/10.1038/s41467-021-27162-2.
Bakken TE, Jorstad NL, Hu Q, Lake BB, Tian W, Kalmbach BE, et al. Comparative cellular analysis of motor cortex in human, marmoset and mouse. Nature. 2021;598:111–9. https://doi.org/10.1038/s41586-021-03465-8.
Liu T, Li J, Yu L, Sun H-X, Li J, Dong G, et al. Cross-species single-cell transcriptomic analysis reveals pre-gastrulation developmental differences among pigs, monkeys, and humans. Cell Discov. 2021;7:8. https://doi.org/10.1038/s41421-020-00238-x.
Kolodziejczyk AA, Kim JK, Svensson V, Marioni JC, Teichmann SA. The technology and biology of single-cell RNA sequencing. Mol Cell. 2015;58:610–20. https://doi.org/10.1016/j.molcel.2015.04.005.
Wu AR, Wang J, Streets AM, Huang Y. Single-cell transcriptional analysis. Annu Rev Anal Chem (Palo Alto Calif). 2017;10:439–62. https://doi.org/10.1146/annurev-anchem-061516-045228.
Islam S, Kjällquist U, Moliner A, Zajac P, Fan J-B, Lönnerberg P, Linnarsson S. Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq. Genome Res. 2011;21:1160–7. https://doi.org/10.1101/gr.110882.110.
Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell. 2015;161:1202–14. https://doi.org/10.1016/j.cell.2015.05.002.
Svensson V, Vento-Tormo R, Teichmann SA. Exponential scaling of single-cell RNA-seq in the past decade. Nat Protoc. 2018;13:599–604. https://doi.org/10.1038/nprot.2017.149.
Khan SJ, Abidi SNF, Tian Y, Skinner A, Smith-Bolton RK. A rapid, gentle and scalable method for dissociation and fluorescent sorting of imaginal disc cells for mRNA sequencing. Fly (Austin). 2016;10:73–80. https://doi.org/10.1080/19336934.2016.1173296.
Hodges GM, Livingston DC, Franks LM. The localization of trypsin in cultured mammalian cells. J Cell Sci. 1973;12:887–902. https://doi.org/10.1242/jcs.12.3.887.
Snow C, Allen A. The release of radioactive nucleic acids and mucoproteins by trypsin and ethylenediaminetetra-acetate treatment of baby-hamster cells in tissue culture. Biochem J. 1970;119:707–14. https://doi.org/10.1042/bj1190707.
Huang H-L, Hsing H-W, Lai T-C, Chen Y-W, Lee T-R, Chan H-T, et al. Trypsin-induced proteome alteration during cell subculture in mammalian cells. J Biomed Sci. 2010;17:36. https://doi.org/10.1186/1423-0127-17-36.
Vrtačnik P, Kos Š, Bustin SA, Marc J, Ostanek B. Influence of trypsinization and alternative procedures for cell preparation before RNA extraction on RNA integrity. Anal Biochem. 2014;463:38–44. https://doi.org/10.1016/j.ab.2014.06.017.
Gehring W, Noethiger R. The imaginal discs of drosophila. In: Waddington CH, Counce-Nicklas S, editors. Developmental systems: insects. London-New York: Academic Press; 1973.
Garcia-Bellido A, Merriam JR. Cell lineage of the imaginal discs in drosophila gynandromorphs. J Exp Zool. 1969;170:61–75. https://doi.org/10.1002/jez.1401700106.
Becker HJ. Ueber Roentgenmosaikflecken Und Defektmutationen Am Auge Von Drosophila Und Die Entwicklungsphysiologie Des Auges. Z Indukt Abstamm Vererbungsl. 1957;88:333–73. https://doi.org/10.1007/bf00309128.
Weasner BP, Kumar JP. The early history of the eye-antennal disc of drosophila melanogaster. Genetics. 2022. https://doi.org/10.1093/genetics/iyac041.
Cohen SM. Imaginal disc development. In: Bate M, Martinez Arias A, editors. The development of drosophila melanogaster. Cold Spring Harbor: Cold Spring Harbor Laboratory Press; 1993.
Domínguez M, Casares F. Organ specification-growth control connection: new in-sights from the drosophila eye-antennal disc. Dev Dyn. 2005;232:673–84. https://doi.org/10.1002/dvdy.20311.
Kumar JP. The fly eye: through the looking glass. Dev Dyn. 2018;247:111–23. https://doi.org/10.1002/dvdy.24585.
Morata G, Lawrence PA. Development of the eye-antenna imaginal disc of drosophila. Dev Biol. 1979;70:355–71. https://doi.org/10.1016/0012-1606(79)90033-2.
Haynie JL, Bryant PJ. Development of the eye-antenna imaginal disc and morphogenesis of the adult head in drosophila melanogaster. J Exp Zool. 1986;237:293–308. https://doi.org/10.1002/jez.1402370302.
Keesey IW, Grabe V, Gruber L, Koerte S, Obiero GF, Bolton G, et al. Inverse resource allocation between vision and olfaction across the genus drosophila. Nat Commun. 2019;10:1162. https://doi.org/10.1038/s41467-019-09087-z.
Posnien N, Hopfen C, Hilbrant M, Ramos-Womack M, Murat S, Schönauer A, et al. Evolution of eye morphology and rhodopsin expression in the drosophila melanogaster species subgroup. PLoS ONE. 2012;7: e37346. https://doi.org/10.1371/journal.pone.0037346.
Hämmerle B, Ferrús A. Expression of enhancers is altered in drosophila melanogaster hybrids. Evol Dev. 2003;5:221–30. https://doi.org/10.1046/j.1525-142x.2003.03030.x.
Norry FM, Gomez FH. Quantitative trait loci and antagonistic associations for two developmentally related traits in the drosophila head. J Insect Sci. 2017. https://doi.org/10.1093/jisesa/iew115.
Gaspar P, Arif S, Sumner-Rooney L, Kittelmann M, Bodey AJ, Stern DL, et al. Characterization of the genetic architecture underlying eye size variation within drosophila melanogaster and drosophila simulans. G3 (Bethesda). 2020;10:1005–18. https://doi.org/10.1534/g3.119.400877.
Ramaekers A, Claeys A, Kapun M, Mouchel-Vielh E, Potier D, Weinberger S, et al. Altering the temporal regulation of one transcription factor drives evolutionary trade-offs between head sensory organs. Dev Cell. 2019;50:780-792.e7. https://doi.org/10.1016/j.devcel.2019.07.027.
Reis M, Wiegleb G, Claude J, Lata R, Horchler B, Ha N-T, et al. Multiple loci linked to inversions are associated with eye size variation in species of the Drosophila virilis phylad. Sci Rep. 2020;10:12832. https://doi.org/10.1038/s41598-020-69719-z.
Norry FM, Vilardi JC, Hasson E. Negative genetic correlation between traits of the Drosophila head, and interspecific divergence in head shape. Heredity (Edinb). 2000;85(Pt 2):177–83. https://doi.org/10.1046/j.1365-2540.2000.00735.x.
Casares F, McGregor AP. The evolution and development of eye size in flies. Wiley Interdiscip Rev Dev Biol. 2021;10: e380. https://doi.org/10.1002/wdev.380.
Buchberger E, Bilen A, Ayaz S, Salamanca D, Matas de Las Heras C, Niksic A, et al. Variation in pleiotropic hub gene expression is associated with interspecific differences in head shape and eye size in drosophila. Mol Biol Evol. 2021;38:1924–42. https://doi.org/10.1093/molbev/msaa335.
Arif S, Hilbrant M, Hopfen C, Almudi I, Nunes MDS, Posnien N, et al. Genetic and developmental analysis of differences in eye and face morphology between drosophila simulans and drosophila mauritiana. Evol Dev. 2013;15:257–67. https://doi.org/10.1111/ede.12027.
Ariss MM, Islam ABMMK, Critcher M, Zappia MP, Frolov MV. Single cell RNA-sequencing identifies a metabolic aspect of apoptosis in Rbf mutant. Nat Commun. 2018;9:5024. https://doi.org/10.1038/s41467-018-07540-z.
Bravo González-Blas C, Quan X-J, Duran-Romaña R, Taskiran II, Koldere D, Davie K, et al. Identification of genomic enhancers through spatial integration of single-cell transcriptomics and epigenomics. Mol Syst Biol. 2020;16: e9438. https://doi.org/10.15252/msb.20209438.
Pappenheimer AM. Experimental studies upon lymphocytes: I. The reactions of lymphocytes under various experimental conditions. J Exp Med. 1917;25:633–50. https://doi.org/10.1084/jem.25.5.633.
Chan LL-Y, Kuksin D, Laverty DJ, Saldi S, Qiu J. Morphological observation and analysis using automated image cytometry for the comparison of trypan blue and fluorescence-based viability detection method. Cytotechnology. 2015;67:461–73. https://doi.org/10.1007/s10616-014-9704-5.
Rabani M, Levin JZ, Fan L, Adiconis X, Raychowdhury R, Garber M, et al. Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells. Nat Biotechnol. 2011;29:436–42. https://doi.org/10.1038/nbt.1861.
Rabani M, Raychowdhury R, Jovanovic M, Rooney M, Stumpo DJ, Pauli A, et al. High-resolution sequencing and modeling identifies distinct dynamic RNA regulatory strategies. Cell. 2014;159:1698–710. https://doi.org/10.1016/j.cell.2014.11.015.
McInnes L, Healy J, Melville J. UMAP: uniform manifold approximation and projection for dimension reduction. arXiv 2020. https://doi.org/10.48550/arXiv.1802.03426.
Dong PDS, Dicks JS, Panganiban G. Distal-less and homothorax regulate multiple targets to pattern the drosophila antenna. Development. 2002;129:1967–74. https://doi.org/10.1242/dev.129.8.1967.
Blochlinger K, Jan LY, Jan YN. Postembryonic patterns of expression of cut, a locus regulating sensory organ identity in drosophila. Development. 1993;117:441–50. https://doi.org/10.1242/dev.117.2.441.
Luo LQ, Le M-M, White K. Identification, secretion, and neural expression of APPL a drosophila protein similar to human amyloid protein precursor. J Neurosci. 1990;10:3849–61. https://doi.org/10.1523/JNEUROSCI.10-12-03849.1990.
Roark M, Sturtevant MA, Emery J, Vaessin H, Grell E, Bier E. scratch, a pan-neural gene encoding a zinc finger protein related to snail, promotes neuronal development. Genes Dev. 1995;9:2384–98. https://doi.org/10.1101/gad.9.19.2384.
Wurmbach E, Wech I, Preiss A. The enhancer of split complex of drosophila melanogaster harbors three classes of notch responsive genes. Mech Dev. 1999;80:171–80. https://doi.org/10.1016/S0925-4773(98)00212-3.
Wolff T, Ready DF. Pattern formation in the drosophila retina. In: Bate M, Martinez-Arias A, editors. The development of drosophila melanogaster. II. Cold Spring Harbor: Cold Spring Harbor Laboratory Press; 1993. p. 1277–325.
Martin PF. Direct determination of the growth rate of drosophila imaginal discs. J Exp Zool. 1982;222:97–102. https://doi.org/10.1002/jez.1402220113.
Wu H, Kirita Y, Donnelly EL, Humphreys BD. Advantages of single-nucleus over single-cell RNA sequencing of adult kidney: rare cell types and novel cell states revealed in fibrosis. J Am Soc Nephrol. 2019;30:23–32. https://doi.org/10.1681/ASN.2018090912.
Yim AKY, Wang PL, Bermingham JR, Hackett A, Strickland A, Miller TM, et al. Disentangling glial diversity in peripheral nerves at single-nuclei resolution. Nat Neurosci. 2022;25:238–51. https://doi.org/10.1038/s41593-021-01005-1.
Lake BB, Codeluppi S, Yung YC, Gao D, Chun J, Kharchenko PV, et al. A comparative strategy for single-nucleus and single-cell transcriptomes confirms accuracy in predicted cell-type expression from nuclear RNA. Sci Rep. 2017;7:6031. https://doi.org/10.1038/s41598-017-04426-w.
Grindberg RV, Yee-Greenbaum JL, McConnell MJ, Novotny M, O’Shaughnessy AL, Lambert GM, et al. RNA-sequencing from single nuclei. Proc Natl Acad Sci U S A. 2013;110:19802–7. https://doi.org/10.1073/pnas.1319700110.
10X Genomics. Nuclei isolation from cell suspensions & tissues for single cell RNA sequencing: CG000124 Rev F. 2021. https://assets.ctfassets.net/an68im79xiti/2HNFgXau0ntv1BhS4ffn6n/71a29daf18e5f7c30cf06b5b4f829e44/CG000124_Demonstrated_Protocol_Nuclei_isolation_RevF.pdf. Accessed 28 Apr 2022.
Litvinukova M, Lindberg E, Maatz H, Zhang H, Radke M, Gotthardt M, et al. Single cell and single nuclei analysis human heart tissue. Protocols. 2018. https://doi.org/10.17504/protocols.io.veae3ae.
Tosti L, Hang Y, Debnath O, Tiesmeyer S, Trefzer T, Steiger K, et al. Single-nucleus and in situ RNA-sequencing reveal cell topographies in the human pancreas. Gastroenterology. 2021;160:1330-1344.e11. https://doi.org/10.1053/j.gastro.2020.11.010.
Butler SJ, Ray S, Hiromi Y. Klingon, a novel member of the drosophila immunoglobulin superfamily, is required for the development of the R7 photoreceptor neuron. Development. 1997;124:781–92. https://doi.org/10.1242/dev.124.4.781.
St Pierre SE, Galindo MI, Couso JP, Thor S. Control of drosophila imaginal disc development by rotund and roughened eye: differentially expressed transcripts of the same gene encoding functionally distinct zinc finger proteins. Development. 2002;129:1273–81. https://doi.org/10.1242/dev.129.5.1273.
Lake BB, Chen S, Sos BC, Fan J, Kaeser GE, Yung YC, et al. Integrative single-cell analysis of transcriptional and epigenetic states in the human adult brain. Nat Biotechnol. 2018;36:70–80. https://doi.org/10.1038/nbt.4038.
Bakken TE, Hodge RD, Miller JA, Yao Z, Nguyen TN, Aevermann B, et al. Single-nucleus and single-cell transcriptomes compared in matched cortical cell types. PLoS ONE. 2018;13: e0209648. https://doi.org/10.1371/journal.pone.0209648.
Estella C, Mann RS. Non-redundant selector and growth-promoting functions of two sister genes, buttonhead and Sp1, in drosophila leg development. PLoS Genet. 2010;6: e1001001. https://doi.org/10.1371/journal.pgen.1001001.
Grubbs N, Leach M, Su X, Petrisko T, Rosario JB, Mahaffey JW. New components of drosophila leg development identified through genome wide association studies. PLoS ONE. 2013;8: e60261. https://doi.org/10.1371/journal.pone.0060261.
Gorfinkiel N, Morata G, Guerrero I. The homeobox gene distal-less induces ventral appendage development in drosophila. Genes Dev. 1997;11:2259–71. https://doi.org/10.1101/gad.11.17.2259.
Duncan DM, Burgess EA, Duncan I. Control of distal antennal identity and tarsal development in drosophila by spineless-aristapedia, a homolog of the mammalian dioxin receptor. Genes Dev. 1998;12:1290–303. https://doi.org/10.1101/gad.12.9.1290.
Seimiya M, Gehring WJ. The drosophila homeobox gene optix is capable of inducing ectopic eyes by an eyeless-independent mechanism. Development. 2000;127:1879–86. https://doi.org/10.1242/dev.127.9.1879.
Czerny T, Halder G, Kloter U, Souabni A, Gehring WJ, Busslinger M. Twin of eyeless, a second pax-6 gene of drosophila, acts upstream of eyeless in the control of eye development. Mol Cell. 1999;3:297–307. https://doi.org/10.1016/s1097-2765(00)80457-8.
Bernardo-Garcia FJ, Fritsch C, Sprecher SG. The transcription factor glass links eye field specification with photoreceptor differentiation in drosophila. Development. 2016;143:1413–23. https://doi.org/10.1242/dev.128801.
Morrison CA, Chen H, Cook T, Brown S, Treisman JE. Glass promotes the differentiation of neuronal and non-neuronal cell types in the drosophila eye. PLoS Genet. 2018;14: e1007173. https://doi.org/10.1371/journal.pgen.1007173.
Moses K, Ellis MC, Rubin GM. The glass gene encodes a zinc-finger protein required by drosophila photoreceptor cells. Nature. 1989;340:531–6. https://doi.org/10.1038/340531a0.
Lai Z. Negative control of photoreceptor development in drosophila by the product of the yan gene, an ETS domain protein. Cell. 1992;70:609–20. https://doi.org/10.1016/0092-8674(92)90430-K.
Weber U, Pataki C, Mihaly J, Mlodzik M. Combinatorial signaling by the frizzled/PCP and Egfr pathways during planar cell polarity establishment in the drosophila eye. Dev Biol. 2008;316:110–23. https://doi.org/10.1016/j.ydbio.2008.01.016.
Rebay I, Rubin GM. Yan functions as a general inhibitor of differentiation and is negatively regulated by activation of the Ras1/MAPK pathway. Cell. 1995;81:857–66. https://doi.org/10.1016/0092-8674(95)90006-3.
Andrews TS, Atif J, Liu JC, Perciani CT, Ma X-Z, Thoeni C, et al. Single-cell, single-nucleus, and spatial RNA sequencing of the human liver identifies cholangiocyte and mesenchymal heterogeneity. Hepatol Commun. 2022;6:821–40. https://doi.org/10.1002/hep4.1854.
Atkins M, Mardon G. Signaling in the third dimension: the peripodial epithelium in eye disc development. Dev Dyn. 2009;238:2139–48. https://doi.org/10.1002/dvdy.22034.
Treisman JE. Retinal differentiation in drosophila. Wiley Interdiscip Rev Dev Biol. 2013;2:545–57. https://doi.org/10.1002/wdev.100.
Denisenko E, Guo BB, Jones M, Hou R, de Kock L, Lassmann T, et al. Systematic assessment of tissue dissociation and storage biases in single-cell and single-nucleus RNA-seq workflows. Genome Biol. 2020;21:130. https://doi.org/10.1186/s13059-020-02048-6.
Ding J, Adiconis X, Simmons SK, Kowalczyk MS, Hession CC, Marjanovic ND, et al. Systematic comparison of single-cell and single-nucleus RNA-sequencing methods. Nat Biotechnol. 2020;38:737–46. https://doi.org/10.1038/s41587-020-0465-8.
Krishnaswami SR, Grindberg RV, Novotny M, Venepally P, Lacar B, Bhutani K, et al. Using single nuclei for RNA-seq to capture the transcriptome of postmortem neurons. Nat Protoc. 2016;11:499–524. https://doi.org/10.1038/nprot.2016.015.
Habib N, Li Y, Heidenreich M, Swiech L, Avraham-Davidi I, Trombetta JJ, et al. Div-seq: single-nucleus RNA-Seq reveals dynamics of rare adult newborn neurons. Science. 2016;353:925–8. https://doi.org/10.1126/science.aad7038.
Slyper M, Porter CBM, Ashenberg O, Waldman J, Drokhlyansky E, Wakiro I, et al. A single-cell and single-nucleus RNA-seq toolbox for fresh and frozen human tumors. Nat Med. 2020;26:792–802. https://doi.org/10.1038/s41591-020-0844-1.
Darmanis S, Sloan SA, Zhang Y, Enge M, Caneda C, Shuer LM, et al. A survey of human brain transcriptome diversity at the single cell level. Proc Natl Acad Sci U S A. 2015;112:7285–90. https://doi.org/10.1073/pnas.1507125112.
Palazzo AF, Lee ES. Sequence determinants for nuclear retention and cytoplasmic export of mRNAs and lncRNAs. Front Genet. 2018;9:440. https://doi.org/10.3389/fgene.2018.00440.
Wickramasinghe VO, Laskey RA. Control of mammalian gene expression by selective mRNA export. Nat Rev Mol Cell Biol. 2015;16:431–42. https://doi.org/10.1038/nrm4010.
Gehring NH, Roignant J-Y. Anything but ordinary—emerging splicing mechanisms in eukaryotic gene regulation. Trends Genet. 2021;37:355–72. https://doi.org/10.1016/j.tig.2020.10.008.
Galloway A, Cowling VH. mRNA cap regulation in mammalian cell function and fate. Biochim Biophys Acta Gene Regul Mech. 2019;1862:270–9. https://doi.org/10.1016/j.bbagrm.2018.09.011.
Mittleman BE, Pott S, Warland S, Zeng T, Mu Z, Kaur M, et al. Alternative polyadenylation mediates genetic regulation of gene expression. Elife. 2020. https://doi.org/10.7554/eLife.57492.
Zheng GXY, Terry JM, Belgrader P, Ryvkin P, Bent ZW, Wilson R, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017;8:14049. https://doi.org/10.1038/ncomms14049.
Torres-Oliva M, Almudi I, McGregor AP, Posnien N. A robust (re-)annotation approach to generate unbiased mapping references for RNA-seq-based analyses of differential expression across closely related species. BMC Genomics. 2016;17:392. https://doi.org/10.1186/s12864-016-2646-x.
Shumate A, Salzberg SL. Liftoff: accurate mapping of gene annotations. Bioinformatics. 2020. https://doi.org/10.1093/bioinformatics/btaa1016.
Hao Y, Hao S, Andersen-Nissen E, Mauck WM, Zheng S, Butler A, et al. Integrated analysis of multimodal single-cell data. Cell. 2021;184:3573-3587.e29. https://doi.org/10.1016/j.cell.2021.04.048.
Hafemeister C, Satija R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biol. 2019;20:296. https://doi.org/10.1186/s13059-019-1874-1.
Waltman L, van Eck NJ. A smart local moving algorithm for large-scale modularity-based community detection. Eur Phys J B. 2013. https://doi.org/10.1140/epjb/e2013-40829-0.
Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H, Vilo J. g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 2019;47:W191–8. https://doi.org/10.1093/nar/gkz369.
Kolberg L, Raudvere U, Kuzmin I, Vilo J, Peterson H. gprofiler2 an R package for gene list functional enrichment analysis and namespace conversion toolset g:Profiler. Fresidency. 2020;9:109. https://doi.org/10.12688/f1000research.24956.2.
Ge SX, Jung D, Yao R. ShinyGO: a graphical gene-set enrichment tool for animals and plants. Bioinformatics. 2020;36:2628–9. https://doi.org/10.1093/bioinformatics/btz931.
Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, et al. Comprehensive integration of single-cell data. Cell. 2019;177:1888-1902.e21. https://doi.org/10.1016/j.cell.2019.05.031.
Bunis DG, Andrews J, Fragiadakis GK, Burt TD, Sirota M. dittoSeq: universal user-friendly single-cell and bulk RNA sequencing visualization toolkit. Bioinformatics. 2020. https://doi.org/10.1093/bioinformatics/btaa1011.
Zappia L, Oshlack A. Clustering trees: a visualization for evaluating clusterings at multiple resolutions. Gigascience. 2018. https://doi.org/10.1093/gigascience/giy083.
Linlin Y. ggvenn: draw venn diagram by ‘ggplot2’: Linlin Yan; 2021.
Brunson JC. ggalluvial: alluvial plots in ‘ggplot2’: Brunson JC; Read QD; 2020.
Many thanks to the CMCB Flow Cytometry Core Facility in Dresden and the Core Facility Cell-Sorting at the UMG in Göttingen for their support with FACS sorting. Thanks to the group of Holger Bastians for providing us with equipment for flow cytometry. We also want to thank Xingbo Xu at the UMG, as well as Ana Veloso at the MDC Berlin for advice on single-nuclei RNA sequencing. Furthermore, we appreciate helpful general input by Anna Schönauer, Jordi Solana and Argyris Papantonis and the critical assessment of the manuscript by Juliana Roscito. We thank the two anonymous reviewers for valuable suggestions which improved the manuscript considerably. We acknowledge support by the Open Access Publication Funds of the Göttingen University.
Open Access funding enabled and organized by Projekt DEAL. This work is supported by a Research Grant of the German Research Foundation (PO 1648/4-1) to NP.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
FACS-plot of D. melanogaster eye-antennal disc cells after live/dead cell staining. (A) Counterstaining of Propidium Iodide to label dead cells (y-axis, Q1) and Calcein violet to label live cells (x-axis, Q4). Double positive signals might indicate dying cells or incompletely separated cells (Q2). This method allows removing debris (Q3) efficiently. (B) The cell population P4 (i.e. 16,208 living cells) were isolated and used for the scRNAseq run using 10x Genomics. Fig. S2. Quality and quantity of cDNA after reverse transcription of mRNA fraction (polyA-based enrichment) and full-length cDNA amplification from cell lysate of 30 cells sorted from cell suspension of eye-antennal discs run on Fragment Analyzer (Agilent). Size distribution of all fragments shows little impact on degradation (almost no cDNA detectable below 400 bp). Fig. S3. Contribution of mitochondrial gene expression to scRNAseq and snRNAseq datasets. (A) Total amount of genes (features) over percentage of mitochondrial reads, per cell each. The dashed line indicates a threshold of 10% of reads attributed to mitochondrial genes. In scRNAseq data, approximately 14% of cells show a high (>10%) proportion of mitochondrial gene reads on the total number of reads. (B) Total amount of genes (features) over Percentage of mitochondrial reads, per cell each in snRNAseq data. The dashed line indicates a threshold of 10% of reads attributed to mitochondrial genes. In most nuclei, only a low percentage of reads is attributed to mitochondrial genes. Fig. S4. Clustering and cluster annotation for scRNAseq data. (A) The heatmap shows the score for each potential cell type (Y-axis) in each cluster (X-axis). The cell types are annotated based on the highest scoring identity in the heatmap. The clusters are grouped based on their transcriptional similarity to each other. For clusters which express an equal number of marker genes for two different identities both identities were assigned (e.g. cluster 3:Antenna|EAB). Clusters with unresolved identities (i.e. more than two equal assignments) are called “Other”. The colors of the cluster names correspond to the colors in UMAP in (B). The marker score is calculated using a matrix of published marker genes (see Additional file Table S2). (B) UMAP of scRNAseq data. The clusters were annotated based on the heatmap in (A). This UMAP is identical to the UMAP with combined cluster annotation shown in Fig. 3A. Note that the color code is not comparable to the one used in Additional file Figs. S7 and S10. Fig. S5. Fluorescence intensity curves from Bioanalyzer for fresh- and cryopreserved nuclei obtained by different nuclei extraction protocols. (A) The RNA was extracted directly from a fresh sample (36 eye-antennal discs), which was dissociated using the 10x Genomics protocol with 0.1% IGEPAL as a detergent. (B) RNA isolated from a cryopreserved sample, which was dissociated using a protocol based on Triton X-100 as a detergent and a variety of RNAse inhibitors . Note that the sample was thawed for 3.5h before being frozen again. Both curves are close to the expectation of RNA isolated from D. melanogaster . Fig. S6. BioAnalyzer results comparing different nuclei isolation protocols for frozen samples. Samples 1, 2, 4 and 5 were dissociated using the protocol based on Triton X-100 as a detergent and a variety of RNAse inhibitors . Samples 1 and 2 were dissociated by pipetting up and down and samples 4 and 5 were dissociated using a Dounce homogenizer. Samples 7 and 8 were dissociated using the protocol “10x Genomics® Isolation of Nuclei for Single Cell RNA Sequencing” . Samples 7 and 8 were dissociated using only citric acid buffer and samples 10 und 11 were dissociated using only a detergent. Note that for D. melanogaster, intense bands are expected at about 40s and weaker bands at 25s and 45-50s . Each run was repeated once. Fig. S7. Clustering and cluster annotation for snRNAseq data. (A) The heatmap shows the score for each potential cell type (Y-axis) in each cluster (X-axis). The cell types are annotated based on the highest scoring identity in the heatmap. The clusters are grouped based on their transcriptional similarity to each other. For clusters which express an equal number of marker genes for two different identities both identities were assigned (e.g. cluster 15:MFurrow|SMW). Clusters with unresolved identities (i.e. more than two equal assignments) are called “Other”. The colors of the cluster names correspond to the colors in UMAP in (B). The marker score is calculated using a matrix of published marker genes (see Additional file Table S2). (B) UMAP of snRNAseq data. The clusters were annotated based on the heatmap in (A). This UMAP is identical to the UMAP with combined cluster annotation shown in Fig. 3C. Note that the color code is not comparable to the one used in Additional file Figs. S4 and S10. Fig. S8. Gene ontology enrichment analysis for genes with most variable expression. (A) Top 3000 genes unique to scRNAseq (i.e. 1520 genes). (B) Top 3000 genes unique to snRNAseq (i.e. 1520 genes). (C) Top 3000 genes shared between scRNAseq and snRNAseq (i.e. 1480 genes). See also Supplementary Table S9 for a full list of enriched GO terms. Fig. S9. Comparison of dataset specific and shared differentially expressed genes for each cell type. The radar plot shows for each cell type the percentage of cluster specific differentially expressed genes unique for the scRNAseq and snRNAseq data, respectively (red and blue lines), as well as the percentage of differentially expressed genes shared between both datasets (black line). The total number of genes fulfilling the differential expression criteria (FDR 0.05 and log2-fold change > 0.25) for each cell type is shown in brackets. Fig. S10. Clustering and cluster annotation of integrated scRNAseq and snRNAseq dataset. (A) The heatmap shows the score for each potential cell type (Y-axis) in each cluster (X-axis). The cell types are annotated based on the highest scoring identity in the heatmap. The clusters are grouped based on their transcriptional similarity to each other. For clusters which express an equal number of marker genes for two different identities both identities were assigned (e.g. cluster 14:Dorsal|Ocelli). Clusters with unresolved identities (i.e. more than two equal assignments) are called “Other”. The colors of the cluster names correspond to the colors in the UMAP in (B). The marker score is calculated using a matrix of published marker genes (see Additional file Table S2). (B) UMAP of integrated scRNAseq and snRNAseq data. Cells are colored by clusters identified based on the (A). Note that the color code in A and B is not comparable to the one used in Additional file Figs. S4 and S7. Table S1. Overview of different dissociation conditions. Samples within blocks (highlighted in grey and white) were prepared in parallel. The Flow Cytometer only provides percentages of survival because it stops after a defined number of events (i.e. ~50,000 cells) and therefore absolute numbers are not meaningful. “Pipetting” refers to the number of strokes during and after incubation. The cells obtained by experiment/block 12 were subjected to 10X Genomics scRNAseq.
Score Matrix used to annotate cell types and list of references for individual marker genes used for cluster/cell type annotation.
Results of automatic cluster annotation and assignment to combined clusters.
List of marker genes for each cluster in scRNAseq analysis.
GO analysis of cell clusters identified in scRNAseq.
List of marker genes for each cluster in snRNAseq analysis.
GO analysis of cell clusters identified in snRNAseq.
List of top 3000 variable genes for scRNAseq and snRNAseq data, respectively.
Top 3000 variable genes and GO enrichment results for scRNAseq, snRNAseq and shared.
Comparison of differentially expressed genes for each comparable cluster between scRNAseq and snRNAseq data.
About this article
Cite this article
Wiegleb, G., Reinhardt, S., Dahl, A. et al. Tissue dissociation for single-cell and single-nuclei RNA sequencing for low amounts of input material. Front Zool 19, 27 (2022). https://doi.org/10.1186/s12983-022-00472-x
- Single-cell RNAseq
- Single-nuclei RNAseq
- Drosophila melanogaster
- Eye-antennal disc