1
|
Yang H, Wang G, Li Z, Li H, Zheng J, Hu Y, Cao X, Liao C, Ye H, Tian Q. Artificial intelligence for neuro MRI acquisition: a review. MAGMA (NEW YORK, N.Y.) 2024; 37:383-396. [PMID: 38922525 DOI: 10.1007/s10334-024-01182-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 06/11/2024] [Accepted: 06/14/2024] [Indexed: 06/27/2024]
Abstract
OBJECT To review recent advances of artificial intelligence (AI) in enhancing the efficiency and throughput of the MRI acquisition workflow in neuroimaging, including planning, sequence design, and correction of acquisition artifacts. MATERIALS AND METHODS A comprehensive analysis was conducted on recent AI-based methods in neuro MRI acquisition. The study focused on key technological advances, their impact on clinical practice, and potential risks associated with these methods. RESULTS The findings indicate that AI-based algorithms have a substantial positive impact on the MRI acquisition process, improving both efficiency and throughput. Specific algorithms were identified as particularly effective in optimizing acquisition steps, with reported improvements in workflow efficiency. DISCUSSION The review highlights the transformative potential of AI in neuro MRI acquisition, emphasizing the technological advances and clinical benefits. However, it also discusses potential risks and challenges, suggesting areas for future research to mitigate these concerns and further enhance AI integration in MRI acquisition.
Collapse
Affiliation(s)
- Hongjia Yang
- School of Biomedical Engineering, Tsinghua University, Beijing, China
| | - Guanhua Wang
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA
| | - Ziyu Li
- Wellcome Centre for Integrative Neuroimaging, FMRIB, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, UK
| | - Haoxiang Li
- School of Biomedical Engineering, Tsinghua University, Beijing, China
| | - Jialan Zheng
- School of Biomedical Engineering, Tsinghua University, Beijing, China
| | - Yuxin Hu
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Xiaozhi Cao
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Congyu Liao
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Huihui Ye
- State Key Laboratory of Extreme Photonics and Instrumentation, College of Optical Science and Engineering, Zhejiang University, Hangzhou, China
| | - Qiyuan Tian
- School of Biomedical Engineering, Tsinghua University, Beijing, China.
- Tsinghua Laboratory of Brain and Intelligence, Tsinghua University, Beijing, China.
| |
Collapse
|
2
|
Kilic T, Liebig P, Demirel OB, Herrler J, Nagel AM, Ugurbil K, Akçakaya M. Unsupervised deep learning with convolutional neural networks for static parallel transmit design: A retrospective study. Magn Reson Med 2024; 91:2498-2507. [PMID: 38247050 PMCID: PMC10997461 DOI: 10.1002/mrm.30014] [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: 04/07/2023] [Revised: 12/28/2023] [Accepted: 12/29/2023] [Indexed: 01/23/2024]
Abstract
PURPOSE To mitigateB 1 + $$ {B}_1^{+} $$ inhomogeneity at 7T for multi-channel transmit arrays using unsupervised deep learning with convolutional neural networks (CNNs). METHODS Deep learning parallel transmit (pTx) pulse design has received attention, but such methods have relied on supervised training and did not use CNNs for multi-channelB 1 + $$ {B}_1^{+} $$ maps. In this work, we introduce an alternative approach that facilitates the use of CNNs with multi-channelB 1 + $$ {B}_1^{+} $$ maps while performing unsupervised training. The multi-channelB 1 + $$ {B}_1^{+} $$ maps are concatenated along the spatial dimension to enable shift-equivariant processing amenable to CNNs. Training is performed in an unsupervised manner using a physics-driven loss function that minimizes the discrepancy of the Bloch simulation with the target magnetization, which eliminates the calculation of reference transmit RF weights. The training database comprises 3824 2D sagittal, multi-channelB 1 + $$ {B}_1^{+} $$ maps of the healthy human brain from 143 subjects.B 1 + $$ {B}_1^{+} $$ data were acquired at 7T using an 8Tx/32Rx head coil. The proposed method is compared to the unregularized magnitude least-squares (MLS) solution for the target magnetization in static pTx design. RESULTS The proposed method outperformed the unregularized MLS solution for RMS error and coefficient-of-variation and had comparable energy consumption. Additionally, the proposed method did not show local phase singularities leading to distinct holes in the resulting magnetization unlike the unregularized MLS solution. CONCLUSION Proposed unsupervised deep learning with CNNs performs better than unregularized MLS in static pTx for speed and robustness.
Collapse
Affiliation(s)
- Toygan Kilic
- Electrical and Computer Engineering, University of Minnesota, Minneapolis, Minnesota, USA
- Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, Minnesota, USA
| | | | - Omer Burak Demirel
- Electrical and Computer Engineering, University of Minnesota, Minneapolis, Minnesota, USA
- Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, Minnesota, USA
| | | | - Armin M Nagel
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
- Division of Medical Physics in Radiology, German Cancer Research Centre (DKFZ), Heidelberg, Germany
| | - Kamil Ugurbil
- Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, Minnesota, USA
| | - Mehmet Akçakaya
- Electrical and Computer Engineering, University of Minnesota, Minneapolis, Minnesota, USA
- Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, Minnesota, USA
| |
Collapse
|
3
|
Runderkamp BA, Roos T, van der Zwaag W, Strijkers GJ, Caan MWA, Nederveen AJ. Whole-liver flip-angle shimming at 7 T using parallel-transmit k T -point pulses and Fourier phase-encoded DREAM B 1 + mapping. Magn Reson Med 2024; 91:75-90. [PMID: 37799015 DOI: 10.1002/mrm.29819] [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: 03/16/2023] [Revised: 06/18/2023] [Accepted: 07/13/2023] [Indexed: 10/07/2023]
Abstract
PURPOSE To obtain homogeneous signal throughout the human liver at 7 T. Flip angle (FA) shimming in 7T whole-liver imaging was performed through parallel-transmit kT -point pulses based on subject-specific multichannel absoluteB 1 + $$ {\mathrm{B}}_1^{+} $$ maps from Fourier phase-encoded dual refocusing echo acquisition mode (PE-DREAM). METHODS The optimal number of Fourier phase-encoding steps for PE-DREAMB 1 + $$ {\mathrm{B}}_1^{+} $$ mapping was determined for a 7T eight-channel parallel-transmission system. FA shimming experiments were performed in the liver of 7 healthy subjects with varying body mass index. In these subjects, firstB 0 $$ {\mathrm{B}}_0 $$ shimming and Fourier PE-DREAMB 1 + $$ {\mathrm{B}}_1^{+} $$ mapping were performed. Subsequently, three small-flip-angle 3D gradient-echo scans were acquired, comparing a circularly polarized (CP) mode, a phase shim, and a kT -point pulse. Resulting homogeneity was assessed and compared with estimated FA maps and distributions. RESULTS Fourier PE-DREAM with 13 phase-encoding steps resulted in a good tradeoff betweenB 1 + $$ {\mathrm{B}}_1^{+} $$ accuracy and scan time. Lower coefficient of variation values (average [min-max] across subjects) of the estimated FA in the volume of interest were observed using kT -points (7.4 [6.6%-8.0%]), compared with phase shimming (18.8 [12.9%-23.4%], p < 0.001) and CP (43.2 [39.4%-47.1%], p < 0.001). kT -points delivered whole-liver images with the nominal FA and the highest degree of homogeneity. CP and phase shimming resulted in either inaccurate or imprecise FA distributions. Here, locations having suboptimal FA in the estimated FA maps corresponded to liver areas suffering from inconsistent signal intensity and T1 -weighting in the gradient-echo scans. CONCLUSION Homogeneous whole-liver 3D gradient-echo acquisitions at 7 T can be obtained with eight-channel kT -point pulses calculated based on subject-specific multichannel absolute Fourier PE-DREAMB 1 + $$ {\mathrm{B}}_1^{+} $$ maps.
Collapse
Affiliation(s)
- Bobby A Runderkamp
- Department of Radiology and Nuclear Medicine, Amsterdam University Medical Centers, University of Amsterdam, Amsterdam, the Netherlands
| | - Thomas Roos
- Spinoza Centre for Neuroimaging, Royal Netherlands Academy for Arts and Sciences (KNAW), Amsterdam, the Netherlands
- High-Field Research Group, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Wietske van der Zwaag
- Spinoza Centre for Neuroimaging, Royal Netherlands Academy for Arts and Sciences (KNAW), Amsterdam, the Netherlands
- Computational and Cognitive Neuroscience and Neuroimaging, Netherlands Institute for Neuroscience, KNAW, Amsterdam, the Netherlands
| | - Gustav J Strijkers
- Department of Biomedical Engineering and Physics, Amsterdam University Medical Centers, University of Amsterdam, Amsterdam Cardiovascular Sciences, Amsterdam, the Netherlands
| | - Matthan W A Caan
- Department of Biomedical Engineering and Physics, Amsterdam University Medical Centers, University of Amsterdam, Amsterdam Cardiovascular Sciences, Amsterdam, the Netherlands
| | - Aart J Nederveen
- Department of Radiology and Nuclear Medicine, Amsterdam University Medical Centers, University of Amsterdam, Amsterdam, the Netherlands
| |
Collapse
|
4
|
Herrler J, Williams SN, Liebig P, Ding B, McElhinney P, Allwood-Spiers S, Meixner CR, Gunamony S, Maier A, Dörfler A, Gumbrecht R, Porter DA, Nagel AM. The effects of RF coils and SAR supervision strategies for clinically applicable nonselective parallel-transmit pulses at 7 T. Magn Reson Med 2023; 89:1888-1900. [PMID: 36622945 DOI: 10.1002/mrm.29569] [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/02/2022] [Revised: 12/13/2022] [Accepted: 12/14/2022] [Indexed: 01/10/2023]
Abstract
PURPOSE To investigate the effects of using different parallel-transmit (pTx) head coils and specific absorption rate (SAR) supervision strategies on pTx pulse design for ultrahigh-field MRI using a 3D-MPRAGE sequence. METHODS The PTx universal pulses (UPs) and fast online-customized (FOCUS) pulses were designed with pre-acquired data sets (B0 , B1 + maps, specific absorption rate [SAR] supervision data) from two different 8 transmit/32 receive head coils on two 7T whole-body MR systems. For one coil, the SAR supervision model consisted of per-channel RF power limits. In the other coil, SAR estimations were done with both per-channel RF power limits as well as virtual observation points (VOPs) derived from electromagnetic field (EMF) simulations using three virtual human body models at three different positions. All pulses were made for nonselective excitation and inversion and evaluated on 132 B0 , B1 + , and SAR supervision datasets obtained with one coil and 12 from the other. At both sites, 3 subjects were examined using MPRAGE sequences that used UP/FOCUS pulses generated for both coils. RESULTS For some subjects, the UPs underperformed when simulated on a different coil from which they were derived, whereas FOCUS pulses still showed acceptable performance in that case. FOCUS inversion pulses outperformed adiabatic pulses when scaled to the same local SAR level. For the self-built coil, the use of VOPs showed reliable overestimation compared with the ground-truth EMF simulations, predicting about 52% lower local SAR for inversion pulses compared with per-channel power limits. CONCLUSION FOCUS inversion pulses offer a low-SAR alternative to adiabatic pulses and benefit from using EMF-based VOPs for SAR estimation.
Collapse
Affiliation(s)
- Jürgen Herrler
- Department of Neuroradiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany.,Siemens Healthcare, Erlangen, Germany
| | | | | | | | - Paul McElhinney
- Imaging Center of Excellence, University of Glasgow, Glasgow, UK
| | | | - Christian R Meixner
- Siemens Healthcare, Erlangen, Germany.,Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Shajan Gunamony
- Imaging Center of Excellence, University of Glasgow, Glasgow, UK.,MR CoilTech, Glasgow, UK
| | - Andreas Maier
- Department of Computer Science, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Arnd Dörfler
- Department of Neuroradiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | | | - David A Porter
- Imaging Center of Excellence, University of Glasgow, Glasgow, UK
| | - Armin M Nagel
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany.,Division of Medical Physics in Radiology, German Cancer Research Center, Heidelberg, Germany
| |
Collapse
|
5
|
Williams SN, McElhinney P, Gunamony S. Ultra-high field MRI: parallel-transmit arrays and RF pulse design. Phys Med Biol 2023; 68. [PMID: 36410046 DOI: 10.1088/1361-6560/aca4b7] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 11/21/2022] [Indexed: 11/22/2022]
Abstract
This paper reviews the field of multiple or parallel radiofrequency (RF) transmission for magnetic resonance imaging (MRI). Currently the use of ultra-high field (UHF) MRI at 7 tesla and above is gaining popularity, yet faces challenges with non-uniformity of the RF field and higher RF power deposition. Since its introduction in the early 2000s, parallel transmission (pTx) has been recognized as a powerful tool for accelerating spatially selective RF pulses and combating the challenges associated with RF inhomogeneity at UHF. We provide a survey of the types of dedicated RF coils used commonly for pTx and the important modeling of the coil behavior by electromagnetic (EM) field simulations. We also discuss the additional safety considerations involved with pTx such as the specific absorption rate (SAR) and how to manage them. We then describe the application of pTx with RF pulse design, including a practical guide to popular methods. Finally, we conclude with a description of the current and future prospects for pTx, particularly its potential for routine clinical use.
Collapse
Affiliation(s)
- Sydney N Williams
- Imaging Centre of Excellence, University of Glasgow, Glasgow, United Kingdom
| | - Paul McElhinney
- Imaging Centre of Excellence, University of Glasgow, Glasgow, United Kingdom
| | - Shajan Gunamony
- Imaging Centre of Excellence, University of Glasgow, Glasgow, United Kingdom.,MR CoilTech Limited, Glasgow, United Kingdom
| |
Collapse
|
6
|
Yetisir F, Poser BA, Grant PE, Adalsteinsson E, Wald LL, Guerin B. Parallel transmission 2D RARE imaging at 7T with transmit field inhomogeneity mitigation and local SAR control. Magn Reson Imaging 2022; 93:87-96. [PMID: 35940379 PMCID: PMC9789791 DOI: 10.1016/j.mri.2022.08.006] [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/04/2022] [Revised: 07/15/2022] [Accepted: 08/02/2022] [Indexed: 12/26/2022]
Abstract
PURPOSE We develop and test a parallel transmit (pTx) pulse design framework to mitigate transmit field inhomogeneity with control of local specific absorption rate (SAR) in 2D rapid acquisition with relaxation enhancement (RARE) imaging at 7T. METHODS We design large flip angle RF pulses with explicit local SAR constraints by numerical simulation of the Bloch equations. Parallel computation and analytical expressions for the Jacobian and the Hessian matrices are employed to reduce pulse design time. The refocusing-excitation "spokes" pulse pairs are designed to satisfy the Carr-Purcell-Meiboom-Gill (CPMG) condition using a combined magnitude least squares-least squares approach. RESULTS In a simulated dataset, the proposed approach reduced peak local SAR by up to 56% for the same level of refocusing uniformity error and reduced refocusing uniformity error by up to 59% (from 32% to 7%) for the same level of peak local SAR compared to the circularly polarized birdcage mode of the pTx array. Using explicit local SAR constraints also reduced peak local SAR by up to 46% compared to an RF peak power constrained design. The excitation and refocusing uniformity error were reduced from 20%-33% to 4%-6% in single slice phantom experiments. Phantom experiments demonstrated good agreement between the simulated excitation and refocusing uniformity profiles and experimental image shading. CONCLUSION PTx-designed excitation and refocusing CPMG pulse pairs can mitigate transmit field inhomogeneity in the 2D RARE sequence. Moreover, local SAR can be decreased significantly using pTx, potentially leading to better slice coverage, enabling larger flip angles or faster imaging.
Collapse
Affiliation(s)
- Filiz Yetisir
- Fetal-Neonatal Neuroimaging & Developmental Science Center, Boston Children's Hospital, Harvard Medical School, 300 Longwood Avenue, Boston, MA 02115, USA.
| | - Benedikt A Poser
- Department of Cognitive Neuroscience, Faculty of Psychology and Neuroscience, Maastricht University, Universiteitssingel 40, 6229 ER Maastricht, the Netherlands
| | - P Ellen Grant
- Fetal-Neonatal Neuroimaging & Developmental Science Center, Boston Children's Hospital, Harvard Medical School, 300 Longwood Avenue, Boston, MA 02115, USA; Department of Pediatrics, Harvard Medical School, 25 Shattuck Street, Boston, MA 02115, USA; Department of Radiology, Harvard Medical School, 25 Shattuck Street, Boston, MA 02115, USA
| | - Elfar Adalsteinsson
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, 77 MA Avenue, Cambridge, MA 02139, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, 77 MA Avenue, Cambridge, MA 02139, USA; Institute for Medical Engineering and Science, Massachusetts Institute of Technology, 77 MA Avenue, Cambridge, MA 02139, USA
| | - Lawrence L Wald
- Department of Radiology, Harvard Medical School, 25 Shattuck Street, Boston, MA 02115, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, 77 MA Avenue, Cambridge, MA 02139, USA; Athinoula A. Martinos Center for Biomedical Imaging, MA General Hospital, 149 13th Street, Charlestown, MA 02129, USA
| | - Bastien Guerin
- Department of Radiology, Harvard Medical School, 25 Shattuck Street, Boston, MA 02115, USA; Athinoula A. Martinos Center for Biomedical Imaging, MA General Hospital, 149 13th Street, Charlestown, MA 02129, USA
| |
Collapse
|
7
|
Mooiweer R, Clark IA, Maguire EA, Callaghan MF, Hajnal JV, Malik SJ. Universal pulses for homogeneous excitation using single channel coils. Magn Reson Imaging 2022; 92:180-186. [PMID: 35820546 DOI: 10.1016/j.mri.2022.07.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Revised: 06/24/2022] [Accepted: 07/07/2022] [Indexed: 11/29/2022]
Abstract
PURPOSE Universal Pulses (UPs) are excitation pulses that reduce the flip angle inhomogeneity in high field MRI systems without subject-specific optimization, originally developed for parallel transmit (PTX) systems at 7 T. We investigated the potential benefits of UPs for single channel (SC) transmit systems at 3 T, which are widely used for clinical and research imaging, and for which flip angle inhomogeneity can still be problematic. METHODS SC-UPs were designed using a spiral nonselective k-space trajectory for brain imaging at 3 T using transmit field maps (B1+) and off-resonance maps (B0) acquired on two different scanner types: a 'standard' single channel transmit system and a system with a PTX body coil. The effect of training group size was investigated using data (200 subjects) from the standard system. The PTX system was used to compare SC-UPs to PTX-UPs (15 subjects). In two additional subjects, prospective imaging using SC-UP was studied. RESULTS Average flip angle homogeneity error fell from 9.5 ± 0.5 % for 'default' excitation to 3.0 ± 0.6 % using SC-UPs trained over 50 subjects. Performance of the UPs was found to steadily improve as training group size increased, but stabilized after ~15 subjects. On the PTX-enabled system, SC-UPs again outperformed default excitation in simulations (4.8 ± 0.6 % error versus 10.6 ± 0.8 % respectively) though greater homogenization could be achieved with PTX-UPs (3.9 ± 0.6 %) and personalized pulses (SC-PP 3.6 ± 1.0 %, PTX-PP 2.9 ± 0.6 %). MP-RAGE imaging using SC-UP resulted in greater separation between grey and white matter signal intensities than default excitation. CONCLUSIONS SC-UPs can improve excitation homogeneity in standard 3 T systems without further calibration and could be used instead of a default excitation pulse for nonselective neuroimaging at 3 T.
Collapse
Affiliation(s)
- Ronald Mooiweer
- School of Biomedical Engineering and Imaging Sciences, Faculty of Life Sciences and Medicine, King's College London, London, United Kingdom; Center for the Developing Brain, School of Biomedical Engineering and Imaging Sciences, King's College London, St. Thomas' Hospital, London, United Kingdom; MR Research Collaborations, Siemens Healthcare Limited, Camberley, United Kingdom
| | - Ian A Clark
- Wellcome Centre for Human Neuroimaging, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom
| | - Eleanor A Maguire
- Wellcome Centre for Human Neuroimaging, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom
| | - Martina F Callaghan
- Wellcome Centre for Human Neuroimaging, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom
| | - Joseph V Hajnal
- School of Biomedical Engineering and Imaging Sciences, Faculty of Life Sciences and Medicine, King's College London, London, United Kingdom; Center for the Developing Brain, School of Biomedical Engineering and Imaging Sciences, King's College London, St. Thomas' Hospital, London, United Kingdom
| | - Shaihan J Malik
- School of Biomedical Engineering and Imaging Sciences, Faculty of Life Sciences and Medicine, King's College London, London, United Kingdom; Center for the Developing Brain, School of Biomedical Engineering and Imaging Sciences, King's College London, St. Thomas' Hospital, London, United Kingdom.
| |
Collapse
|
8
|
Okada T, Fujimoto K, Fushimi Y, Akasaka T, Thuy DHD, Shima A, Sawamoto N, Oishi N, Zhang Z, Funaki T, Nakamoto Y, Murai T, Miyamoto S, Takahashi R, Isa T. Neuroimaging at 7 Tesla: a pictorial narrative review. Quant Imaging Med Surg 2022; 12:3406-3435. [PMID: 35655840 PMCID: PMC9131333 DOI: 10.21037/qims-21-969] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Accepted: 02/05/2022] [Indexed: 01/26/2024]
Abstract
Neuroimaging using the 7-Tesla (7T) human magnetic resonance (MR) system is rapidly gaining popularity after being approved for clinical use in the European Union and the USA. This trend is the same for functional MR imaging (MRI). The primary advantages of 7T over lower magnetic fields are its higher signal-to-noise and contrast-to-noise ratios, which provide high-resolution acquisitions and better contrast, making it easier to detect lesions and structural changes in brain disorders. Another advantage is the capability to measure a greater number of neurochemicals by virtue of the increased spectral resolution. Many structural and functional studies using 7T have been conducted to visualize details in the white matter and layers of the cortex and hippocampus, the subnucleus or regions of the putamen, the globus pallidus, thalamus and substantia nigra, and in small structures, such as the subthalamic nucleus, habenula, perforating arteries, and the perivascular space, that are difficult to observe at lower magnetic field strengths. The target disorders for 7T neuroimaging range from tumoral diseases to vascular, neurodegenerative, and psychiatric disorders, including Alzheimer's disease, Parkinson's disease, multiple sclerosis, epilepsy, major depressive disorder, and schizophrenia. MR spectroscopy has also been used for research because of its increased chemical shift that separates overlapping peaks and resolves neurochemicals more effectively at 7T than a lower magnetic field. This paper presents a narrative review of these topics and an illustrative presentation of images obtained at 7T. We expect 7T neuroimaging to provide a new imaging biomarker of various brain disorders.
Collapse
Affiliation(s)
- Tomohisa Okada
- Human Brain Research Center, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Koji Fujimoto
- Department of Real World Data Research and Development, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Yasutaka Fushimi
- Department of Diagnostic Imaging and Nuclear Medicine, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Thai Akasaka
- Human Brain Research Center, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Dinh H. D. Thuy
- Human Brain Research Center, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Atsushi Shima
- Human Brain Research Center, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Nobukatsu Sawamoto
- Department of Human Health Sciences, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Naoya Oishi
- Medial Innovation Center, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Zhilin Zhang
- Department of Psychiatry, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Takeshi Funaki
- Department of Neurosurgery, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Yuji Nakamoto
- Department of Diagnostic Imaging and Nuclear Medicine, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Toshiya Murai
- Department of Psychiatry, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Susumu Miyamoto
- Department of Neurosurgery, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Ryosuke Takahashi
- Department of Neurology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Tadashi Isa
- Human Brain Research Center, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| |
Collapse
|
9
|
Plumley A, Watkins L, Treder M, Liebig P, Murphy K, Kopanoglu E. Rigid motion-resolved B1+ prediction using deep learning for real-time parallel-transmission pulse design. Magn Reson Med 2022; 87:2254-2270. [PMID: 34958134 PMCID: PMC7613077 DOI: 10.1002/mrm.29132] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2021] [Revised: 11/30/2021] [Accepted: 12/03/2021] [Indexed: 11/19/2022]
Abstract
PURPOSE Tailored parallel-transmit (pTx) pulses produce uniform excitation profiles at 7 T, but are sensitive to head motion. A potential solution is real-time pulse redesign. A deep learning framework is proposed to estimate pTx B 1 + distributions following within-slice motion, which can then be used for tailored pTx pulse redesign. METHODS Using simulated data, conditional generative adversarial networks were trained to predict B 1 + distributions in the head following a displacement. Predictions were made for two virtual body models that were not included in training. Predicted maps were compared with ground-truth (simulated, following motion) B1 maps. Tailored pTx pulses were designed using B1 maps at the original position (simulated, no motion) and evaluated using simulated B1 maps at displaced position (ground-truth maps) to quantify motion-related excitation error. A second pulse was designed using predicted maps (also evaluated on ground-truth maps) to investigate improvement offered by the proposed method. RESULTS Predicted B 1 + maps corresponded well with ground-truth maps. Error in predicted maps was lower than motion-related error in 99% and 67% of magnitude and phase evaluations, respectively. Worst-case flip-angle normalized RMS error due to motion (76% of target flip angle) was reduced by 59% when pulses were redesigned using predicted maps. CONCLUSION We propose a framework for predicting B 1 + maps online with deep neural networks. Predicted maps can then be used for real-time tailored pulse redesign, helping to overcome head motion-related error in pTx.
Collapse
Affiliation(s)
- Alix Plumley
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom
| | - Luke Watkins
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom
- School of Physics & Astronomy, Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff, United Kingdom
| | - Matthias Treder
- School of Computer Science and Informatics, Cardiff University, Cardiff, United Kingdom
| | | | - Kevin Murphy
- School of Physics & Astronomy, Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff, United Kingdom
| | - Emre Kopanoglu
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom
| |
Collapse
|
10
|
Ma X, Uğurbil K, Wu X. Mitigating transmit‐B
1
artifacts by predicting parallel transmission images with deep learning: A feasibility study using high‐resolution whole‐brain diffusion at 7 Tesla. Magn Reson Med 2022; 88:727-741. [PMID: 35403237 PMCID: PMC9324974 DOI: 10.1002/mrm.29238] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Revised: 02/24/2022] [Accepted: 02/28/2022] [Indexed: 11/12/2022]
Abstract
Purpose To propose a novel deep learning (DL) approach to transmit‐B1 (B1+)‐artifact mitigation without direct use of parallel transmission (pTx), by predicting pTx images from single‐channel transmission (sTx) images. Methods A deep encoder–decoder convolutional neural network was constructed and trained to learn the mapping from sTx to pTx images. The feasibility was demonstrated using 7 T Human‐Connectome Project (HCP)‐style diffusion MRI. The training dataset comprised images acquired on 5 healthy subjects using commercial Nova RF coils. Relevant hyperparameters were tuned with a nested cross‐validation, and the generalization performance evaluated using a regular cross‐validation. Results Our DL method effectively improved the image quality for sTx images by restoring the signal dropout, with quality measures (including normalized root‐mean‐square error, peak SNR, and structural similarity index measure) improved in most brain regions. The improved image quality was translated into improved performances for diffusion tensor imaging analysis; our method improved accuracy for fractional anisotropy and mean diffusivity estimations, reduced the angular errors of principal eigenvectors, and improved the fiber orientation delineation relative to sTx images. Moreover, the final DL model trained on data of all 5 subjects was successfully used to predict pTx images for unseen new subjects (randomly selected from the 7 T HCP database), effectively recovering the signal dropout and improving color‐coded fractional anisotropy maps with largely reduced noise levels. Conclusion The proposed DL method has potential to provide images with reduced B1+ artifacts in healthy subjects even when pTx resources are inaccessible on the user side.
Collapse
Affiliation(s)
- Xiaodong Ma
- Center for Magnetic Resonance Research, Radiology, Medical School University of Minnesota Minneapolis Minnesota USA
| | - Kâmil Uğurbil
- Center for Magnetic Resonance Research, Radiology, Medical School University of Minnesota Minneapolis Minnesota USA
| | - Xiaoping Wu
- Center for Magnetic Resonance Research, Radiology, Medical School University of Minnesota Minneapolis Minnesota USA
| |
Collapse
|
11
|
Eberhardt B, Poser BA, Shah NJ, Felder J. B1 field map synthesis with generative deep learning used in the design of parallel-transmit RF pulses for ultra-high field MRI. Z Med Phys 2022; 32:334-345. [PMID: 35144850 PMCID: PMC9948838 DOI: 10.1016/j.zemedi.2021.12.003] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2021] [Revised: 11/27/2021] [Accepted: 12/23/2021] [Indexed: 10/19/2022]
Abstract
Spoke trajectory parallel transmit (pTX) excitation in ultra-high field MRI enables B1+ inhomogeneities arising from the shortened RF wavelength in biological tissue to be mitigated. To this end, current RF excitation pulse design algorithms either employ the acquisition of field maps with subsequent non-linear optimization or a universal approach applying robust pre-computed pulses. We suggest and evaluate an intermediate method that uses a subset of acquired field maps combined with generative machine learning models to reduce the pulse calibration time while offering more tailored excitation than robust pulses (RP). The possibility of employing image-to-image translation and semantic image synthesis machine learning models based on generative adversarial networks (GANs) to deduce the missing field maps is examined. Additionally, an RF pulse design that employs a predictive machine learning model to find solutions for the non-linear (two-spokes) pulse design problem is investigated. As a proof of concept, we present simulation results obtained with the suggested machine learning approaches that were trained on a limited data-set, acquired in vivo. The achieved excitation homogeneity based on a subset of half of the B1+ maps acquired in the calibration scans and half of the B1+ maps synthesized with GANs is comparable with state of the art pulse design methods when using the full set of calibration data while halving the total calibration time. By employing RP dictionaries or machine-learning RF pulse predictions, the total calibration time can be reduced significantly as these methods take only seconds or milliseconds per slice, respectively.
Collapse
Affiliation(s)
- Boris Eberhardt
- Institute of Neuroscience and Medicine 4, Forschungszentrum Jülich, Jüich, Germany; RWTH Aachen University, Aachen, Germany.
| | - Benedikt A. Poser
- Department of Cognitive Neuroscience, Faculty of Psychology and Neuroscience, Maastricht University, Maastricht, Netherlands
| | - N. Jon Shah
- Institute of Neuroscience and Medicine 4, Forschungszentrum Jülich, Jüich, Germany,Institute of Neuroscience and Medicine 11, Forschungszentrum Jülich, Jülich, Germany,Department of Neurology, RWTH Aachen University, Aachen, Germany,JARA-BRAIN, Translational Medicine, Aachen, Germany
| | - Jörg Felder
- Institute of Neuroscience and Medicine 4, Forschungszentrum Jülich, Jüich, Germany; RWTH Aachen University, Aachen, Germany.
| |
Collapse
|
12
|
Le Ster C, Mauconduit F, Massire A, Boulant N, Gras V. Standardized universal pulse: A fast RF calibration approach to improve flip angle accuracy in parallel transmission. Magn Reson Med 2022; 87:2839-2850. [PMID: 35122302 DOI: 10.1002/mrm.29180] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2021] [Revised: 01/10/2022] [Accepted: 01/11/2022] [Indexed: 11/09/2022]
Abstract
PURPOSE In parallel transmission (pTX), subject-tailored RF pulses allow achieving excellent flip angle (FA) accuracy but often require computationally extensive online optimizations, precise characterization of the static field ( Δ B 0 ), and the transmit RF field ( B 1 + ) distributions. This costs time and requires expertise from the MR user. Universal pulses (UPs) have been proposed to reduce this burden, yet, with a penalty in FA accuracy. This study introduces the concept of standardized universal pulses (SUPs), where pulses are designed offline and adjusted to the subject through a fast online calibration scan. METHODS A SUP is designed offline using a so-called standardized database, wherein each B 1 + map has been normalized to a reference transmit RF field distribution. When scanning a new subject, a 3-slice B 1 + acquisition (scan time < 10 s) is performed and used to adjust the SUP to the subject through a linear transform. SUP performance was assessed at 7T with simulations by computing the FA-normalized root mean square error (FA-NRMSE) and the FA pattern stability as measured by the average and coefficient of variation of the FA across 15 control subjects, along with in vivo experiments using an MP2RAGE sequence implementing the SUP variant for the FLASH readout. RESULTS Adjusted SUP improved the FA-NRMSE (8.8 % for UP vs. 7.1 % for adjusted SUP). Experimentally in vivo, this translated in an improved signal homogeneity and more accurate T 1 quantification using MP2RAGE. CONCLUSION The proposed SUP approach improves excitation accuracy (FA-NRMSE) while preserving the same offline pulse design principle as offered by UPs.
Collapse
Affiliation(s)
- Caroline Le Ster
- NeuroSpin, CEA, CNRS, BAOBAB, Université Paris-Saclay, Gif-Sur-Yvette, France
| | - Franck Mauconduit
- NeuroSpin, CEA, CNRS, BAOBAB, Université Paris-Saclay, Gif-Sur-Yvette, France
| | | | - Nicolas Boulant
- NeuroSpin, CEA, CNRS, BAOBAB, Université Paris-Saclay, Gif-Sur-Yvette, France
| | - Vincent Gras
- NeuroSpin, CEA, CNRS, BAOBAB, Université Paris-Saclay, Gif-Sur-Yvette, France
| |
Collapse
|
13
|
Berrington A, Považan M, Mirfin C, Bawden S, Park YW, Marsh DC, Bowtell R, Gowland PA. Calibration-free regional RF shims for MRS. Magn Reson Med 2021; 86:611-624. [PMID: 33749010 DOI: 10.1002/mrm.28749] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2020] [Revised: 01/07/2021] [Accepted: 02/04/2021] [Indexed: 12/27/2022]
Abstract
PURPOSE Achieving a desired RF transmit field ( B 1 + ) in small regions of interest is critical for single-voxel MRS at ultrahigh field. Radio-frequency (RF) shimming, using parallel transmission, requires B 1 + mapping and optimization, which limits its ease of use. This work aimed to generate calibration-free RF shims for predefined target regions of interest, which can be applied to any participant, to produce a desired absolute magnitude B 1 + (| B 1 + |). METHODS The RF shims were found offline by joint optimization on a database comprising B 1 + maps from 11 subjects, considering regions of interest in occipital cortex, hippocampus and posterior cingulate, as well as whole brain. The | B 1 + | achieved was compared with a tailored shimming approach, and MR spectra were acquired using tailored and calibration-free shims in 4 participants. Global and local 10g specific-absorption-rate deposition were estimated using Duke and Ella dielectric models. RESULTS There was no difference in the mean | B 1 + | produced using calibration-free versus tailored RF shimming in the occipital cortex (p = .15), hippocampus (p = .5), or posterior cingulate (p = .98), although differences were observed in the RMS error | B 1 + |. Spectra acquired using calibration-free shims had similar SNR and low residual water signal. Under identical power settings, specific-absorption-rate deposition was lower compared with operating in quadrature mode. For example, the total head specific absorption rate was around 35% less for the occipital cortex. CONCLUSION This work demonstrates that static RF shims, optimized offline for small regions, avoid the need for B 1 + mapping and optimization for each region of interest and participant. Furthermore, power settings may be increased when using calibration-free shims, to better take advantage of RF shimming.
Collapse
Affiliation(s)
- Adam Berrington
- Sir Peter Mansfield Imaging Centre, School of Physics and Astronomy, University of Nottingham, Nottingham, UK
| | - Michal Považan
- Russell H. Morgan Department of Radiology, Johns Hopkins University School of Medicine, Baltimore, MD, USA.,F. M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, USA
| | - Christopher Mirfin
- Sir Peter Mansfield Imaging Centre, School of Physics and Astronomy, University of Nottingham, Nottingham, UK
| | - Stephen Bawden
- Sir Peter Mansfield Imaging Centre, School of Physics and Astronomy, University of Nottingham, Nottingham, UK
| | - Young Woo Park
- Center for Magnetic Resonance Research, Department of Radiology, University of Minnesota Medical School, Minneapolis, MN, USA
| | - Daniel C Marsh
- Sir Peter Mansfield Imaging Centre, School of Physics and Astronomy, University of Nottingham, Nottingham, UK
| | - Richard Bowtell
- Sir Peter Mansfield Imaging Centre, School of Physics and Astronomy, University of Nottingham, Nottingham, UK
| | - Penny A Gowland
- Sir Peter Mansfield Imaging Centre, School of Physics and Astronomy, University of Nottingham, Nottingham, UK
| |
Collapse
|
14
|
Vinding MS, Aigner CS, Schmitter S, Lund TE. DeepControl: 2DRF pulses facilitating B 1 + inhomogeneity and B 0 off-resonance compensation in vivo at 7 T. Magn Reson Med 2021; 85:3308-3317. [PMID: 33480029 DOI: 10.1002/mrm.28667] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2020] [Revised: 12/09/2020] [Accepted: 12/09/2020] [Indexed: 01/26/2023]
Abstract
PURPOSE Rapid 2DRF pulse design with subject-specific B 1 + inhomogeneity and B0 off-resonance compensation at 7 T predicted from convolutional neural networks is presented. METHODS The convolution neural network was trained on half a million single-channel transmit 2DRF pulses optimized with an optimal control method using artificial 2D targets, B 1 + and B0 maps. Predicted pulses were tested in a phantom and in vivo at 7 T with measured B 1 + and B0 maps from a high-resolution gradient echo sequence. RESULTS Pulse prediction by the trained convolutional neural network was done on the fly during the MR session in approximately 9 ms for multiple hand-drawn regions of interest and the measured B 1 + and B0 maps. Compensation of B 1 + inhomogeneity and B0 off-resonances has been confirmed in the phantom and in vivo experiments. The reconstructed image data agree well with the simulations using the acquired B 1 + and B0 maps, and the 2DRF pulse predicted by the convolutional neural networks is as good as the conventional RF pulse obtained by optimal control. CONCLUSION The proposed convolutional neural network-based 2DRF pulse design method predicts 2DRF pulses with an excellent excitation pattern and compensated B 1 + and B0 variations at 7 T. The rapid 2DRF pulse prediction (9 ms) enables subject-specific high-quality 2DRF pulses without the need to run lengthy optimizations.
Collapse
Affiliation(s)
- Mads Sloth Vinding
- Center of Functionally Integrative Neuroscience (CFIN), Department of Clinical Medicine, Faculty of Health, Aarhus University, Aarhus N, Denmark
| | | | - Sebastian Schmitter
- Physikalisch-Technische Bundesanstalt, Braunschweig and Berlin, Germany.,Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, Minnesota, USA
| | - Torben Ellegaard Lund
- Center of Functionally Integrative Neuroscience (CFIN), Department of Clinical Medicine, Faculty of Health, Aarhus University, Aarhus N, Denmark
| |
Collapse
|
15
|
Herrler J, Liebig P, Gumbrecht R, Ritter D, Schmitter S, Maier A, Schmidt M, Uder M, Doerfler A, Nagel AM. Fast online-customized (FOCUS) parallel transmission pulses: A combination of universal pulses and individual optimization. Magn Reson Med 2021; 85:3140-3153. [PMID: 33400302 DOI: 10.1002/mrm.28643] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Revised: 11/03/2020] [Accepted: 11/23/2020] [Indexed: 12/12/2022]
Abstract
PURPOSE To mitigate spatial flip angle (FA) variations under strict specific absorption rate (SAR) constraints for ultra-high field MRI using a combination of universal parallel transmit (pTx) pulses and fast subject-specific optimization. METHODS Data sets consisting of B0 , B 1 + maps, and virtual observation point (VOP) data were acquired from 72 subjects (study groups of 48/12 healthy Europeans/Asians and 12 Europeans with pathological or incidental findings) using an 8Tx/32Rx head coil on a 7T whole-body MR system. Combined optimization values (COV) were defined as combination of spiral-nonselective (SPINS) trajectory parameters and an energy regularization weight. A set of COV was optimized universally by simulating the individual RF pulse optimizations of 12 training data sets (healthy Europeans). Subsequently, corresponding universal pulses (UPs) were calculated. Using COV and UPs, individually optimized pulses (IOPs) were calculated during the sequence preparation phase (maximum 15 s). Two different UPs and IOPs were evaluated by calculating their normalized root-mean-square error (NRMSE) of the FA and SAR in simulations of all data sets. Seven additional subjects were examined using an MPRAGE sequence that uses the designed pTx excitation pulses and a conventional adiabatic inversion. RESULTS All pTx pulses resulted in decreased mean NRMSE compared to a circularly polarized (CP) pulse (CP = ~28%, UPs = ~17%, and IOPs = ~12%). UPs and IOPs improved homogeneity for all subjects. Differences in NRMSE between study groups were much lower than differences between different pulse types. CONCLUSION UPs can be used to generate fast online-customized (FOCUS) pulses gaining lower NRMSE and/or lower SAR values.
Collapse
Affiliation(s)
- Jürgen Herrler
- Department of Neuroradiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | | | | | | | - Sebastian Schmitter
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig and Berlin, Germany
| | - Andreas Maier
- Department of Computer Science, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Germany
| | - Manuel Schmidt
- Department of Neuroradiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Michael Uder
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Arnd Doerfler
- Department of Neuroradiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Armin M Nagel
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany.,Institute of Medical Physics, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany.,Division of Medical Physics in Radiology, German Cancer Research Centre (DKFZ), Heidelberg, Germany
| |
Collapse
|