Challenges and Insights in Patch-Clamp Studies: From Cell-Attached to Whole-Cell Configurations
Sheng-Nan Wu, Ya-Jean Wang, Rasa Liutkevičienė

TL;DR
This paper discusses how patch-clamp configurations affect electrophysiological recordings and how drugs influence potassium channels in cells.
Contribution
The paper provides insights into drug effects on specific potassium channels and the transition to AI in patch-clamp experiments.
Findings
Switching from cell-attached to whole-cell patch-clamp increases current amplitude and abolishes action currents.
Mitoxantrone suppresses Kir channels without altering single-channel conductance in RAW 264.7 cells.
GAL-021 inhibits BKCa channels by shifting activation curves, acting as a gating modifier.
Abstract
The patch-clamp technique is widely regarded as the gold standard in cellular electrophysiology and can be applied in several configurations. In the cell-attached (C-A) mode, it enables the recording of single-channel currents, whereas the whole-cell (W-C) mode allows for the measurement of macroscopic currents, representing the collective activity of many channels. When the recording configuration was switched from C-A to W-C on the same cell, the current amplitude increased dramatically, while action currents (ACs) were completely abolished, indicating a profound alteration in the cell’s electrophysiological response under the new setup. In excitable cells, the occurrence of ACs, representing propagated action potentials, can interfere with C-A single-channel recordings. To address this, a high-K+ solution is typically applied to the bath to suppress the ACs. The inwardly rectifying…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8- —National Science and Technology Council
- —An Nan Hospital
- —the Lithuania-Latvia-Taiwan collaborative project
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsIon channel regulation and function · Cardiac electrophysiology and arrhythmias · Ion Channels and Receptors
1. Introduction
The patch-clamp technique is widely recognized as the gold standard in electrophysiology, owing to its unparalleled ability to provide high-resolution insights into the electrical activity of cells and their ion channels. Its strength lies in the versatility of its experimental configurations, which allow researchers to precisely manipulate and control the conditions on both sides of the cell membrane [1]. The commonly used patch-clamp configurations are cell-attached (C-A), inside-out, outside-out, and whole-cell (W-C) modes [1,2]. Of these, the cell-attached (C-A) and whole-cell (W-C) modes are the two most frequently used. Figure 1 is a schematic cartoon that briefly illustrates the main differences between the two configurations. However, their operation often involves specific challenges, such as the direction of ionic flow across the membrane and the associated rectification properties. Rectification in ion channel dynamics describes the phenomenon where ions pass more readily in one direction than in the other, a behavior analogous to that observed in semiconductor devices. This article aims to present essential information and clear key concepts commonly encountered during operation, serving as a useful reference.
2. Ionic Currents Evoked by Ramp or Rectangular Pulse Recorded Under Cell-Attached (C-A) and Whole-Cell (W-C) Voltage-Clamp Configurations
As illustrated in Figure 2A, during these measurements, INS-1 pancreatic β-cells were maintained in standard Tyrode’s solution, and recordings were obtained using pipettes filled with a K^+^-rich internal solution. INS-1 pancreatic β-cells are a rat insulinoma-derived cell line widely used as a model to study pancreatic β cell physiology, glucose metabolism, and insulin secretion [3,4]. In the C-A configuration, after achieving a gigaohm seal, ramp pulses of 1 s duration evoked only relatively small ionic currents. Notably, a repetitive inward current was observed during the upsloping phase of the ramp, with a threshold potential near −22 mV. The amplitude of this inward current (shown in blue) gradually diminished and broadened as the membrane depolarized.
Upon transition to the W-C configuration following patch rupture, the ramp-evoked currents (shown in red) increased dramatically—approximately 20-fold compared to C-A recordings. In the same cells, membrane resistance between –40 and +40 mV decreased sharply from 2.21 ± 0.03 to 0.089 ± 0.01 GΩ (n = 12, p < 0.05). Interestingly, the trains of action currents (ACs) evident in C-A mode were abolished under W-C conditions, leaving only a small inward current around 0 mV. Furthermore, both inwardly rectifying K^+^ currents of modest amplitude and prominent outwardly rectifying K^+^ currents were detected in response to the ramp protocol.
To better characterize the biophysical properties of ionic currents in INS-1 cells, command voltage was switched to a rectangular pulse protocol. As shown in Figure 2B, under the C-A mode, when the cell was held at –50 mV and depolarized to +50 mV for 1 s, a small outward current (blue trace) was observed together with a train of inward currents (red arrowhead, labeled b). These inward currents exhibited spike-frequency adaptation and are referred to as ACs, corresponding to the generation of propagated action potentials (APs) [5,6,7,8,9]. Spike-frequency adaptation in excitable cells refers to the gradual reduction in firing rate during sustained constant current injection, caused by intrinsic ionic mechanisms that limit excitability over time. Notably, once the patch membrane was ruptured, the current amplitude was increased substantially. Transitioning from C-A to W-C mode resulted in a marked reduction in membrane resistance during depolarization from –50 to +50 mV, decreasing from 3.48 ± 0.09 to 0.35 ± 0.02 GΩ (n = 12, p < 0.05). In the same cells, both a transient voltage-gated Na^+^ current, characterized by rapid activation and inactivation (red arrowhead a), and a delayed-rectifier K^+^ current with slow inactivation were simultaneously recorded. However, while the ACs were observable in the C-A configuration, they were no longer discernible upon transition to the W-C mode. Overall, under these experimental conditions, approximately 70% of INS-1 cells exhibited repetitive AC firing.
3. Measurement of Resting Membrane Potential Under C-A Mode
Determining the resting membrane potential can be difficult to some extent in the C-A mode. Figure 3 provides a suitable example for illustration. When using C-A single-channel recordings to measure inwardly rectifying K^+^ (Kir) channels, if the bath solution is normal Tyrode’s solution (containing ~5.4 mM K^+^) and the recording pipette is filled with a K^+^-rich solution (~140 mM), the reversal potential—that is, the potential at which single-channel current neither flows inward nor outward—becomes positive as a function of Δvoltage (Figure 3B). This occurs because the patch potential is the sum of the resting potential and the pipette potential, while Δvoltage is defined as the potential relative to the bath [10,11,12]. In contrast, under the W-C configuration, the resting membrane potential of the cell is approximately –70 mV (Figure 2).
Figure 3A shows Kir-channel activity recorded from RAW 264.7 osteoclast precursors, which is sensitive to suppression by mitoxantrone. An osteoclast is a large, multinucleated cell responsible for breaking down bone tissue through a process called bone resorption. Mitoxantrone (Novantrone^®^, 1,4-dihydroxy-5,8-bis[2-(2-hydroxyethylamino)ethylamino]anthracene-9,10-dione), a synthetic anthracenedione chemotherapeutic agent used in cancer treatment and multiple sclerosis [14], has been observed to suppress Kir channels [13]. Notably, the channel activity was observed in the C-A mode, when outward currents were recorded at a potential of +50 mV relative to the bath (pipette potential), which, from the perspective of the cell membrane, actually represent inward K^+^ current. This voltage (i.e., patch potential) corresponds to approximately 120 mV relative to the resting membrane potential (50 plus 70 mV). Of note, this form of single-channel current exhibited a relatively long mean open time, making it readily distinguishable from other ion channels present in the same cells [10,13,15].
Figure 3B illustrates the current-voltage (I-V) relationships of Kir channels in lipopolysaccharide-differentiated RAW 264.7 cells, both in the absence and presence of mitoxantrone. Mitoxantrone treatment did not significantly alter the single-channel conductance, as reflected by the unchanged linear slope of the I-V curve. Importantly, the recorded currents are negative, which correspond to inward currents from the cellular perspective—indicating that K^+^ ions move from the extracellular space into the cytoplasm, driven by their electrochemical gradient [13,15].
The Kir2.x subfamily has been identified as the predominant inward rectifier expressed in differentiated RAW 264.7 cells [15]. However, it remains unclear whether mitoxantrone-induced inhibition of Kir channels represents an intrinsic property of neurons or glial cells, including microglia [10,15]. This mechanism is of particular interest, as modulation of Kir-channel activity may hold therapeutic potential in multiple sclerosis [14,16]. Multiple sclerosis is a chronic autoimmune demyelinating disorder in which the immune system attacks the central nervous system, leading to damage of the protective myelin sheath surrounding axons and subsequent impairment of neuronal signaling [17].
4. Differences Between ATP-Sensitive K+ (KATP) Channels and Large-Conductance Ca2+-Activated K+ (BKCa) Channels in the C-A Configuration
When using the C-A single-channel recordings, the holding potential displayed on the instrument (i.e., patch amplifier) is negative, but for the cell being examined, it is actually positive. Therefore, what appears as a downward deflection representing inward current on the recording is, in fact, an upward deflection representing outward current from the cell’s perspective, as demonstrated in Figure 4. Because of this discrepancy, some early studies appeared to mistakenly report K_ATP_ channel activity when they were actually observing the activity of BK_Ca_ channels.
The large-conductance Ca^2+^-activated K^+^ channels—also known as BK, BK_Ca_, or maxi-K channels—are a unique type of channel. They open in response to membrane depolarization and elevated intracellular Ca^2+^, permitting large outward-rectifying K^+^ currents to flow [8,16,18,19]. The single-channel conductance of BK_Ca_ channels in symmetrical K^+^ solution (~140 mM) is 150–250 pS. However, when the bath solution contains low K^+^ (i.e., normal Tyrode’s solution), their conductance can be around 100–150 pS, a value which is close to that of K_ATP_ channel.
The ATP-sensitive K^+^ channels (K_ATP_ channels) are specialized K^+^ channels that couple a cell’s metabolic state to its electrical activity. When intracellular ATP levels are high, these channels close; when ATP levels decline, they open. Through this mechanism, K_ATP_ channels play a critical role in regulate physiological processes, such as hormonal release (i.e., insulin release), vascular tome, and cellular protection during metabolic stress [20].
Both of these channels exhibit mean open times on the order of a few milliseconds, whereas classical Kir channels (such as Kir2.1), as described above in RAW 264.7 cells, typically remain open for several hundred milliseconds (Figure 3), making them much easier to distinguish. It is thus important to note that K_ATP_ channels display bursting kinetics characterized by a mean intraburst open time of ~1 ms, whereas BK_Ca_ channels typically lack this behavior [19]. In addition, rectification properties provide an important means of distinguishing between the two: K_ATP_ channels display mild inward rectification, while BK_Ca_ channels exhibit strong outward rectification [16,18,19]. Figure 4 illustrates a schematic representation of K_ATP_ channel activity.
5. Interference of Single-Channel Activity by Action Currents (ACs) in C-A Current Recordings
When using C-A single-channel recordings, if the bath solution is normal Tyrode’s solution (containing ~5.4 mM K^+^), the recording electrode is filled with K^+^-rich solution, and the holding potential is very negative (e.g., below −80 mV), corresponding to a Δvoltage of about more than +10 mV relative to the bath, the appearance of inward currents observed from the pipette’s perspective may indicate the presence of BK_Ca_-channel currents. However, such currents could also be action currents (ACs), if the cell is highly electrically excitable. In this situation, the single-channel current traces are easily disturbed and distorted [5,6,8,9,19,20,21,22,23,24,25,26]. Consequently, analyzing the single-channel kinetics of BK_Ca_ channels (such as mean open time or mean closed time, derived from open-time or closed-time histograms) becomes severely compromised due to contamination by the concurrent activity of ACs. To minimize interference from ACs, one effective strategy is to replace the bath solution with a high-K^+^ solution (about 140 mM). This adjustment shifts the reversal potential for K^+^ ions toward 0 mV, thereby suppressing AC activity through prolonged depolarization of the cells. It needs to be noted, however, that prolonged exposure to high-K^+^ solution may diminish overall cellular activity.
When using C-A single-channel recordings in electrically excitable cells, ACs often appear (which correspond to the action potentials [APs] observed in W-C current-clamp recordings) [2,5,6,7,8,9]. Typically, their deflections in the same direction as those of BK_Ca_ channels, although their directions are biphasic to some extent. Moreover, the duration of ACs is relatively short, whereas the duration of BK_Ca_ channel openings can vary, sometimes brief and sometimes prolonged. However, when the activity of BK_Ca_ channels is relatively low, their opening duration is also short, making them indistinguishable from ACs.
In W-C recordings, when using the voltage-clamp configuration combined with variable voltage-clamp protocols enables the observation of many distinct types of macroscopic ionic currents. Switching to the current-clamp configuration with the holding current set to 0 pA allows for measurement of the cell’s the resting membrane potential. Furthermore, in electrically excitable cells, appropriate current injection can elicit action potentials (APs), and spontaneous APs may occasionally be observed [2,8,16]. Sustained AP firing can also occur when a constant holding current is applied, thereby depolarizing the membrane potential [6,8].
In the early days, commonly used microelectrode amplifiers such as the Axoclamp-2B (Molecular Devices, San Jose, CA, USA) made switching from single-electrode voltage-clamp mode to current-clamp mode rather cumbersome and difficult to operate. In contrast, many modern patch amplifiers—such as Axopatch 1D, Axopatch 200B (Molecular Devices), EPC-7 (HEKA Elektronik, Lambrecht, Germany), or RK-400 (Bio-Logic, Claix, France)—allow for much easier switching between voltage-clamp and current-clamp modes [2,8,12]. Nevertheless, the observed phenomena are a consequence of the patch-clamp method, independent of the specific amplifier model used.
In C-A voltage-clamp recordings, single-channel events are frequently detectable; however, action currents (ACs), which correspond to action potentials (APs) observed in whole cell current-clamp recordings, may also be recorded [8,9,12,24,25,27,28]. Moreover, a shift from C-A voltage-clamp to current-clamp recording mode often results in observable changes in membrane potential. Yet, under this circumstance, due to differences in the position of the holding current, the true resting membrane potential cannot be accurately measured, particular when the bath solution contains low K^+^ (approximately 5.4 mM K^+^). In addition, under C-A current-clamp mode, spontaneous APs can easily be confused with subthreshold depolarizations or synaptic potentials accompanied by spiking firing. Subthreshold oscillations refer to rhythmic fluctuations in a neuron’s membrane potential that occur below the threshold required to trigger an AP [7,8,24,26,27,28,29,30].
6. Modification of BKCa-Channel Activity by GAL-021 in GH3 Cells
In this study, we further demonstrated that inhibitory effect of GAL-021 on the channel opening of BK_Ca_ channels. The compound GAL-021 (also known as ENA-001, 6-(methoxy(methyl)amino)-N^2^,N^4^-di-propyl-1,3,5-triazine-2,4-diamine) is a small-molecule drug developed as a novel breathing control modulator and a respiratory stimulant. As demonstrated in Figure 5, under high-K^+^ conditions, the activities of these channels can be readily observed in an excised inside-out patch. Notably, in the inside-out configuration, the activity of ACs observed in C-A mode is absent. When the detached patch was exposed to GAL-021, the activity of channel openings was progressively reduced [18,31,32,33]. The open-state probability of the channel measured at the level of +60 mV under control condition (no GAL-021) was found to b1 0.182 ± 0.012 (n = 8). The addition of GAL-21 (10 μM) to the bath medium decreased the channel activity to 0.022 ± 0.002 (n = 8, p < 0.01). However, there was no significant difference in the amplitude of the unitary outward current between the absence and presence of GAL-021. It is clear that the presence of GAL-021 can decrease the opening probability of BK_Ca_ channels in these cells, but not that of K_ATP_ channels [19,33].
The study also investigated the effect of GAL-021 on the single-channel conductance of BK_Ca_ channels using the voltage-clamp technique on inside-out patches. This method allowed for the construction of single-channel amplitude versus membrane potential (I-V) relationships by applying voltage ramp pulses from 0 to +120 mV for 1 s over 1 s. Figure 6A depicts the single-channel amplitude versus membrane potential (I-V) relationships of BK_Ca_ channels with or without the addition of GAL-021 (3 μM). The single-channel conductance of BK_Ca_ channels achieved from the linear I-V relationship in control (no GAL-021) was 165 ± 7 pS (n = 11) with a reversal potential of 0 ± 1 mV (n = 11). The value was not found to differ significantly from that (164 ± 6 pS; n = 11, p > 0.05) measured in the presence of GAL-021 (3 μM). Therefore, GAL-021 produced no significant change in the single-channel conductance of BK_Ca_ channels, but it suppressed the channel activity in these cells, along with a prolongation in mean closed time of the channel [18]. The results suggest a mechanism where GAL-021 stabilizes the channel in a non-conducting (closed) state without physically blocking the open pore.
Figure 6B illustrates the activation curve of BK_Ca_ channels taken from the absence or presence of GAL-021 (3 μM). To determine the effect of GAL-021 on the activation curve of BK_Ca_ channels, the upsloping ramp pulses from +20 to +120 mV with a duration of 1 s were designed and creased from pCLAMP 10.7, and they were then applied to the membrane patch through digital-to-analog conversion. During the single-channel current recordings, the linear ramp pulses were delivered from +20 to +120 mV for 1 s. The activation curve obtained during ramp pulses were calculated by averaging the current traces from 20 voltage ramps, then normalizing each mean current value to the corresponding single-channel amplitude at that potential, after correcting for the leakage component. The plot of the relative open probability as a function of membrane potential applied were thereafter constructed and then fitted with a modified Boltzmann function. In control, V_1/2_ = +66.8 ± 2.3 mV and q = 4.3 ± 0.2 e (n = 7), whereas in the presence of GAL-021 (3 μM), V_1/2_ = +79.3 ± 2.6 mV and q = 4.1 ± 0.2 e (n = 7). Consequently, the presence of GAL-021 (3 μM) caused a 45% decrease in the maximal open probability of BK_Ca_ channels. However, there was no significant effects on the gating charge of the curve occurring during exposure to GAL-021. This study led us to indicate that GAL-021 could suppress the activity of these BK_Ca_ channels in a voltage-dependent fashion in GH_3_ cells [18]. GAL-021 may possess a distinctive structural property that enables functional interaction with the a-subunit of the BK_Ca_ (KCNMA1) channel, thereby reducing the amplitude of whole-cell Ca^2+^-activated K^+^ currents [18].
7. Brain Slice Patch-Clamp Recordings
A brain slice is a thin section of brain tissue, usually prepared with a vibratome, that is kept alive in artificial cerebrospinal fluid. Researchers use brain slices because they preserve much of the brain’s local circuitry while allowing for direct experimental access [9,33,34]. Brain slice patch-clamp recording is a commonly used electrophysiological method. However, because the patch electrode has a relatively large tip diameter, only sufficiently large neurons can be recorded appropriately. In addition, since the cell membrane is very fragile, it is often difficult during the procedure to distinguish whether the recording configuration in W-C mode or C-A mode, as demonstrated in Figure 7. When the electrode tip is too small (e.g., less than 1 μm), the membrane may reseal after membrane rupture in W-C mode, thereby reverting a C-A configuration. Furthermore, interference from excitatory or inhibitory synaptic currents adds to the complexity [2,7,9,12,24,25,26,27,32,33,34,35]. As a result, interpretation of the recordings becomes more challenging, and the reproducibility of experimental outcomes is greatly affected [2,7,24,32,34]. Therefore, in experiments such as brain slice patch clamping, the size of the cell is an important factor for the operator to consider when selecting cells for electrophysiological measurements (Figure 7).
8. Conclusions
In the cell-attached (C-A) configuration, many highly excitable cells readily display spontaneous or stimulation-evoked action currents (ACs) [5,8,9,25,26,27,28]. These ACs often correspond to the action potential firing observed under the whole-cell (W-C) configuration. In addition, the C-A mode allows for the detection of single ion channel activity in various cells. However, channels such as large-conductance Ca^2+^-activated K^+^ (BK_Ca_) channels are influenced by the presence of ACs, which complicates the analysis of their kinetics [6,7,21,22]. Finally, it should to be noted that in brain slice patch-clamp recordings, when the target cells are very small, distinguishing between W-C and C-A modes becomes difficult. Consequently, the assessment of AC activity is also challenging [26,28,31,32,33].
In earlier electrophysiological experiments, amplifiers were typically connected to oscilloscopes to observe the dynamic changes in voltage and current. Subsequently, through a data acquisition board, analog-to-digital conversion was performed to transform analog signals into digital form, allowing for effective storage on personal computers [2,8,36]. Nowadays, digital-to-analog conversion capability has also been incorporated, enabling various digitized voltage profiles—such as specific ramp voltage shapes or realistic action potential waveforms—to be reconverted into analog signals. These analog signals are subsequently employed to precisely regulate the patch amplifier and, in turn, the electrical environment of the cells, as described from the pCLAMP^TM^ 11 software suite or the HEKA system (https://www.heka.com/). Moreover, while strip chart recorders were once used to print experimental results, personal computers are now employed for analysis and output.
As a result, these developments have made it difficult for senior electrophysiologists to update or replace their instruments or devices in line with advances in technology and artificial intelligence. Meanwhile, younger scholars tend to focus more on computer operation, often neglecting the practical skills and techniques of patch-clamp recording. Additionally, automated patch-clamp system (such as QPatch^®^; Sophion Bioscience; Ballerup, Denmark, and the Port-a-Patch^®^ system; Nanion Technologies GmbH; Munich, Germany), as well as HEKA’s manual patch-clamp amplifier (HEKA Technologies GmbH), has led to significant breakthroughs in the study of various voltage-gated and ligand-gated ion currents. Although automated patch clamp systems efficiently replace manual methods for isolated cells, brain slice patch-clamp recordings remain a manual technique. Effective robotic platforms used for complex tissue slices awaits further technological developments [28,31,33,34,37,38,39,40]. For example, artificial intelligence (AI)-driven image recognition can identify suitable cells for recordings under the microscope. Machine learning algorithms can guide micropipette positioning with high precision, reducing human error. It is hoped that AI can help determine when to transition from C-A to W-C mode by analyzing resistance and capacitance changes automatically in the near future. Nevertheless, the scientific findings rather than technical background need to emphasized.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Hamill O.P. Marty A. Neher E. Sakmann B. Sigworth F.J. Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches Pflugers Arch.19813918510010.1007/BF 006569976270629 · doi ↗ · pubmed ↗
- 2Wilson J.R. Clark R.B. Banderali U. Giles W.R. Measurement of the membrane potential in small cells using patch clamp methods Channels 2011553053710.4161/chan.5.6.1748421829090 PMC 3265801 · doi ↗ · pubmed ↗
- 3Kusunoki M. Hirota K. Shoji T. Uba T. Matsuo Y. Hayashi M. Dexmedetomidine suppresses glucose-stimulated insulin secretion in pancreatic β-cells FEBS Open Bio 20251563464610.1002/2211-5463.13960 PMC 1196139139704036 · doi ↗ · pubmed ↗
- 4Parker M. Yadavalli N.S. Peranidze K. Boland E. Reukov V. Minko S. Insulin-producing INS-1 cell cultures on biomimetic 3D scaffolds J. Mater. Chem. B 202513121341214510.1039/D 5TB 00519 A 40899143 · doi ↗ · pubmed ↗
- 5O’Neill R.J. Tung L. Cell-attached patch clamp study of the electropermeabilization of amphibian cardiac cells Biophys. J.1991591028103910.1016/S 0006-3495(91)82318-91907865 PMC 1281338 · doi ↗ · pubmed ↗
- 6Haug T.M. Hafting T. Sand O. Inhibition of BK channels contributes to the second phase of the response to TRH in clonal rat anterior pituitary cells Acta Physiol. Scand.200418034735710.1111/j.1365-201X.2004.01266.x 15030376 · doi ↗ · pubmed ↗
- 7Mason M.J. Simpson A.K. Mahaut-Smith M.P. Robinson H.P. The interpretation of current-clamp recordings in the cell-attached patch-clamp configuration Biophys. J.20058873975010.1529/biophysj.104.04986615516522 PMC 1305049 · doi ↗ · pubmed ↗
- 8Wu S.N. Wu Y.H. Chen B.S. Lo Y.C. Liu Y.C. Underlying mechanism of actions of tefluthrin, a pyrethroid insecticide, on voltage-gated ion currents and on action currents in pituitary tumor (GH 3) cells and Gn RH-secreting (GT 1-7) neurons Toxicology 2009258707710.1016/j.tox.2009.01.00919378468 · doi ↗ · pubmed ↗
