Abstract
Single-cell transcriptomic approaches are revolutionizing neuroscience. Integrating this wealth of data with morphology and physiology, for the comprehensive study of neuronal biology, requires multiplexing gene expression data with complementary techniques. To meet this need, multiple groups in parallel have developed “Patch-seq,” a modification of whole-cell patch-clamp protocols that enables mRNA sequencing of cell contents after electrophysiological recordings from individual neurons and morphologic reconstruction of the same cells. In this review, we first outline the critical technical developments that enabled robust Patch-seq experimental efforts and analytical solutions to interpret the rich multimodal data generated. We then review recent applications of Patch-seq that address novel and long-standing questions in neuroscience. These include the following: (1) targeted study of specific neuronal populations based on their anatomic location, functional properties, lineage, or a combination of these factors; (2) the compilation and integration of multimodal cell type atlases; and (3) the investigation of the molecular basis of morphologic and functional diversity. Finally, we highlight potential opportunities for further technical development and lines of research that may benefit from implementing the Patch-seq technique. As a multimodal approach at the intersection of molecular neurobiology and physiology, Patch-seq is uniquely positioned to directly link gene expression to brain function.
Introduction
Single-cell approaches have been a cornerstone of neuroscience research for over a century, beginning with the Golgi staining (Golgi, 1873) and the exquisite descriptions of neuronal morphologies by Ramón y Cajal (1909; De Carlos and Borrell, 2007). Later on, the development of electrophysiological recording techniques allowed the study of individual neuronal responses (Hodgkin and Huxley, 1952; Hubel and Wiesel, 1959). Among these, the patch-clamp technique enabled isolated access to an individual neuron to study the ionic currents traversing the plasma membrane (Neher and Sakmann, 1976; Sakmann and Neher, 1984). Initial studies of the gene expression patterns of specific neuronal types relied on the purification of genetically labeled cells, followed by mRNA preparation and expression microarrays or RNA sequencing (Sugino et al., 2006; Okaty et al., 2009, 2011; Fishell and Rudy, 2011). These experiments looked at bulk samples of hundreds to millions cells, and so were limited by both the availability of known cell type-specific markers and the population average of gene expression levels. More recently, single-cell methods have entered the transcriptomic domain and revolutionized neuroscience, providing unparalleled insight into numerous neuronal features, such as cell type diversity (Zeisel et al., 2015, 2018; Gokce et al., 2016; Tasic et al., 2016, 2018; Stanley et al., 2020), developmental and differentiation trajectories (Habib et al., 2016; Kee et al., 2017; Nowakowski et al., 2017; Mi et al., 2018; Raj et al., 2018), and activity state (Chevée et al., 2018; Hrvatin et al., 2018; Jaeger et al., 2018; Sathyamurthy et al., 2018), all revealed by genome-wide gene expression patterns. However, the overall complexity and diversity of neurons often requires a combination of methodological approaches to address experimental questions. For example, recent developments in spatial transcriptomics have started profiling the expression of hundreds or thousands of genes in neurons, while preserving tissue architecture (Chen et al., 2015; Wang et al., 2018; Eng et al., 2019; Rodriques et al., 2019; Soldatov et al., 2019), and implementations of well established techniques for the combined evaluation of morphologic features and electrophysiological properties in the same individual neurons have been fundamental for improving our understanding of the cell types that are the building blocks of neuronal circuits (Jiang et al., 2015; Markram et al., 2015; Zeng and Sanes, 2017; Gouwens et al., 2019).
The study of genome-wide gene expression patterns concomitantly with other neuronal features has only recently become possible. Initial steps toward establishing direct links between electrophysiological properties and gene expression involved the profiling of a small number of genes using RT-PCR (Eberwine et al., 1992; Lambolez et al., 1992; Sucher and Deitcher, 1995; Chiang, 1998; Sucher et al., 2000; Toledo-Rodriguez et al., 2004; Pfeffer et al., 2013; Toledo-Rodriguez and Markram, 2014) or higher-throughput expression microarrays (Subkhankulova et al., 2010) on cytosolic aspirates sampled using the patch pipette. Though groundbreaking, these experiments were limited to small numbers of genes and/or cells. More recently, a number of groups have simultaneously developed Patch-seq, a multimodal method that combines electrophysiological recordings with single-cell whole-transcriptome profiling in the same individual cells (Fig. 1). The first attempt at single-cell RNA sequencing (scRNAseq) from neurons collected with a patch pipette profiled a very small number of samples (Table 1) and yielded sequencing data of low quality (Qiu et al., 2012). Improvements in scRNAseq (Picelli et al., 2014b) and adaptations of the patch-clamp technique facilitated the successful and reliable implementation of Patch-seq (Bardy et al., 2016; Cadwell et al., 2016; Chen et al., 2016; Földy et al., 2016; Fuzik et al., 2016). Following these initial studies, the Patch-seq technique has been further optimized and applied to a number of model systems (Table 1). Moreover, these new types of experiments necessitated the development of innovative data analysis tools suitable for handling multimodal datasets (Kobak et al., 2018; Gala et al., 2019; Bernaerts et al., 2020).
Summary of published Patch-seq datasets
Schematic workflows of Patch-seq using various experimental model systems. A, On acute brain slices, electrophysiological recordings are performed on identified neurons, and the cell processes are filled with dye for subsequent morphologic reconstruction. During sample collection, most of the cytoplasmic contents are aspirated (generally including the nucleus) and transferred to an individual tube containing lysis buffer. The single-cell transcriptomes are then obtained using scRNAseq protocols. B, Electrophysiological recordings can also be performed on neuronal cultures (e.g., iPSC-derived neurons). Simultaneous, live imaging of the recorded neuron is possible if the recording pipette contains a fluorescent dye. After recording and imaging, the whole cell is collected using the patch pipette and the transcriptome obtained via scRNAseq. C, In vivo patch-clamp recordings are used to study the responses of individual neurons to sensory stimuli. Using the patch pipette, each neuron can then be collected and processed for scRNAseq and for in vivo Patch-seq.
In this review, we outline both experimental and analytical strategies for generating and interpreting high-quality Patch-seq data. We provide examples of new biological insights that have resulted from this transformative technique and further discuss outstanding questions especially well suited to the rich datasets produced by this multimodal approach.
The Patch-seq experimental workflow
Patch-seq experiments profile the electrophysiological properties and transcriptome of the same individual neurons, with the goal of identifying the underlying relationships between gene expression and neuronal function. Additionally, recorded neurons can be backfilled with appropriate dies to evaluate cell morphology.
Sample collection
Patch-seq experiments begin with whole-cell patch-clamp recordings from individual neurons to characterize their electrophysiological properties. There are a number of published protocols for Patch-seq (Table 1) that have been optimized for collecting mRNA from neurons recorded in acute brain slices (Cadwell et al., 2017; Lee et al., 2020; Lipovsek et al., 2020; Mahfooz and Ellender, 2021), or in vitro cell cultures (van den Hurk et al., 2018), and allowing subsequent morphologic reconstruction (Cadwell et al., 2017; van den Hurk et al., 2018; Scala et al., 2019, 2020; Gouwens et al., 2020; Que et al., 2020). Below, we discuss commonalities and differences between these protocols.
Given the ubiquitous nature of mRNA-degrading RNases in the environment and the small amount of mRNA present in each individual cell, Patch-seq protocols emphasize the importance of maintaining strict RNase-free conditions during sample collection and before reverse transcription (RT). This includes meticulously cleaning any surfaces that may come into contact with the mRNA sample or internal solution using an RNase decontamination solution and maintaining strict RNase-free conditions while preparing all pre-RT solutions.
Standard electrophysiological protocols can be performed in Patch-seq experiments to quantify the physiological properties of interest. The diameter of the recording pipette tip may need to be adjusted (typically, slightly larger tips are used; but see the study by Lipovsek et al., 2020) to facilitate the aspiration of cell contents into the pipette. Additionally, smaller than normal volumes of internal solution are used to avoid dilution of the RNA and reduce the surface area of glass in contact with the sample, which can lead to sample loss. Finally, the internal solution can be modified to further promote RNA preservation by adding calcium chelators, RNA carriers, and/or RNase inhibitors (Cadwell et al., 2017; van den Hurk et al., 2018; Lee et al., 2020; Lipovsek et al., 2020; Mahfooz and Ellender, 2021). In light of these modifications, the recordings may be more challenging compared with standard patch-clamp conditions and the overall recording time may need to be shortened to ensure minimal RNA degradation. In addition, the use of calcium chelators in the internal solution to reduce RNase activity (Sucher and Deitcher, 1995; Sucher et al., 2000) may also reduce the amplitude of calcium-dependent currents. With these caveats, there is no limit to the types of electrophysiological experiments that can be performed before aspiration of the cell contents for transcriptomic analysis.
After electrophysiological recordings are completed, the cytosol is aspirated into the recording pipette by applying negative pressure. The main objective is to collect the cell contents while maintaining a stable electrical seal between the pipette and the cell membrane, thus avoiding contamination from the extracellular solution. In recent protocols, it has been demonstrated that extracting both the cytosol and the nucleus, for example, using the nucleated patch technique, can greatly improve RNA yield and transcriptomic data quality (Gouwens et al., 2020; Lee et al., 2020; Lipovsek et al., 2020; Scala et al., 2020). After collection, the cell contents are deposited in individual tubes containing lysis and/or RT buffer.
During recording, neurons can be filled with biocytin, enabling subsequent recovery of detailed cellular morphology (Cadwell et al., 2017; Fig. 1A). While this was not realized by the earliest applications of Patch-seq (Cadwell et al., 2016; Fuzik et al., 2016), recent adaptations to the internal and external solutions have substantially improved cell filling (Cadwell et al., 2017; Gouwens et al., 2020; Lee et al., 2020; Scala et al., 2020). Subsequent staining and morphologic reconstruction are further enhanced by the slow and careful retraction of the pipette during the sample collection step after electrophysiological recordings (Gouwens et al., 2020; Lee et al., 2020). Alternatively, if rhodamine is added to the intracellular solution, fluorescent images of the live neuron can be acquired before sample collection (Bardy et al., 2016; van den Hurk et al., 2018; Fig. 1B).
The potential applications of Patch-seq approaches span model systems, including traditional rodent models (Cadwell et al., 2016; Földy et al., 2016; Fuzik et al., 2016, 2019; Muñoz-Manchado et al., 2018; Ellender et al., 2019; Laboissonniere et al., 2019; Oberst et al., 2019; Scala et al., 2019, 2020; Wang et al., 2019; Winterer et al., 2019; Cadwell et al., 2020; Gouwens et al., 2020; He et al., 2020; Liu et al., 2020; Oláh et al., 2020; Que et al., 2020), neurons derived from induced pluripotent stem cells (iPSCs; Bardy et al., 2016; Chen et al., 2016), and less genetically tractable models such as nonhuman primates (Bakken et al., 2020) and primary human tissue (Berg et al., 2020; Kim et al., 2020b). Patch-seq protocols have also been applied to non-neuronal pancreatic islet cells (Camunas-Soler et al., 2020). The basic procedure is very similar regardless of the model system, although the patching and sample collection techniques are tailored to the specific features of the tissue and cell type at hand (e.g., cell size, number of projections, age, myelinization). For example, the negative pressure applied to aspirate the cell contents needs to be carefully calibrated. In some instances (e.g., nonhuman primates and human tissue), cell viability is more variable, and so the extraction pressure needs to be adapted to each individual cell, with typically higher pressure applied for larger cell sizes. Also, if the neuron being extracted has a large number of projections embedded within the tissue, a careful balance between high and low negative pressure is needed to extract the cell contents without tearing the delicate processes.
The application of Patch-seq to the study of human iPSC-derived neurons in vitro was developed independently (Bardy et al., 2016; Chen et al., 2016; van den Hurk et al., 2018) but shares features with Patch-seq methods tailored for acute brain slices (van den Hurk et al., 2018). The main advantages of the cell culture-based approach reside in the ability to perform high-resolution imaging of the live cells for morphology reconstructions and an enhanced capability to collect RNA from the entire neuron (Fig. 1B). The latter provides a more accurate representation of the complete transcriptional program of the cell, including mRNAs localized to synapses or dendrites (Cajigas et al., 2012), otherwise hidden in most other single-cell analyses that are restricted to the mRNA localized to the soma or the nucleus. Extraction of the complete cell (including soma and neurites) may prove more or less technically challenging depending on the stiffness of extracellular matrix, the size of the dendritic arborization, and neuronal adherence. In summary, Patch-seq is a versatile methodology, and the choice of protocol and model system depends on the experimental question at hand.
cDNA synthesis, library preparation, and sequencing
Having collected the contents from individually recorded neurons, subsequent steps aim to generate quantitative gene expression information using scRNAseq workflows. These involve cDNA synthesis and preparation of sequencing libraries. The protocol to be used is chosen before sample collection, as this will dictate the composition of the collection buffer. The cell contents aspirated using the patch recording pipette can be deposited into lysis buffer or on a buffer containing components for the RT reaction. Depending on the collection buffer and experimental workflow, RT reactions are performed on the day of recording or, following short-term storage at −80°C, RT reactions can be performed in batches after sample collection is completed.
Patch-seq is not compatible with high-throughput droplet-based scRNAseq methods, such as 10× Genomics or Drop-seq (Macosko et al., 2015), but can be combined with tube/plate-based methods for the preparation of cDNA and sequencing libraries. Synthetic mRNA spike-ins [e.g., ERCC (External RNA Controls Consortium) spike-ins; Jiang et al., 2011] can be included at this point, or in the internal solution used for electrophysiological recordings, to help control for technical variability. The most commonly used cDNA synthesis protocols are commercial (e.g., SMART-seq version 4, Clontech) or “off-the-shelf” versions of the popular Smart-Seq2 protocol (Picelli et al., 2014b). These are based on full-length cDNA synthesis and amplification using a template-switching oligonucleotide, followed by amplification of full-length cDNA. Other protocols allow for the incorporation of unique molecular identifiers (UMIs) followed by isolation of 5′-end fragments (Zeisel et al., 2015). The choice between full-length cDNA versus UMI-based approaches will depend on the experimental questions and data analysis workflow. The former is more sensitive and provides full-length coverage, allowing for analysis of splice variants (Ziegenhain et al., 2017), while the latter provides absolute numbers of molecule counts, correcting for library amplification bias. The newly developed Smart-Seq3 protocol combines full cDNA coverage and UMIs, greatly improving mRNA quantification and identification of isoforms (Hagemann-Jensen et al., 2020), and is compatible with Patch-seq workflows.
The preparation of sequencing libraries is usually performed using tagmentation-based methods, with either commercially available kits (e.g., Illumina Nextera XT) or off-the-shelf reagents (Picelli et al., 2014a,b). Both the cDNA and sequencing libraries undergo standard quality control steps. These may include capillary electrophoresis (e.g., Bioanalyzer or TapeStation) and quantification (e.g., Qubit) of amplified cDNA from individual cells (Cadwell et al., 2017; van den Hurk et al., 2018), as well as quantitative RT-PCR for housekeeping or marker genes and/or technical spike-ins (if included in the preparation; van den Hurk et al., 2018). Low-quality samples can be excluded from sequencing to conserve resources. Finally, short-read sequencing is typically performed on Illumina platforms. cDNA/library preparation and sequencing amount to most of the costs of running Patch-seq experiments. The choice between using commercially available kits or off-the-shelf reagents (see above) also greatly impacts overall expenses. For example, the cost of the SMART-seq version 4 kit (Takara) ranges from $40 to $90 per cell (excluding sequencing and equipment costs), whereas using off-the-shelf reagents can reduce the cost of library preparation to ∼$20 per cell (Cadwell et al., 2017; Lee et al., 2020).
Data processing and analysis
Bioinformatics pipelines for low-level processing of Patch-seq data are typically identical to those currently in use for dissociated single-cell and single-nucleus RNA-seq data (Lun et al., 2016; Luecken and Theis, 2019; Vieth et al., 2019). In brief, sequenced reads are first aligned to the reference genome and then quantified, typically at the gene level. In Patch-seq, up to ∼50% of reads come from unspliced transcripts (presumed to originate from the cell nucleus; Cadwell et al., 2017; Scala et al., 2020), similar to single-nucleus RNA-seq (Bakken et al., 2018), and so it can be beneficial to use both exonic and intronic read counts to quantify gene expression (Gouwens et al., 2020; Lee et al., 2020; Scala et al., 2020).
Following mapping and read quantification, there are a number of sample-level quality metrics that should be considered for Patch-seq. These include standard metrics commonly used for scRNAseq, such as the number of reads and genes detected per cell, the fraction of reads that map to the genome, the proportion of reads aligning to mitochondrial genes and spike-ins (if included), and the distribution of reads across intergenic, exonic, and intronic regions (Lun et al., 2016; Luecken and Theis, 2019; Vieth et al., 2019). These assessments are especially important when laboratories are in the early stages of implementing Patch-seq, as the harvesting and sequencing of degraded mRNA often result in low-quality data. Additional quality control specific to Patch-seq involves quantifying mRNA contamination because of inadvertent sampling of mRNA from off-target sources, such as surrounding cells and neuropil (Tripathy et al., 2018; Gouwens et al., 2020; Lee et al., 2020; Scala et al., 2020). Using high-quality dissociation-based scRNAseq as a reference, noticeable glial contamination (e.g., microglial markers) has often been observed in Patch-seq samples obtained from acute brain slices. Computational approaches have been developed for quantifying and addressing the effect of these confounds in the context of downstream analyses (Tripathy et al., 2018; Lee et al., 2020).
Data analysis workflows will depend on the specific research questions being asked. When the objective is to study cell types, a common first step is to “map” Patch-seq transcriptomes to an existing transcriptomic reference atlas (Scala et al., 2019, 2020; Cadwell et al., 2020; Gouwens et al., 2020; Que et al., 2020). This is a key analysis step for classifying Patch-seq samples, as dissociation-based scRNAseq datasets typically encompass several orders of magnitude more cells and therefore provide more robust transcriptomic identification of cell types. Despite differences in expression profiles between Patch-seq and reference scRNAseq data, accurate mapping is typically possible using correlation-based methods (Fuzik et al., 2016; Stuart et al., 2019; Berg et al., 2020; Gouwens et al., 2020; Scala et al., 2020). This allows confident assignment of reference transcriptomic types to each Patch-seq cell and their positioning onto a reduced dimension (tSNE or UMAP) reference visualization (Kobak and Berens, 2019). Patch-seq cells with poor or inconsistent mapping to the reference atlas can then be discarded in an additional quality control step (Gouwens et al., 2020; Scala et al., 2020).
Finally, given the multimodal nature of Patch-seq data, there is strong interest in dedicated statistical tools for multimodal data exploration and visualization. Several methods relating gene expression to electrophysiological properties and allowing for low-dimensional visualizations have been developed, based on reduced-rank regression (Kobak et al., 2018), coupled autoencoders (Gala et al., 2019), and bottleneck networks (Bernaerts et al., 2020). It remains a challenge to include morphologic data in these statistical frameworks because of lower sample sizes and difficulties with adequately representing neural morphologies as feature vectors.
Using Patch-seq to address novel biological questions
The advantages of targeted cell sampling
Though relatively low throughput, the Patch-seq sample collection presents one crucial advantage over standard dissociation-based scRNAseq approaches: the experimenter knows exactly from which cell the transcriptomic sample is being collected (Fig. 2A). This knowledge has been exploited to address multiple questions by pairing transcriptomes with spatial context, morphology, connectivity, and functional phenotypes (to name a few). For example, sample collection using the Patch-seq approach identified correlations between transcriptomes and neuronal position (Pfeffer and Beltramo, 2017; Gouwens et al., 2020; Scala et al., 2020), projection pattern, sensory response type (Pfeffer and Beltramo, 2017), and long-range axonal connectivity (Luo et al., 2019). The transcriptomes of individual cells identified solely by their functional response properties, such as thermosensitive neurons of the preoptic area of the hypothalamus (Wang et al., 2019), light-sensitive and insensitive neurons of the visual cortex (Liu et al., 2020), and glucose-inhibited/excited neurons of the ventromedial hypothalamic nucleus (He et al., 2020) have been contrasted, via Patch-seq, for the first time, revealing new molecular markers and genes involved in differential functional responses.
Diagrams of typical analytic workflows for Patch-seq experiments. A, Patch-seq allows the direct targeting of identified cells of interest (e.g., genetically labeled, retrogradely labeled by targeting the projection site or by their unique responses to a stimulus). B, The gene expression profile of Patch-seq samples can be used to map cells onto a reference scRNAseq dataset, thus assigning a multimodal phenotype to transcriptomically identified cell clusters. C, Multimodal Patch-seq data can be used to train a machine learning algorithm that can be used to predict the functional state of cells collected in independent scRNAseq datasets. D, Patch-seq data can be used to identify gene expression patterns that correlate with morphoelectric properties and other phenotypes.
Targeted Patch-seq experiments also allow studying neurons with a shared developmental lineage. For example, Patch-seq revealed that cortical excitatory neurons derived from a common early progenitor are transcriptomically diverse and no more likely to belong to the same transcriptomic cell type than randomly selected neurons in similar anatomic locations (Cadwell et al., 2020). However, Patch-seq performed on excitatory neurons derived from different progenitor pools at later developmental stages revealed that apical intermediate progenitors contribute to a transcriptionally distinct layer 2/3 neuronal subtype (Ellender et al., 2019). Finally, in combination with elegant transplantation experiments, sample collection via Patch-seq identified a prominent role for Wnt signaling in the temporal specification of apical neuronal progenitors (Oberst et al., 2019).
A Rosetta stone for integrating multimodal atlases of cell type diversity
The latest and largest scRNAseq studies of the mammalian brain describe hundreds of different transcriptomically defined cell types (La Manno et al., 2020; Yao et al., 2020), identified by clustering analysis. However, for many of these, the morphoelectric identity remains unknown. Patch-seq approaches are bridging this gap by directly profiling morphologic and electrophysiological phenotypes, alongside the transcriptomes, followed by mapping these multimodally defined cells onto comprehensive scRNAseq atlases (Fig. 2B). The first multimodal Patch-seq studies identified morphoelectrically and transcriptomically distinct neurons in layers 1 and 2 of the cortex (Cadwell et al., 2016; Fuzik et al., 2016). The latest multimodal classification efforts have profiled thousands of neurons in the mouse visual and motor cortices, comprising the majority of known transcriptomically defined cell types (Gouwens et al., 2020; Lee et al., 2020; Scala et al., 2020; Table 1).
The identification of distinct cell types is not always straightforward. In the striatum, transcriptomic classification of spiny projection neurons identified continuous as well as discrete variation between cell types (Gokce et al., 2016; Stanley et al., 2020). Application of Patch-seq to striatum interneurons uncovered a continuum of transcriptomic gradients that correlate with a continuum of electrophysiological phenotypes (Muñoz-Manchado et al., 2018), indicating that graded phenotypes may be a defining feature of neuronal diversity, at least in the striatum. In the hippocampus, Patch-seq of morphologically distinct parvalbumin-positive interneurons found them to be biophysically and transcriptionally homogeneous, but identified a developmental change in transcriptomic state that is characterized by the unexpected onset of hemoglobin expression (Que et al., 2020). Further exploiting the potential of Patch-seq data, multimodal analysis refined the classification of layer 4 interneurons of the mouse visual and somatosensory cortex, suggesting that transcriptomically similar neurons can differ in terms of morphology and electrophysiology, depending on their cortical area (Scala et al., 2019).
The picture emerging from cross-modality studies is that no single feature is sufficient to unequivocally identify cell types. A study in the motor cortex suggested that, while broad families of transcriptomic cell types have distinct and essentially nonoverlapping morphoelectric phenotypes, closely related transcriptomic cell types are not well separated in the morphoelectric space (Scala et al., 2020). Also, corticocortical projection neurons of the visual cortex, identified as a single transcriptomically defined cell cluster, correspond to distinct cell types, defined by their anatomic projections and connectivity patterns (Kim et al., 2020a). Finally, continuous variation of morphologic, electrophysiological, and/or transcriptomic features between related cell types can limit the ability to predict one from the other, although using more than one modality can enhance prediction (Gouwens et al., 2020). The examples discussed above highlight the requirement for multimodal profiling to better determine cell identity and have crucial implications for the conceptual and empirical definitions of cell types, subtypes, and states (Zeng and Sanes, 2017).
Uncovering the molecular basis of morphologic and functional diversity
The unbiased nature of single-cell transcriptomics makes Patch-seq ideal for generating novel hypotheses about the molecular underpinnings of morphologic and functional diversity. By targeting recordings to previously known functionally distinct cell types, the transcriptomic basis of such differences can now be potentially elucidated using Patch-seq (Wang et al., 2019; He et al., 2020; Liu et al., 2020). The predictive power of the Patch-seq approach to reveal features of neuronal function within local circuitry was illustrated in CCK+ interneurons. Here, transcriptomic data identified differential expression of 5-HT3 receptor subunits between interneuron types. At the functional level, this resulted in differential responses to serotonin and the further demonstration that, in interneurons that lack 5-HT3 receptors mRNA, responses to serotonin are driven by gap junction coupling to neighboring cells (Fuzik et al., 2016). Similarly, Patch-seq identified novel differentially expressed genes including cell adhesion and synaptic regulatory molecules, that may underlie differences in the connectivity of distinct subtypes of layer 1 interneurons (Cadwell et al., 2016).
Multimodal experiments have identified cell types with distinct functional or morphologic features that show highly similar transcriptomic profiles (Scala et al., 2019; Que et al., 2020; Kim et al., 2020a). In-depth analysis of one such example in the rat hippocampus found that CCK+ interneurons of the CA3 area that display different excitability properties express distinct isoforms of voltage-gated potassium channel auxiliary subunits (Oláh et al., 2020). This demonstrates that subtle differences in transcriptomic profile can have profound consequences for neuronal function.
The potential for Patch-seq to predict functional states has also been demonstrated in human iPSC-derived neuronal cultures (Bardy et al., 2016). Such cultures often comprise neuronal populations of heterogeneous electrophysiological states or neurodevelopmental maturity. Applied to a Patch-seq dataset of iPSC-derived neurons, a machine-learning classifier efficiently predicted the electrophysiological state of individual neurons from their transcriptional profile (Fig. 2C). Once trained on Patch-seq data, such a classifier was applied to thousands of neurons profiled by high-throughput scRNAseq, allowing the study of functional heterogeneity in iPSC-derived cultures (Bardy et al., 2016). Similar machine learning approaches for predicting aspects of neuronal physiology from genome-wide gene expression have also been implemented on layer 1 cortical interneurons (Cadwell et al., 2016). Patch-seq can further be used to identify specific genes whose expression correlates strongly with electrophysiological properties (Fig. 2D). For example, mRNA level of Kcna1, which codes for subunits of the Kv1.1 voltage-gated potassium channel, correlates strongly with cell-to-cell variability in action potential width among inhibitory cell types in the hippocampus and striatum (Bomkamp et al., 2019).
Finally, Patch-seq methods are not restricted to neurons. An extensive study of human pancreas samples combined electrophysiological measurements of exocytosis and channel activity with scRNAseq of pancreatic β-cells (Camunas-Soler et al., 2020). Subsequent analyses of cell function and gene expression networks identified a gene set associated with functional heterogeneity that can be used to predict endocrine physiology. The authors go further to demonstrate that Patch-seq can be used to study cryopreserved tissue (Camunas-Soler et al., 2020). The application of these novel approaches to patient samples holds great promise for future research and therapeutic avenues for the treatment of diabetes.
Outstanding questions and future directions
Patch-seq is a multimodal approach that can simultaneously profile multiple aspects of neuronal biology, namely their transcriptome, electrophysiological properties, anatomic location, and morphology. As such, it holds immense potential to address a multitude of outstanding questions in neuroscience. Approaches using Patch-seq have already made significant contributions to the identification of cell types and subtypes (Muñoz-Manchado et al., 2018; Gouwens et al., 2020; Scala et al., 2020), providing an essential backdrop for current and future neuroscience research.
The Patch-seq methodology presents three significant challenges. First, it is relatively low throughput, as the number of cells that can be profiled in a given experiment is governed by the practical considerations of patch-clamp electrophysiology (as opposed to the high throughput of droplet-based scRNAseq), made more cumbersome by the specific requirements of Patch-seq sample collection. Second, Patch-seq is a laborious technique to learn and master. However, recent protocol improvements have substantially improved the success rates of critical mRNA sample collection and morphology recovery steps (Cadwell et al., 2017; Gouwens et al., 2020; Lee et al., 2020; Lipovsek et al., 2020; Scala et al., 2020). Last, Patch-seq requires interdisciplinarity, as it yields multiple types of data, each calling for their own skillset to analyze and interpret. Higher-throughput methods can individually profile the modalities targeted by Patch-seq. Using dedicated data analysis tools, these can be combined post hoc to start exploring their relationships (Tripathy et al., 2017; Butler et al., 2018). However, Patch-seq is unique in the sense that all data modalities are acquired from each individual cell simultaneously, providing unique and integrative opportunities for cross-modal data exploration.
Overall, Patch-seq is a versatile methodology and can be applied to any sample suitable for patch-clamp experiments. Decades of electrophysiological research provide virtually endless tried experimental model systems. Going beyond traditional acute mouse or rat brain slices, Patch-seq has already been implemented on neurosurgically resected human tissues to identify increased glutamatergic neuron diversity in the human cortex (Berg et al., 2020) and in a comparative study of mouse, marmoset, and human motor cortex (Bakken et al., 2018).
Improvement, refinement, and modifications of patch-clamp recordings and sequencing methodologies will undoubtedly extend the reach of Patch-seq methods. To date, Patch-seq has already been applied to in vivo recordings (Cadwell et al., 2016; Pfeffer and Beltramo, 2017; Liu et al., 2020), thus opening up the opportunity to co-opt invaluable experimental paradigms (e.g., learning and memory, sensory processing, and motor control) into multimodal workflows to ultimately study the relationships between gene expression and behavior. Furthermore, incorporating additional steps after sample collection to separate the mRNA from the rest of the cell contents (e.g., G&T protocol for simultaneous genome and transcriptome profiling; Macaulay et al., 2016; Lipovsek et al., 2020) could provide the opportunity to study additional modalities. For example, chromatin accessibility or DNA methylation patterns can be profiled on the genomic DNA (Smallwood et al., 2014; Buenrostro et al., 2015) to investigate the epigenomic correlates of functional phenotypes, and advances in single-cell proteomics will enable assaying protein abundance levels (Labib and Kelley, 2020). Finally, Patch-seq could be combined with CRISPR-based gene expression perturbation experiments (Gomez et al., 2019) to simultaneously study the transcriptomic and functional effects of specific mutations at the single-cell level.
Future experimental questions will undoubtedly drive protocol innovations. Placed at the intersection between the functional and molecular domains, Patch-seq approaches are ideally suited to explore open questions on cell state dynamics. For a defined cell type, how is the functional state affected by external input? Which are the key genes that define a functional state? How does a change in functional state relate to changes in transcriptomic state? Addressing these questions will contribute to identifying the molecular players shaping the landscape of possible states a given cell can adopt during development, aging, and disease (Trapnell, 2015).
A yet untapped area of investigation for Patch-seq are the transcriptomic mechanisms of synapse function and plasticity. These experiments would require improvements in both recording conditions and scRNAseq methods to mitigate the effects on mRNA stability of extended recording periods and to capture transcriptomic changes related to mRNA transport and localization. Additionally, Patch-seq approaches can be a starting point for studying circuit assembly. For example, a description of a cell type-specific code of synaptic and cell adhesion molecules and its relationship with a functional phenotype could be used to predict circuit connectivity from single cell-level transcriptomic data (Földy et al., 2016). With the expansion of Patch-seq datasets, this approach can first be validated on known connected neurons and subsequently exploited to predict hitherto unknown connections. Finally, multi-Patch-seq experiments, in which pairs (or groups) of neurons are simultaneously recorded to determine their connectivity and subsequently collected for sequencing, would allow a direct correlation between gene expression and connectivity of transcriptomically and morphoelectrically defined neurons.
Another exciting area of investigation is using Patch-seq to disentangle the causality between molecules and function. For example, while many genes might appear correlated with a quantitative functional feature, such as an electrophysiological property (Cadwell et al., 2016; Földy et al., 2016), pinpointing which of these genes is likely to drive the causal relationship is a major goal of multimodal experimental approaches (Tripathy et al., 2017; Bomkamp et al., 2019). Addressing these questions will require further experimental manipulations alongside the development of new data analysis frameworks to distill understanding from complex multimodal datasets.
Overall, multimodal approaches to study cell types and function are becoming the norm, as new experimental and analytical tools are developed, allowing the exploration of long-standing and novel questions in neuroscience. Many of these approaches are focused on increasing the experimental throughput, aiming at profiling the highest number of cells possible. In parallel, lower-throughput approaches that are tailored to addressing specific questions can be equally or more powerful. Though Patch-seq is a laborious approach, by directly interrogating the connections between organizational levels it is already proving invaluable for advancing our understanding of neuronal function.
Footnotes
The authors declare no competing financial interests.
M.L. was supported by a Consolidator Grant from the European Research Council (Grant 725729; FUNCOPLAN) to Matthew Grubb. C.B. was supported by the Ian Potter Foundation, the Flinders Foundation, Impact Philanthropy, Cancer Australia, and the Australian Department of Health. K.H. was supported by the Allen Institute for Brain Sciences. D.K. was supported by the Deutsche Forschungsgemeinschaft (Grant BE5601/4-1 and the Cluster of Excellence “Machine Learning—New Perspectives for Science”, EXC 2064, Project 390727645), the Federal Ministry of Education and Research (Grants FKZ 01GQ1601 and 01IS18039A), and the National Institute of Mental Health of the National Institutes of Health under Award Number U19MH114830. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. S.J.T. was supported by the CAMH Discovery Fund, the Krembil Foundation, the Kavli Foundation, the Canadian Institutes of Health Research (Grant NGN-171423), and the Natural Sciences and Engineering Research Council of Canada (Grants RGPIN-2020-05834 and DGECR-2020-00048). We thank Matthew Grubb for helpful comments on the manuscript and for kickstarting the organization of the mini-symposium.
- Correspondence should be addressed to Marcela Lipovsek at marcela.lipovsek{at}kcl.ac.uk or Shreejoy J. Tripathy at shreejoy.tripathy{at}utoronto.ca
