1
|
Friesen E, Hari K, Sheft M, Thiessen JD, Martin M. Magnetic resonance metrics for identification of cuprizone-induced demyelination in the mouse model of neurodegeneration: a review. MAGMA (NEW YORK, N.Y.) 2024:10.1007/s10334-024-01160-z. [PMID: 38635150 DOI: 10.1007/s10334-024-01160-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Revised: 03/17/2024] [Accepted: 03/26/2024] [Indexed: 04/19/2024]
Abstract
Neurodegenerative disorders, including Multiple Sclerosis (MS), are heterogenous disorders which affect the myelin sheath of the central nervous system (CNS). Magnetic Resonance Imaging (MRI) provides a non-invasive method for studying, diagnosing, and monitoring disease progression. As an emerging research area, many studies have attempted to connect MR metrics to underlying pathophysiological presentations of heterogenous neurodegeneration. Most commonly, small animal models are used, including Experimental Autoimmune Encephalomyelitis (EAE), Theiler's Murine Encephalomyelitis (TMEV), and toxin models including cuprizone (CPZ), lysolecithin, and ethidium bromide (EtBr). A contrast and comparison of these models is presented, with focus on the cuprizone model, followed by a review of literature studying neurodegeneration using MRI and the cuprizone model. Conventional MRI methods including T1 Weighted (T1W) and T2 Weighted (T2W) Imaging are mentioned. Quantitative MRI methods which are sensitive to diffusion, magnetization transfer, susceptibility, relaxation, and chemical composition are discussed in relation to studying the CPZ model. Overall, additional studies are needed to improve both the sensitivity and specificity of MRI metrics for underlying pathophysiology of neurodegeneration and the relationships in attempts to clear the clinico-radiological paradox. We therefore propose a multiparametric approach for the investigation of MR metrics for underlying pathophysiology.
Collapse
Affiliation(s)
- Emma Friesen
- Chemistry, University of Winnipeg, Winnipeg, Canada.
| | - Kamya Hari
- Physics, University of Winnipeg, Winnipeg, Canada
- Electronics and Communication Engineering, SSN College of Engineering, Chennai, India
| | - Maxina Sheft
- Physics, University of Winnipeg, Winnipeg, Canada
- Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, USA
| | - Jonathan D Thiessen
- Imaging Program, Lawson Health Research Institute, London, Canada
- Medical Biophysics, Western University, London, Canada
- Medical Imaging, Western University, London, Canada
| | | |
Collapse
|
2
|
Rowley CD, Campbell JSW, Leppert IR, Nelson MC, Pike GB, Tardif CL. Optimization of acquisition parameters for cortical inhomogeneous magnetization transfer (ihMT) imaging using a rapid gradient echo readout. Magn Reson Med 2023; 90:1762-1775. [PMID: 37332194 DOI: 10.1002/mrm.29754] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 04/25/2023] [Accepted: 05/22/2023] [Indexed: 06/20/2023]
Abstract
PURPOSE Imaging biomarkers with increased myelin specificity are needed to better understand the complex progression of neurological disorders. Inhomogeneous magnetization transfer (ihMT) imaging is an emergent technique that has a high degree of specificity for myelin content but suffers from low signal to-noise ratio (SNR). This study used simulations to determine optimal sequence parameters for ihMT imaging for use in high-resolution cortical mapping. METHODS MT-weighted cortical image intensity and ihMT SNR were simulated using modified Bloch equations for a range of sequence parameters. The acquisition time was limited to 4.5 min/volume. A custom MT-weighted RAGE sequence with center-out k-space encoding was used to enhance SNR at 3 T. Pulsed MT imaging was studied over a range of saturation parameters, and the impact of the turbo factor on the effective ihMT resolution was investigated. 1 mm isotropic ihMTsat maps were generated in 25 healthy adults. RESULTS Greater SNR was observed for larger number of bursts consisting of 6-8 saturation pulses each, combined with a high readout turbo factor. However, that protocol suffered from a point spread function that was more than twice the nominal resolution. For high-resolution cortical imaging, we selected a protocol with a higher effective resolution at the cost of a lower SNR. We present the first group-average ihMTsat whole-brain map at 1 mm isotropic resolution. CONCLUSION This study presents the impact of saturation and excitation parameters on ihMTsat SNR and resolution. We demonstrate the feasibility of high-resolution cortical myelin imaging using ihMTsat in less than 20 min.
Collapse
Affiliation(s)
- Christopher D Rowley
- McConnell Brain Imaging Centre, Montreal Neurological Institute and Hospital, McGill University, Montreal, Québec, Canada
- Department of Neurology and Neurosurgery, McGill University, Montreal, Québec, Canada
| | - Jennifer S W Campbell
- McConnell Brain Imaging Centre, Montreal Neurological Institute and Hospital, McGill University, Montreal, Québec, Canada
| | - Ilana R Leppert
- McConnell Brain Imaging Centre, Montreal Neurological Institute and Hospital, McGill University, Montreal, Québec, Canada
| | - Mark C Nelson
- McConnell Brain Imaging Centre, Montreal Neurological Institute and Hospital, McGill University, Montreal, Québec, Canada
- Department of Neurology and Neurosurgery, McGill University, Montreal, Québec, Canada
| | - G Bruce Pike
- Hotchkiss Brain Institute and Departments of Radiology and Clinical Neuroscience, Cumming School of Medicine, University of Calgary, Calgary, Alberta, Canada
| | - Christine L Tardif
- McConnell Brain Imaging Centre, Montreal Neurological Institute and Hospital, McGill University, Montreal, Québec, Canada
- Department of Neurology and Neurosurgery, McGill University, Montreal, Québec, Canada
- Department of Biomedical Engineering, McGill University, Montreal, Québec, Canada
| |
Collapse
|
3
|
Wiggermann V, Endmayr V, Hernández‐Torres E, Höftberger R, Kasprian G, Hametner S, Rauscher A. Quantitative magnetic resonance imaging reflects different levels of histologically determined myelin densities in multiple sclerosis, including remyelination in inactive multiple sclerosis lesions. Brain Pathol 2023; 33:e13150. [PMID: 36720269 PMCID: PMC10580011 DOI: 10.1111/bpa.13150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2022] [Accepted: 11/16/2022] [Indexed: 02/02/2023] Open
Abstract
Magnetic resonance imaging (MRI) of focal or diffuse myelin damage or remyelination may provide important insights into disease progression and potential treatment efficacy in multiple sclerosis (MS). We performed post-mortem MRI and histopathological myelin measurements in seven progressive MS cases to evaluate the ability of three myelin-sensitive MRI scans to distinguish different stages of MS pathology, particularly chronic demyelinated and remyelinated lesions. At 3 Tesla, we acquired two different myelin water imaging (MWI) scans and magnetisation transfer ratio (MTR) data. Histopathology included histochemical stainings for myelin phospholipids (LFB) and iron as well as immunohistochemistry for myelin proteolipid protein (PLP), CD68 (phagocytosing microglia/macrophages) and BCAS1 (remyelinating oligodendrocytes). Mixed-effects modelling determined which histopathological metric best predicted MWF and MTR in normal-appearing and diffusely abnormal white matter, active/inactive, inactive, remyelinated and ischemic lesions. Both MWI measures correlated well with each other and histology across regions, reflecting the different stages of MS pathology. MTR data showed a considerable influence of components other than myelin and a strong dependency on tissue storage duration. Both MRI and histology revealed increased myelin densities in inactive compared with active/inactive lesions. Chronic inactive lesions harboured single scattered myelin fibres indicative of low-level remyelination. Mixed-effects modelling showed that smaller differences between white matter areas were linked to PLP densities and only to a small extent confounded by iron. MWI reflects differences in myelin lipids and proteins across various levels of myelin densities encountered in MS, including low-level remyelination in chronic inactive lesions.
Collapse
Affiliation(s)
- Vanessa Wiggermann
- Department of Physics and AstronomyUniversity of British ColumbiaVancouverBritish ColumbiaCanada
- Department of PediatricsUniversity of British ColumbiaVancouverBritish ColumbiaCanada
- Danish Research Centre for Magnetic ResonanceCopenhagen University Hospital Amager & HvidovreCopenhagenDenmark
| | - Verena Endmayr
- Division of Neuropathology and Neurochemistry, Department of NeurologyMedical University of ViennaViennaAustria
- Centre for Brain ResearchMedical University of ViennaViennaAustria
| | - Enedino Hernández‐Torres
- Danish Research Centre for Magnetic ResonanceCopenhagen University Hospital Amager & HvidovreCopenhagenDenmark
- Faculty of Medicine (Division Neurology)University of British ColumbiaVancouverBritish ColumbiaCanada
| | - Romana Höftberger
- Division of Neuropathology and Neurochemistry, Department of NeurologyMedical University of ViennaViennaAustria
| | - Gregor Kasprian
- Department of Biomedical Imaging and Image‐Guided TherapyMedical University of ViennaViennaAustria
| | - Simon Hametner
- Division of Neuropathology and Neurochemistry, Department of NeurologyMedical University of ViennaViennaAustria
- Centre for Brain ResearchMedical University of ViennaViennaAustria
| | - Alexander Rauscher
- Department of Physics and AstronomyUniversity of British ColumbiaVancouverBritish ColumbiaCanada
- Department of PediatricsUniversity of British ColumbiaVancouverBritish ColumbiaCanada
- Department of RadiologyUniversity of British ColumbiaVancouverBritish ColumbiaCanada
- BC Children's Hospital Research InstituteUniversity of British ColumbiaVancouverBritish ColumbiaCanada
| |
Collapse
|
4
|
Alsop DC, Ercan E, Girard OM, Mackay AL, Michal CA, Varma G, Vinogradov E, Duhamel G. Inhomogeneous magnetization transfer imaging: Concepts and directions for further development. NMR IN BIOMEDICINE 2023; 36:e4808. [PMID: 35916067 DOI: 10.1002/nbm.4808] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/18/2022] [Revised: 07/27/2022] [Accepted: 07/28/2022] [Indexed: 05/23/2023]
Abstract
Off-resonance radio frequency irradiation can induce the ordering of proton spins in the dipolar fields of their neighbors, in molecules with restricted mobility. This dipolar order decays with a characteristic relaxation time, T1D , that is very different from the T1 and T2 relaxation of the nuclear alignment with the main magnetic field. Inhomogeneous magnetization transfer (ihMT) imaging is a refinement of magnetization transfer (MT) imaging that isolates the MT signal dependence on dipolar order relaxation times within motion-constrained molecules. Because T1D relaxation is a unique contrast mechanism, ihMT may enable improved characterization of tissue. Initial work has stressed the high correlation between ihMT signal and myelin density. Dipolar order relaxation appears to be much longer in membrane lipids than other molecules. Recent work has shown, however, that ihMT acquisitions may also be adjusted to emphasize different ranges of T1D . These newer approaches may be sensitive to other microstructural components of tissue. Here, we review the concepts and history of ihMT and outline the requirements for further development to realize its full potential.
Collapse
Affiliation(s)
- David C Alsop
- Department of Radiology, Beth Israel Deaconess Medical Center and Harvard Medical School, Boston, Massachusetts, USA
| | - Ece Ercan
- MR Clinical Science, Philips, Best, The Netherlands
| | | | - Alex L Mackay
- Department of Physics and Astronomy, University of British Columbia, Vancouver, British Columbia, Canada
- Department of Radiology, University of British Columbia, Vancouver, British Columbia, Canada
| | - Carl A Michal
- Department of Physics and Astronomy, University of British Columbia, Vancouver, British Columbia, Canada
| | - Gopal Varma
- Department of Radiology, Beth Israel Deaconess Medical Center and Harvard Medical School, Boston, Massachusetts, USA
| | - Elena Vinogradov
- Department of Radiology and Advanced Imaging Research Center, UT Southwestern Medical Center, Dallas, Texas, USA
| | | |
Collapse
|
5
|
Soustelle L, Troalen T, Hertanu A, Ranjeva JP, Guye M, Varma G, Alsop DC, Duhamel G, Girard OM. Quantitative magnetization transfer MRI unbiased by on-resonance saturation and dipolar order contributions. Magn Reson Med 2023. [PMID: 37154400 DOI: 10.1002/mrm.29678] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 03/22/2023] [Accepted: 04/01/2023] [Indexed: 05/10/2023]
Abstract
PURPOSE To demonstrate the bias in quantitative MT (qMT) measures introduced by the presence of dipolar order and on-resonance saturation (ONRS) effects using magnetization transfer (MT) spoiled gradient-recalled (SPGR) acquisitions, and propose changes to the acquisition and analysis strategies to remove these biases. METHODS The proposed framework consists of SPGR sequences prepared with simultaneous dual-offset frequency-saturation pulses to cancel out dipolar order and associated relaxation (T1D ) effects in Z-spectrum acquisitions, and a matched quantitative MT (qMT) mathematical model that includes ONRS effects of readout pulses. Variable flip angle and MT data were fitted jointly to simultaneously estimate qMT parameters (macromolecular proton fraction [MPF], T2,f , T2,b , R, and free pool T1 ). This framework is compared with standard qMT and investigated in terms of reproducibility, and then further developed to follow a joint single-point qMT methodology for combined estimation of MPF and T1 . RESULTS Bland-Altman analyses demonstrated a systematic underestimation of MPF (-2.5% and -1.3%, on average, in white and gray matter, respectively) and overestimation of T1 (47.1 ms and 38.6 ms, on average, in white and gray matter, respectively) if both ONRS and dipolar order effects are ignored. Reproducibility of the proposed framework is excellent (ΔMPF = -0.03% and ΔT1 = -19.0 ms). The single-point methodology yielded consistent MPF and T1 values with respective maximum relative average bias of -0.15% and -3.5 ms found in white matter. CONCLUSION The influence of acquisition strategy and matched mathematical model with regard to ONRS and dipolar order effects in qMT-SPGR frameworks has been investigated. The proposed framework holds promise for improved accuracy with reproducibility.
Collapse
Affiliation(s)
- Lucas Soustelle
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | | | - Andreea Hertanu
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Jean-Philippe Ranjeva
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Maxime Guye
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Gopal Varma
- Division of MR Research, Department of Radiology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts, USA
| | - David C Alsop
- Division of MR Research, Department of Radiology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts, USA
| | - Guillaume Duhamel
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Olivier M Girard
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| |
Collapse
|
6
|
High-resolution magnetization-transfer imaging of post-mortem marmoset brain: Comparisons with relaxometry and histology. Neuroimage 2023; 268:119860. [PMID: 36610679 DOI: 10.1016/j.neuroimage.2023.119860] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2022] [Revised: 12/28/2022] [Accepted: 01/03/2023] [Indexed: 01/07/2023] Open
Abstract
Cell membranes and macromolecules or paramagnetic compounds interact with water proton spins, which modulates magnetic resonance imaging (MRI) contrast providing information on tissue composition. For a further investigation, quantitative magnetization transfer (qMT) parameters (at 3T), including the ratio of the macromolecular and water proton pools, F, and the exchange-rate constant as well as the (observed) longitudinal and the effective transverse relaxation rates (at 3T and 7T), R1obs and R2*, respectively, were measured at high spatial resolution (200 µm) in a slice of fixed marmoset brain and compared to histology results obtained with Gallyas' myelin stain and Perls' iron stain. R1obs and R2* were linearly correlated with the iron content for the entire slice, whereas distinct differences were obtained between gray and white matter for correlations of relaxometry and qMT parameters with myelin content. The combined results suggest that the macromolecular pool interacting with water consists of myelin and (less efficient) non-myelin contributions. Despite strong correlation of F and R1obs, none of these parameters was uniquely specific to myelination. Due to additional sensitivity to iron stores, R1obs and R2* were more sensitive for depicting microstructural differences between cortical layers than F.
Collapse
|
7
|
Inhomogeneous Magnetization Transfer (ihMT) imaging in the acute cuprizone mouse model of demyelination/remyelination. Neuroimage 2023; 265:119785. [PMID: 36464096 DOI: 10.1016/j.neuroimage.2022.119785] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2022] [Revised: 10/31/2022] [Accepted: 12/01/2022] [Indexed: 12/05/2022] Open
Abstract
BACKGROUND To investigate the association of ihMT (inhom signals with the demyelination and remyelination phases of the acute cuprizone mouse model in comparison with histology, and to assess the extent of tissue damage and repair from MRI data. METHODS Acute demyelination by feeding 0.2% cuprizone for five weeks, followed by a four-week remyelination period was applied on genetically modified plp-GFP mice. Animals were scanned at different time points of the demyelination and remyelination phases of the cuprizone model using a multimodal MRI protocol, including ihMT T1D-filters, MPF (Macromolecular Proton Fraction) and R1 (longitudinal relaxation rate). For histology, plp-GFP (proteolipid protein - Green Fluorescent Protein) microscopy and LFB (Luxol Fast Blue) staining were employed as references for the myelin content. Comparison of MRI with histology was performed in the medial corpus callosum (mCC) and cerebral cortex (CTX) at two brain levels whereas ROI-wise and voxel-based analyses of the MRI metrics allowed investigating in vivo the spatial extent of myelin alterations. RESULTS IhMT high-pass T1D-filters, targeted toward long T1D components, showed significant temporal variations in the mCC consistent with the effects induced by the cuprizone toxin. In addition, the corresponding signals correlated strongly and significantly with the myelin content assessed by GFP fluorescence and LFB staining over the demyelination and the remyelination phases. The signal of the band-pass T1D-filter, which isolates short T1D components, showed changes over time that were poorly correlated with histology, hence suggesting a sensitivity to pathological processes possibly not related to myelin. Although MPF was also highly correlated to histology, ihMT high-pass T1D-filters showed better capability to characterize the spatial-temporal patterns during the demyelination and remyelination phases of the acute cuprizone model (e.g., rostro-caudal gradient of demyelination in the mCC previously described in the literature). CONCLUSIONS IhMT sequences selective for long T1D components are specific and sensitive in vivo markers of demyelination and remyelination and have successfully captured the spatially heterogeneous pattern of the demyelination and remyelination mechanisms in the cuprizone model. Interestingly, differences in signal variations between the ihMT high-pass and band-pass T1D-filter, suggest a sensitivity of the ihMT sequences targeted to short T1Ds to alterations other than those of myelin. Future studies will need to further address these differences by examining more closely the origin of the short T1D components and the variation of each T1D component in pathology.
Collapse
|
8
|
Hertanu A, Soustelle L, Le Troter A, Buron J, Le Priellec J, Carvalho VND, Cayre M, Durbec P, Varma G, Alsop DC, Girard OM, Duhamel G. T 1D -weighted ihMT imaging - Part I. Isolation of long- and short-T 1D components by T 1D -filtering. Magn Reson Med 2022; 87:2313-2328. [PMID: 35037302 DOI: 10.1002/mrm.29139] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Revised: 10/29/2021] [Accepted: 12/12/2021] [Indexed: 01/15/2023]
Abstract
PURPOSE To identify T1D -filtering methods, which can specifically isolate various ranges of T1D components as they may be sensitive to different microstructural properties. METHODS Modified Bloch-Provotorov equations describing a bi-T1D component biophysical model were used to simulate the inhomogeneous magnetization transfer (ihMT) signal from ihMTRAGE sequences at high RF power and low duty-cycle with different switching time values for the dual saturation experiment: Δt = 0.0, 0.8, 1.6, and 3.2 ms. Simulations were compared with experimental signals on the brain gray and white matter tissues of healthy mice at 7T. RESULTS The lengthening of Δt created ihMT high-pass T1D -filters, which efficiently eliminated the signal from T1D components shorter than 1 ms, while partially attenuating that of longer components (≥ 1 ms). Subtraction of ihMTR images obtained with Δt = 0.0 ms and Δt = 0.8 ms generated a new ihMT band-pass T1D -filter isolating short-T1D components in the 100-µs to 1-ms range. Simulated ihMTR values in central nervous system tissues were confirmed experimentally. CONCLUSION Long- and short-T1D components were successfully isolated with high RF power and low duty-cycle ihMT filters in the healthy mouse brain. Future studies should investigate the various T1D -range microstructural correlations in in vivo tissues.
Collapse
Affiliation(s)
- Andreea Hertanu
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France.,APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Lucas Soustelle
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France.,APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Arnaud Le Troter
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France.,APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Julie Buron
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France.,APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France.,Aix Marseille Univ, CNRS, IBDM, Marseille, France
| | | | - Victor N D Carvalho
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France.,APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France.,Aix Marseille Univ, CNRS, ICR, Marseille, France
| | - Myriam Cayre
- Aix Marseille Univ, CNRS, IBDM, Marseille, France
| | | | - Gopal Varma
- Division of MR Research, Radiology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts, USA
| | - David C Alsop
- Division of MR Research, Radiology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts, USA
| | - Olivier M Girard
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France.,APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Guillaume Duhamel
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France.,APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| |
Collapse
|