Abstract
Biological networks are constructed of repeated simplified patterns, or modules, called network motifs. Network motifs can be found in a variety of organisms including bacteria, plants, and animals, as well as intracellular transcription networks for gene expression and signal transduction processes in neuronal circuits. Standard models of signal transduction events for synaptic plasticity and learning often fail to capture the complexity and cooperativity of the molecular interactions underlying these processes. Here, we apply network motifs to a model for signal transduction during an in vitro form of eyeblink classical conditioning that reveals an underlying organization of these molecular pathways. Experimental evidence suggests there are two stages of synaptic AMPA receptor (AMPAR) trafficking during conditioning. Synaptic incorporation of GluR1-containing AMPARs occurs early to activate silent synapses conveying the auditory conditioned stimulus and this initial step is followed by delivery of GluR4 subunits that supports acquisition of learned conditioned responses (CRs). Overall, the network design of the two stages of synaptic AMPAR delivery during conditioning describes a coherent feed-forward loop (C1-FFL) with AND logic. The combined inputs of GluR1 synaptic delivery AND the sustained activation of 3-phosphoinositide-dependent protein-kinase-1 (PDK-1) results in synaptic incorporation of GluR4-containing AMPARs and the gradual acquisition of CRs. The network architecture described here for conditioning is postulated to act generally as a sign-sensitive delay element that is consistent with the non-linearity of the conditioning process. Interestingly, this FFL structure also performs coincidence detection. A motif-based approach to modeling signal transduction can be used as a new tool for understanding molecular mechanisms underlying synaptic plasticity and learning and for comparing findings across forms of learning and model systems.
Introduction
Motor learning is fundamental to neural circuits controlling voluntary movements in order to appropriately plan and adapt them to environmental situations. In recent years a great deal of effort has gone into elucidating the cellular and molecular mechanisms that underlie learning processes in the brain. To understand the signal transduction events involved in these processes a relatively simple form of associative learning called classical conditioning has been extensively studied. In this form of learning, a predictive relationship is acquired between a neutral stimulus that ordinarily produces no behavior, like the ringing of a bell, and a noxious stimulus such as an airpuff to the eye that produces a defensive motor behavior, in this case an eyeblink. After only a few pairings of the bell and airpuff animals and humans learn to blink in response to the bell alone. The generation of motor responses to novel stimuli is necessary for survival and therefore this form of learning is rapid and robust. Voluntary movements are controlled by motor signals generated by widespread brain regions acting in parallel, including the motor cortex, basal ganglia, and cerebellum, and the recognition of such large-scale network activity in brain function has lead to the development of the distributed processing module (DPM) model of behavior and cognition (Houk, ). How the DPM model applies to classical conditioning through the activity of distributed brain networks that coordinate training input signals with adaptive motor output is outlined by Houk (). Using this as a foundation, we are constructing a multilevel model of classical conditioning starting here with the molecular events that underlie an in vitro form of eyeblink classical conditioning initially described several years ago (Keifer et al., ) and recently reviewed in Keifer and Zheng ().
Biological networks are constructed of repeated simplified patterns called network motifs (Alon, ). These building block networks have their own limited information processing capacity and combine to build larger more complex systems. Network motifs can be found in a variety of organisms including bacteria, yeast, plants, and animals, as well as intracellular transcription networks for gene expression and signal transduction processes in neuronal circuits of interest here. That biological systems follow general rules is not new. However, as applied to studies of learning, models for signal transduction processes are typically designed as sequential and parallel pathways that can be daunting in their complexity. The application of network motifs to such models reveals simplifying recurring patterns, or modules, in the design structure of these signal transduction pathways. Using these simplified principals, the structure of potential unknown signaling pathways may also be predicted. As a first step, in this report we apply common biological network motifs mainly using feed-forward loops (FFLs) to model signal transduction during in vitro eyeblink classical conditioning.
Signal Transduction Events in in vitro Classical Conditioning
A neural analog of eyeblink classical conditioning can be studied using an isolated brainstem preparation from turtles (Keifer and Zheng, ). This preparation, illustrated in Figure 1A, is unique because turtle brain tissue is highly resistant to hypoxia such that the entire brain and brainstem can be maintained for hours or days in a dish allowing for in vitro studies of learning processes. Considerable experimental advantage is achieved since neuronal circuits under study are accessible to recording electrodes and application of pharmacological compounds or small interfering RNAs (siRNAs) using incubation procedures. Instead of using a tone or airpuff as in behaving animals, paired stimulation of the auditory nerve (the “tone” conditioned stimulus, CS) with the trigeminal nerve (the “airpuff” unconditioned stimulus, US) generates neural discharge in the abducens nerve, which controls blinking in this species, that is characteristic of conditioned eyeblink responses (Figure 1B, arrow). Once CS–US pairing is initiated, conditioned responses (CRs) are acquired rapidly in about 1 h, or during the second pairing session, compared to controls receiving unpaired stimuli that show no CRs (Figure 1C). Each pairing session during the training consists of 50 paired CS–US stimuli (lasting 25 min in duration) followed by a 30-min rest period in which there is no stimulation. Our studies indicate that trafficking of AMPA type glutamate receptors (AMPARs) containing GluR1 and GluR4 subunits support conditioning in this preparation. Based on data describing the action of multiple signal transduction elements, a two-stage model of synaptic AMPAR delivery in the abducens motor neurons during acquisition of eyeblink conditioning was developed and is illustrated in Figure 1D (Zheng and Keifer, ; Keifer and Zheng, ). First, GluR1-containing AMPAR synaptic incorporation occurs early in conditioning to activate silent synapses containing only NMDA receptors (NMDARs) that convey the auditory CS (Figure 1D, Early). This step is initiated within 15 min of training with the phosphorylation of protein kinase A (PKA) and the calcium/calmodulin-dependent protein kinases (CaMK) II and IV (Zheng and Keifer, ). This leads to activation of the transcription factor cAMP response element-binding (CREB) protein which ultimately results in production of brain-derived neurotrophic factor (BDNF) that is required for AMPAR delivery and CR acquisition (Li and Keifer, , ). BDNF is hypothesized to activate extracellular signal-regulated kinase (ERK) by signaling through the BDNF receptor tropomyosin-related kinase B (TrkB) that induces synaptic delivery of GluR1-containing AMPARs by translocation of existing receptors into the synapse. The incorporation of GluR1 into auditory nerve synapses early in conditioning is an essential step that serves to activate NMDARs, thereby allowing post-synaptic intracellular calcium (Ca2+) entry that triggers the second stage of AMPAR trafficking for conditioning involving GluR4 subunits (Mokin et al., ). Therefore, the first stage of AMPAR delivery is NMDAR-independent while the second stage requires NMDAR activation. In the second step of AMPAR trafficking, the delivery of GluR1 is followed by replacement of those subunits with synaptic incorporation of newly synthesized GluR4-containing AMPARs that underlies the acquisition of CRs (Figure 1D, Late; Mokin et al., ; Zheng and Keifer, ; Keifer and Zheng, ). Pharmacological data indicate that this step requires the coordinated actions of Ca2+-dependent and independent isoforms of PKC, as well as CaMKII and ERK (Zheng and Keifer, , ).
Figure 1
The multistage physiologically based model shown in Figure 1D necessarily illustrates signal transduction events during conditioning as a series of sequential steps. However, it fails to convey the multiple interactions among the elements of these molecular pathways and the cooperativity in their actions. Here, as a way of improving our understanding of these events, we present a new motif-based model that is founded on generalized design principles derived from a broad selection of biological systems (Alon, ). Our initial goals are: (1) to construct a biologically realistic model of signal transduction for classical conditioning based on network motifs, and (2) to predict the potential structure of unknown signaling pathways to guide hypothesis construction for experimental examination. The motif-based model describes the two stages of synaptic AMPAR trafficking during conditioning in Parts A and B below. Potential interactions of these networks form a more generalized network motif pattern that is presented in Part C. The model helps to clarify potential interactions between signal transduction elements during acquisition of conditioning, timing of these synaptic processes, and guides hypotheses for as yet undiscovered signaling events that may occur during coincidence detection of the CS and US. The application of network motif design principals to signal transduction mechanisms shown here should have broad applications to studies of learning.
Analysis
Repeated simplified patterns, or network motifs, form the building blocks of information processing systems. Some common ones are shown in Figure 2 and are discussed extensively by Alon (). The FFL is a network pattern having three nodes, X, Y, and Z, in which X has a direct path to Z and an indirect path to Z through Y (Figure 2). Coherent FFLs are those in which the direct path has the same overall sign as the indirect path, while the indirect path in the incoherent FFL has the opposite sign as the direct one. The complexity of the FFL in particular can be enhanced by duplication of one of its nodes as illustrated in Figure 3. The basic FFL is shown in Figure 3 (left) as is the network structure resulting from a single duplication of X which forms a multi-input FFL (right), in this case a two-input FFL. Therefore, complexity of the information processing capacity of the network can be increased by simple replication. The FFL is a prominent network motif in transcriptional (Shen-Orr et al., ) and neuronal networks (White, ; Milo et al., ). It’s structure and function are reviewed by Mangan and Alon (). FFLs have also been proposed to support the function of large neuronal arrays in the DPM model of brain function (Houk et al., ; Houk, , ), for example, as anticipatory commands for sensorimotor control.
Figure 2
Figure 3

The basic feed-forward loop (FFL) has a direct pathway from X to Z and an indirect pathway to Z through Y (left). Complexity of the FFL is increased by a simple replication of one of its nodes, in this case X, which forms a multi-input FFL.
Part A. Multi-Input FFL with Diamond Motif for GLUR1 Delivery
The FFL network is prominent in signal transduction cascades engaged by classical conditioning. Figure 4 illustrates the network motif design for the early signaling events leading to synaptic GluR1 AMPAR incorporation during in vitro classical conditioning. In the initial stages of conditioning, PKA and CaMKs II and IV are phosphorylated within 15 min of CS–US pairing (Zheng and Keifer,
Figure 4

A multi-input FFL with diamond motif describes the initial stage of signal transduction for synaptic GluR1-containing AMPAR delivery during in vitro classical conditioning. This step results in synaptic incorporation of GluR1 subunits to unsilence auditory nerve synapses that convey the CS by activation of NMDARs. See text for details.
The design architecture of signal transduction events leading to synaptic GluR1 delivery during conditioning can be described as a multi-input, in this case two-input, FFL as shown in Figure 4 (see also Figure 3, right). Two inputs from PKA and the CaMKs (X) directly, and indirectly through a CREB–BDNF–ERK pathway (Y), regulate GluR1 trafficking (Z). Embedded in the Y pathway is a diamond network motif for CREB–BDNF regulation. The multi-input FFL and diamond are strong network motifs found in neuronal networks and signal transduction processes of some relatively simple systems (Alon,
Part B. Multi-Input FFL for GluR4 Delivery
The motif-based model for the known signal transduction events leading to synaptic delivery of GluR4-containing AMPARs in the second stage of conditioning is illustrated in Figure 5. Delivery of GluR4 AMPAR subunits is an NMDAR-dependent step. Synaptic incorporation of GluR1 AMPAR subunits into NMDAR-containing silent auditory nerve synapses in Part A (GluR1/NMDAR in Figure 5) unsilences those synapses allowing for post-synaptic Ca2+ entry into abducens motor neurons to activate Ca2+-sensitive conventional PKC (PKCc, namely α and β). Activation of PKC is required for GluR4-containing AMPAR trafficking and acquisition of CRs but pharmacological inhibition of PKC does not affect synaptic GluR1 delivery (Zheng and Keifer,
Figure 5

A multi-input FFL describes synaptic delivery of GluR4-containing AMPARs in the second stage of conditioning. AMPARs containing GluR4 replace GluR1 subunits and underlie the acquisition of CRs. See text for details.
The signal transduction events leading to synaptic incorporation of GluR4-containing AMPARs are not yet as well described as for GluR1 during conditioning in this system. Nevertheless, it is instructive to model the known interactions in terms of the multi-input FFL-like architecture shown in Figure 5. This model helps to clarify two previous experimental findings. First, application of either chelerythrine (an inhibitor of PKCc, i.e., PKCαβ but not ζ) or the PKCζ pseudosubstrate peptide inhibitor ZIP (an inhibitor of PKCζ and downstream ERK but not PKCαβ) suppresses GluR4 delivery and conditioning (Zheng and Keifer,
Part C. Network Interactions between the FFLs
Viewing these signal transduction pathways from a more global perspective, the FFL of Part A is hypothesized to directly interact with the FFL of Part B in a feed-forward manner in two distinct ways (Figure 6; blue arrows). First, the synaptic insertion of AMPARs containing GluR1 subunits into NMDAR-containing silent auditory nerve synapses activates them and results in post-synaptic intracellular Ca2+ entry required to initiate Part B (Figure 6; blue arrow in the middle). Second, PKA and CaMKII activated in Part A are likely to regulate GluR4 trafficking perhaps by direct phosphorylation of that subunit (Figure 6; outer two blue arrows). While this has not yet been directly demonstrated experimentally, this assertion is supported by the observation that once activated after onset of the conditioning stimuli both PKA and CaMKII are maintained in a phosphorylated state for hours allowing them to affect later AMPAR trafficking events (Zheng and Keifer,
Figure 6

Network interactions between the two stages of AMPAR trafficking in conditioning. The later stage of synaptic incorporation of GluR4 AMPAR subunits (Part B) depends on the initial signal transduction events leading to synaptic delivery of GluR1 (Part A) and activation of PDK-1. Approximate activation time after the onset of CS–US pairing is indicated. See text for details.
The overall structure of the signal transduction pathways for conditioning shown in Figure 6 describes a coherent type-1 FFL (C1-FFL) with an AND gate as discussed by Alon (
Figure 7

(A) General network structure of a C1-FFL with AND logic. (B) This network motif as applied to signal transduction during in vitro classical conditioning detailed in Figure 3 and described in the text.
General Discussion and Hypotheses
The C1-FFL motif in classical conditioning
Feed-forward loops have been applied to many biological systems to simplify and explain their function. The most common FFL in biological networks has been proposed to be the C1-FFL (Alon,
Experimental evidence for the function of network motifs in biological systems is derived mainly from organisms such as bacteria. Mangan et al. (
Coincidence detection
The two-input FFL motif embedded in the initial stages of conditioning shown in Figure 4 and discussed in Part A is particularly relevant to models of classical conditioning because it can perform coincidence detection. Two different inputs may have a cooperative effect if they occur close enough together in time. The dynamics of this integrative effect depend on the strength of the two inputs, the threshold of the output to be achieved (activation threshold), and the rate function of activation or decay. In classical conditioning, molecular events leading to coincidence detection are usually applied to the initial association of a CS with a coincident or briefly delayed US. In order to form a learned association of the two stimuli they must occur close enough together in time to have a cooperative effect. Many ideas for coincidence detection have been founded on Donald Hebb’s original hypothesis that synaptic modifications underlying learning occur in response to simultaneous pre- and post-synaptic activity (Hebb,
Motif-based models as a new tool for understanding signal transduction during synaptic plasticity and learning
Current efforts in depicting the complex signaling events that underlie different forms of synaptic plasticity and learning are undoubtedly helpful in conveying the basic physiological events underlying these processes (e.g., Derkach et al.,
Statements
Acknowledgments
Supported by NIH grants NS051187 and P20 RR015567 which is designated as a Center of Biomedical Research Excellence (COBRE) to Joyce Keifer, and P01 NS44383 to James C. Houk. We thank Dr. Clive Bramham for helpful discussions and review of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AlonU. (2007). An Introduction to Systems Biology. Boca Raton, FL: Chapman & Hall/CRC.
2
BramhamC. R.MessaoudiE. (2005). BDNF function in adult synaptic plasticity: the synaptic consolidation hypothesis. Prog. Neurobiol.76, 99–125.10.1016/j.pneurobio.2005.06.003
3
CarvalhoA. L.KameyamaK.HuganirR. L. (1999). Characterization of phosphorylation sites on the glutamate receptor 4 subunit of the AMPA receptors. J. Neurosci.19, 4748–4754.
4
CasamayorA.MorriceN. A.AlessiD. R. (1999). Phosphorylation of Ser-241 is essential for the activity of 3-phosphoinositide-dependent protein kinase-1: identification of five sites of phosphorylation in vivo. Biochem. J.342, 287–292.10.1042/0264-6021:3420287
5
DerkachV. A.OhM. C.GuireE. S.SoderlingT. R. (2007). Regulatory mechanisms of AMPA receptors in synaptic plasticity. Nat. Rev. Neurosci.8, 101–113.10.1038/nrm2101
6
DessauerC. W. (2009). Adenylyl cyclase-A-kinase anchoring protein complexes: the next dimension in cAMP signaling. Mol. Pharmacol.76, 935–941.10.1124/mol.109.059345
7
GervasiN.TchenioP.PreatT. (2010). PKA dynamics in a Drosophila learning center: coincidence detection by rutabaga adenylyl cyclase and spatial regulation by dunce phosphodiesterase. Neuron65, 516–529.10.1016/j.neuron.2010.01.014
8
HebbD. O. (1949). The Organization of Behavior. New York, NY: John Wiley.
9
HoukJ. C. (2005). Agents of the mind. Biol. Cybern.92, 427–437.10.1007/s00422-005-0569-8
10
HoukJ. C. (2010). “Voluntary movement: control, learning and memory,” in Encyclopedia of Behavioral Neuroscience, eds KoobG. F.Le MoalM.ThompsonR. F. (Oxford: Academic Press), 3, 455–458.
11
HoukJ. C.KeiferJ.BartoA. G. (1993). Distributed motor commands in the limb premotor network. Trends Neurosci.16, 27–33.10.1016/0166-2236(93)90049-R
12
KeiferJ.ArmstrongK. E.HoukJ. C. (1995). In vitro classical conditioning of abducens nerve discharge. J. Neurosci.15, 5036–5048.
13
KeiferJ.SabirzhanovB. E.ZhengZ.LiW.ClarkT. G. (2009). Cleavage of proBDNF to BDNF by a tolloid-like metalloproteinase is required for acquisition of in vitro eyeblink classical conditioning. J. Neurosci.29, 14956–14964.10.1523/JNEUROSCI.3649-09.2009
14
KeiferJ.ZhengZ. (2010). AMPA receptor trafficking and learning. Eur. J. Neurosci.32, 269–277.10.1111/j.1460-9568.2010.07339.x
15
LiW.KeiferJ. (2008). Coordinate action of pre- and postsynaptic brain-derived neurotrophic factor is required for AMPAR trafficking and acquisition of in vitro classical conditioning. Neuroscience155, 686–697.10.1016/j.neuroscience.2008.06.043
16
LiW.KeiferJ. (2009). BDNF-induced synaptic delivery of AMPAR subunits is differentially dependent on NMDA receptors and requires ERK. Neurobiol. Learn. Mem.91, 243–249.10.1016/j.nlm.2008.10.002
17
LorenzettiF. D.BaxterD. A.ByrneJ. H. (2008). Molecular mechanisms underlying a cellular analog of operant reward learning. Neuron59, 815–828.10.1016/j.neuron.2008.07.019
18
ManganS.AlonU. (2003). Structure and function of the feed-forward loop network motif. Proc. Natl. Acad. Sci. U.S.A.100, 11980–11985.10.1073/pnas.2133841100
19
ManganS.ZaslaverA.AlonU. (2003). The coherent feedforward loop serves as a sign-sensitive delay element in transcription networks. J. Mol. Biol.334, 197–204.10.1016/j.jmb.2003.09.049
20
MiguesP. V.HardtO.WuD. C.GamacheK.SacktorT. C.WangY. T.NaderK. (2010). PKMζ maintains memories by regulating GluR2-dependent AMPA receptor trafficking. Nat. Neurosci.13, 630–634.10.1038/nn.2531
21
MiloR.Shen-OrrS.ItzkovitzS.KashtanN.ChklovskiiD.AlonU. (2002). Network motifs: simple building blocks of complex networks. Science298, 824–827.10.1126/science.298.5594.824
22
MokinM.ZhengZ.KeiferJ. (2007). Conversion of silent synapses into the active pool by selective GluR1-3 and GluR4 AMPAR trafficking during in vitro classical conditioning. J. Neurophysiol.98, 1278–1286.10.1152/jn.00212.2007
23
MullasserilP.DosemeciA.LismanJ. E.GriffithL. C. (2007). A structural mechanism for maintaining the “on-state” of the CaMKII memory switch in the post-synaptic density. J. Neurochem.103, 357–364.
24
NewtonA. C. (2003). Regulation of the ABC kinases by phosphorylation: protein kinase C as a paradigm. Biochem. J.370, 361–371.10.1042/BJ20021626
25
PangP. T.TengH. K.ZaitsevE.WooN. T.SakataK.ZhenS.TengK. K.YungW.-H.HempsteadB. L.LuB. (2004). Cleavage of proBDNF by tPA/plasmin is essential for long-term hippocampal plasticity. Science306, 487–491.10.1126/science.1100135
26
PattersonS. L.PittengerC.MorozovA.MartinK. C.ScanlinH.DrakeC.KandelE. R. (2001). Some forms of cAMP-mediated long-lasting potentiation are associated with release of BDNF and nuclear translocation of phospho-MAP kinase. Neuron32, 123–140.10.1016/S0896-6273(01)00443-3
27
SabirzhanovB. E.KeiferJ.ClarkT. G. (2007). Characterization of a novel reptilian tolloid-like gene in the pond turtle, Pseudemys scripta elegans. Brain Res.1154, 22–30.10.1016/j.brainres.2007.04.003
28
SatoS.FujitaN.TsuruoT. (2004). Involvement of 3-phosphoinositide-dependent protein kinase-1 in the MEK/MAPK signal transduction pathway. J. Biol. Chem.279, 33759–33767.10.1074/jbc.M314332200
29
SharmaS. K.SherffC. M.StoughS.HsuanV.CarewT. J. (2006). A tropomyosin-related kinase B ligand is required for ERK activation, long-term synaptic facilitation, and long-term memory in Aplysia. Proc. Natl. Acad. Sci. U.S.A.103, 14206–14210.10.1073/pnas.0603412103
30
Shen-OrrS.MiloR.ManganS.AlonU. (2002). Network motifs in the transcriptional regulation network of Escherichia coli. Nat. Genet.31, 64–68.10.1038/ng881
31
ThomasG. M.LinD.-T.NuriyaM.HuganirR. L. (2008). Rapid and bi-directional regulation of AMPA receptor phosphorylation and trafficking by JNK. EMBO J.27, 361–372.10.1038/sj.emboj.7601969
32
TsienJ. Z. (2000). Linking Hebb’s coincidence-detection to memory formation. Curr. Opin. Neurobiol.10, 266–273.10.1016/S0959-4388(00)00070-2
33
WhiteJ. G. (1985). Neuronal connectivity in Caenorhabditis elegans. Trends Neurosci.8, 277–283.10.1016/0166-2236(85)90124-9
34
YangJ.SiaoC.-J.NagappanG.MarinicT.JingD.McGrathK.ChenZ.-Y.MarkW.TessarolloL.LeeF. S.LuB.HempsteadB. L. (2009). Neuronal release of proBDNF. Nat. Neurosci.12, 113–115.10.1038/nn.2244
35
ZhengZ.KeiferJ. (2008). Protein kinase C-dependent and independent signaling pathways regulate synaptic GluR1 and GluR4 AMPAR subunits during in vitro classical conditioning. Neuroscience156, 872–884.10.1016/j.neuroscience.2008.08.042
36
ZhengZ.KeiferJ. (2009). PKA has a critical role in synaptic delivery of GluR1- and GluR4-containing AMPARs during initial stages of acquisition of in vitro classical conditioning. J. Neurophysiol.101, 2539–2549.10.1152/jn.91282.2008
37
ZhuJ. J.QinY.ZhaoM.van AelstL.MalinowR. (2002). Ras and Rap control AMPA receptor trafficking during synaptic plasticity. Cell110, 443–455.10.1016/S0092-8674(02)00897-8
Summary
Keywords
classical conditioning, AMPA receptor trafficking, network motifs, model, signal transduction, eyeblink, in vitro, feed-forward loops
Citation
Keifer J and Houk JC (2011) Modeling Signal Transduction in Classical Conditioning with Network Motifs. Front. Mol. Neurosci. 4:9. doi: 10.3389/fnmol.2011.00009
Received
22 December 2010
Accepted
22 June 2011
Published
07 July 2011
Volume
4 - 2011
Edited by
Alistair N. Garratt, Max Delbrück Center for Molecular Medicine, Germany
Reviewed by
Verena Tretter, Medical University Vienna, Austria; Alistair N. Garratt, Max Delbrück Center for Molecular Medicine, Germany
Copyright
© 2011 Keifer and Houk.
This is an open-access article subject to a non-exclusive license between the authors and Frontiers Media SA, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and other Frontiers conditions are complied with.
*Correspondence: Joyce Keifer, Neuroscience Group, Division of Basic Biomedical Sciences, University of South Dakota School of Medicine, 414 East Clark Street, Vermillion, SD 57069, USA. e-mail: jkeifer@usd.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.