How Can Transcranial Magnetic Stimulation Be Used to Modulate Episodic Memory?: A Systematic Review and Meta-Analysis

A systematic review and meta-analysis were conducted to synthesize the existing literature on how transcranial magnetic stimulation (TMS) has been used to modulate episodic memory. Given the numerous parameters of TMS protocols and experimental design characteristics that can be manipulated, a mechanistic understanding of how changes in the combination of parameters (e.g., frequency, timing, intensity, targeted brain region, memory task) modulate episodic memory is needed. To address this, we reviewed 59 studies and conducted a meta-analysis on 245 effect sizes from 37 articles on healthy younger adults (N = 1,061). Analyses revealed generally more beneficial effects of 1-Hz rTMS vs. other frequencies on episodic memory. Moderation analyses revealed complex interactions as online 20-Hz rTMS protocols led to negative effects, while offline 20-Hz rTMS led to enhancing effects. There was also an interaction between stimulation intensity and frequency as 20-Hz rTMS had more negative effects when applied below- vs. at-motor threshold. Conversely, 1-Hz rTMS had more beneficial effects than other frequencies when applied below- vs. at- or above-motor threshold. No reliable aggregate or hypothesized interactions were found when assessing stimulation site (frontal vs. parietal cortex, left vs. right hemisphere), stimulated memory process (during encoding vs. retrieval), the type of retrieval (associative/recollection vs. item/familiarity), or the type of control comparison (active vs. sham or no TMS) on episodic memory. However, there is insufficient data to make strong inference based on the lack of aggregate or two-way interactions between these factors, or to assess more complex (e.g., 3-way) interactions. We reviewed the effects on other populations (healthy older adults and clinical populations), but systematic comparison of parameters was also prevented due to insufficient data. A database of parameters and effects sizes is available as an open source repository so that data from studies can be continuously accumulated in order to facilitate future meta-analysis. In conclusion, modulating episodic memory relies on complex interactions among the numerous moderator variables that can be manipulated. Therefore, rigorous, systematic comparisons need to be further investigated as the body of literature grows in order to fully understand the combination of parameters that lead to enhancing, detrimental or null effects on episodic memory.

A systematic review and meta-analysis were conducted to synthesize the existing literature on how transcranial magnetic stimulation (TMS) has been used to modulate episodic memory. Given the numerous parameters of TMS protocols and experimental design characteristics that can be manipulated, a mechanistic understanding of how changes in the combination of parameters (e.g., frequency, timing, intensity, targeted brain region, memory task) modulate episodic memory is needed. To address this, we reviewed 59 studies and conducted a meta-analysis on 245 effect sizes from 37 articles on healthy younger adults (N = 1,061). Analyses revealed generally more beneficial effects of 1-Hz rTMS vs. other frequencies on episodic memory. Moderation analyses revealed complex interactions as online 20-Hz rTMS protocols led to negative effects, while offline 20-Hz rTMS led to enhancing effects. There was also an interaction between stimulation intensity and frequency as 20-Hz rTMS had more negative effects when applied below-vs. at-motor threshold. Conversely, 1-Hz rTMS had more beneficial effects than other frequencies when applied below-vs. at-or above-motor threshold. No reliable aggregate or hypothesized interactions were found when assessing stimulation site (frontal vs. parietal cortex, left vs. right hemisphere), stimulated memory process (during encoding vs. retrieval), the type of retrieval (associative/recollection vs. item/familiarity), or the type of control comparison (active vs. sham or no TMS) on episodic memory. However, there is insufficient data to make strong inference based on the lack of aggregate or two-way interactions between these factors, or to assess more complex (e.g., 3-way) interactions. We reviewed the effects on other populations (healthy older adults and clinical populations), but systematic comparison of parameters was also prevented due to insufficient data. A database of parameters and effects sizes is available as an open source repository so that data from studies can be continuously accumulated in order to facilitate future meta-analysis. In conclusion, modulating episodic memory relies on complex interactions among the numerous moderator variables that can be manipulated. Therefore, rigorous, systematic comparisons need to be further investigated as the body of literature grows in order to fully understand the combination of parameters that lead to enhancing, detrimental or null effects on episodic memory.
Over the past 30 years, transcranial magnetic stimulation (TMS) has enabled researchers to move beyond correlational research and address causal relations between brain and behavior in humans. While powerful, the mechanisms of TMS and the resulting outcomes on higher order cognitive functioning such as episodic memory are complex and mired by differences across studies. As the body of research on TMS continues to emerge, rigorous, and systematic comparisons need to be investigated to better understand the moderating effects of the various parameters of TMS that can be manipulated. To this aim, we first reviewed studies that have used repetitive TMS (rTMS) in attempt to modulate episodic memory. Next, we conducted a meta-analysis to bridge the gap in understanding of aspects of rTMS protocols that lead to memory enhancement or impairment. Lastly, we reviewed and compared the limited number of studies that have used rTMS to modulate episodic memory in healthy older adults and those who suffer from clinical disorders.
In a systematic review, we focus on examining the effects of rTMS on episodic memory, which entails vividly remembering past events along with their spatial, temporal, and/or source details (Tulving, 2002). A full understanding of the mechanisms underlying cognitive processes in episodic memory is critical for developing effective interventions to target memory deficits that occur in healthy aging and in those with clinical disorders. However, the neural substrates of episodic memory are difficult to directly observe and manipulate and the field has generally relied on neuroimaging techniques that have provided important, albeit correlational, findings linking neural activity in specific brain regions with episodic memory processes. The introduction of TMS techniques has made it possible to make more causal claims about the role of specific brain regions (and functionally coupled regions) in specific processes by temporarily and reversibly activating neuronal firing in targeted cortical regions. Such exogenous stimulation can induce the activation of specific neural circuits and/or perturb or disrupt endogenous neural activity in targeted and functionally connected areas, and then allows for investigation of the resulting effects on behavior, as well as on associated neural signals when TMS is coupled with simultaneous neuroimaging. Therefore, it is critical to both understand the ways in which varying rTMS parameters can reveal the role of neurocognitive processes in supporting episodic memory, and how rTMS may be used to enhance or impair episodic memory functioning (for early reviews, see Grafman and Wassermann, 1998;Manenti et al., 2012). Until now, there has been no systematic review or meta-analysis on the effects of rTMS on episodic memory.

PRINCIPLES OF rTMS AND EPISODIC MEMORY THAT ARE IMPORTANT FOR RESEARCHERS TO UNDERSTAND
When conducted properly, TMS is a relatively safe and painless technique that utilizes electromagnetism to induce current in the brain that can cause neurons to fire (Rossi et al., 2009). During TMS, an electrical current passes through a coil of insulated wiring and produces brief magnetic "pulses" to induce neural firing. These transitory currents cause rapid depolarization (i.e., action potentials) of neurons beneath the coil and can modulate cortical excitability with a relatively high level of spatial and temporal precision (Wagner et al., 2009). Importantly, through polysynaptic connections, action potentials can propagate to distal brain regions that are functionally connected at the time of stimulation.
The application of TMS protocols can be roughly classified with respect to the amount of time between pulses. For example, single and paired pulse procedures vs. repetitive TMS reflect differences in the time between delivering pulses. rTMS has been implemented to modulate cortical excitability with the intent to either enhance or attenuate subsequent cognitive and behavioral outcomes (Pascual-Leone et al., 1998). Stimulation parameters such as frequency, intensity, and stimulation timing have been shown to be critical in achieving the desired response. Multimodal techniques (i.e., TMS-EEG, TMS-fMRI) and machine learning analyses (i.e., multivariate pattern analysis) have revealed the importance of understanding the relationship between the anatomical and functional brain regions being targeted and the cognitive representations and processes that are coded in the neural signals and may be manipulated by stimulation. In other words, stimulation parameters interact with ongoing patterns of neural activity to drive changes in cognition and behavior (Silvanto et al., 2008;Romei et al., 2016). TMS techniques have grown in popularity due the feasibility of non-invasively modulating brain activity to study the role of specific brain regions in cognitive processes and modify behavior in awake, behaving humans. The techniques have far reaching implications within the field of memory for the development of treatments for people suffering from memory deficits due to normal aging or neurodegenerative disorders. As reviewed below, findings are often quite mixed and variable between individuals; thus, increased understanding of the most proficient rTMS protocols for effectively enhancing episodic memory is critical in order to achieve translational applications in clinical settings.

Episodic Memory
Memory is not an isomorphic construct and, of the various forms of memory, episodic memory is the form that is absolutely critical for maintaining one's self-identity. According to Tulving (2002), episodic memory is the type of declarative memory that represents knowledge of particular events that one has experienced as well as their contextual acquisition (e.g., autobiographical memory), such as when and where an event occurred (e.g., who you met at a dinner party last week). It is measured by retrieval accuracy on declarative memory tests such as recall and recognition. Episodic memory requires and extends a related type of declarative memory, i.e., semantic memory or one's general knowledge of the world. Semantic memory involves retrieving knowledge (i.e., facts, concepts, meaning) without the accompanying contextual information (e.g., knowing how many planets are in the solar system, but without remembering in vivid detail when or where this was learned). The critical role of the hippocampus and surrounding structures of the medial temporal lobes (MTL) in the encoding, consolidation, and retrieval of episodic memories has been highlighted by striking cases of amnesia following either surgical resection (as in patient HM), disease (as in Clive Wearing), or brain trauma (as in patient KC), and have been elucidated by decades of animal/behavioral, neuropsychological, neuroimaging, and neurostimulation research (for reviews see, Moscovitch et al., 2016;Eichenbaum, 2017).
The hippocampus alone does not support episodic memory, however. In order for an event to be remembered, it must first be encoded through some combination of perceptual (e.g., visual, acoustic), affective, contextual, and semantic processes that integrate an event with its contextual associations to form episodic memories. Additionally, executive control processes are thought to facilitate the multi-modal binding of features and associations for episodic memory representations and are thought to be mediated by regions of the prefrontal cortex (PFC) and parietal cortex (PC; Buckner et al., 1999). Encoding and retrieval processes are inherently intertwined. Once encoded, successful retrieval demonstrates access to a memory, which also serves as a measure for successful encoding. For example, the encoding specificity principle (Tulving and Thomson, 1973) and related concept of "transfer-appropriate-processing" (Morris et al., 1977;Blaxton, 1989) emphasize the role of a reinstatement of the encoding context during retrieval for episodic memory. At the neural level, context reinstatement is operationalized by the degree of overlap in neural activity elicited during encoding and retrieval, and functional neuroimaging studies have revealed the involvement of regions in the prefrontal (Prince, 2005), medial temporal (Rugg and Vilberg, 2013), and parietal cortices (Wagner et al., 2005) (for review, see Craik and Rose, 2012).

The Role of Encoding and Retrieval
Although encoding and retrieval are intertwined, there are important differences in the neural correlates of successful episodic encoding and retrieval. For example, there is often a hemispheric asymmetry in left vs. right PFC activation during encoding vs. retrieval, depending on the task stimuli, familiarity, and the cognitive processes engaged (see below). With regards to rTMS, it is important to note that the neural substrates of episodic memory are impacted by the underlying state of activation in a stimulated brain region that is driven by the task demands and stimuli (e.g., words vs. pictures; Miniussi et al., 2013;Romei et al., 2016).

Verbal Stimuli
It is generally theorized that, in younger adults, the left PFC is critical for encoding, while the right PFC is to be preferentially involved in retrieval (Habib et al., 2003). This is known as the hemispheric encoding-retrieval asymmetry (HERA) model, which has been supported by some evidence from rTMS findings . However, the role of the left and right PFC in episodic memory may also vary due to numerous factors, such as task stimuli. For example, in two studies, left DLPFC stimulation impaired memory for verbal information (Rami et al., 2003), while right DLPFC stimulation impaired memory for non-verbal stimuli (Epstein et al., 2002;Floel et al., 2004). However, differences between the parameters of stimulation such as frequency and intensity across studies complicate direct comparison.

Stimulus Familiarity
The effects of rTMS on episodic memory has also been shown to vary with an individual's familiarity with the type of information that is to be remembered. For example, online 20 Hz rTMS to the left DLPFC during retrieval impaired novel compared to familiar information (Sandrini et al., 2003). Analogously, rTMS has been shown to preferentially affect memory for contextual information of studied associations. For example, online 10 Hz rTMS to the left inferior frontal gyrus (IFG) during encoding of novel faces paired with a "context" stimulus (e.g., the word "lawyer") impaired subsequent recognition of faces compared to faces studied without a context stimulus (i.e., faces studied with "no context"; Feurra et al., 2010). Thus, effects of rTMS on episodic memory likely depends not only on the specific rTMS parameters (e.g., frequency, intensity, stimulation timing, or target site), but also on the cognitive processes that are tapped by the type of information presented or the way in which it is to be remembered.

Cognitive Processes Engaged
In a similar vein, the effects of rTMS on episodic memory likely also depends on the cognitive strategies involved during encoding and/or retrieval (Manenti et al., 2010bHawco et al., 2013). For example, Innocenti et al. (2010) had individuals engage in "deep" (semantic) or "shallow" (perceptual) levels of processing while encoding words and applying online 10 Hz rTMS to either the left or right DLPFC. There were dissociable effects on recognition of words that were deeply or shallowly encoded, with left DLPFC stimulation abolishing the typical benefit of deep processing on memory (Innocenti et al., 2010). Relatedly, in a face-name pair encoding task, participants received fMRI-guided 10 Hz rTMS to either the left or right DLPFC during retrieval. After subdividing performance for those who reported using a retrieval strategy from those who reported using no retrieval strategy, it was concluded that right DLPFC stimulation disrupted retrieval for strategy users, while left DLPFC stimulation disrupted retrieval for non-strategy users (Manenti et al., 2010b).
Taken together, the effects of rTMS on episodic memory likely depend not only on the specific parameters of the protocol, but also on the brain state and cognitive processes engaged at the time of stimulation. This brain-state dependence is consistent with research that has used TMS in a wide variety of domains that has shown that the effects of TMS to a targeted region interacts with the ongoing, endogenous neural activity in that region and functionally coupled regions at the time of stimulation. Thus, while patterns are complicated across conditions, this principle enables researchers to investigate the causal role of the targeted and functionally coupled regions in specific neurocognitive mechanisms through the systematic manipulation of different task conditions (for review see, Romei et al., 2016).

The Role of Retrieval Processes
According to dual-process models of episodic memory, the cognitive processes engaged during retrieval can be supported by two separate processes (i.e., familiarity and recollection). A central theme of dual-process models is that previously studied information supports recognition by being generally more "familiar" than new/non-studied information. A distinct retrieval process, "recollection, " can further facilitate recognition of previously studied information through the retrieval of associated spatial, temporal, or contextual information (Yonelinas, 2002). The remember/know procedure is commonly implemented to examine differences in the nature or subjective experience of retrieval. Specifically, during a recognition task, participants indicate if the decision to respond that a presented item was previously seen (i.e., old) was based on recollection of the encoding details or based on a feeling of familiarity; they are instructed to respond with "remember" or "know, " respectively, (Tulving, 1985). rTMS enables the ability to causally investigate the underlying neural substrates of familiarity and recollection processes. For example, Turriziani et al. (2008) applied online 20 Hz rTMS to the left or right DLPFC during either encoding or retrieval (i.e., remember/know paradigm) of faces. They found that recollection (remember responses) was impaired following right compared to left DLPFC stimulation during encoding. However, familiarity was impaired following left and right DLPFC stimulation during encoding. These findings suggest that both the left and right DLPFC support familiarity and recollective retrieval processes, but the left vs. right DLPFC may contribute differentially depending on the type of retrieval process engaged.

The Role of Cortical Regions
The traditional view of episodic memory has focused on the hippocampus and its role in binding specific details of an experience (such as what events have occurred) with when (temporal), where (spatial), or how (source) the events occurred (Tulving, 1985(Tulving, , 2002. A large body of work on patient H.M., who had extensive surgical resection of the bilateral hippocampi and surrounding structures of the MTL in young adulthood, revealed the importance of the hippocampus in the formation of new episodic memories (Squire, 2009). However, inferences from brain lesion studies are often confounded by compensatory reorganization mechanisms (e.g., plasticity) that occur following brain injury and/or damage that is not isolated to specific brain regions. Alternatively, TMS can be used to causally modulate neural activity in healthy brain regions and networks at the precise time that a hypothesized cognitive operation is thought to occur; it affords a direct test of how specific brain regions contribute to episodic memory. However, the limited spatial ability to target deep regions in the MTL has led to substantial focus on targeting neocortical regions. The combination of TMS with neuroimaging (e.g., fMRI, EEG) and functional connectivity analyses provides a promising route to understanding brain network interactions, and has enabled further understanding of how TMS related changes in local neural activity propagate through polysynaptic connections to distal regions that are functionally connected at the time of stimulation (Shafi et al., 2012). For example, several studies have shown effects on episodic memory from targeting nodes in neocortical networks in prefrontal cortex or parietal cortex that are functionally connected with the hippocampus (e.g., Wang et al., 2014).

Prefrontal Cortex
The PFC in particular has received substantial attention as a stimulation site in rTMS protocols. For example, a fMRIlocalized rTMS study targeted the left inferior prefrontal cortex (LIPFC) online during a word-pair encoding task with 7 Hz rTMS. In comparison to two active control conditions (i.e., RIPFC, left superior partial cortex; LPC), stimulation of the LIPFC enhanced subsequent memory on an immediate word recognition test. The authors suggested that this causal link between LIPFC and episodic memory may have been driven by LIPFC stimulation triggering elaborative processing of word pairs that led to more distinctive memory cues (Köhler et al., 2004). One important design aspect of this study that we will briefly return to later is the fMRI-localization targeting procedure that accounts for individual differences in the structure and function of stimulated brain regions for more accurate targeting of regions of interest in a specific task.

Parietal Cortex
While PFC-MTL interactions are a cornerstone of neurocognitive models of episodic memory, it is also clear that there is extensive, often task-dependent recruitment of distributed networks of areas that process sensory, semantic, or emotional information to support various aspects of episodic memory (Ciaramelli et al., 2008;Cabeza and Moscovitch, 2013). Thus, research has examined the posterior parietal cortex due to its hypothesized role in attending to memory representations. It has been suggested that regions of parietal cortex are critical for episodic retrieval because of their role in either the formation of memory representations or shifts in attention for decision making processes, or both (Dobbins et al., 2003).
For example, Manenti et al. (2010a) implemented an rTMS-fMRI design to investigate the causal functioning of both the parietal cortex and DLPFC during encoding and retrieval. The authors applied online 10 Hz rTMS to either the left or right DLPFC or individualized regions of the left or right parietal cortex (either supramarginal or angular gyrus) during either encoding or retrieval. The authors demonstrated word recognition decrements (i.e., slower reaction times) following rTMS to the left parietal cortex (compared to sham). Furthermore, those with larger rTMS-related memory decrements were those who demonstrated a larger degree of right parietal activation during retrieval in the initial fMRI-localization task (Manenti et al., 2010a). One implication from these findings is that different aspects of episodic memory performance (e.g., accuracy, reaction times, subjective confidence) may be supported differentially by the PFC vs. PC.
Additionally, the null memory results reported (Manenti et al., 2010a) contrast with previously reported findings of Innocenti et al. (2010), despite the use of many of the same parameters (i.e., 10 Hz, left and right DLPFC, encoding, and retrieval of words). One possible explanation may be attributable to differences in other parameters, such as stimulation intensity (i.e., 100 vs. 90%). Similarly, the findings of Manenti et al. (2010a) also differed from Köhler et al. (2004) despite implementing many of the same parameters (i.e., online stimulation at the same intensity to the PFC and PC during word encoding). However, differences in the stimulation frequency (i.e., 7 vs. 10 Hz) might account for the different results.
It should be noted that studies often select a stimulation site to have a selective (or, more likely, stronger) effect on one cognitive domain, but it must be acknowledged that stimulated brain regions may be implicated in numerous related cognitive processes. For example, a meta-analysis revealed that rTMS to the DLPFC improved working memory performance (e.g., accuracy, reaction times; Brunoni and Vanderhasselt, 2014). The studies included in the Brunoni and Vanderhasselt (2014) metaanalysis targeted the DLPFC because of its role in short-term (STM)/working memory (WM). Thus, subsequent assessments of the rTMS effects on long-term episodic memory should be interpreted based on possible modulation of STM/WM processes on encoding and maintenance processes that may subsequently result in modulation of long-term memory. Indeed, a number of studies have examined the effects of rTMS to the DLPFC during STM/WM and has assessed the effects of rTMS on subsequent long-term memory (Nilakantan et al., 2017;Marin et al., 2018). For example, Marin et al. (2018) found that continuous theta burst (cTBS) to the right DLPFC prior to the encoding of objects and their spatial location was impaired for subsequent longterm memory tasks (i.e., object recognition and spatial location recall), but short-term memory performance (i.e., short-term spatial location recall) was left unaffected compared to controls.
Taken together, the findings reviewed thus far build upon previous correlational neuroimaging findings that the PFC, PC, and MTL are important for episodic memory, but systematic examination of encoding, consolidation, and retrieval processes, as well as potential interactions among neocortical regions, is needed. For example, the functional roles of these regions in recollection and familiarity (Yonelinas, 2005) remains to be fully understood.

The Role of Stimulation Frequency, Intensity, and Timing
How rTMS can enhance or impair memory (or leave it unaffected) depends on numerous factors regarding the ways in which the rTMS protocols are implemented. A detailed review of all technical considerations regarding rTMS is beyond the scope of the current paper (for details, see Kammer et al., 2001;Rossi et al., 2009;Sauvé and Crowther, 2014). The main stimulation characteristics considered here are stimulation frequency, intensity, and timing. A common approach in selecting a stimulation frequency generally involves applying a protocol hypothesized to either enhance cortical excitability, such as a high frequency (>5 Hz) or intermittent theta burst (iTBS) protocol, or suppress cortical excitability, such as a low frequency (<1 Hz) or cTBS protocol (Huang et al., 2005;Hallett, 2007). However, these hypotheses stem largely from physiological findings that have been investigated in the motor cortex with the assumption that this remains true for other cortical regions.
Another way in which different stimulation frequencies may modulate neural activity is through altering neural oscillatory activity through the resetting of natural oscillations, and in some cases, driving neural entrainment that may influence communication between brain regions and cognitive performance (Luber and Lisanby, 2014). Importantly, this modulation of neural activity is assumed to be specific to functionally relevant frequencies. Because relevant oscillations may vary in a continuous fashion between tasks, across phases of a task, and/or between individuals, considerable research is needed to fully understand how rTMS frequency can modulate episodic memory processes. A major aim of the current meta-analysis is to help elucidate the aggregate and complex interactions associated with the effects of different rTMS frequencies on episodic memory.
One of the many factors that will influence the effectiveness of rTMS is stimulation intensity. Generally, researchers will adjust stimulation intensity based upon an individual's resting or active motor threshold (RMT, AMT; Ngomo et al., 2012), which is defined as the minimal stimulator intensity needed to induce a muscle evoked potential in 5 out of 10 occurrences in either a relaxed (RMT) or actively engaged muscle (AMT). Prior work has found evidence that low (below) and high (above) stimulation intensities mirror facilitatory vs. inhibitory neural activity that is also associated with enhancements vs. decrements in behavioral outcomes (Silvanto et al., 2018). In addition, recent findings have shown that the effects of stimulation intensity on cortical excitability also may interact with specific frequencies, which in turn can lead to divergent effects on cognitive performance (Chung et al., 2018). Furthermore, the effects of varying stimulation intensity may covary with numerous factors such as the "brain state" or the cortical excitability of targeted and functionally coupled regions at the time of stimulation (Romei et al., 2016). Thus, a systematic investigation of the effects of stimulation intensity and how it interacts with frequency is needed. Elucidating the effects of varying stimulation intensity on episodic memory was another aim of the current study.
When stimulation is applied (i.e., stimulation timing parameter) is also an important factor that enables researchers to investigate the time course of episodic memory processes. The timing of the stimulation protocol reflects differences in whether stimulation occurs simultaneously with a cognitive process (online; e.g., rTMS during picture encoding) or prior to or following a cognitive process (offline; e.g., rTMS prior to picture encoding). One reason that online and offline stimulation may affect episodic memory differently is that they may affect different mechanisms, with online stimulation having an immediate influence on brain activity while offline stimulation often results in after-effects up to an hour after stimulation (Bergmann et al., 2016). In other words, online approaches can interfere or enhance neural activity through immediate depolarization of targeted neurons, while offline approaches can modulate long-term potentiation-or long-term depression-like plasticity mechanisms to either enhance or inhibit episodic memory processes. In a previous meta-analysis of different forms of non-invasive brain stimulation on various cognitive measures, differences in stimulation timing were found, with larger enhancing effects occurring for online stimulation (Hsu et al., 2015). Importantly, stimulation timing may interact with stimulation intensity, which could lead to facilitating or impairing effects on cognition (Silvanto and Cattaneo, 2017).

The Importance of Proper Control Conditions
The implementation of a proper control conditions is important to note given the variety of different control comparisons that are available. A proper control comparison is critical in order to assess whether the observed effects on performance are due specifically to the effects of applying TMS to the targeted region as opposed to some other variable, such as non-specific effects associated with the sensations of stimulation (Sandrini et al., 2011). Active stimulation control conditions that involve stimulating some other brain region that is not hypothesized to be as involved in the specific process of interest as the experimentally targeted region at the same intensity and frequency in both conditions are ideal because they most closely replicate the same aspects of the experimental/treatment TMS procedures. Sham stimulation controls (e.g., sham coil or tilting the coil away from the scalp) are implemented to mimic the non-specific effects of TMS (i.e., the discharging sound and scalp contact of the coil). "No TMS" conditions provide the least amount of control for non-specific effects of rTMS. In order to test if the type of control condition moderated the size of rTMS effects on episodic memory we examined whether there were systematic differences between effect sizes associated with comparisons to active, active vertex, sham, or no TMS control conditions.

THE PRESENT STUDY
A clear understanding of when episodic memory is enhanced, impaired, or left unaffected should depend on a wide variety of factors regarding the ways in which the rTMS protocols are implemented, and a subset of these factors can be roughly classified into specific rTMS parameters (e.g., frequency, intensity, timing, stimulation site, and control comparisons) and experimental design characteristics (e.g., memory processes: encoding vs. consolidation vs. retrieval, and retrieval type: familiarity vs. recollection). Contrasting the effects of different parameters allows researchers to test hypotheses about the role of specific neurocognitive processes that support different aspects of episodic memory. For example, rTMS could be used to influence many different processes associated with the acquisition of an event (encoding), its storage and preservation over time (consolidation processes), and/or the ability to remember the event (retrieval processes). Collectively, these different stages may engage both, similar, or distinct brain regions and, thus, the targeted brain regions and the nature of stimulation should result in variable effects. With this in mind, we conducted our survey of the literature to examine how rTMS has been shown to modulate episodic memory with a specific focus on how, when, and where rTMS was applied, as well as how episodic memory was measured.
A significant hurdle with understanding how rTMS impacts episodic memory is that most studies differed in a number of ways, such as in stimulation parameters (e.g., frequency, intensity, timing, and stimulation site), or with respect to how episodic memory was measured. Such differences make it difficult to pinpoint the specific protocols that will enhance or impair episodic memory. In the current study, we conducted a comprehensive systematic review of the literature with the aim of identifying the aggregate effects of rTMS on episodic memory and interactions among various factors. The focus of the meta-analysis was on healthy younger adults due to the small number of older adult rTMS studies. To the extent of our knowledge, this is the first meta-analysis to investigate this relationship. Such a meta-analysis is important for drawing consistent conclusions regarding the circumstances that lead to rTMS impacting episodic memory.
In order to carry this out, our primary aim was to investigate the effects of potential moderators of effect sizes associated with specific rTMS parameters or experimental design characteristics. Specific rTMS parameters included: stimulation frequency (1, 5, 10, 20, iTBS, cTBS), intensity (above, at, and below 100% MT), timing (online vs. offline), targeted hemisphere (left vs. right), targeted cortical region (frontal cortex vs. parietal cortex), and control conditions (active, active vertex, sham, and no TMS). Experimental design moderators included: memory processes stimulated (encoding vs. retrieval) and retrieval type (recollection vs. familiarity). Relatedly, we attempted to elucidate the complex relationship between rTMS and episodic memory by examining two-way interactions between factors that are hypothesized to interact. This allowed us to begin to bridge the gap in the understanding of the circumstances in which rTMS will lead to memory enhancement, impairment, or leave memory unaffected. Although it is not yet possible to conduct a similar meta-analysis of the effects of rTMS on episodic memory in older adults and some clinical conditions, we provide a brief narrative review of the available studies for comparison.

METHOD
We pursued all possible pertinent articles that reported data (e.g., group means, mean differences, tables, and bar graphs with error information) regarding the relationship between rTMS and episodic memory performance, and coded relevant parameters for potential moderator variables. First, we detail how the literature search was conducted for the systematic review, including the inclusion criteria that were implemented and selection protocol. Next, we discuss the subset used for the meta-analysis, which will focus on rTMS and episodic memory in healthy younger adults given that they account for a majority of the findings reported in the literature. Details for the coding system and how effect sizes were calculated and implemented in a three-level random effects meta-analysis are reported below.

Literature Search
An exhaustive PubMed search was conducted to select studies for the review and meta-analysis with a combination of the FIGURE 1 | A total of 59 articles were included in the systematic review. The meta-analysis focused on healthy younger adults and rTMS in episodic memory resulting in 37 articles with 245 effect sizes.
following search terms: "transcranial magnetic stimulation, " "TMS, " "rTMS, " "Theta burst stimulation, " "TBS, " "memory, " "episodic memory, " "long term memory, " "associative memory, " "recollection, " "mild cognitive impairment." The search for studies ended in July 2018. Our literature search yielded 484 studies from PubMed, and references from these relevant articles were examined for possible inclusion in our analysis. In line with previous protocols, when the title and abstract were insufficient to identify the study for inclusion, the full-text underwent further examination. In cases of uncertainty, articles were examined by the authors and research assistants and removed based on mutual agreement. The most frequent reasons for exclusion were that articles (a) did not apply TMS or (b) did not assess episodic memory, see Figure 1.

Inclusion Criteria and Selected Studies
For inclusion in the systematic review, the studies were limited to (1) single pulse TMS, paired pulse TMS, rTMS, or TBS protocols that measured (2) episodic memory with a free-recall, cuedrecall, or recognition test (3) in healthy younger adults, (4) healthy older adults, or (5) older adults with a clinical condition, (6) and those that were published in English.
A total of 59 studies were incorporated into the systematic review based on the search terms and inclusion criteria. However, inclusion in the meta-analysis was focused on rTMS studies of healthy younger adults. Seven of the remaining studies were excluded from the meta-analysis because there was insufficient reporting of data necessary to calculate effect sizes (e.g., means, SDs, t-, or F-values). Thus, 37 articles examining the effects of rTMS on episodic memory in healthy younger adults were included in the meta-analysis. See Figure 1 for a flow diagram illustrating this process. A total number of 245 effect sizes were extracted from the selected studies. For experiments that pooled participants' data across multiple studies, only participants who were uniquely recruited for each study were analyzed. Details about the calculation of effect sizes are described below in the data extraction section.

Quality Assessment
In order to assess the methodological quality of the selected studies, we were guided by the assessment criteria recommended by Higgins et al. (2008), which is to categorize study characteristics as low risk, unknown risk, or high risk based on whether there was random allocation of participants to conditions, blinding of participants, blinding of experimenter, blinding of experimenter assessing the outcome measure(s), selective reporting of outcome measures, and participant attrition. Characteristics were coded as unknown when insufficient details were provided; studies were categorized as high risk when it was clear that the conditions were not met, such as when it was implausible for the participants and/or experimenters to have been blind to the administration of the experimental vs. control conditions. We also extended these criteria to include additional assessments pertinent to TMS studies, namely the targeting procedure and type of control condition. For targeting procedure, if a study used structural or functional MRI neuronavigation, international 10-20 EEG positioning system, or distance measurements, they were classified as low risk, unknown risk, and high risk, respectively. A variety of control conditions were implemented, such as active, sham, and no TMS conditions, that were classified, as low, unknown, and high risk, respectively.

Coding of the Studies
Moderators that were hypothesized to modulate the size of effects of rTMS on episodic memory were coded for each study. In order to code for moderators, the full-text articles were examined for explicitly stated information regarding the moderators of interest. In situations where it was not clearly stated, the moderators were inferred from the study protocol (e.g., the study did not state if stimulation occurred "online" or "offline, " but, if the task figures showed that stimulation occurred prior to a cognitive task, then we coded stimulation timing as offline). Studies were independently coded by both the first author and research assistants to establish reliable coding and any discrepancies were resolved by the authors.

Coding of rTMS Parameter Moderators
Frequency was coded as a categorical variable with 1, 5, 10, and 20 Hz being coded as "one, " "five, " "ten, " and "twenty, " respectively. Theta burst frequencies were separated and coded as intermittent theta burst (iTBS) or continuous theta burst (cTBS). One study used a modified iTBS protocol with short iTBS (2 s) trains interleaved with stimulus presentation and was classified as online iTBS. None of the obtained studies implemented an intermediate theta burst (imTBS) procedure. A remaining 15 effect sizes across 3 studies implemented stimulation frequencies (i.e., 6.8, 7, 10.7, 17.5, and 18.7 Hz) that did not fall in these categories and, thus, were not used in analyses examining the effects of frequency.
Motor threshold intensity (MTI) can be calculated numerous ways, such as based upon individual RMT or AMT. Although procedures that estimate motor thresholds using AMT vs. RMT or MEP vs. visually-evoked responses result in slightly lower levels of stimulator intensity (Rossini et al., 2015), there were not enough studies that clearly implemented either RMT or AMT with either MEP or visually-determined thresholds to allow for systematic investigation of the role of different thresholding procedures. Therefore, MTI was treated as a categorical variable. Studies that reported using stimulation parameters at, above, or below 100% of motor threshold were coded as "at, " "above, " or "below, " respectively.
Studies were coded as "online" protocols when stimulation occurred simultaneously with the study's experimental task, such as when rTMS co-occurred with stimulus (e.g., words or pictures) presentation during either encoding or retrieval. Although studies differed with regards to the precise timing of the delivery of TMS pulses, for example, at stimulus onset or varying intervals from stimulus onset, there were not enough studies to examine possible differences in the effects of timing differences in the administration of online rTMS protocols at this time. Therefore, we collapsed across different online timing protocols. Studies where stimulation occurred in isolation without any cooccurring task, such as before or after an encoding or retrieval task, were coded as "offline" protocols.
The hemisphere that was stimulated was coded as "left" or "right." Studies that targeted cortical regions in frontal cortex or parietal cortex were coded as "FC" and "PC, " respectively. There were too few studies that targeted regions beyond the FC and PC (e.g., occipital cortex) for valid analysis and inference of aggregated effect sizes at this time. Therefore, a metaanalysis investigating the effects of stimulation that targeted either consolidation processes or regions outside the PC and FC was not possible at this time. Again, these gaps in the literature highlight areas where more research is needed.
Lastly, we coded studies based on the type of control condition used for comparison to the primary "treatment" condition of interest in order to assess the potential differences between different types of control conditions. Classification of an "active" control occurred when experimental stimulation was compared to any active rTMS control site except for the vertex. A "vertex" control comparison was coded separately. When studies implemented either a sham coil or applied sham stimulation by tilting the coil away from the participant's head, they were coded as "sham" controls. In situations where an experimental stimulation group was compared to a condition in which no rTMS was applied, they were coded as "no TMS."

Coding of Experimental Design Moderators
The "stimulated memory processes" (encoding or retrieval) were coded as the main memory processes that were most likely to be primarily affected by stimulation based on when stimulation was applied. A study was coded as targeting "encoding" if the stimulation parameters occurred within seconds to minutes prior to, during, or immediately after (e.g., within 400 ms of stimulus offset) encoding each item. For example, conditions with a train of cTBS administered within minutes prior to encoding all items was coded as primarily targeting "encoding." Similarly, conditions with brief trains of 1 Hz rTMS during picture presentation were coded as targeting "encoding." Situations in which stimulation occurred shortly after picture presentation offset (i.e., within milliseconds), were coded as targeting encoding. While such protocols may affect both encoding and early consolidation processes, studies were categorized using this range in order to maximize the number of effect sizes for reliable analysis. The critical distinction for this analysis was to compare the aggregate effect size of rTMS protocols that clearly impacted retrieval processes vs. encoding/early-consolidation processes. Unfortunately, there was only one study that applied stimulation during consolidation, which prevented its inclusion in the metaanalysis. If stimulation occurred immediately before or during a retrieval task (e.g., cTBS prior to attempting to retrieve all previously encoded words or 1 Hz rTMS during presentation of each picture in a recognition test), the study was coded as primarily targeting "retrieval." To test for a difference in the effects of rTMS on episodic retrieval processes that are hypothesized to be distinct by dual-process theories, effect sizes were extracted for conditions hypothesized to primarily measure either recollection or familiarity. Conditions were aggregated and coded as assessing "recollection" if they assessed free recall, cued recall, "remember" responses in remember/know recognition memory paradigms, and paradigms that reported associativeor source-recognition memory performance. Associative-or source-recognition estimates (i.e., associative-or source-hits, corrected recognition scores, or d ′ ) were extracted from recognition memory tasks that required the retrieval of qualitative information about the study of a specific piece of information (e.g., if participants studied a soccer ball object and were asked to recall its spatial location), as in source memory paradigms. Conversely, conditions were aggregated and coded as assessing "familiarity" if they assessed and reported recognition memory performance (e.g., item memory hit rates, corrected scores, d ′ ) on tests that did not require additional qualitative information to be retrieved. This coding scheme included "know" responses for remember/know paradigms.

Data Extraction
Means (M), standard deviations (SD), and sample sizes (n) were obtained for relevant episodic memory outcome measures for each study. When standard errors (SE) were reported, they were converted to SD. In cases where there was insufficient statistical information provided in the text or tables, they were collected from graphical representations using GetData Graph Digitizer data extraction software (downloaded at http://www.getdatagraph-digitizer.com). In the event that insufficient data could be obtained, attempts were made to contact the corresponding authors for additional data. Descriptive statistics and effect size calculations were done independently by the first author and research assistants to cross-check the calculations and resolve any discrepancies. The table of experiments included in the meta-analysis, their moderator variables, and all of the effect sizes that were extracted from the experiments are presented in an open-access database here: https://osf.io/6kbqe/? view_only=226e600d6e1546748b52f6bbfc593140. For additional information or to contribute additional data, please contact the first author. Our intention for sharing the data is so that effect size calculations can be confirmed and so that additional studies can be continuously added as the literature grows, which will facilitate future meta-analysis and systematic investigation.
Effect sizes for between-subject comparisons were computed with the following formula:

SD Pooled
When studies reported comparisons between gain scores (e.g., Post-experimental-Pre-experimental) the effect size was calculated with the following formula: Effect sizes for within-subject comparisons were calculated with the following formula: This was done to obtain an accurate effect size estimate for within-subject designs (Cumming, 2012) using the commonly recommended approach for meta-analyses (Lakens, 2013). Alternative routes that address dependencies in within-subject designs, such as averaging effect sizes, can lead to a reduction in statistical power and the loss of information can result in less precise estimates and SE's (Hedges and Pigott, 2001;Cheung, 2014). In a similar vein, limitations also arise when calculating Cohen's d rm , which attempts to account for the dependencies by correcting for the correlation in within-subject designs. As correlations between conditions are rarely reported, estimations may be too conservative, especially when correlations are high (Lakens, 2013). Hedge's g and g av effect sizes were calculated, which provides a less biased estimate of the true effect than Cohen's d, especially for studies with small sample sizes (Hedges and Olkin, 1985). In order to do this, we first calculated Cohen's d and then applied a Hedge's g correction (see Cumming, 2012) with the following formula: With the exception of 5 studies, we were able to collect multiple effect sizes per study. Thus, we can no longer assume statistical independence, as within-study effect sizes may be more similar than effect sizes from other studies. Thus, dependent effect sizes may contain redundant information that becomes generally less informative when effect sizes are correlated. This form of dependency may occur when multiple effect sizes can be calculated from multiple outcome measures (Van den Noortgate et al., 2013). For example, in examining the effects of rTMS on episodic memory, a study can test a set of stimulation parameters (e.g., 10 Hz rTMS at 100% AMT to the DLPFC) on multiple memory outcomes, such as recall, recognition and source memory. A second dependency occurs when effect sizes are correlated in situations where an experimental condition is compared to multiple control groups or multiple experimental conditions are compared to one control group. For instance, many of the studies compared experimental stimulation to combinations of active control, sham, or no stimulation conditions, while other effect sizes were comparisons between multiple experimental stimulation conditions to a single control condition. This can be problematic in meta-analyses and increase type 1 errors if dependent effect sizes are assumed to contain unique information. Therefore, to address multiple dependency concerns we conducted a three-level random effects metaanalysis with Hedge's g and g av corrections, which is detailed below (for additional information on dependencies and threelevel models see, Van den Noortgate et al., 2013).

Three Level Modeling
We conducted a three-level random effects model using the Metafor package in R (Viechtbauer, 2010). The typical random effects approach allows for effect sizes to vary due to sampling variation and due to differences between studies. In order to account for studies that contain multiple effect sizes, we can extend this model to a three-level random effects model to take into account the correlation of within-studies effect sizes. This approach does not assume independence among effect sizes as is done with other meta-analysis techniques. Multilevel modeling allows for an alternative for dealing with the dependency of effect sizes when correlations between outcomes are not reported or there is insufficient data to estimate the dependency (for additional details, see Assink and Wibbelink, 2016). The hierarchical structure of the data allows effect sizes to vary between participants (level 1; sampling variance), outcomes (level 2; within studies), and studies (level 3; between studies). In this case, individual studies (denoted by the subscript k) can contain one or multiple outcomes (denoted by the subscript j). The studies (i.e., Study ID) and individual effect sizes (i.e., Effect Size ID) were entered into the model as random effects. The simplest model with no moderators and three residuals (v k, u jk, e jk ) that are assumed to be normally distributed with a mean of zero is: where an observed effect size for outcome j within study k is represented by ES jk. The overall population effect size mean across all outcomes and studies is represented by β 0. The v k element is the random mean deviation effect in study k from the overall effect. The random deviation of jth population effect in study k from the mean effect in study k is denoted by u jk. Lastly, e jk is the random residual error from the population effect due to sampling variation. Parameter estimates in the meta-analysis include the between-study variance (σ 2 v ) and within-study variance (σ 2 u ), and the sampling variance is not reported because it is calculated from reported data prior to the analysis. Therefore, observed significant between-study variance differences reflect that effect sizes systematically vary between studies to a greater degree than would be expected due to random differences in sampling variance (due to chance). Meanwhile significant within-study variance above zero signifies that the observed effects within a study vary across individuals to a greater degree than chance.
In order to implement a parsimonious model, we examined if our three-level model provided a better model fit over more simplistic models with multiple ANOVAs. There was a significant difference between a single and two-level model (within-study variance), likelihood ratio test (LRT) = 41.75, p < 0.01. There was also a significant difference between two and three-level models (between-study variance), LRT = 41.83, p < 0.01. Thus, the threelevel modeling of the dependent effect sizes provided a better model fit than treating the effect sizes as independent.

Publication Bias
A prevalent issue in systematic reviews and meta-analyses is how to address missing data due to publication bias. This problem occurs when significant results are more likely to be published than non-significant findings. If left unaccounted for, publication bias may lead to spurious results. In this instance, publication bias may lead to the true effect of rTMS on episodic memory being inflated. We used two methods to assess the degree of possible publication bias in the data: funnel plots and an extension of Egger's regression method (Egger et al., 1997)

Funnel Plots
One approach to identify and estimate missing data is to examine funnels plots for asymmetry (Sterne et al., 2011). In order to do this, a scatter plot is created with the effect size on the x-axis and the sampling variance (or some variation taking into account the sample size) on the y-axis. Publication bias is less likely to be a concern when effect sizes are distributed symmetrically around the mean effect and are clustered around the mean effect size. Conversely, support for possible publication bias manifests visually in these plots with an asymmetrical distribution of effect sizes around the mean effect size, often with larger sampling variances producing larger effect sizes.
Contour-enhanced plots offer an extension of assessing publication bias by drawing a reference line at zero with contour color changes at different significance levels (0.01 and 0.05) as the effect sizes are plotted against precision (1/SE, inverse of the standard error). In this case, publication bias can be examined with missing effect sizes observed in nonsignificant regions of the plot. However, there are other possible causes of asymmetry besides publication bias. Asymmetry due to heterogeneity of potential moderators is another factor that can influence asymmetry and can be assessed with residual plots. Here, residuals are plotted against their standard errors at different levels of the moderators (e.g., online vs. offline protocols). Similar to the contour-enhanced plots, significant levels can be added, and asymmetry of the moderators can be examined in relation to the predicted effect sizes. In an attempt to address both publication bias and other potential forms of publication bias we examined contour enhanced and residuals plots.

Extension of Egger's Test
Additional traditional approaches for formally testing for publication bias, such as trim and fill and Egger's test, are complicated with three-level modeling. Although correcting (e.g., trim and fill) for publication bias is important (Copas, 1999), obtaining more precise estimates through correction methods may only occur if the bias is large (Hedges and Vevea, 1996). Therefore, we decided to implement an extension of Egger's test for more complex models by entering the standard error into the model as a moderator. Essentially, this provides an extension of Egger's regression method which involves estimating funnel plot asymmetry by measuring how much the regression line (effect size vs. precision, 1/SE) deviates from zero. However, future research is needed to determine how to properly assess missing data in three level meta-analyses, as the numerous available methods to measure and correct for publication bias each come with limitations and, to our knowledge, have not been evaluated in the three-level meta-analysis approach.

RESULTS
The systematic review incorporated 1,532 participants across 59 articles, see Figure 1. This resulted in an average age of 33.12 (18-80) with 61% females and 39% males in the studies. Of the 59 articles, 44 contained rTMS procedures focusing on younger adults, while 12 articles included healthy older adults and those with memory impairments (i.e., aMCI, Alzheimer's disease) or other disorders (i.e., depression, fibromyalgia, alcohol dependency). The small number of studies focusing on older adults made a meta-analysis unfeasible for comparison with younger adults. Therefore, our approach was 2-fold. First, we conducted a three-level meta-analysis focusing on the effects of rTMS on episodic memory in younger adults. Second, we     synthesized the results of studies using rTMS to modulate episodic memory in older adults in the form of a narrative review. A total of 1,016 participants were included in the metaanalysis of the 37 articles that were identified. The average age of participants across the studies was 25.70 years old with 39% male and 61% female. The average sample size was 24.19 (range 10-69). For a breakdown of the study designs and protocols, see Table 1.

Quality Assessment
The proportion of studies classified as low, medium, and high risk for the different categories of the quality assessment analysis are presented in Figure 2. Overall, many of the studies implemented appropriate methodological protocols for random allocation of participants to conditions, blinding of the experimenters assessing effects on the outcome measure(s), and non-selective reporting of outcome measures. However, several studies either did not implement proper blinding of participants and/or experimenters, or did not provide enough information to confirm proper blinding procedures. There was low concern for the degree of participant attrition that was reported. For our additional quality assessment measure of control type, we found that the majority of studies implemented either active or sham controls. With regards to the rTMS targeting procedure, we found that the majority of studies implemented either low risk neuronavigational procedures or used the 10-20 system for targeting cortical regions of interest.

Evaluation of Publication Bias
Several approaches were used to examine possible publication bias. First, we examined the enhanced-contour plot (i.e., effect size vs. 1/standard error). The reference line is set to g = −0.06 with the contour shading representing p-value significance. The enhanced-contour plot revealed some asymmetry, but no systematic patterns of missing non-significant effect sizes, see Figure 3. Given the numerous effect sizes collected, the larger effect sizes observed in the funnel plots are expected. Importantly, they appear to be symmetrically distributed. An examination of influential points revealed no Cook's d values above 0.5.
Next, we examined the residual plots (residuals vs. standard error) for each moderator variable, see Figure 4. The reference line is set to average effect size with contour shading once again representing p-value significance. The residual plots revealed some asymmetry across the moderator variables.
Interpreting funnel plots is inherently difficult (Terrin et al., 2005) and may be further complicated in three-level models (e.g., observed asymmetry may not be problematic). For example, clustered data points due to the three-level structure could potentially be misinterpreted as bias. Therefore, in a more formal test of asymmetry we entered standard error into the model. This revealed no significant moderating effect of standard error on effect size, F (1, 244) = 1.76, p = 0.18, indicating that the funnel plot was not significantly asymmetrical. Collectively, this suggests that publication bias is unlikely to be problematic in this analysis. However, future research is needed to properly address how to handle detecting and correcting for missing data (e.g., publication bias) in three-level models.

Overall Relationship Between rTMS and Episodic Memory
We initially entered individual effect sizes (i.e., Effect Size ID) and studies (i.e., Study ID) as random effects. This model revealed the average effect size (g = −0.06, SE = 0.08) was not significantly different from zero, t (244) = −0.75, p = 0.45, see Figure 5. Thus, we do not see an overall effect of rTMS on episodic memory performance, which is unsurprising given that most studies selected and implemented protocol parameters to either enhance or impair episodic memory functioning.

Heterogeneity of Effect Sizes
We tested the heterogeneity of effect sizes to examine if differences between effect sizes were systematic and not due to random sampling variance. If a substantial amount of betweenor within-study effect size variance is found, potential moderators can be added into the model to explain this variability. There was significant within-study variance (level 2), estimate = 0.19, p < 0.01 and between-study variance (level 3), estimate = 0.12, p < 0.01, demonstrating that effect sizes varied both within and between studies. Consequently, 44.58% of the variance is accounted for by differences in effect sizes within studies, 27.25% by variance between studies, and 28.07% by random sampling variance. The significant variability across both levels highlights the substantial differences in effect sizes between studies and within the same study. Therefore, we next conducted a moderation analysis in order to assess whether the observed differences in effect sizes were moderated by other factors.

Effect of Frequency and Motor Threshold Intensity
Although we did not find support for MTI moderating episodic memory performance overall, we anticipated that MTI may interact with frequency. Recent findings have shown that varying MTI's during iTBS follow an inverted U-shaped pattern in modulating working memory performance (Chung et al., 2018). An investigation of MTI (above, at, or below) and frequency (1,  Tables 3, 4. In order to examine what was driving this interaction we considered how differences between MTI (above vs. at vs. below) at each frequency affected episodic memory. However, the number of comparisons were restricted due to the limited number of studies that included FIGURE 5 | Overview of effect sizes included in meta-analysis. Reference line set to the average effect size (g = −0.06). Legend represents study sample size ranges: 10-19, 20-29, 30-39, 40-49, 50-59, 60-69. frequencies at different MTIs. Follow-up tests revealed that the advantage of 1 Hz rTMS compared to 10, 20, iTBS, and cTBS protocols in producing significantly stronger enhancing effects of experimental/treatment rTMS vs. control/sham rTMS on memory performance was specific to intensities below MTI. There was also a significant effect of 20 Hz rTMS "below" MTI (β = −0.38, SE = 0.17) in that it produced a larger negative effect of experimental/treatment rTMS vs. control/sham rTMS compared to 20 Hz "at" MTI (β = 0.40, SE = 0.35), F (1, 202) = 3.98, p = 0.04. There were no other differences between levels of MTI for any other frequency, see Tables 3-5 1 .

Effect of Frequency and Stimulation Timing
There was a significant interaction of stimulation timing (online vs. offline) and frequency (1, 5, 10, 20, iTBS, cTBS) on episodic memory effects, F (9, 221) = 4.16, p < 0.01. We first examined differences between frequencies within online and offline stimulation protocols. Subsequently, we looked at differential outcomes between online vs. offline protocols for specific frequencies. However, the number of comparisons we could make was limited due to the small number of studies that included online or offline protocols at each frequency.
Following offline stimulation protocols, 1 Hz rTMS (β = 0.58, SE = 0.21) led to significantly larger enhancing effects compared to both iTBS (β = −1.84, SE = 0.53) and cTBS (β = −0.06, SE = 0.13). In addition, offline 20 Hz rTMS (β = 0.94, SE = 0.45) led to significantly larger enhancing effects compared to both iTBS and cTBS. Interestingly, offline iTBS led to significantly larger negative effects compared to cTBS. No other comparisons were significant, although a full decomposition of the interaction was not possible because the number of studies with offline stimulation for each type of frequency were limited, see Table 7.
Additional follow-up tests revealed that stimulation timing significantly influenced the directionality of the effects of 20 Hz rTMS. Specifically, online 20 Hz rTMS led to larger negative effects, while offline 20 Hz rTMS led to larger enhancing effects of experimental/treatment rTMS vs. control/sham rTMS on memory performance, F (1, 221) = 7.33, p < 0.01. There were no other significant differences.

Memory Processes and Hemispheric Stimulation
Findings of individual studies have reported that there may be different effects of rTMS when stimulating the left or right hemisphere during either encoding or retrieval. To assess potential aggregate effects across studies, we examined the interaction between stimulated hemisphere (left vs. right hemisphere) and memory process stimulated (encoding vs. retrieval processes). However, there was no significant interaction  Level 2 is the variance attributed to all the moderators that varies within studies. Level 3 is the variance attributed to between study differences. * p < 0.05, ** p < 0.01, *** p < 0.001.

Retrieval Type and Cortical Regions
Although it is plausible that the effects of stimulating either frontal or parietal cortex may differ for familiarity-vs. recollection-based retrieval processes, we found no significant interaction between effect sizes of the stimulated cortical region (FC vs. PC) and retrieval type (familiarity vs. recollection), F (4, 233) = 0.45, p = 0.77, although, again, inferences are currently limited by the relatively small numbers of studies that could be included in such a moderation-based meta-analysis.

DISCUSSION
A systematic review of the literature spanning the past 30 years yielded 59 articles that focused on how TMS can modulate episodic memory in younger and older adults. To better understand the circumstances of when rTMS will enhance or impair memory, we conducted a three-level random effects meta-analysis of 245 effect sizes aggregated from 37 articles focused on young adults. Because the effects of some studies used parameters designed to impair memory to probe the necessity of a stimulated brain region in a particular process relevant for episodic memory, whereas others implemented designs to facilitate episodic memory, the mean effect size, collapsing across all parameters is relatively uninformative.

Enhancing Effects of 1 Hz rTMS
Moderation analyses revealed that stimulation frequency modulated the effect of rTMS on episodic memory in younger adults, with larger facilitatory effects of 1 Hz (g = 0.59) rTMS compared to 10 (g = −0.28), 20 (g = −0.22), iTBS (g = −0.42), and cTBS (g = −0.05) protocols. The facilitatory effect of 1 Hz rTMS may appear counter-intuitive because low frequency (e.g., < 1 Hz) rTMS is typically thought to lead to inhibitory effects on neural processing of targeted regions, whereas, many high frequency protocols are generally thought to lead to enhancing cortical excitability (Walsh and Cowey, 2000;. Relatedly, the numerically larger negative effects of iTBS compared to cTBS paradigms is also surprising given that cTBS is typically assumed to have longterm-depression-like inhibitory effects (LTD) that can suppress cortical excitability in targeted regions whereas iTBS is assumed Only available frequency comparisons are listed that contained observable data for at least one of the motor threshold levels. Empty cells represent no observed effect sizes for comparison *p < 0.05, **p < 0.01.
to lead to long-term potentiation-like effects (LTP) that can enhance cortical excitability (Huang et al., 2005). However, these findings stem from physiological studies that modulated cortical activity in the motor cortex (Jelić et al., 2015), while the studies in the present meta-analysis generally targeted prefrontal and parietal regions. One possible explanation is that the effects of TMS may differ according to the structural and functional architecture of the targeted region (e.g., motor cortex vs. prefrontal/parietal cortex). However, strong inferences based on these findings are currently limited due to the numerous differences between studies in the specific rTMS parameters and experimental designs that were used, which may interact to affect episodic memory.

Interactions Between rTMS Frequency and Stimulation Intensity
Interaction analyses revealed facilitatory effects of 1 Hz rTMS for both online and offline stimulation. However, the effects of 1 Hz rTMS was driven by stimulating below MTI. While the Means, standard errors, and confidence intervals. The number of studies and effect sizes are denoted by k and n, respectively. Given the available studies not all possible combinations were represented.
underlying mechanisms of 1 Hz rTMS on episodic memory are poorly understood, one possible explanation comes from the principle of stochastic resonance, whereby neural communication is enhanced by low levels of noise that can push a neural signal past threshold (Schwarzkopf et al., 2011;Silvanto and Cattaneo, 2017). Thus, 1 Hz rTMS may enhance episodic memory by interacting with intensity to enhance the signal-to-noise ratio in the stimulated network of neurons. However, this process is affected by a variety of other factors that further complicate the process. For example, the effects of MTI are known to be sensitive to "brain states" and changes in brain activity can influence if below-or above-MTI leads to facilitatory or inhibitory effects (Romei et al., 2016). Stimulation intensity also interacted with 20 Hz rTMS with attenuating effects of below-MTI compared to at-MTI. These findings are not consistent with stochastic resonance, where below MTI stimulation should lead to enhancing effects. However, the situations in which different levels of intensity will lead to either enhancing or inhibitory effects have also been known to depend on brain states. One reason for the observed effects may be due to differences in experimental designs, as task stimuli or cognitive processes engaged likely resulted in different activation states in the targeted brain regions. Relatedly, previous work has shown that cortical excitability varies with different MTI levels following iTBS (Chung et al., 2018) protocols. Thus, it is possible that the enhancing or inhibitory effects of rTMS may vary with MTI, which may also differ depending on the frequency of stimulation.

Interaction Between rTMS Frequency and Stimulation Timing
Additional analyses revealed that the effects of episodic memory performance in young adults depended on frequency and stimulation timing. Specifically, the enhancing effects of 1 Hz compared to 20 Hz was driven by online protocols. On the other hand, both 1 Hz and 20 Hz rTMS protocols had enhancing effects compared to either iTBS or cTBS paradigms when stimulation occurred offline. Interestingly, online 20 Hz rTMS protocols had inhibitory effects on episodic memory; however, offline 20 Hz rTMS had enhancing effects on episodic memory. These findings suggest that 20 Hz rTMS modulates episodic memory through different mechanisms when applied online vs. offline. For example, online paradigms may disrupt ongoing neural activity while offline paradigms may potentiate episodic memory processes through LTPlike effects. Taken together, while the meta-analysis suggests novel perspectives on the effects of interacting factors on rTMS studies of episodic memory, the limited numbers of studies with combinations of factors spanning the entire parameter space may have hindered our ability to reliably estimate the role of many factors (e.g., targeted hemisphere, cortical region, memory process, and retrieval type) that are hypothesized to moderate the effects of rTMS on episodic memory. Moreover, the current meta-analysis did not find a moderating effect of the type of control condition (i.e., active, active vertex, sham, or no TMS) that was implemented. Importantly, these findings do not suggest that there are no differences between using active, vertex, sham, or no TMS control conditions. There are substantial gains in the quality of inference that can be obtained with comparison to proper control conditions. Inferences about the role of a particular brain region in a specific process are strengthened when effects are contrasted against conditions in which active stimulation at the same level of intensity and frequency is applied to either the targeted brain region at a different time point or to a different brain region that is not hypothesized to be involved in the process of interest at the time of stimulation. Therefore, researchers should still take careful consideration in the type of control condition that is most appropriate and feasible for specific experimental designs (for further information on control conditions see, Davis et al., 2013;Duecker and Sack, 2015).
In the next section we discuss how the meta-analysis on younger adults relates to the pattern of observed effects of rTMS effects on episodic memory functioning in healthy older adults and those with clinical disorders. This systematic review of the effects of rTMS on older adults revealed several similar findings that were observed in the meta-analysis on younger adults, thereby providing some sense of cross-validation of the findings in different populations.

COMPARISONS TO rTMS EFFECTS ON OLDER ADULTS AND CLINICAL POPULATIONS Enhancing Effects of 1 Hz rTMS When Applied Below MTI
A similar pattern of effects of 1 Hz rTMS on episodic memory that was revealed in the meta-analysis on young adults was found in 1 Hz rTMS studies on older adults with below-MTI stimulation enhancing performance (Turriziani et al., 2012) while above-MTI had no effect (Davis et al., 2017). Turriziani et al. (2012) administered offline 1 Hz rTMS at 90% MT to the left or right DLPFC of patients with mild cognitive impairment following encoding of faces, buildings, and words and found enhanced recognition memory. In contrast, Davis et al. (2017) administered offline 1 or 5 Hz rTMS at 120% MT to the left DLPFC of healthy older adults prior to an encoding task and observed null effects on associative memory (although there were differential effects on brain network activity related to successful memory performance, for details see Davis et al., 2017). Although several parameters differed between these studies, the difference in outcomes of 1 Hz rTMS as a function of MTI is consistent with the findings of studies included in the meta-analysis irrespective of differences among those studies. However, there are not enough studies available to compare all combinations of parameters of offline 1 Hz rTMS below-, at-, or above-MTI stimulation.

Impairments From Online 20 Hz rTMS Below MTI
Twenty Hz rTMS findings in healthy older adults found a similar pattern that was revealed in the meta-analysis with impairments to episodic memory when 20 Hz rTMS was administered online or below MTI.  reported attenuating effects of online 20 Hz rTMS at 90% MT to the left DLPFC during picture retrieval. In addition, Manenti et al. (2011) found attenuating effects on episodic memory when online 20 Hz rTMS at 90% RMT was applied to left vs. right DLPFC during either encoding or retrieval. However, this effect was largest in low performing adults when stimulation was applied to left compared to right DLPFC during encoding. Thus, it is possible that the effects of 20 Hz rTMS on episodic memory may rely on more complex interactions that were not possible to examine in the current meta-analysis due to the limited number of studies that have investigated the combinations of parameters.

Enhancing Effects of Offline 20 Hz rTMS
Two studies on older adults with memory disorders found a similar pattern that was revealed in the meta-analysis with enhancing effects of multi-session administration of offline 20 Hz rTMS (Wang et al., 2014;Nilakantan et al., 2017). In prodromal Alzheimer's patients, Koch et al. (2018) reported that 10 sessions (5 sessions/week) with 2-week separation between cross over stimulation (active vs. sham) of offline 20 Hz rTMS at 100% RMT to the precuneus enhanced delayed recall (Koch et al., 2018). Similarly, in an individual with mild cognitive impairment, Cotelli et al. (2012) found enhanced associative memory compared to healthy controls following 10 sessions (5 sessions/week) of offline 20 Hz rTMS at 100% MT to the left parietal cortex .
Importantly, the enhancing effect of offline 20 Hz contrasts with an overall impairment of online 20 Hz rTMS on episodic memory (see Table 2). However, offline 20 Hz studies implemented multi-session protocols (e.g., 5 sessions/week for 2 weeks) while all nine online 20 Hz studies implemented single session protocols (see Table 1), so this confounds interpretation of the potential difference in the direction of effects of 20 Hz rTMS when applied offline vs. online. This also reveals a gap in the literature as few studies have implemented multi-session protocols with younger adults. Collectively, these findings are encouraging with regards to the potential for applying multi-session offline 20 Hz rTMS protocols to enhance episodic memory in older adults, even those with MCI or AD.

Enhancing Effects of Offline 5 Hz rTMS Below MTI
Three studies implemented 5 Hz rTMS and found enhanced behavioral outcomes (Solé-Padullés et al., 2006;Peña-Gomez et al., 2012) or modulation of neural activity (Davis et al., 2017). In low performing older adults, Solé-Padullés et al. (2006) reported enhancing effects when offline 5 Hz rTMS at 80% MT was applied prior to encoding face-name associations. In addition, in carriers of a major genetic risk factor for Alzheimer's disease (Apolipoprotein ε4 allele; APOE), enhancing effects were reported following offline 5 Hz rTMS at 80% MT prior to a face-name encoding task (Peña-Gomez et al., 2012). In healthy older adults, Davis et al. (2017) reported modulated neural activity, but no differences in memory performance following offline 5 Hz rTMS at 120% MT prior to an associative encoding task. There are no studies on young adults that administered offline 5 Hz rTMS. This prevents making comparisons and highlights the need for future work to address this gap.

Enhancing Effects of Offline 10 Hz rTMS
Five studies applied offline 10 Hz rTMS to older adults or clinical populations that included a wide age range and all but two found enhancing effects on episodic memory (see Table 1). The offline 10 Hz studies that found null effects used a different stimulation intensity than the others (Rektorova et al., 2005) or applied rTMS to motor cortex in fibromyalgia patients (Baudic et al., 2013). The studies that found enhancing effects of offline 10 Hz used stimulation intensities below-or above-MTI (Drumond Marra et al., 2015;Qiao et al., 2016;Kavanaugh et al., 2018). These findings diverged from the results of the meta-analysis that revealed a trend for impairments of online 10 Hz rTMS in younger adults. Differences in the stimulation timing between studies is one potential reason for the different effects; however, direct comparison is complicated because no studies to date have applied offline 10 Hz in young adults.

LIMITATIONS OF CONDUCTING A META-ANALYSIS ON AVAILABLE rTMS RESEARCH
Although the systematic review and meta-analysis revealed some promising findings, a considerable amount of research is needed to more fully understand the roles of both specific aspects of the rTMS parameters (frequency, intensity, timing, hemisphere, cortical region, and control conditions) and design considerations (targeted memory process and retrieval type) on the effects on episodic memory functioning, as well as how these numerous factors may interact. Because there are not published studies for every pairwise combination of the rTMS parameters and design aspects of interest, there were many missing cells in the meta-analysis. Relatedly, because of missing data cells and small numbers of observed effects sizes for each possible comparison, we weren't able to investigate higherorder interactions.

Effects on Objective vs. Subjective Aspects of Episodic Memory
It should also be noted that the current meta-analysis focused on accuracy as the outcome measure of episodic memory performance. We did not address the potential moderating effects of rTMS and study design parameters on other measures of episodic memory performance such as reaction times and confidence ratings. These are important issues to investigate in future research as the body of literature accumulates.
For example, while the meta-analysis on accuracy measures failed to reveal a main effect or interaction with the cortical region that was targeted or the memory process that was stimulated, based on the systematic review of the published studies, it seemed to be the case that several studies found that applying rTMS to regions of parietal cortex affected the subjective experience of remembering as evidenced by modulation of confidence ratings even in the absence of effects on objective measures of accuracy. For example, Ye et al. (2018) recorded participants playing an action-adventure video game and then applied offline 1 Hz rTMS at 110% AMT to the precuneus prior to snapshot images from the video game and then tested subjects' temporal order memory and confidence ratings. The authors found decreases in confidence ratings during retrieval following precuneus compared to vertex stimulation, despite no differences in memory performance (Ye et al., 2018).
In a similar fashion, Bonnì et al. (2015) applied offline cTBS at 100% MT to the precuneus, posterior parietal cortex, or vertex prior to associative (color-context) memory tests of studied object images (e.g., apple outlined in red). While there were no differences in hit rates among the conditions, there was a significant decrease in the number of source memory errors (e.g., correct judgment that the apple was outlined in red not green) that were made following precuneus stimulation (Bonnì et al., 2015). Relatedly, offline 20 Hz rTMS at 100% MT to the left angular gyrus prior to retrieval of previously encoded pictures led to decreased recognition relative to superior parietal lobe stimulation, and a decrease in recollective details (by altering response bias) compared to either superior parietal or vertex stimulation (Sestieri et al., 2013). Moreover, offline cTBS at 70% RMT to the left angular gyrus prior to encoding auditory stimuli reduced subjective confidence ratings for source memory while leaving free-and cued-recall of word pairs unaffected (Yazar et al., 2014). Taken together, it seems that rTMS to parietal regions tends to affect the subjective experience of remembering even when there are null effects on episodic memory accuracy. However, too few studies have been published to allow for meta-analysis of the effects of moderators of the size of effects on episodic memory. The same was true for other important measures of episodic memory performance such as reaction times. Future research should address these limitations.

Effects on Encoding, Consolidation, and/or Retrieval
Another limitation has to do with the complication of interpreting the effects of rTMS on episodic memory due to the design of many studies, particularly studies with multiple rTMS sessions or offline rTMS, because it is difficult to know if the effects of rTMS are modulating encoding, consolidation, and/or retrieval processes. Relatedly, to our knowledge, only one study implemented rTMS to definitively target consolidationbased processes. In a pair of studies, Sandrini et al. (2013) investigated the role of right DLPFC in the reconsolidation of memories for studied words of objects by administering offline 1 Hz rTMS at 100% MT. The first experiment spanned 3 days. On day one, participants learned a list of words with contextual reminders (e.g., specific room and colored bag containing the words). On day two, participants were assigned to receive either rTMS to right DLPFC during a pure consolidation phase without reactivating the memories or rTMS to right DLPFC or the vertex before reactivating the memories in a contextual reminder reconsolidation phase. Finally, on day three, participants underwent a free recall task in which participants recalled as many words as possible. The findings revealed that rTMS prior to memory reactivation led to memory enhancements compared to the no reactivation and vertex control conditions. To examine if rTMS effects were specific to reconsolidation and rather than either encoding or retrieval processes, a follow up experiment used a similar procedure with the exception that memory was assessed 1 h after stimulation instead of the next day. Given that reconsolidation processes take several hours, finding no memory differences across the conditions after 1 h supported the conclusion that the DLPFC was specifically involved in reconsolidation processes (Sandrini et al., 2013). The use of multiple active control conditions further helped the researchers to isolate the effects to reconsolidation processes (as opposed to initial encoding or retrieval). This study demonstrates the unique feasibility of using rTMS to systematically manipulate (re)consolidation processes. However, because there is so little research that has investigated the effects of rTMS on consolidation processes, the current metaanalysis was restricted to studies that have investigated effects on encoding and retrieval stages. This highlights an area where future TMS research is needed.

Variability in the Quality of Targeting Procedures
Another limitation of the current meta-analysis is the possibility that we may have underestimated the true effect of rTMS on episodic memory performance across conditions. The quality assessment revealed that a large number of studies targeted cortical regions using the international 10-20 electrode positioning system. This approach is imprecise and does not take into account inter-individual differences in the location of specific brain structures. Structurally or functionally guided TMS targeting using a neuronavigation system with each subject's structural and/or functional MRI can more precisely stimulate the targeted brain region of interest. Modeling work that has contrasted the variety of targeting methods has shown that larger TMS effects are obtained with individual fMRI-guided TMS neuronavigation when targeting parietal regions than targeting with the 10-20 EEG positioning system, which yielded smaller effect sizes (Sack et al., 2009). Power analyses revealed that a sample size of 47 was required to detect a significant effect when the 10-20 EEG positioning system is used for TMS targeting (Sack et al., 2009). Given that a large number of studies implemented the 10-20 targeting system and the selected studies had an average sample size of 24.19 in the meta-analysis, it is possible that the existing effects sizes available in the published literature underestimate the true effect of rTMS on episodic memory. Another important aspect to note is that only published effect sizes could be included in the meta-analysis and, therefore, excluding any unpublished null effects would obscure and overestimate the true effect of rTMS on episodic memory. Reassuringly, the steps we took to investigate a potential publication bias suggests that this was unlikely to be of major concern. However, it should be noted that absence of smallstudy biases does not prove the absence of publication bias, and further research is needed to properly address publication bias in multi-level models.

Personalized rTMS
Below we highlight and discuss potentially important factors that need to be systematically investigated over the next 30 years of TMS research and beyond. One pervasive issue in the field has been the highly variable findings across both behavioral and neural outcomes. Certain individuals may respond to a set of rTMS parameters, while other individuals may have null or negative effects to the same rTMS protocol. One possible solution to this heterogeneity in the response to rTMS is to implement more personalized rTMS parameters titrated to match each individual's physiological characteristics that likely underlie the variability between individuals. We have already noted the need to use structural and functional neuronavigated targeting procedures. Other characteristics that are likely to be important concern the titration of the frequency of stimulation to the precise frequency of endogenous oscillations of interest in targeted brain regions.
For example, in a pair of studies, Waldhauser et al. (2016) implemented a personalized rTMS-EEG design to examine if reactivation of sensory information is necessary for recognition of studied items and/or their studied context/source. In experiment one, participants intentionally encoded visual stimuli (i.e., objects) presented laterally (left and right) to a fixation cross; at test, participants made old/new judgments and source judgments about whether the item initially appeared on the left or right side of the screen while undergoing concurrent EEG recording. Identification of a time window (33.5-204.5 ms), frequency (17.5 Hz), and source (MNI coordinates ± 40, −78, 0) of oscillatory behavior for successful source memory judgments were obtained by examining the lateralization of alpha/beta activity during encoding and retrieval. In a second experiment using the same encoding/retrieval procedure they applied rTMS (17.5 Hz) to the occipital lobe (MNI coordinates ± 40, −78, 0) during retrieval from 33.5 to 204.5 ms after cue presentation (Waldhauser et al., 2016). Importantly, when rTMS was applied using the specific temporal, spatial, and frequency parameters related to the task demands for successful episodic retrieval, there was a decrease in source memory performance relative to control.
In a similar study, Hanslmayr et al. (2014) compared online rTMS-EEG at 90% RMT in the beta frequency range (18.7 Hz) to 10.7, 6.8 Hz, and sham stimulation to the left IFG during word encoding. On a subsequent free recall task, participants showed memory impairment following 18.7 Hz stimulation compared to sham. Critically, despite targeting the same task-dependent region, there were no differences in memory performance for 10.7 or 6.8 Hz stimulation compared to sham. At the neural level, the degree of beta entrainment following 18.7 Hz stimulation was correlated with memory impairment, with higher levels of beta entrainment leading to greater decrements in memory performance (Hanslmayr et al., 2014). The findings that artificial beta synchronization leads to memory impairment supports electrophysiology findings that alpha/beta desynchronization is associated with memory enhancement . Thus, the implementation of more nuanced rTMS frequencies affords the possibility to more effectively modulate episodic memory.

Aspects of Episodic Memory
In this review, many of the studies that assessed the effects of rTMS on episodic memory involved spatial locations/associations with stimuli (e.g., object-location associations). While spatial context information is often one important element of episodic memory, there are other important aspects of episodic memories (e.g., object, temporal, social, and emotional information) that may result in dissociable effects of rTMS to different brain regions. Furthermore, it is possible that the type of episodic information (e.g., visual, auditory, multi-modal) may interact with different aspects of episodic memory (e.g., object, temporal, and spatial). Future research should consider these potentially important moderating variables when investigating the effects of rTMS on various aspects of episodic memory.

Closed Loop rTMS
The majority of rTMS findings have focused on researching the brain as a "black box" in an open-loop manner, whereby a priori input stimulation parameters remain constant and one observes the corresponding output. Emerging evidence from multi-modal approaches has begun to clarify how stimulation parameters such as timing and frequency modulate brain-behavior relationships. This has also led to a shift in focus from mapping segregated brain regions to a more interactive-systems approach by attempting to understand how fast-changing patterns of neural activity interact with local and distal brain regions to give rise to human behavior. Understanding the frequency and temporally specific manner that different cortical regions interact has helped uncover how neural processing is linked with behavior. However, the optimal stimulation parameters that lead to episodic memory enhancement or detriment have not been clearly identified and the sparseness of available findings await replication and extension with systematic manipulations to map the parameter space. Moreover, there has been little consideration paid to how memory outcomes may be dependent on interactions between stimulation protocols and temporal patterns of initial brain state activity and ongoing neural activity during a task. The results of the current meta-analysis reveal that it is of paramount importance that future studies take this into consideration in order to determine optimal stimulation parameters for episodic memory (for additional information see, Bergmann et al., 2016) While this "black box" approach has proven fruitful, the brain is a dynamic system capable of receiving, modulating, and generating both further inputs (modified representations, novel associations) as well as outputs (Berkman and Lieberman, 2011). A wave of technological advancements has made it possible to make stimulation adjustments based on realtime neural signals. This has given rise to an informed closed-loop paradigm where neural activity from EEG/MEG provide amplitude, phase, or power of oscillatory activity markers that allow for stimulation in a temporally specific fashion. This has also paved the way for fully adaptive closed-loops that aim to use stimulated neural markers to not only trigger stimulation but also to use neural activity to adaptively alter stimulation protocols in real time .

Multi-Modal rTMS
Given the wealth of knowledge that has been obtained using multimodal approaches, extending these typically two-way combinations (e.g., TMS-EEG and TMS-fMRI) to three-way TMS-EEG-fMRI tools may provide a more holistic view of brain functioning. Such a combination with TMS enables novel stimulation findings based on the temporal precision of EEG to be complemented with fMRI to examine changes in functionally related neural activity in cortical and subcortical structures with high spatial resolution. The feasibility of combining TMS-EEG-fMRI has already been established to be safe and provide reliable signal recordings (Peters et al., 2013).

CONCLUSION
The findings of the meta-analysis on younger adults revealed that frequency interacts with MTI as well as with stimulation timing. Specifically, both online and offline 1 Hz rTMS led to enhancing effects, which was driven by below-MTI stimulation. In addition, offline 20 Hz rTMS had enhancing effects whereas, online 20 Hz rTMS and 20 Hz rTMS at below-MTI led to impairing effects on episodic memory. A systematic review of the older adults and those with clinical disorders revealed a similar pattern of enhancing and attenuating effects of rTMS on episodic memory performance. However, important differences did arise with older adults and clinical populations implementing both offline 5 and 10 Hz rTMS protocols that had enhancing effects on episodic memory.
In sum, these findings highlight the importance of the contextual aspects of stimulation to reveal brain-behavior relationships at a more causal level than permitted by methods that are inherently correlational. Important contextual aspects to consider include the brain state of the targeted region and functionally connected regions at the time of stimulation. Furthermore, researchers must take into consideration the time-dependent processes engaged during specific encoding, consolidation, and retrieval phases that are coded by the brain region of interest, as well as variability both between and within individuals in the precise timing and magnitude of activation of such relevant brain areas. Multi-modal and closed-loop approaches offer promise for addressing these likely pernicious contributors to the variability that is seen across the literature, which may be necessary to parse out the role of exogenous stimulation of brain regions and functionally connected regions for modulating episodic memory. Although still in the initial stages of development, the future of closed-loop paradigms offers the ability to analyze stimulation induced brain activity to dynamically fine-tune stimulation parameters to facilitate rTMS effects on cognition. It is exciting to ponder what the next 30 years will elucidate regarding the use of non-invasive brain stimulation for causal modulation of brain-behavior relationships and, ideally, the enhancement of episodic memory in both healthy and clinical populations.

AUTHOR CONTRIBUTIONS
NY conducted the review and meta-analysis and wrote the initial drafts. NR helped design the study, provide edits, conceptual and theoretical feedback, and write the final draft.