- Open Access
Longitudinal diffusion tensor magnetic resonance imaging analysis at the cohort level reveals disturbed cortical and callosal microstructure with spared corticospinal tract in the TDP-43G298S ALS mouse model
Translational Neurodegenerationvolume 8, Article number: 27 (2019)
In vivo diffusion tensor imaging (DTI) of the mouse brain was used to identify TDP-43 associated alterations in a mouse model for amyotrophic lateral sclerosis (ALS).
Ten mice with TDP-43G298S overexpression under control of the Thy1.2 promoter and 10 wild type (wt) underwent longitudinal DTI scans at 11.7 T, including one baseline and one follow-up scan with an interval of about 5 months. Whole brain-based spatial statistics (WBSS) of DTI-based parameter maps was used to identify longitudinal alterations of TDP-43G298S mice compared to wt at the cohort level. Results were supplemented by tractwise fractional anisotropy statistics (TFAS) and histological evaluation of motor cortex for signs of neuronal loss.
Alterations at the cohort level in TDP-43G298S mice were observed cross-sectionally and longitudinally in motor areas M1/M2 and in transcallosal fibers but not in the corticospinal tract. Neuronal loss in layer V of motor cortex was detected in TDP-43G298S at the later (but not at the earlier) timepoint compared to wt.
DTI mapping of TDP-43G298S mice demonstrated progression in motor areas M1/M2. WBSS and TFAS are useful techniques to localize TDP-43G298S associated alterations over time in this ALS mouse model, as a biological marker.
Diffusion tensor imaging (DTI) has become an important tool to study the anatomy of the mouse brain in vivo [1, 2, 13, 30], and its non-invasive nature enables longitudinal studies of transgenic murine disease models [12, 34]. Ultra-high fields (11.7 T) and dedicated resonators (cryogenic cooled resonator-CCR) allow the recording of high-resolution DTI datasets with in-plane resolutions down to 156 μm × 250 μm with an axial slice thickness of 250 μm  while the development of fast DTI protocols has reduced acquisition time to 35 min, enabling the logistics for the monitoring of larger cohorts . The combination of these approaches enables the study of the large-scale structural integrity of white matter and axonal tracts as well as of the microstructure of brain structures over time . Moreover, DTI allows the exploration of (micro) anatomical integrity and changes in a brain-wide, unbiased and operator-independent manner.
The large-scale structure of myelinated, long-range axonal projections has recently emerged as imaging biomarker of disease progression in Amyotrophic Lateral Sclerosis (ALS): DTI investigation of the integrity of cortico-efferent axonal tracts has revealed the progressive involvement of cortico-spinal, cortico-rubral and cortico-striatal subpopulations of projection neurons in human patients , in agreement with the pattern of ALS-related pathobiochemistry (i.e., tar DNA binding protein 43 (TDP-43) inclusions) detected in neuropathological studies . Moreover, DTI has been instrumental in comparing the patterns of anatomical involvement in ALS clinical variants .
Nevertheless, much less is known about the large-scale involvement of axonal projections in ALS murine models and in particular about the anatomical patterns and the overall degree of structural involvement in the brain. A comparative study of retrograde-AAV projection tracing in the SOD1G93A mouse model has revealed substantial similarities between the connectivity remodeling in human patients and murine models . Likewise, the loss of axons in the spinal cord has been detected in vivo by DTI approaches . Very little is known about non-SOD1 models, which may recapitulate different aspects of the ALS/FTD disease spectrum. In fact, mutant TDP-43 transgenic mice develop axonal loss and progressive cognitive phenotypes [3, 21], together with biochemical signs of TDP-43 aggregation . However, the extent of brain involvement in TDP-43G298S transgenic mice has not been fully investigated because histological techniques cannot be easily implemented on every brain structure and may suffer from significant biases in quantification. Thus, the extent of involvement of the motor cortex as well as any pathological change in extra-motor areas and their progression over time remain unclear. We focused on the TDP-43G298S transgenic mouse model because, in contrast to the mutant SOD1 model, it may more faithfully mimic the pathogenic cascades and the neuropathological hallmarks observed in human patients affected by sporadic ALS, frontotemporal dementia (FTD), and related conditions [3, 4, 9]. Here, whole brain-based spatial statistics (WBSS) of DTI data  was applied to the TDP-43G298S transgenic mouse model of ALS in order to establish the large-scale pattern of abnormalities in the cerebral microstructure in a longitudinal experimental design.
Materials and methods
Mice expressing full-length mutant (G298S) human TARDBP (hereafter referred to as TDP-43G298S mice) were provided by Phil Wong (line available at the Jackson Laboratory, Stock No: 017589, Thy1.2-TDP43*G298S line S97). Mice were maintained at 22 °C with a 12/12 h light/dark cycle and had food and water ad libitum. All animal experiments were performed in accordance with institutional guidelines of Ulm University and were approved by the according regulatory authority (Regierungspräsidium Tübingen, Germany; animal permission no. 1242). Twice a week, male mice were subjected to weighing and disease scoring. Ten adult wild type (wt) mice and ten heterozygous TDP-43G298S mice (mean age 9.3 months at baseline scan) underwent the whole-brain DTI-MRI protocol. Heterozygous TDP-43G298S mice develop a robust motor neuron disease phenotype with tremor and amyotrophic paralysis of the limbs progressing to flaccid paralysis at around 2 years of age (data not shown; [39, 43]).
Data acquisition was performed under isoflurane anesthesia (5% for induction and 1.5% for maintenance). The animals were placed in a stereotactic head support (Bruker Biospin, Ettlingen, Germany) to immobilize the head. Body temperature was controlled by an integrated water-based heating device. The body temperature of the mouse was monitored by a rectal temperature probe and respiration was monitored by a respiratory pillow positioned under the abdomen of the mouse. The breathing frequency was maintained at 75–80 cycles per minute. The mice rapidly recovered (< 5 min) after the termination of anesthesia at the end of the MRI procedure.
Imaging was performed with an 11.7 T small bore animal scanner (Biospec 117/16, Bruker, Ettlingen, Germany). A two-element transmit/receive 1H mouse cryogenic surface coil (Cryo-Probe, Bruker BioSpin) was used for data acquisition. Imaging parameters of the optimized rapid diffusion prepared spin echo EPI imaging protocol were as: TE/TR 36.0 ms / 5000 ms, matrix 180 × 40, in-plane resolution 102 μm × 102 μm, 70 slices with a slice thickness of 250 μm. Thirty diffusion directions with b = 1000 s/mm2 and 5 unweighted b = 0 volumes (standard gradient scheme as provided by the Bruker software), 1 signal average, were acquired, resulting in a total acquisition time of 22 min.
The sequence was respiratory gated on a slice by slice basis. The mice underwent two scans, one baseline and one follow-up scan after a time interval of approximately 5 months. If possible, an additional scan after a further follow-up interval of 5 months was acquired. These additional scans were used to replace the follow-up scans in order to select the optimum data quality for baseline and follow-up scan. One wt mouse and one TDP-43G298S mouse had only one scan useful for analysis, the remaining scans of these mice did not pass quality control due to motion artefacts. Thus, in summary 9 TDP-43G298S mice and 9 wt mice, each with a baseline and one follow-up scan were used for analysis.
Data processing was performed with the Tensor Imaging and Fiber Tracking (TIFT) software package  which has been successfully applied to animal DTI group studies (e.g. [26, 27]). A schematic overview of the analysis cascade is provided in Fig. 1.
Recorded data were transformed into a 50 μm isogrid (nearest neighbor interpolation) in order to minimize partial volume effects. The slice thickness to in-plane resolution ratio of 1.6 as well as the recorded brain grid were in the same order as in human DTI studies since the transformation to an iso-grid of 50 μm corresponds to an iso-grid of 1 mm in human studies. Spatial normalization to a stereotaxic standard space (Fig. 1a) was performed using a study-specific b0-template and an FA-template  (Fig. 1c). The whole normalization process was iterative: scanner- and sequence specific b0- and fractional anisotropy (FA)-templates were created in a first step by arithmetically averaging data sets of all mice after linear transformation according to manually set landmarks identified using a stereotaxic mouse atlas . After this first iteration, data were non-linearily normalized during a second iteration step in order to further optimize the normalization matrices. This process was iteratively repeated until the correlation between the individual FA-maps and the FA-template was > 0.7 which was achieved after two iterations.
Motion artefacts were eliminated in each data set separately by a dedicated quality check procedure  (Fig. 1b). This quality check procedure was also used to select two high-quality scans of each mouse that were compared longitudinally.
DTI-based maps (FA, axial diffusivity – AD, radial diffusivity – RD, mean diffusivity - MD) were calculated from these MNI normalized data sets. DTI metrics maps (FA, AD, RD, MD) of each data set were calculated and were smoothed with a Gaussian filter of 200 μm full-width-at-half-maximum (Fig. 1b). The filter size which is about 2–3 times the recording voxel size provides a good balance between sensitivity and specificity.
Whole brain-based spatial statistics
Whole-brain-based spatial statistics (WBSS – ) was performed for cross-sectional comparison of TDP-43G298S mice vs wt mice at baseline and at follow-up, using an FA threshold of 0.2 . Statistical comparisons of DTI metrics maps for TDP-43G298S vs wt were performed voxelwise by means of the Student’s t-test. Statistical results were corrected for multiple comparisons with the false-discovery-rate (FDR) algorithm at p < 0.05 . Further reduction of the alpha error was performed by a spatial correction algorithm in the size range of the smoothing kernel leading to a cluster size threshold of 256 voxels.
WBSS of longitudinal DTI metrics’ map differences (ΔDM) was performed by calculating voxelwise differences between FA maps of baseline and follow-up scans for TDP-43G298S mice and wt mice; differences were then linearly normalized to an identical time-interval prior to statistical comparison:
t1 and t2 are the date of baseline and follow-up scans, DM is the respective DTI metric (AD, RD, MD). Statistical comparisons of ΔDM were performed voxelwise by means of the Student’s t-test and results were corrected for multiple comparisons and also by clustering with a cluster size threshold of 256 voxels (Fig. 1b).
In order to test for hemispherical symmetry, a hypothesis-guided approach was performed: spherical regions-of-interest (ROIs) were placed in defined anatomical regions, in this study in the motor areas, in the retrosplenial cortex, in the agranular insula, and as a reference in the visual cortex. Mean DM values within the respective ROI were calculated from each single DM map.
Tractwise fractional anisotropy statistics
An averaged data set was created from all contributing data sets, and fiber tracking (FT) was performed by deterministic streamline tracking technique  at an FA threshold of 0.2. In order to obtain a quantitative access to the tractography results, tractwise fractional anisotropy statistics (TFAS)  was applied for a hypothesis-guided analysis of FT bundles (Fig. 1d).
For the study of histological counterparts of DTI findings, wt or TDP-43G298S mice were sacrificed at 6–8 months (N = 3 and N = 3, respectively) or at 14–17 months of age (N = 3 and N = 4, respectively) by perfusion fixation (performed as previously reported; ). The brain was then dissected out, post-fixed in 4% PFA in PBS for 24 h at 4 °C, cryoprotected in 30% Sucrose in PBS, embedded in OCT (TissueTek, Fischer Scientifics) and sectioned at 40 μm thickness in a Leica Cryotome (Leica Cryotom AG Protect, CM 1950). For each brain, 6 sections were considered for further processing, starting 0.6 frontal of bregma and every 40 μm thereafter for 500 μm; all sections were processed together to ensure the homogeneity of the staining procedure. Immunostaining was performed as previously reported : briefly, brain sections were blocked and permeabilized in 3%BSA/2%FBS/0.1% Triton X-100 for 1 h at 24 °C and then incubated for 48 h a 4 °C with anti NeuN antibody (diluted 1:100 in blocking buffer); thereafter, the sections were washed in PBS (3 × 30 min, 24 °C), exposed to secondary antibody (donkey anti mouse conjugated with Alexa647 fluorochrome, 1:1000 in blocking buffer, 2 h at 24 °C), washed once again (3 × 30 min in PBS) and mounted with ProLong Gold (ThermoFisher). Each brain section was imaged in its entirety using a Leica microscope (Leica Live Cell Imaging Microscope DMI 6000B equipped with a 10X objective) with 12-bits image depth; imaging parameters were set to avoid over- or undersaturated areas and were kept constant when imaging the different sections.
For the quantification, each image was imported in ImageJ, subject to background subtraction (rolling-ball mode), thresholded at 130 (grey value in arbitrary units ranging between 0 and 4095) and contrast-enhanced (range 130–1300). A rectangular region of interest (approx. Size 2.5*105 μm2) was manually located in correspondence of the layer V of primary motor cortex, identified based on the anatomical landmarks according to the Allen Reference Brain Atlas (available at https://mouse.brain-map.org/static/atlas). Within the region of interest, neurons were manually counted by an experimenter blind to the genotype. Neuronal counts were obtained for left and right motor cortex; sites containing processing or staining artifacts were excluded. For each mouse, average neuronal density was calculated from at least 4 artifact-free sections. Statistical analysis was performed using two-way ANOVA (age and genotype) with Tukey-corrected post-hoc comparison. Statistical significance was set at p < 0.05 after correction for multiple comparisons.
Whole brain-based spatial statistics show loss of FA in motor cortex of TDP-43 mice
The experimental design included two sets of pre-specified analysis: first, the cross-sectional comparison of wt and TDP-43G298S cohorts at baseline and at follow-up and second, the longitudinal comparison of wt and TDP-43G298S mice compared to their baseline scans. Based on previously published data on longitudinal progression of the TDP-43G298S mice, we set the timepoint for the baseline at approx. 8–9 months of age, when grip strength is still comparable in TDP-43G298S and control mice and the timepoint for the follow-up at on average 14 months of age, when TDP-43G298S mice display a significant decrease in grip strength ; within this timeframe the performance on the running wheel is significantly reduced in TDP-43G298S mice compared to wt mice at both timepoints but does not decline between the timepoints and mice do not yet show weight loss or visible amyotrophy. Thus, the experimental design is set to identify abnormalities that may correlated with the onset of disturbances in the grip strength in TDP-43G298S mice.
The FA maps for wt mice displayed remarkable stability over time. No significant clusters in the comparison of baseline and follow-up FA maps of wt mice were detected.
At baseline at 9 months of age TDP-43G298S and wt did not display any significant regional differences across the whole brain volume. However, at the follow-up at 14 months of age the comparison of TDP-43G298S vs. wt mice revealed a significant reduction in FA in a cluster of voxel spanning the boundary between the primary motor cortex (M1) and the adjacent secondary motor cortex (M2) in the right hemisphere (Table 1).
Notably, when baseline and follow-up scans of the TDP-43G298S mice were compared (so to obtain an intra-subject comparison and further decrease variability), a significant decline in FA was identified not only in the M1/M2 area but also in clusters located in the cingulate cortex, in the retrosplenial granular cortex, and in the posterior thalamic nucleus (Table 1).
Longitudinal voxelwise differences of TDP-43G298S mice were compared to longitudinal voxelwise differences of wt mice. ΔFA was linearly normalized to an identical time-interval prior to statistical comparison. The WBSS of ΔFA for TDP-43G298S vs wt mice not only confirmed the loss of FA in the M1/M2 cluster, but also demonstrated an additional cluster in the agranular insular cortex (Fig. 2, Table 1).
Region-of-interest analyses of the lateralization of the FA loss in motor cortex
Since the WBSS data detected a hemispheric asymmetry in the decrease in FA in TDP-43G298S mice, which may be a statistical artifact of the whole-brain comparison, a hypothesis-guided ROI analysis was performed with spherical ROIs located either in the right M1 and M2 or in the left M1 and M2 areas. Confirming the WBSS data, a significant FA reduction was observed at follow-up for the TDP-43G298S mice in the right motor areas, whereas the left motor areas showed a trend (p = 0.1) for FA reduction (Fig. 3).
Diffusion metrics reveals a selective increase in mean and radial diffusivity in motor cortex of TDP-43G298S mice
Although FA changes could be reliably identified in the motor cortex of TDP-43G298S mice, these findings alone do not provide information on the underlying the histopathological substrate. We then explored whether the FA loss could be due to changes in the microstructure of the motor cortex. To this aim, we considered the diffusivity metrics of the motor cortex in wt and TDP-43G298S mice. While FA is sensitive to microstructural changes, it does not indicate a specific type of lesion; on the other hand, among the FA metrics, MD is an inverse measure of the membrane density, AD tends to be strongly affected by axonal injury whereas RD is sensitive to white matter damage due to demyelination and less to changes in the axonal density or size [19, 37, 41, 44].
Although no difference was detected in the comparison of wt and TDP-43G298S mice at baseline (in agreement with the FA values), the whole-brain voxelwise comparison of AD, MD and RD in wt vs TDP-43 mice at follow-up revealed a significant increase of RD in the motor cortex, whereas clusters on decreased AD were detected in correspondence of the retrosplenial cortex (Fig. 4).
Thus, these findings supported the hypothesis that FA alterations may be due to selective disturbances of the microstructure in the motor and extra-motor cortex of TDP-43G298S mice.
Tractwise fractional anisotropy reveals no damage to corticospinal and transcallosal fibers
Primary and secondary motor cortex give origin to multiple fiber tracts, connecting them either to sub-cortical and spinal targets (i.e., corticostriatal and corticospinal tracts – CST) or to other cortical areas (e.g., contralateral homonymous regions; ). Since degeneration of the CST is a hallmark of ALS in human subjects (and detectable by DTI; [16, 17]), we investigated whether the regional FA decrease observed in motor areas was specifically correlated with the disturbance of a subset of efferent fibers from motor areas. To this aim, TFAS was performed for the CST and in the transcallosal projections between ipsi- and contralateral motor cortices. Despite the detection of FA loss in the motor cortex, values of MD, RD and AD in the CST remained comparable in wt and TDP-43G298S mice at all timepoints, indicating that the disturbances observed in the whole-brain analysis are not due to CST degeneration. On the other hand, a distinct change in MD and RD values was detected in trans-callosal fibers (Fig. 5). Thus, changes in cortical FA metrics are correlated with alterations in a subset of projection neurons not including those giving rise to the CST.
Histological analysis of neuronal density reveals neuronal loss corresponding to FA loss in primary motor cortex
The anatomical basis underlying the microstructural changes (detected by FA) was investigated by evaluating the extent of neuronal loss in M1 in two samples of mice at two timepoints (independent of those subject to MRI imaging); the earlier timepoint was at 6–8 months of age (wt (N = 3) and TDP-43G298S (N = 3)) and the later timepoint was at 13–17 months of age (wt (N = 3) and TDP-43G298S (N = 4)). The density of NeuN+ cells in layer V of M1 (where large projection neurons, such as corticospinal neurons, are located) was assessed. Two-way ANOVA identified a significant effect of time (F (1,9) = 46.90, p < 0.001) and genotype (F (1,9) = 7.90, p = 0.020) and a significant interaction between the factors (F (1,9) = 14.75, p = 0.004); post-hoc analysis (Tukey correction for multiple comparisons) revealed comparable neuronal density in wt at the earlier and later timepoint (122.3 ± 2.7 NeuN+/105μm2 at 6 months and 114.9 ± 1.9 NeuN+/105μm2 at 14 months, p > 0.05), whereas TDP-43G298S mice were comparable to wt mice at 6 months (124.8 ± 5.5 NeuN+/105μm2, p > 0.05) but displayed a significant decline over time at 14 months (98.5 ± 5.5 NeuN+/105μm2, p = 0.001) (Fig. 6). However, due to the small sample size for histological analysis, a test for a lateralization effect was not possible.
In the present work we have demonstrated (i) the technical possibilities for long-term longitudinal DTI analyses in murine models of neurodegenerative diseases, specifically in an ALS model using ultrahigh field (11.7 T) with 30 gradient directions, (ii) the appearance over time of FA abnormalities in the motor cortex as well cingulate, insular and retrosplenial cortices of mice with mutant TDP-43 overexpression in neurons (Thy1.2 promoter) mice, and (iii) the dissociation between pathologic values of diffusivity indexes in motor cortex and normal values in the CST.
Technical considerations: reproducible longitudinal data in murine DTI imaging
Although previous studies have demonstrated the broad feasibility of repeated DTI in murine imaging [6, 22], they were limited in spatial resolution and by the number of gradient directions; thus, previous analysis used segmentation techniques that did not allow a full-brain comparison and were restricted to identify changes in spinal cord . In this study, the feasibility of microstructural analysis in a whole-brain spatial statistics analysis and in vivo fiber tracking of the DTI-based analysis of the mouse brain at the cohort level was demonstrated. Technically, this study has shown that, by spatial rescaling of the DTI data generated in high-resolution ultrahigh field with 30 gradient directions, the multiparametric analysis cascade as it has been established in human DTI data could be transferred to an analysis cascade for murine DTI data (in a whole-brain spatial statistics analysis, not previously reported in murine ALS models). Moreover, due to scanning at the ultrahigh field of 11.7 T, murine DTI data could be recorded with a high signal to noise ratio (SNR). Furthermore, rodent scans could be repeated frequently, that way offering the opportunity to obtain high quality repeated scans.
TDP-43G298S mice show disturbed cortical microstructure with spared corticospinal tract
Several studies in SOD1 transgenic ALS murine models have revealed DTI disturbances in the spinal cord, findings interpreted in terms of axonal loss affecting the cortico-spinal tract and other long white matter tracts stretching between the brain and the spinal cord [6, 22, 40]. However, despite the well-accepted observation of DTI abnormalities in ALS patients [5, 17], technical limitations have prevented so far the analysis of DTI patterns in the brain of murine models over time. Having overcome these limitations, our findings suggest that TDP-43G298S mice display involvement of the motor cortex and other cortical regions and this is seemingly not associated with involvement of the corticospinal tract (whose degeneration is an hallmark of ALS – ) but is rather linked to the involvement of local projection neurons (such as those projecting through the corpus callosum) and to intracortical microstructural change (highlighted by AD changes in the retrosplenial cortex). Importantly, these alterations are not detectable at a stage in which TDP-43G298S mice show normal grip strength but appear once grip strength is strongly affected , suggesting that they may provide an imaging correlate of neurodegenerative changes. In fact, a neuronal loss in M1 was shown where DTI demonstrates significant changes in FA, indicating that degeneration of neurons (and of dendritic and axonal structures associated) may be one (not necessarily the only one) of the processes generating the abnormal FA signal. Thus, FA may be used as a non-invasive readout of the neurodegenerative process in TDP-43G298S mice.
The pattern of abnormalities in the TDP-43G298S mice, as detected at this stage of the disease progression, shows closer similarities with FTD patients rather than with ALS patients: whereas in the latter the involvement of motor cortex and CST is prominent early in disease progression [5, 17], in the former a subset of patients with intact CST can be identified . FTD patients also show DTI changes in frontal callosal fibers  and in the thalamus . It is important to note that the dissociation between grip-strength abnormalities, FA changes and normal CST metrics may be limited to the stage of progression under study (about 450 days of age); whether CST abnormalities might appear at later stages in the TDP-43G298S mice, with the progression of the involvement of neuronal subpopulations, remains object of active investigation.
Furthermore, the cortical pattern of involvement observed in the TDP-43G298S mice, involving primary and secondary motor cortices as well as agranular insula and retrosplenial cortices shows similarities with the involvement of frontal cortex, insular and posterior cingualate areas in FTD and the related Semantic Dementia [35, 36]. In agreement with this interpretation, FTD-like behavioural abnormalities (such as hyperphagic behavior) have been reported in unrelated TDP-43 transgenic model . Our findings include, intriguingly, an asymmetry in the involvement of the primary and secondary motor cortices (with statistically significant changes appearing in the right but not in the left hemisphere) in the TDP-43G298S mice; it is worth noting that asymmetric cortical involvement is a well-known feature of FTD and ALS [14, 35] and that mice display a significant degree of lateralization in gene expression patterns and neurochemical properties (e.g., ). Analysis of larger cohorts may reveal if the asymmetry observed in this study is due to statistical noise or reflects an increased vulnerability of the right hemisphere in the TDP-43G298S mouse model.
FA metrics can be sensitive to multiple pathological changes and, although they may offer insights into the ongoing processes, they cannot be unequivocally attributed to a single histological alteration . The histological analysis of mouse brains suggested that at least one process revealed by FA analysis is the loss of neurons (and related structures) in the motor cortex. In addition, and more speculatively, the differential effect on RD (but not on AD) observed in motor cortex, and the simultaneous involvement of RD and MD in transcallosal fibers, may point toward the disproportionate involvement of myelin and myelinated fibers. Indeed, involvement of white matter in the form of white matter hyperintensities has been shown in FTD patients using FLAIR and T2 MRI imaging [33, 38] and, in a case for which pathology is available, has been related to activation of microglia rather than vascular changes  and may result from primarily axonal damage. On the other hand, the selective AD changes in retrosplenial cortex appear to point toward a selective disruption of axonal microstructure; the simultaneous involvement of retrosplenial and M1/M2 may suggest that the reciprocal circuit linking these two areas  may be selectively involved.
The study was not without limitations. The number of mice was comparatively small, especially in contrast to previous human studies with high subject numbers longitudinally , but are similar to other small animal studies. Thus, the significance levels are limited in a cohort of mice with limited sample size. A further limitation was that follow-up scans could not be acquired in all animals since single wt mice and TDP-43G298S mice died during the course of the study. Moreover, mutant TDP-43 protein is expressed under control of the Thy1.2 promoter in this model, thus lacking expression in non-neuronal cells, while in the human situation TDP-43 is ubiquitously expressed, including glial cells; thus, our study may have revealed the cell-autonomous effects of mutant TDP-43 expression, and more severe phenotypes may occur in case of neuronal and non-neuronal expression of the transgene.
This study demonstrated the technical feasibility of high-resolution, longitudinal DTI for preclinical cohort studies in murine models at 11.7 T. DTI-based analysis of the brain of TDP-43G298S mice showed alterations of motor areas M1/M2 (showing a temporal correlation with neuronal loss) which increase over time as well as a broader involvement of retrosplenial cortex and agranular insula. Taken together, the DTI data suggest that, within the ALS-FTD spectrum, the TDP-43G298S mice show similarities to FTD-ALS syndrome rather than pure motoneuron disease. In this context, degeneration of transcallosal and intracortical projections (with notable sparing of the CST) appears to be the most prominent pathological process ongoing in TDP-43G298S mice.
Availability of data and materials
The datasets generated during and/or analysed during the current study are not publicly available as the acquired DTI data are property of the University of Ulm, but are available from the corresponding author on reasonable request.
Amyotrophic lateral sclerosis
Diffusion tensor imaging
Standard error of the mean
Signal to noise ratio
Tractwise fractional anisotropy statistics
Tensor Imaging and Fiber Tracking
Whole brain-based spatial statistics
- wt :
Aggarwal M, Zhang J, Miller MI, Sidman RL, Mori S. Magnetic resonance imaging and micro-computed tomography combined atlas of developing and adult mouse brains for stereotaxic surgery. Neuroscience. 2009;162:1339–50.
Aggarwal M, Mori S, Shimogori T, Blackshaw S, Zhang J. Three-dimensional diffusion tensor microimaging for anatomical characterization of the mouse brain. Magn Reson Med. 2010;64:249–61.
Bargsted L, Medinas DB, Martínez Traub F, Rozas P, Muñoz N, Nassif M, Jerez C, Catenaccio A, Court FA, Hetz C, Matus S. Disulfide cross-linked multimers of TDP-43 and spinal motoneuron loss in a TDP-43A315T ALS/FTD mouse model. Sci Rep. 2017;7:14266.
Braak H, Brettschneider J, Ludolph AC, Lee VM, Trojanowski JQ, Del Tredici K. Amyotrophic lateral sclerosis - a model of corticofugal axonal spread. Nat rev Neurol. 2013;9:708–14.
Canu E, Agosta F, Riva N, Sala S, Prelle A, Caputo D, Perini M, Comi G, Filippi M. The topography of brain microstructural damage in amyotrophic lateral sclerosis assessed using diffusion tensor MR imaging. AJNR Am J Neuroradiol. 2011;32:1307–14.
Caron I, Micotti E, Paladini A, Merlino G, Plebani L, Forloni G, Modo M, Bendotti C. Comparative magnetic resonance imaging and histopathological correlates in two SOD1 transgenic mouse models of amyotrophic lateral sclerosis. PLoS One. 2015;10:e0132159.
Commisso B, Ding L, Varadi K, Gorges M, Bayer D, Boeckers TM, Ludolph AC, Kassubek J, Mueller O, Roselli F. Stage-dependent remodeling of projections to motor cortex in ALS mouse model revealed by a new variant retrograde-AAV9. Elife. 2018;7:e36892.
Daianu M, Mendez MF, Baboyan VG, Jin Y, Melrose RJ, Jimenez EE, Thompson PM. An advanced white matter tract analysis in frontotemporal dementia and early-onset Alzheimer's disease. Brain Imaging Behav. 2016;10:1038–53.
Ebstein SY, Yagudayeva I, Shneider NA. Mutant TDP-43 causes early-stage dose-dependent motor neuron degeneration in a TARDBP Knockin mouse model of ALS. Cell Rep. 2019;26:364–373.e4.
Genovese CR, Lazar NA, Nichols T. Thresholding of statistical maps in functional neuroimaging using the false discovery rate. Neuroimage. 2002;15:870–8.
Grabrucker S, Haderspeck JC, Sauer AK, Kittelberger N, Asoglu H, Abaei A, Rasche V, Schön M, Boeckers TM, Grabrucker AM. Brain lateralization in mice is associated with zinc signaling and altered in prenatal zinc deficient mice that display features of autism Spectrum disorder. Front Mol Neurosci. 2018;10:450.
Haber M, Hutchinson EB, Sadeghi N, Cheng WH, Namjoshi D, Cripton P, Irfanoglu MO, Wellington C, Diaz-Arrastia R, Pierpaoli C. Defining an analytic framework to evaluate quantitative MRI markers of traumatic axonal injury: preliminary results in a mouse closed head injury model. eNeuro. 2017;4(5). ENEURO 0164-172017.
Harsan LA, Paul D, Schnell S, Kreher BW, Hennig J, Staiger JF, von Elverfeldt D. In vivo diffusion tensor magnetic resonance imaging and fiber tracking of the mouse brain. NMR Biomed. 2010;23:884–96.
Irwin DJ, McMillan CT, Xie SX, Rascovsky K, Van Deerlin VM, Coslett HB, Hamilton R, Aguirre GK, Lee EB, Lee VMY, Trojanowski JQ, Grossman M. Asymmetry of post-mortem neuropathology in behavioural-variant frontotemporal dementia. Brain. 2018;141:288–301.
Jakabek D, Power BD, Macfarlane MD, Walterfang M, Velakoulis D, van Westen D, Lätt J, Nilsson M, Looi JCL, Santillo AF. Regional structural hypo- and hyperconnectivity of frontal-striatal and frontal-thalamic pathways in behavioral variant frontotemporal dementia. Hum Brain Mapp. 2018;39:4083–93.
Kassubek J, Müller HP, Del Tredici K, Brettschneider J, Pinkhardt EH, Lulé D, Böhm S, Braak H, Ludolph AC. Diffusion tensor imaging analysis of sequential spreading of disease in amyotrophic lateral sclerosis confirms patterns of TDP-43 pathology. Brain. 2014;137:1733–40.
Kassubek J, Müller H-P, Del Tredici K, Lulé D, Gorges M, Braak H, Ludolph AC. Imaging the pathoanatomy of amyotrophic lateral sclerosis in vivo: targeting a propagation-based biological marker. J Neurol Neurosurg Psychiatry. 2018a;89:374–81.
Kassubek J, Müller HP, Del Tredici K, Hornberger M, Schroeter ML, Müller K, Anderl-Straub S, Uttner I, Grossman M, Braak H, Hodges JR, Piguet O, Otto M, Ludolph AC. Longitudinal diffusion tensor imaging resembles patterns of pathology progression in behavioral variant frontotemporal dementia (bvFTD). Front Aging Neurosci. 2018b;10:47.
Kim JH, Wu TH, Budde MD, Lee JM, Song SK. Noninvasive detection of brainstem and spinal cord axonal degeneration in an amyotrophic lateral sclerosis mouse model. NMR Biomed. 2011;24:163–9.
Kunimatsu A, Aoki S, Masutani Y, Abe O, Hayashi N, Mori H, Masumoto T, Ohtomo K. The optimal trackability threshold of fractional anisotropy for diffusion tensor tractography of the corticospinal tract. Magn Reson Med Sci. 2004;3:11–7.
Lee JD, Levin SC, Willis EF, Li R, Woodruff TM, Noakes PG. Complement components are upregulated and correlate with disease progression in the TDP-43Q331K mouse model of amyotrophic lateral sclerosis. J Neuroinflammation. 2018;15:171.
Marcuzzo S, Bonanno S, Figini M, Scotti A, Zucca I, Minati L, Riva N, Domi T, Fossaghi A, Quattrini A, Galbardi B, D'Alessandro S, Bruzzone MG, García-Verdugo JM, Moreno-Manzano V, Mantegazza R, Bernasconi P. A longitudinal DTI and histological study of the spinal cord reveals early pathological alterations in G93A-SOD1 mouse model of amyotrophic lateral sclerosis. Exp Neurol. 2017;293:43–52.
Müller HP, Unrath A, Ludolph AC, Kassubek J. Preservation of diffusion tensor properties during spatial normalization by use of tensor imaging and fibre tracking on a normal brain database. Phys Med Biol. 2007a;52:N99–109.
Müller HP, Unrath A, Sperfeld AD, Ludolph AC, Riecker A, Kassubek J. Diffusion tensor imaging and tractwise fractional anisotropy statistics: quantitative analysis in white matter pathology. Biomed Eng Online. 2007b;6:42.
Müller HP, Süßmuth SD, Landwehrmeyer GB, Ludolph AC, Tabrizi SJ, Klöppel S, Kassubek J. Stability effects on results of diffusion tensor imaging analysis by reduction of the number of gradient directions due to motion artifacts: an application to presymptomatic Huntington’s disease. PLoS Curr. 2011;3:RRN1292.
Müller HP, Unrath A, Huppertz HJ, Ludolph AC, Kassubek J. Neuroanatomical patterns of cerebral white matter involvement in different motor neuron diseases as studied by diffusion tensor imaging analysis. Amyotroph Lateral Scler. 2012a;13:254–64.
Müller H-P, Vernikouskaya I, Ludolph AC, Kassubek J, Rasche V. Fast diffusion tensor magnetic resonance imaging of the mouse brain at ultrahigh-field: aiming at cohort studies. PLoS One. 2012b;7:e53389.
Müller H-P, Kassubek J, Vernikouskaya I, Ludolph AC, Stiller D, Rasche V. Diffusion tensor magnetic resonance imaging of the brain in APP transgenic mice: a cohort study. PLoS One. 2013;8:e67630.
Müller HP, Gorges M, Kassubek R, Dorst J, Ludolph AC, Kassubek J. Identical patterns of cortico-efferent tract involvement in primary lateral sclerosis and amyotrophic lateral sclerosis: a tract of interest-based MRI study. Neuroimage Clin. 2018;18:762–9.
Nouls JC, Badea A, Anderson RBJ, Cofer GP, Allan JG. Diffusion tensor imaging using multiple coils for mouse brain connectomics. NMR Biomed. 2018;31:e3921.
Oh SW, Harris JA, Ng L, Winslow B, Cain N, Mihalas S, Wang Q, Lau C, Kuan L, Henry AM, Mortrud MT, Ouellette B, Nguyen TN, Sorensen SA, Slaughterbeck CR, Wakeman W, Li Y, Feng D, Ho A, Nicholas E, Hirokawa KE, Bohn P, Joines KM, Peng H, Hawrylycz MJ, Phillips JW, Hohmann JG, Wohnoutka P, Gerfen CR, Koch C, Bernard A, Dang C, Jones AR, Zeng H. A mesoscale connectome of the mouse brain. Nature. 2014;508:207–14.
Paxinos G, Franklin KBJ. The mouse brain in stereotaxic coordinates. 3rd ed: Academic. London: Elsevier; 2007.
Paternicó D, Premi E, Gazzina S, Cosseddu M, Alberici A, Archetti S, Cotelli MS, Micheli A, Turla M, Gasparotti R, Padovani A, Borroni B. White matter hyperintensities characterize monogenic frontotemporal dementia with granulin mutations. Neurobiol Aging. 2016;38:176–80.
Petrella LI, Castelhano JM, Ribeiro M, Sereno JV, Gonçalves SI, Laço MN, Hayden MR, Rego AC, Castelo-Branco M. A whole brain longitudinal study in the YAC128 mouse model of Huntington's disease shows distinct trajectories of neurochemical, structural connectivity and volumetric changes. Hum Mol Genet. 2018;27:2125–37.
Rosen HJ, Gorno-Tempini ML, Goldman WP, Perry RJ, Schuff N, Weiner M, Feiwell R, Kramer JH, Miller BL. Patterns of brain atrophy in frontotemporal dementia and semantic dementia. Neurology. 2002;58:198–208.
Scheltens NME, van der Weijden K, Adriaanse SM, van Assema D, Oomen PP, Krudop WA, Lammertsma AA, Barkhof F, Koene T, Teunissen CE, Scheltens P, van der Flier WM, Pijnenburg YAL, Yaqub M, Ossenkoppele R, van Berckel BNM. Hypometabolism of the posterior cingulate cortex is not restricted to Alzheimer's disease. Neuroimage Clin. 2018;19:625–32.
Song S-K, Sun S-W, Ramsbottom MJ, Chang C, Russell J, Cross AH. Dysmyelination revealed through MRI as increased radial (but unchanged axial) diffusion of water. Neuroimage. 2002;17:1429–36.
Sudre CH, Bocchetta M, Cash D, Thomas DL, Woollacott I, Dick KM, van Swieten J, Borroni B, Galimberti D, Masellis M, Tartaglia MC, Rowe JB, Graff C, Tagliavini F, Frisoni G, Laforce R Jr, Finger E, de Mendonça A, Sorbi S, Ourselin S, Cardoso MJ, Rohrer JD. Genetic FTD initiative, GENFI White matter hyperintensities are seen only in GRN mutation carriers in the GENFI cohort. Neuroimage Clin. 2017;15:171–80.
Tsao W, Jeong YH, Lin S, Ling J, Price DL, Chiang PM, Wong PC. Rodent models of TDP-43: recent advances. Brain Res. 2012;1462:26–39.
Underwood CK, Kurniawan ND, Butler TJ, Cowin GJ, Wallace RH. Non-invasive diffusion tensor imaging detects white matter degeneration in the spinal cord of a mouse model of amyotrophic lateral sclerosis. Neuroimage. 2011;55:455–61.
Wang Y, Sun P, Wang Q, Trinkaus K, Schmidt RE, Naismith RT, Cross AH, Song SK. Differentiation and quantification of inflammation, demyelination and axon injury or loss in multiple sclerosis. Brain. 2015;138:1223–38.
White MA, Kim E, Duffy A, Adalbert R, Phillips BU, Peters OM, Stephenson J, Yang S, Massenzio F, Lin Z, Andrews S, Segonds-Pichon A, Metterville J, Saksida LM, Mead R, Ribchester RR, Barhomi Y, Serre T, Coleman MP, Fallon JR, Bussey TJ, Brown RH Jr, Sreedharan J. TDP-43 gains function due to perturbed autoregulation in a Tardbp knock-in mouse model of ALS-FTD. Nat Neurosci. 2018;21:552–63.
Wiesner D, Tar L, Linkus B, Chandrasekar A, Olde Heuvel F, Dupuis L, Tsao W, Wong PC, Ludolph A, Roselli F. Reversible induction of TDP-43 granules in cortical neurons after traumatic injury. Exp Neurol. 2018;299:15–25.
Winklewski PJ, Sabisz A, Naumczyk P, Jodzio K, Szurowska E, Szarmach A. Understanding the physiopathology behind axial and radial diffusivity changes-what do we know? Front Neurol. 2018;9:92.
Yamawaki N, Radulovic J, Shepherd GM. A Corticocortical circuit directly links Retrosplenial cortex to M2 in the mouse. J Neurosci. 2016;36:9365–74.
Woollacott IOC, Bocchetta M, Sudre CH, Ridha BH, Strand C, Courtney R, Ourselin S, Cardoso MJ, Warren JD, Rossor MN, Revesz T, Fox NC, Holton JL, Lashley T, Rohrer JD. Pathological correlates of white matter hyperintensities in a case of progranulin mutation associated frontotemporal dementia. Neurocase. 2018;16:1–9.
The authors would like to thank the Ulm University Center for Translational Imaging MoMAN for its support. Johannes Ulmer, Clara Bruno, Kirsten Sieverding are thankfully acknowledged for housing the mice. FR was supported by the Ulm University Medical School Baustein program, by Deutsche ForschungsGemeinschaft (DFG) as part of the SFB1149 and by the Thierry-Latran Foundation. This study was supported by the German Center for Neurodegenerative Diseases (Deutsches Zentrum für Neurodegenerative Erkrankungen – DZNE) a foundation of the Helmholtz-Gemeinschaft. Tobias M. Boeckers, Institute for Anatomy and Cell Biology, Ulm University, Ulm, Germany, is thankfully acknowledged for providing the Leica Live Cell Imaging Microscope.
Ethics approval and consent to participate
All animal experiments were performed in accordance with institutional guidelines of Ulm University and were approved by the regulatory authority (Regierungspräsidium Tübingen, Germany; animal permission no. 1242).
Consent for publication
The authors declare that they have no competing interests.