Original Research ARTICLE
Functional and Structural Neuroplasticity Induced by Short-Term Tactile Training Based on Braille Reading
- 1Laboratory of Neuroplasticity, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Warsaw, Poland
- 2CNS Lab, Nalecz Institute of Biocybernetics and Biomedical Engineering, Polish Academy of Sciences, Warsaw, Poland
- 3Bioimaging Research Center, World Hearing Center, The Institute of Physiology and Pathology of Hearing, Warsaw, Poland
- 4Laboratory of Psychophysiology, Nencki Institute of Experimental Biology, Warsaw, Poland
- 5Department of Psychology, University of Social Sciences and Humanities, Warsaw, Poland
- 6Department of Psychology, Jagiellonian University, Cracow, Poland
Neuroplastic changes induced by sensory learning have been recognized within the cortices of specific modalities as well as within higher ordered multimodal areas. The interplay between these areas is not fully understood, particularly in the case of somatosensory learning. Here we examined functional and structural changes induced by short-term tactile training based of Braille reading, a task that requires both significant tactile expertise and mapping of tactile input onto multimodal representations. Subjects with normal vision were trained for 3 weeks to read Braille exclusively by touch and scanned before and after training, while performing a same-different discrimination task on Braille characters and meaningless characters. Functional and diffusion-weighted magnetic resonance imaging sequences were used to assess resulting changes. The strongest training-induced effect was found in the primary somatosensory cortex (SI), where we observed bilateral augmentation in activity accompanied by an increase in fractional anisotropy (FA) within the contralateral SI. Increases of white matter fractional anisotropy were also observed in the secondary somatosensory area (SII) and the thalamus. Outside of somatosensory system, changes in both structure and function were found in i.e., the fusiform gyrus, the medial frontal gyri and the inferior parietal lobule. Our results provide evidence for functional remodeling of the somatosensory pathway and higher ordered multimodal brain areas occurring as a result of short-lasting tactile learning, and add to them a novel picture of extensive white matter plasticity.
Neuroplastic changes accompanying the processes of new skill acquisition or associative learning have been studied in both humans and animals. It is well established that following sensory learning, the cortex of the modality specific to the task shows modifications in electrophysiological properties, neuronal responsiveness, synaptic, and functional connectivity, gray matter volume, and the structure of white matter (for review see: de Villers-Sidani and Merzenich, 2011; Kilgard, 2012; Zatorre et al., 2012; Lövdén et al., 2013; Sur et al., 2013; Takeuchi and Izumi, 2013). In the somatosensory domain, tactile sensory discrimination was reported to alter receptive fields and excitability of neurons of the primary somatosensory cortex, and induce representational plasticity in both humans and animals (Recanzone et al., 1992; Siucinska and Kossut, 1996 Pleger et al., 2003; Gebel et al., 2013; Ladda et al., 2014; Andrew et al., 2015). Rapid recruitment of the primary visual cortex for touch after short-term visual deprivation (Merabet et al., 2008) and alterations in activity of supramodal cortical regions have also been shown (Pleger et al., 2003; Saito et al., 2006; Eckhoff et al., 2008; Groussard et al., 2010; Sathian et al., 2013).
Sensory learning was also found to increase neuronal excitability in animal learning models (Saar and Barkai, 2003; Ohl and Scheich, 2005; Matthews and Disterhoft, 2009; Bekisz et al., 2010) and stronger activation may induce increases in myelination (Sampaio-Baptista et al., 2013; Gibson et al., 2014; McKenzie et al., 2014) or axonal sprouting (Beaulieu, 2002; Boele et al., 2013), which may be reflected by its fractional anisotropy (FA). Moreover, only a few functional activation studies were combined with investigation of white matter integrity (Scholz et al., 2009; Taubert et al., 2010; Loui et al., 2011; Gebauer et al., 2012; Schlegel et al., 2012; Lövdén et al., 2013; Draganski et al., 2014; Chavan et al., 2015) and none of them focused on tactile learning.
To address this issue, we investigated the patterns of functional and structural reorganization induced by tactile learning. The learning task was to read Braille exclusively by touch. Our subjects were individuals with normal vision who had no prior experience with the tactile Braille alphabet. Braille reading consists of several distinct processes including simple finger movement, the perception of series of variously arranged raised dots, pattern recognition, and semantic decoding/lexical processing. Learning to read Braille combines mastering the tactile discrimination task with higher cognitive functions. Since the somatosensory modality plays a crucial role in this process, we expected to find changes along the somatosensory pathway including the primary (SI) and secondary somatosensory (SII) cortices, as well as other areas of the ventral somatosensory pathway engaged in processing of tactile information (Burton and Sinclair, 2000; Pleger et al., 2003; Reed et al., 2004).
While use-dependent plasticity within the primary somatosensory representation is well documented, much less attention has been paid to the secondary somatosensory area. The SII is a higher-order, somatotopically organized cortical area known for its participation in tactile and learning processes, decision making, and susceptibility to the influence of attention (Johansen-Berg et al., 2000; Fujiwara et al., 2002; Romo et al., 2002; Sadato et al., 2002; Pleger et al., 2003). It is also considered to be a multimodal area that integrates information from both sides of the body (Huttunen et al., 1996), but its exact role is still unknown. We previously reported bilateral plasticity of the vibrissae representation within the SII of mice resulting from a short-term sensory training, which was the first demonstration of functional changes to this area induced by associative learning in rodents (Debowska et al., 2011). Unilateral structural (voxel based morphometry, VBM) change of the SII cortex associated with tool-use learning was reported in primates by Quallo et al. (2009). Therefore, the second objective of our study was to understand the relative contributions of the somatosensory cortices to learning by determining whether the secondary somatosensory cortex in humans also undergoes neuroplastic changes induced by sensory training.
Besides somatosensory areas, we expected to observe training-induced changes within multisensory cortical regions known to be involved in tactile recognition: the fusiform gyrus/lateral occipital cortex and associative parietal areas (Amedi et al., 2001, 2002; Sadato et al., 2002; Kassuba et al., 2011; Kim and Zatorre, 2011; Marangon et al., 2016).
In the present study, a group of sighted individuals underwent two scanning sessions (fMRI and DTI) before and after they learned to read Braille exclusively by touch, without any visual deprivation (subjects also kept their eyes open during trainings and MRI sessions).
21 right-handed and sighted subjects (11 women, 10 men, min. age = 22, max. age = 26, mean age = 24.1, SD = 1.4) with no history of neurological or psychiatric disorders participated in the study (Oldfield, 1971). All participants gave their written informed consent prior to the start of the experiment and received financial compensation for completed participation. None of the subjects had prior experience with Braille reading. The study was approved by the Ethics Committee of the Warsaw Medical University.
Experimental Procedure—Braille Training
Training consisted of 15 Braille reading lessons (45 min per day) conducted over the course of 3 weeks. The Braille teaching procedure was planned under the supervision of the Polish Association of the Blind. A Braille 36-chapter ABC-book printed on plastic with standard Braille dimensions (Marburg Medium) was used for the purpose of learning to read Braille (Andraszewski, 2006). As the ABC-book is designed for use by visually-impaired people, it enables us to explore and learn subsequent Braille characters exclusively by touch. The new skill of Braille character recognition was acquired by tactile association with embossed (same-size dots) letters of the Latin alphabet. The first couple of chapters introduced the Braille alphabet with the subsequent ones providing exercises and readings solely in Braille. The six-dot cell (rectangles) used in a further part of the experiment for non-Braille tactile recognition (control) did not appear during the training, so this particular spatial dot-pattern remained untrained and ascribed no particular meaning. The greatest emphasis of the training was put on practice, such that the first lesson included reading a series of up to six-letter words. Subjects were trained to use only their index finger of the right hand while learning, while the left hand was used to guide the right hand position on the text. Participants kept their eyes open during the trainings and the Braille text was concealed by a partition. Additionally, they were asked not to look up the visual Braille displays. Trainings were conducted in small groups (3–4 people) by the same teacher. Great care was taken to teach each subject in the same manner and to a similar skill level. The first nine chapters introduced letters, so during trainings subjects repeatedly returned to them throughout the entire course. After completing the basic lessons, subjects were trained in reading short words and then short sentences that comprised short stories. Since reading skill was balanced among the group, subjects were able to progress up to the 30th chapter over 15 sessions.
The first brain imaging session was conducted 1 day before beginning the trainings and the second session 1 day after the last training.
fMRI—Braille Character Discrimination Task
During the functional imaging sessions, the subjects were asked to actively explore and compare two simultaneously presented Braille characters or non-Braille signs (six-dot rectangles; Figure 1A), or simply to move their finger as though they were touching the presented stimuli, but without touching the matrix (control task). The stimuli were perceived by an active touch of the “trained” finger. A total of 32 pairs of Braille characters (chosen from all of the letters of the alphabet, half of which were the same) and 32 pairs of either parallel or shifted positioned rectangles were presented pseudorandomly alternating with 32 simple movement trials, with a constraint of no more than three consecutive trials of the same type in a row. The functional data were obtained using an MRI-compatible Braille Character Stimulator, fully computer-controlled, phneumatically driven, and capable of delivery up to three standard Braille characters simultaneously (2 × 3, Marburg Medium) with addition of delivery also diacritics and meaningless characters i.e., rectangles (for detailes see: Debowska et al., 2013). The stimulator was placed on the subject's thigh, so the arm laid naturally along the body and the “reading” finger was positioned as when reading regular text. Subjects were requested to look at a fixation cross (“+”) located in the center of a black screen during the whole session. Change of the fixation point color was used as a cue. The sequence of colors in a single experimental trial was as follows: starting with white (basic color), a yellow “+” is presented for 2 s as a preparatory cue. The “+” then turns green and the stimuli are presented for 3.5 s, after which it turns back to white. After a few seconds it turns red to signal time for a response and then reverts back to white. In the “movement” condition trials, conducted in order to identify brain activations associated with finger movement per se, the preparatory cue was presented in blue. Intra- and inter-trial intervals varied across tasks and were administered randomly for 4 to 6 s and 4 to 10 s, respectively. Examples of tactile stimuli and a single trial sequence schema are presented in Figures 1A,B respectively. The delay between the stimulation and the response cues was introduced in order to obtain “pure” activation from unilateral stimulation. Subjects submitted their answers by pushing the corresponding button on a ResponsePad® (SMIT-LAB, www.smit-lab.eu) with their left hand. The software package PRESENTATION® (Neurobehavioural Systems, Albany, CA) was used to present stimuli and register subjects' responses. The stimuli application procedure used for the active tactile discrimination task in this study guaranteed precise (Braille) stimuli size, resolution, and presentation timing.
Figure 1. fMRI procedure. Examples of same-different Braille characters and rectangles used as tactile stimuli during functional run (A). The schema of the visual presentation used as a task-cue for subjects. The change of the fixation cross (“+”) was coupled with tactile (Braille/rectangles) stimuli presentation and submission of answers (B).
All images were acquired on a 3T MRI scanner (Magnetom Trio TIM VB17, Siemens Healthcare, Germany) using a 12-channel head matrix coil. All subjects underwent two identical scanning sessions: before and after the Braille training. Detailed anatomical, functional, and diffusion-weighted images were acquired during both sessions.
fMRI Data Acquisition and Analysis
Anatomical data (T1-weighted) were obtained with 1 mm isotropic resolution (TR = 1900 ms/TE = 2.2 ms/FA = 9 deg./FOV = 192 mm/IPAT = 2/TA = 22:11 min, 47 axial slices perpendicular to AC-PC line). In the functional run a total of 440 volumes (plus 5 dummy scans) were obtained during Echo-Planar Imaging sequences with the following parameters: TR = 3000 ms/TE = 30 ms/voxel: 2 × 2 × 3 mm/FA = 90 deg. SPM12 (http://www.fil.ion.ucl.ac.uk/spm/) was used for data preprocessing and analysis. First, functional images were motion corrected. Then, structural images from single subjects were co-registered to the mean functional image. In the next step, T1 images were warped to the MNI (Montreal Neurological Institute) T1 image template (voxel size of 2 × 2 × 2 mm), and normalization parameters were applied to the functional images, wherein these data were smoothed using a 5-mm isotropic Gaussian kernel. A high-pass filter with a 128 s cutoff was applied to remove low-frequency fluctuations in the BOLD response. For each subject and for each time period separately (before and after training), the onsets and stimuli durations were entered into the design matrix and modeled in a general linear model according to different event types. Regressors were convolved with the canonical hemodynamic response function (Friston et al., 1995). Specific condition effects were assessed by the application of linear contrasts, where parameter estimates for events (i.e., discrimination of Braille characters) were compared to the movement condition. Overall, there were four planned contrasts of interest at the first level of analyses: Braille vs. movement and Braille vs. rectangles each before and after the training. At the second level (whole-brain), paired t-tests were used: Braille > movement after > before and Braille > rectangles after > before. The statistical threshold was set at p < 0.001 and corrected to p < 0.05 for multiple comparisons (Family Wise Error) at the cluster-level using cluster size. Detailed results of statistical analysis with MNI coordinates for peaks, clusters size, and exact values used as an extent thresholds are presented in Table 1. Statistical parametric maps for appropriate comparisons superimposed on the gray matter template (GM tissue probability map) are presented in Figures 2, 3. Percentage of signal change was calculated individually within a spherical volume of interest with a diameter of 10 mm using Marsbar (Brett et al., 2002) and then compared using paired t-test.
Table 1. Regions showing significant increases in activation as a result of the Braille training (Braille discrimination vs. control tasks × session interaction).
Figure 2. fMRI results. Statistical parametric maps of the functional plasticity induced by Braille training revealed by Braille > movement after vs. before training comparison, FWE corrected at p < 0.05 (A); percentage of BOLD signal change (B). Percentage of signal change was calculated individually within a spherical volume of interest with a diameter of 10 mm (*p < 0.05, paired t-test). Abbreviations: AnG—angular gyrus, FuG—fusiform gyrus, Ins—insula, IPL—inferior parietal lobule, MFG—medial frontal gyrus, PCC—posterior cingulate cortex, Pcun—precuneus, SI—primary somatosensory cortex, L—left, R—right.
Figure 3. fMRI results. Statistical parametric maps of functional plasticity induced by Braille training. Results from Braille > rectangles after > before comparisons are presented on the axial slices of the gray matter template (MNI). FWE corrected at p < 0.05. Abbreviations: AnG—angular gyrus, FuG—fusiform gyrus, IFG—inferior frontal gyrus, Ins—insula, LinG—lingual gyrus, MFG—medial frontal gyrus, PCC—posterior cingulate cortex, Pcun—precuneus, SFG—superior frontal gyrus, SI—primary somatosensory cortex.
DTI Protocol and Analysis
Spin-echo diffusion weighted echo planar imaging (DW-EPI) sequence was performed with an isotropic (2 × 2 × 2 mm) resolution, b-value of 1000 s/mm2, 64 gradient directions and repeated twice in order to increase signal-to-noise ratio. The eddy distortion correction and diffusion tensor parameter estimation were performed using ExploreDTI (Leemans et al., 2009). Fractional anisotropies (FA) were calculated for each subject before and after training. SPM12 was used for motion correction, normalization procedures including co-registration to the FA template, smoothing using 3-mm isotropic Gaussian kernel algorithm, and statistical analysis employing paired t-tests (Abe et al., 2010). The statistical threshold was set at p < 0.001 and corrected to p < 0.05 for multiple comparisons (Family Wise Error, FWE) at the cluster-level using cluster size. Detailed results of statistical analysis with MNI coordinates for peaks, clusters size, percentage of FA change, and the exact values used as extent thresholds are presented in Table 2. Statistical parametric maps superimposed on the white matter template (WM tissue probability map) are presented in Figure 5. Anatomical localization of cortical and white matter regions was verified on the basis of the MNI Space Utility and confirmed manually using the Atlas of the Human Brain and the MRI Atlas of Human White Matter (Mai et al., 2007; Mori et al., 2010).
Braille Discrimination (Same-Different) Task
Task performance (percentage of correct responses) differed significantly from the first to the second imaging session. A paired t-test (t = 31.46, df = 20) revealed a significant improvement from 37.5% (min. = 34.3, max. = 37.5, SD = 4.4) to 91.85% (min. = 90.6, max. = 96.8, SD = 6.3) in the Braille character discrimination task performance (p < 0.001).
Rectangles Discrimination Task
Comparisons of non-Braille discrimination task performance before and after a posteriori did not reveal significant differences (t = 2.13, df = 20). In the first session the mean percent of correct responses was 56.9% (min. = 53.1, max. = 68.7, SD = 7.6), while in the second session it reached 61.4% (min. = 46.8, max. = 75, SD = 8.2).
We analyzed brain activation patterns accompanying discrimination of Braille characters performed by active touch, before and after a short-term intensive tactile training of Braille reading in sighted subjects (N = 21). Statistical parametric maps of task-related brain activity patterns before and after the Braille-based tactile training separately are shown in Supplementary Materials. The Braille discrimination vs. finger movement comparison showed a similar pattern of activity across the brain with stronger activations after the training (Figure S1A); comparing Braille vs. rectangles revealed a massive reorganization of cortical activation within the left hemisphere (Figure S1B). Analysis of the interaction between the two imaging sessions (after vs. before) and the experimental conditions revealed a set of brain areas that changed their activation as result of the training (Figures 2, 3). A detailed description is given below and the full set of obtained results along with statistics is given in Table 1.
The Effects of Braille-Based Tactile Training
Braille Characters Discrimination Task > Finger Movement after > before
By contrasting the brain activity pattern involved in the Braille discrimination task and motion activity without touching (movement), we aimed to depict the cortical underpinnings of Braille-related somatosensory information processing. Comparing the Braille character discrimination task-related activity to the simple finger movement condition showed a wide network of cortical regions involved in this specific activity before and after the training. Our analysis indicated several structures within this network presenting more pronounced responses while performing the Braille character discrimination task after the training. Notably, a bilateral increase in activity of the primary somatosensory representation (SI) was observed. Stronger activations were also found in the precuneus and the angular gyrus located within the left hemisphere, and in the bilateral inferior parietal lobule. Other areas of increased activation were located within posterior cingulate cortex, the middle frontal gyrus, and the fusiform gyrus bilaterally, as well as in the right insula (Figure 2A). The strength of the BOLD signal change among structures revealed by this comparison is shown in Figure 2B.
Braille Characters Discrimination Task > Rectangles after > before
The Braille vs. rectangles comparison was constructed to isolate the specific effect of Braille character processing as a higher order cognitive task, by subtracting the activation by the same finger movement but crucially without the letter-related connotations. Results depicted a stronger response bilaterally within the precuneus, the posterior, and anterior cingulate cortex, as well as the middle and inferior frontal gyri. Left-lateralized changes in BOLD signal were found within the fusiform and angular gyrus, as well as in primary somatosensory finger representation within the postcentral gyrus (Figure 3).
ROI Analysis of Parietal Operculum/SII
Since the whole-brain analysis did not reveal functional changes within the secondary somatosensory cortex, we decided to conduct additional ROI analyses. Two ROI's were created, for each left and right SII, using group data from the Braille discrimination task > finger movement contrast from both sessions (Figure 4A) and coordinates from our previous study (Figure 4B) revealed by the Braille discrimination task > rest comparison (Debowska et al., 2013). Parameter estimates for Braille and rectangles conditions were calculated individually within a spherical region-of-interest with a diameter of 6 mm. Interaction effects were tested using random-effects analysis, with results showing no significant differences in SII activity before and after the training in both the Braille or rectangles discrimination tasks (Figure 4). Based on these results, we conclude that there is no evidence for functional plasticity within the SII attributed to short-term Braille learning.
Figure 4. ROI's for parietal operculum/SII. Parameter estimates for all of the experimental conditions before and after the training in the secondary somatosensory cortex bilaterally. Coordinates for creating 6 mm spheres were chosen on the basis of our current (A) and previous study (B). p < 0.05, corrected for multiple comparisons.
Voxel-wise comparison of the FA maps before and after the Braille training revealed a set of the white matter regions in which FA changed as a result of tactile learning. FA was increased among the crucial areas of the somatosensory pathway: the postcentral gyrus (SI), the parietal operculum (SII), and the lateral posterior/pulvinar thalamic nuclei of the left hemisphere (Figure 5A). Other changes were located within frontal and parietal-occipital areas: the right inferior and medial frontal gyri, the superior frontal gyrus, fusiform gyrus, and lingual gyrus bilaterally, and the left parahippocampal gyrus (Figure 5B). The opposite comparisons (before > after) of the FA maps did not reach statistical significance. The full set of obtained results along with the statistics and percentage of FA change is given in Table 2.
Figure 5. FA results. White matter areas showing an increase of fractional anisotropy as a result of the Braille training, within somatosensory (A) and the ventral stream (B) regions. Results revealed by paired t-test between the FA maps before and after the training. Statistics presented on the sagittal (x), coronal (y), and axial (z) slices of the white matter tissue template (MNI). FWE corrected at p < 0.05, k = 26.
Relatively short but intensive tactile training based on Braille reading induced modifications in both primary somatosensory and higher ordered areas. We observed BOLD signal changes likely associated with increased synaptic strength, excitability and expansion in local synaptic connectivity (Saar and Barkai, 2009; Wu and Mel, 2009; Bekisz et al., 2010; Froemke et al., 2013; Jasinska et al., 2013). With DTI, we found increases in the fractional anisotrphy index, likely reflecting strengthening of the intra-cortical connections in response to the demand of increased communication during new-skill acquisition (Blumenfeld-Katzir et al., 2011; Sampaio-Baptista et al., 2004). Changes in fMRI activation were congruent with changes in FA within both somatosensory and inferior occiptiotemporal cortices. The simultaneous use of both methods allows a more precise assessment of task-related plasticity. Its exact impact on relevant brain regions is discussed below.
Tactile Training—The Somatosensory Pathway
SI. Since Braille reading as a tactile activity must engage somatosensory processing, we expected to find alteration in cortical response and white matter structure within the somatosensory pathway. The contrast of Braille character discrimination vs. finger movement revealed areas that responded stronger to swiping fingers over Braille dots relative to swiping fingers without touching the surface, thus, finger movement acted as a motor control. Following training, the strongest activation increase observed for this contrast was found in the SI (Figures 2A,B) bilaterally. The primary somatosensory cortex was previously reported to manifest mainly contralateral responsiveness (Hari et al., 1993; Deuchert et al., 2002; Hlushchuk and Hari, 2006) and several more recent papers have reported its bilateral involvement in somatosensory information processing (Blatow et al., 2007; Tamè et al., 2012; Chung et al., 2014). The current study reports robust bilateral activation of the primary somatosensory representation during a tactile discrimination task before and—even stronger—after the training. We suggest that this may be a result of increased interhemispheric communication between homotopic regions related to the task complexity and top-down attention (Verstynen et al., 2005; Perez and Cohen, 2008).
In parallel with change in the magnitude of the BOLD signal, we found FA increases within the postcentral gyrus in the region corresponding to the functional localization of the contralateral finger representation within the SI. The increase of activation and in white matter integrity within the SI (Figure 5A) provides strong evidence for the SI's specific engagement in Braille reading skill acquisition. These results are in line with previous findings concerning plasticity of the primary somatosensory representation induced by both short and long-term tactile training (Pascual-Leone and Torres, 1993; Pleger et al., 2003; Hodzic et al., 2004). The increase of white matter integrity around the altered functional SI representation documents structural plasticity as a novel aspect of learning-induced SI modification.
Interestingly, our results are contrary to those recently presented by Sathian et al. (2013), where no effect of tactile learning was found within the SI. We propose that this might be caused by key differences in the experimental procedure—our subjects were trained to actively discriminate different spatial patterns with assigned meaning (Braille characters), subjects in Sathian's study were trained in pure perceptual learning (micro-spatial task based on a linear three-dot array), and the touch was passive. There are also differences in training duration (fixed vs. individual performance criterion), stimuli presentation protocols (active in our study vs. passive) and duration (3 s per pair in event-related design vs. 1 s per single stimulus presented in block-design).
Our results are also divergent from those of Siuda-Krzywicka et al. (2016) who trained sighted subjects in whole-word tactile reading for 9 months. After such a long period of training, Siuda-Krzywicka et al. (2016) did not find changes in the somatosensory cortex. Instead, their results showed an increase in the left fusiform gyrus activation (Visual Word Form Area, VWFA), where most likely this divergence is due to the very long learning period (9 months vs. 3 weeks) and the complexity of stimuli used (whole words vs. pairs of Braille characters). Experiments that study learning-related plasticity at multiple time points (Lövdén et al., 2013) suggest that at the initial stage of Braille learning described by Siuda-Krzywicka et al. (2016), the somatosensory cortex might have increased its response to Braille words, and as the effects of early sensory learning consolidated in the somatosensory cortex, the cortical focus of learning might have shifted elsewhere—to the ventral visual stream.
We found no functional changes in either ipsi- or contralateral SII, but a highly significant FA increase within white matter in the parietal operculum (anatomical localization of the SII) contralateral to the reading finger (Eickhoff et al., 2010). The changes in white matter structure clearly suggest functional involvement of the contralateral SII in acquisition of the new tactile ability, whereas the absence of increased functional activation after 3 weeks of training raises questions regarding the precise character of this involvement. The importance of the SII in tactile processing, including Braille characters' discrimination in sighted individuals (Sadato et al., 2002), is well documented and we see its strong activation during the Braille character discrimination task used in the present study (Figure 4). We propose that lack of the increase in activation of the SII following training might be due to the rigid timeframe used for training. In other words, we speculate that if checked at different (shorter) time points, e.g., a few days after onset of training, the SII might demonstrate its plasticity in terms of functional involvement—stronger at the beginning of training and coming back to normal when the skill became easier (Pleger et al., 2003; Vahdat et al., 2014) and when the white matter tracts are already remodeled. An alternative interpretation is that SII is not involved into neuroplasticity indcued by training, but only play its usual role in ‘transmitting’ information to brain regions involved in stimuli processing.
Additionally, the lack of long-term functional change in SII while increased activity was observed in the SI cortex might inform an ongoing discussion about hierarchical vs. parallel information processing between the SI and SII (Chung et al., 2014; Klinger et al., 2015). While early animal studies suggested strict dependencies between the SI and SII cortices by demonstrating the extinction of the SII activity after lesioned within SII (Pons et al., 1988), it was later proposed that SII may process sensory information independently (Zhang et al., 2001). In humans, effective connectivity studies showed information flow relating to texture, from the SI to the parietal operculum (Sathian et al., 2011). Here, if increased activity in the SI does not alter the activity within the SII, it is most likely that the SII reacts independently and is not a simple reflection of the SI. Another possible explanation is that although the SII is necessary in tactile stimuli processing, its role in the well-trained tactile pattern discrimination task is no longer needed because of the altered global brain pattern of activity and stronger engagement of regions specified in the object recognition—FuG (Figure 2; Table 1). The FA increase observed within areas of the parietal operculum might then reflect an enhanced intra-cortical communication occurring while mastering the new tactile recognition task.
Outside of the Somatosensory Domain: Fronto-Parietal Cortex and the Ventral Stream
As Braille characters become meaningful to the subjects as a result of training, it should be possible to determine loci where this meaning is extracted. Although the task in the scanner required only a same/different answer, activation of the brain was strongly affected by the meaning of the discriminanda. When we computed the interaction of Braille > rectangles with the difference between the before- and after- training scans, we found activations in the left fusiform cortex, medial and inferior frontal gyri, and the angular gyrus (Figure 3, Table 1). Increased activation of those areas after training might be language related, since they were found at locations typical to those found in readers of a variety of different scripts (Perfetti et al., 2010; Szwed et al., 2014). Increased activity within ventral stream activation suggest that even short-term tactile Braille reading training may result in activation of brain regions involved in categories recognition (Tyler et al., 2013), which might be interpreted as a having semantic conotations. Alternatively, they might be task-related. Unfortunately, we do not have more appropriate behavioral nor imaging data to discriminate between those two.
The increase in activity within the fusiform gyrus (Figure 3, Table 1) is congruent with previously reported involvement of the lateral occipital cortex/fusiform gyrus in tactile object recognition in both blind and sighted subjects (Sadato et al., 2002; Amedi et al., 2003; Striem-Amit et al., 2012) or multisensory object perception (Kassuba et al., 2011, 2014; Schlaffke et al., 2015). Alternatively, it could reflect the emergence of tactile Braille recognition mechanisms similar to those found in the VWFA (Price and Devlin, 2001; Cohen et al., 2002; Dehaene and Cohen, 2011), recently found in sighted whole-word tactile Braille readers (Siuda-Krzywicka et al., 2016). Interestingly, these fMRI results were accompanied in our study by FA increase in the lingual gyrus bilaterally and within the fusiform gyri (Figure 2, Table 2). Our results might thus constitute the first report of white matter integrity changes within higher-ordered and visual areas induced by short-term tactile training in sighted subjects. Apparently even short-term tactile Braille reading training may result in activation of brain regions involved in language in general and specifically—in reading.
Finally, some of our results (FuG, AnG, PCC, Pcun) have much in common with those observed after Morse code deciphering training using auditory stimuli (Schmidt-Wilcke et al., 2010). Our results not only support these findings, but also complement them with a novel aspect, showing significant white matter changes throughout the network involved in acquisition and usage of alphabet signs perceived by a different (somatosensory) sensory modality.
Short-term tactile training based on Braille reading increased SI activation, confirming the presence of learning-induced plasticity in the primary somatosensory cortex. No functional effects were found within secondary somatosensory cortex, however FA increase within the parietal operculum suggest that the SII was involved in acquisition of this new sensory skill, and an fMRI measurement at a third time point might have captured a change of BOLD activity also in the SII. Increased activation of the fusiform gyrus followed by structural alteration in white matter of occipital regions demonstrates that a relatively short tactile training can induce plastic changes in the ventral visual cortex. In summary, our study shows how acquisition and mastering tactile discrimination tasks induce modifications in the pattern and intensity of brain activation and white matter integrity of somatosensory areas, higher-ordered pathways and in the thalamus.
Our training has been arbitrary set for a fixed period of time. It would be extremely interesting and informative to conduct longitudinal study with multiple scanning time points (week by week or even day by day) as well as follow up study to determine how long the observed changes persists. More detailed imaging sequences such as spectroscopy might also brought better insight in to the neurobiological underpinnings of the plastic changes. Moreover, since whole-brain analyses of diffusion indices are of several limitations including Partial Volume Effects and EPI distortions, and requires perfect co-registration (Beaulieu, 2002), other analyses such as tractography combined with i.e., computational anatomy (Draganski et al., 2014) would shed some more light on the exact localization of the reorganized tissue and allow for a better insight into nature of neuroplasticity.
Study conception and design: WD, AN, MS, MK. Acquisition of data: WD, TW, AK. Analysis and interpretation of data: WD, TW, MK. Drafting of manuscript: WD, AN, MK. Critical revision: WD, MS, AN, MK.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Dr. Malgorzata Paplinska (The Maria Grzegorzewska Academy for Special Education in Warsaw) for assistance in setting the Braille training procedure. The study was supported by grant no. 3608 to MK, and no. IP97 to WD by the Polish Ministry of Science and Higher Education, the ETIUDA doctoral scholarship no. 2013/08/T/NZ4/00741 from the National Science Center to WD and by the Foundation for Polish Science ‘Homing Plus’ grant to MS.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnins.2016.00460
Figure S1. Braille training—overall fMRI results. Statistical parametric maps of the Braille vs. movement (A) and Braille vs. rectangles (B) contrasts, separately for before and after the Braille training. Results are presented on a glass-brain in MNI space (view from the top, back, and left of the brain), p < 0.05, FWE.
Amedi, A., Jacobson, G., Hendler, T., Malach, R., and Zohary, E. (2002). Convergence of visual and tactile shape processing in the human lateral occipital complex. Cereb. Cortex 12, 1202–1212. doi: 10.1093/cercor/12.11.1202
Amedi, A., Raz, N., Pianka, P., Malach, R., and Zohary, E. (2003). Early ‘visual’ cortex activation correlates with superior verbal memory performance in the blind. Nat. Neurosci. 6, 758–766. doi: 10.1038/nn1072
Andrew, D., Haavik, H., Dancey, E., Yielder, P., and Murphy, B. (2015). Somatosensory evoked potentials show plastic changes following a novel motor training task with the thumb. Clin. Neurophysiol. 126, 575–580. doi: 10.1016/j.clinph.2014.05.020
Bekisz, M., Garkun, Y., Wabno, J., Hess, G., and Wrobel, A. (2010). Increased excitability of cortical neurons induced by associative learning: an ex vivo study. Eur. J. Neurosci. 32, 1715–1725 doi: 10.1111/j.1460-9568.2010.07453.x
Blatow, M., Nennig, E., Durst, A., Sartor, K., and Stippich, C. (2007). fMRI reflects functional connectivity of human somatosensory cortex. Neuroimage 37, 927–936. doi: 10.1016/j.neuroimage.2007.05.038
Blumenfeld-Katzir, T., Pasternak, O., Dagan, M., and Assaf, Y. (2011). Diffusion MRI of structural brain plasticity induced by a learning and memory task. PLoS ONE 6:e20678. doi: 10.1371/journal.pone.0020678
Boele, H. J., Koekkoek, S. K., De Zeeuw, C. I., and Ruigrok, T. J. (2013). Axonal sprouting and formation of terminals in the adult cerebellum during associative motor learning. J. Neurosci. 33, 17897–17907. doi: 10.1523/JNEUROSCI.0511-13.2013
Brett, M., Anton, J. L., Valabregue, R., and Poline, J. B. (2002). “Region of interest analysis using an SPM toolbox [abstract],” in 8th International Conference on Functional Mapping of the Human Brain, Available on CD-ROM in Neuroimage.
Burton, H., and Sinclair, R. J. (2000). Tactile-spatial and cross-modal attention effects in the primary somatosensory cortical areas 3b and 1-2 of rhesus monkeys. Somatosens Mot. Res. 17, 213–228. doi: 10.1080/08990220050117574
Chavan, C., Mouthon, M., Simonet, M., Hoogewoud, H. M., and Draganski, B. (2015). Sustained enhancements in inhibitory control depend primarily on the reinforcement of fronto-basal anatomical connectivity. Brain Struct. Funct. [Epub ahead of print]. doi: 10.1007/s00429-015-1156-y
Chung, Y. G., Han, S. W., Kim, H. S., Chung, S. C., Park, J. Y., Wallraven, C., et al. (2014). Intra- and inter-hemispheric effective connectivity in the human somatosensory cortex during pressure stimulation. BMC Neurosci. 15:43. doi: 10.1186/1471-2202-15-43
Cohen, L., Lehéricy, S., Chochon, F., Lemer, C., Rivaud, S., and Dehaene, S. (2002). Language-specific tuning of visual cortex? Functional properties of the visual word form area. Brain 125, 1054–1069. doi: 10.1093/brain/awf094
Debowska, W., Liguz-Lecznar, M., and Kossut, M. (2011). Bilateral plasticity of Vibrissae SII representation induced by classical conditioning in mice. J. Neurosci. 31, 5447–5453. doi: 10.1523/JNEUROSCI.5989-10.2011
Debowska, W., Wolak, T., Soluch, P., Orzechowski, M., and Kossut, M. (2013). Design and evaluation of and innovative MRI-compatible Braille stimulator with high spatial and temporal resolution. J. Neurosci. Meth. 213, 32–38. doi: 10.1016/j.jneumeth.2012.12.002
Deuchert, M., Ruben, J., Schwiemann, J., Meyer, R., Thees, S., Krause, T., et al. (2002). Event-related fMRI of the somatosensory system using electrical finger stimulation. Neuroreport 13, 365–369. doi: 10.1097/00001756-200203040-00023
de Villers-Sidani, E., and Merzenich, M. M. (2011). Lifelong plasticity in the rat auditory cortex: basic mechanisms and role of sensory experience. Prog. Brain Res. 191, 119–131. doi: 10.1016/B978-0-444-53752-2.00009-6
Eckhoff, P., Holmes, P., Law, C., Connolly, P. M., and Gold, J. I. (2008). On diffusion processes with variable drift rates as models for decision making during learning. New J. Phys. 10, 015006–015027. doi: 10.1088/1367-2630/10/1/015006
Eickhoff, S. B., Jbabdi, S., Caspers, S., Laird, A. R., Fox, P. T., Zilles, K., et al. (2010). Anatomical and functional connectivity of cytoarchitectonic areas within the human parietal operculum. J. Neurosci. 30, 6409–6421. doi: 10.1523/JNEUROSCI.5664-09.2010
Friston, K. J., Holmes, A. P., Poline, J. B., Grasby, P. J., Williams, S. C., Frackowiak, R. S., et al. (1995). Analysis of fMRI time-series revisited. Neuroimage 2, 45–53. doi: 10.1006/nimg.1995.1007
Froemke, R. C., Carcea, I., Barker, A. J., Yuan, K., Seybold, B. A., Martins, A. R., et al. (2013). Long-term modification of cortical synapses improves sensory perception. Nat. Neurosci. 16, 79–88. doi: 10.1038/nn.3274
Fujiwara, N., Imai, M., Nagamine, T., Mima, T., Oga, T., Takeshita, K., et al. (2002). Second somatosensory area (SII) plays a significant role in selective somatosensory attention. Brain Res. Cogn. Brain Res. 14, 389–397. doi: 10.1016/S0926-6410(02)00141-6
Gebauer, D., Fink, A., Filippini, N., Johansen-Berg, H., Reishofer, G., et al. (2012). Differences in integrity of white matter and changes with training in spelling impaired children: a diffusion tensor imaging study. Brain Struct. Funct. 217, 747–760. doi: 10.1007/s00429-011-0371-4
Gebel, B., Braun, C. H., Kaza, E., Altenmüller, E., and Lotze, M. (2013). Instrument specific brain activation in sensorimotor and auditory representation in musicians. Neuroimage 74, 37–44. doi: 10.1016/j.neuroimage.2013.02.021
Gibson, E. M., Purger, D., Mount, C. W., Goldstein, A. K., Lin, G. L., et al. (2014). Neuronal activity promotes oligodendrogenesis and adaptive myelination in the mammalian brain. Science 344:1252304. doi: 10.1126/science.1252304
Groussard, M., La Joie, R., Rauchs, G., Landeau, B., Chételat, G., Viader, F., et al. (2010). When music and long- term memory interact: effects of musical expertise on functional and structural plasticity in the hippocampus. PLoS ONE 5:e13225. doi: 10.1371/journal.pone.0013225
Hari, R., Karhu, J., Hämäläinen, M., Knuutila, J., Salonen, O., Sams, M., et al. (1993). Functional organization of the human first and second somatosensory cortices: a neuromagnetic study. Eur. J. Neurosci. 5, 724–734. doi: 10.1111/j.1460-9568.1993.tb00536.x
Hlushchuk, Y., and Hari, R. (2006). Transient suppression of ipsilateral primary somatosensory cortex during tactile finger stimulation. J. Neurosci. 26, 5819–5824. doi: 10.1523/JNEUROSCI.5536-05.2006
Hodzic, A., Veit, R., Karim, A. A., Erb, M., and Godde, B. (2004). Improvement and decline in tactile discrimination behavior after cortical plasticity induced by passive tactile coactivation. J. Neurosci. 24, 442–446. doi: 10.1523/JNEUROSCI.3731-03.2004
Huttunen, J., Wikström, H., Korvenoja, A., Seppäläinen, A. M., Aronen, H., and Ilmoniemi, R. J. (1996). Significance of the second somatosensory cortex in sensorimotor integration: enhancement of sensory responses during finger movements. Neuroreport 7, 1009–1012. doi: 10.1097/00001756-199604100-00011
Jasinska, M., Siucinska, E., Jasek, E., Litwin, J. A., Pyza, E., and Kossut, M. (2013). Fear learning increases the number of polyribosomes associated with excitatory and inhibitory synapses in the barrel cortex. PLoS ONE 8:e54301. doi: 10.1371/journal.pone.0054301
Johansen-Berg, H., Christensen, V., Woolrich, M., and Matthews, P. M. (2000). Attention to touch modulates activity in both primary and secondary somatosensory areas. Neuroreport 11, 1237–1241. doi: 10.1097/00001756-200004270-00019
Kassuba, T., Klinge, C., Hölig, C., Menz, M. M., Ptito, M., Röder, B., et al. (2011). The left fusiform gyrus hosts representations of manipulable objects. Neuroimage 56, 1566–1577. doi: 10.1016/j.neuroimage.2011.02.032
Klinger, C. M., Brodoehl, S., Huonker, R., Götz, T., Baumann, L., and Witte, O. W. (2015). Parallel processing of somatosensory information: evidence from dynamic causal modeling of MEG data. Neuroimage 118, 193–198. doi: 10.1016/j.neuroimage.2015.06.028
Ladda, A. M., Pfannmoeller, J. P., Kalisch, T., Roschka, S., Platz, T., Dinse, H. R., et al. (2014). Effects of combining 2 weeks of passive sensory stimulation with active hand motor training in healthy adults. PLoS ONE 9:e84402. doi: 10.1371/journal.pone.0084402
Leemans, A., Jeurissen, B., Sijbers, J., and Jones, D. K. (2009). “ExploreDTI: A graphical toolbox for processing, analyzing, and visualizing diffusion MR data,” in Proceedings of the 17th Scientific Meeting, International Society for Magnetic Resonance in Medicine (Honolulu, HI), 3537.
Lövdén, M., Wenger, E., Mårtensson, J., Lindenberger, U., and Bäckman, L. (2013). Structural brain plasticity in adult learning and development. Neurosci. Biobehav. Rev. 37, 2296–2310. doi: 10.1016/j.neubiorev.2013.02.014
Marangon, M., Kubiak, A., and Króliczak, G. (2016). Haptically guided grasping. fMRI shows right-hemisphere parietal stimulus encoding, and bilateral dorso-ventral parietal gradients of object- and action-related processing during grasp execution. Front. Hum. Neurosci. 9:691. doi: 10.3389/fnhum.2015.00691
Merabet, L. B., Hamilton, R., Schlaug, G., Swisher, J. D., Kiriakopoulos, E. T., Pitskel, N. B., et al. (2008). Rapid and reversible recruitment of early visual cortex for touch. PLoS ONE 3:3046. doi: 10.1371/journal.pone.0003046
Perez, M. A., and Cohen, L. G. (2008). Mechanisms underlying functional changes in the primary motor cortex ipsilateral to an active hand. J. Neurosci. 28, 5631–5640. doi: 10.1523/JNEUROSCI.0093-08.2008
Perfetti, C. A., Nelson, J., Liu, Y., Fiez, J., and Ta, L.-H. (2010). “The neural bases of reading: Universals and writing system variations,” in The Neural Basis of Reading, eds P. Cornelissen, M. Kringelbach, and P. Hansen (Oxford: Oxford University Press), 147–172.
Pleger, B., Foerster, A. F., Ragert, P., Dinse, H. R., Schwenkreis, P., Malin, J. P., et al. (2003). Functional imaging of perceptual learning in human primary and secondary somatosensory cortex. Neuron 40, 643–653. doi: 10.1016/S0896-6273(03)00677-9
Pons, T. P., Garraghty, P. E., and Mishkin, M. (1988). Lesion-induced plasticity in the second somatosensory cortex of adult macaques. Proc. Natl. Acad. Sci. U.S.A. 85, 5279–5281. doi: 10.1073/pnas.85.14.5279
Quallo, M. M., Price, C. J., Ueno, K., Asamizuya, T., and Cheng, K. (2009). Gray and white matter changes associated with tool-use learning in macaque monkeys. Proc. Natl. Acad. Sci. U.S.A. 106, 18379–18384. doi: 10.1073/pnas.0909751106
Recanzone, G. H., Merzenich, M. M., Jenkins, W. M., Grajski, K. A., and Dinse, H. R. (1992). Topographic reorganization of the hand representation in cortical area 3b owl monkeys trained in a frequency-discrimination task. J Neurophysiol. 67, 1031–1056.
Saar, D., and Barkai, E. (2009). Long-lasting maintenance of learning-induced enhanced neuronal excitability: mechanisms and functional significance. Mol. Neurobiol. 39, 171–177. doi: 10.1007/s12035-009-8060-5
Saito, D. N., Okada, T., Honda, M., Yonekura, Y., and Sadato, N. (2006). Practice makes perfect: the neural substrates of tactile discrimination by Mah-Jong experts include the primary visual cortex. BMC Neurosci. 7:79. doi: 10.1186/1471-2202-7-79
Sampaio-Baptista, C., Khrapitchev, A. A., Foxley, S., Schlagheck, T., Scholz, J., Jbabdi, S., et al. (2013). Motor skill learning induces changes in white matter microstructure and myelination. J. Neurosci. 33, 19499–19503. doi: 10.1523/JNEUROSCI.3048-13.2013
Sampaio-Baptista, C., Scholz, J., Jenkinson, M., Thomas, A. G., Filippini, N., Smit, G., et al. (2004). Gray matter volume is associated with rate of subsequent skill learning after a long term training intervention. Neuroimage 96, 158–166. doi: 10.1016/j.neuroimage.2014.03.056
Sathian, K., Deshpande, G., and Stilla, R. (2013). Neural changes with tactile learning reflect decision-level reweighting of perceptual readout. J. Neurosci. 33, 5387–5398. doi: 10.1523/JNEUROSCI.3482-12.2013
Sathian, K., Lacey, S., Stilla, R., Gibson, G. O., Deshpande, G., Hu, X., et al. (2011). Dual pathways for haptic and visual perception of spatial and texture information. Neuroimage 57, 462–475. doi: 10.1016/j.neuroimage.2011.05.001
Schlaffke, L., Rüther, N. N., Heba, S., Haag, L. M., Schultz, T., Rosengarth, K., et al. (2015). From perceptual to lexico-semantic analysis–cortical plasticity enabling new levels of processing. Hum. Brain Mapp. 36, 4512–4528. doi: 10.1002/hbm.22939
Schmidt-Wilcke, T., Rosengarth, K., Luerding, R., Bogdahn, U., and Greenlee, M. W. (2010). Distinct patterns of funcional and structural neuroplasticity associated with learning Morse code. Neuroimage 51, 1234–1241. doi: 10.1016/j.neuroimage.2010.03.042
Siucinska, E., and Kossut, M. (1996). Short-lasting classical conditioning induces reversible changes of representational maps of vibrissae in mouse SI cortex—a 2DG study. Cereb. Cortex 6, 506–513. doi: 10.1093/cercor/6.3.506
Siuda-Krzywicka, K., Bola, Ł., Paplińska, M., Sumera, E., Jednoróg, K., Marchewka, A., et al. (2016). Massive cortical reorganization in sighted Braille readers. Elife 5:e10762. doi: 10.7554/eLife.10762
Striem-Amit, E., Cohen, L., Dehaene, S., and Amedi, A. (2012). Reading with sounds: sensory substitution selectively activates the visual word form area in the blind. Neuron 76, 640–652. doi: 10.1016/j.neuron.2012.08.026
Szwed, M., Qiao, E., Jobert, A., Dehaene, S., and Cohen, L. (2014). Effects of literacy in early visual and occipitotemporal areas of Chinese and French readers. J. Cogn. Neurosci. 26, 459–475. doi: 10.1162/jocn_a_00499
Tamè, L., Braun, C., Lingnau, A., Schwarzbach, J., Demarchi, G., Li Hegner, Y., et al. (2012). The contribution of primary and secondary somatosensory cortices to the representation of body parts and body sides: an fMRI adaptation study. J. Cogn. Neurosci. 24, 2306–2320. doi: 10.1162/jocn_a_00272
Taubert, M., Draganski, B., Anwander, A., Müller, K., Horstmann, A., Villringer, A., et al. (2010). Dynamic properties of human brain structure: learning-related changes in cortical areas and associated fiber connections. J. Neurosci. 30, 11670–11677. doi: 10.1523/JNEUROSCI.2567-10.2010
Tyler, L. K., Chiu, S., Zhuang, J., Randall, B., Devereux, B. J., Wright, P., et al. (2013). Objects and categories: feature statistics and object processing in the ventral stream. J. Cogn. Neurosci. 25, 1723–1735. doi: 10.1162/jocn_a_00419
Vahdat, S., Darainy, M., and Ostry, D. J. (2014). Structure of plasticity in human sensory and motor networks due to perceptual learning. J. Neurosci. 34, 2451–2463. doi: 10.1523/JNEUROSCI.4291-13.2014
Verstynen, T., Diedrichsen, J., Albert, N., Aparicio, P., and Ivry, R. B. (2005). Ipsilateral motor cortex activity during unimanual hand movements relates to task complexity. J. Neurophysiol. 93, 1209–1222. doi: 10.1152/jn.00720.2004
Keywords: neuroplasticity, primary somatosensory cortex, secondary somatosensory cortex, sighted, Braille, fractional anisotropy
Citation: Debowska W, Wolak T, Nowicka A, Kozak A, Szwed M and Kossut M (2016) Functional and Structural Neuroplasticity Induced by Short-Term Tactile Training Based on Braille Reading. Front. Neurosci. 10:460. doi: 10.3389/fnins.2016.00460
Received: 30 April 2016; Accepted: 26 September 2016;
Published: 13 October 2016.
Edited by:Merim Bilalić, Alpen-Adria-Universität Klagenfurt, Austria
Reviewed by:Hidenao Fukuyama, Kyoto University, Japan
Tobias Schmidt-Wilcke, Bergmannsheil and Ruhr-Universität Bochum, Germany
Copyright © 2016 Debowska, Wolak, Nowicka, Kozak, Szwed and Kossut. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Malgorzata Kossut, email@example.com