Short-T2 MRI: Principles and recent advances

Short-T2 MRI: Principles and recent advances

Journal Pre-Proof Short-T2 MRI: Principles and recent advances Markus Weiger, Klaas P. Pruessmann PII: DOI: Reference: S0079-6565(19)30033-0 https://...

3MB Sizes 0 Downloads 68 Views

Journal Pre-Proof Short-T2 MRI: Principles and recent advances Markus Weiger, Klaas P. Pruessmann PII: DOI: Reference:

S0079-6565(19)30033-0 https://doi.org/10.1016/j.pnmrs.2019.07.001 JPNMRS 1483

To appear in:

Progress in Nuclear Magnetic Resonance Spectroscopy

Received Date: Accepted Date:

4 June 2019 26 July 2019

Please cite this article as: M. Weiger, K.P. Pruessmann, Short-T2 MRI: Principles and recent advances, Progress in Nuclear Magnetic Resonance Spectroscopy (2019), doi: https://doi.org/10.1016/j.pnmrs.2019.07.001

This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

© 2019 Elsevier B.V. All rights reserved.

JOURNAL PRE-PROOF Short-T2 MRI: Principles and recent advances Markus Weiger*,[email protected], Klaas P. Pruessmann

Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland *Corresponding

author at: Institute for Biomedical Engineering, ETH Zurich and University of

O

O

F

Zurich, Gloriastrasse 35, CH-8092 Zurich, Switzerland.

PR

Graphical abstract Highlights

Overarching concepts structure the variety of short-T2 MRI techniques

E-

Relaxation in a nutshell

PR

Detailed description of ZTE imaging and its variants

Dedicated hardware requirements for all scanner components

Abstract

AL

Comprehensive literature overview including T2 and T2* values of potential imaging targets

R N

Among current modalities of biomedical and diagnostic imaging, MRI stands out by virtue of its versatile contrast obtained without ionizing radiation. However, in various cases, e.g., water protons in tissues such as bone, tendon, and lung, MRI performance is limited by the rapid decay of resonance

U

signals associated with short transverse relaxation times T2 or T2*. Efforts to address this shortcoming have led to a variety of specialized short-T2 techniques. Recent progress in this field expands the

JO

choice of methods and prompts fresh considerations with regard to instrumentation, data acquisition, and signal processing. In this review, the current status of short-T2 MRI is surveyed. In an attempt to

structure the growing range of techniques, the presentation highlights overarching concepts and basic methodological options. The most frequently used approaches are described in detail, including acquisition strategies, image reconstruction, hardware requirements, means of introducing contrast, sources of artifacts, limitations, and applications. Keywords: Relaxation; Zero echo time (ZTE); Ultra-short echo time (UTE); Single point imaging (SPI); Solid tissues

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

1 Introduction Magnetic resonance imaging (MRI) is a well-established and widely used modality in medical diagnostics due to its ability to depict tissues in three dimensions, with versatile contrast, and without ionizing radiation. It is based on spatial encoding and detection of macroscopic transverse nuclear magnetization. The latter is of limited lifetime and usually decays approximately exponentially with

F

time constants T2 or T2*, depending on the exact nature of the decay. Nuclei of molecules in liquids, in

O

particular water protons in most soft tissues, are characterized by relatively long T2 in the range of tens of milliseconds to seconds, which readily permit extensive spatial encoding by gradient fields and

O

full signal detection.

PR

However, much shorter signal lifetimes can be observed for nuclei in molecules with more restricted mobility. Some tissues, including tendons, ligaments, bone, and teeth, are stringy or solid and owe their short T2 to low water mobility. 1H and other nuclei bound in large molecules such as lipids and

E-

proteins can be even less mobile and thus exhibit yet shorter T2. However, loss of macroscopic transverse magnetization can also stem from dephasing due to B0 inhomogeneity at the microscopic

PR

scale. A special case is lung tissue, in which rapid decay of the water signal associated with strong, microscopic variation in magnetic susceptibility is further influenced by diffusion. Generally, the exact mechanism of coherence loss is relevant when considering methodology. However, in the context of imaging, in particular in vivo, many of the main considerations hold equally for different

AL

types of signal decay. Therefore, unless stated otherwise, all kinds of coherence loss (such as T2 or T2*) will be jointly referred to as T2 decay in the scope of this work.

R N

Rapid signal decay hampers MRI in two ways. In conventional scans, signals with T2 below approximately one millisecond are largely missed during the dead time between MR excitation and detection, causing voids in resulting images. In addition, even if short-T2 signal is still captured, signal decay during data acquisition translates into apodization in k-space, entailing loss of spatial resolution

U

or even artifacts.

JO

To expand the scope of application of MRI to short-T2 tissues, continuing efforts have been made to

image such components by means of dedicated techniques. In fact, the field of short-T2 imaging is

almost as old as MRI itself. In solid-state NMR, short T2 is often addressed by line-narrowing techniques [1-4]. Unfortunately, line-narrowing is barely applicable in vivo since it involves either rapid sample rotation or high levels of radiofrequency (RF) power absorption. Therefore, short-T2 imaging in vivo must primarily rely on sufficient speed of spatial encoding and data acquisition. The chief traditional strategies to this end are single-point imaging (SPI) [5] and ultra-short echo time (UTE) imaging [6], which both employ efficient center-out k-space encoding. Over the past decade, the field has been expanded by zero echo time (ZTE) methods [7], of which an array of

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

implementations is now in use while more variants continue to emerge. Promising results of recent applications have promoted general interest in short-T2 imaging [8] and instrument manufacturers are increasingly following up with commercial implementations of new advances. At present, short-T2 MRI is a diverse field, and it can be a challenge to gain an overview and keep track of ongoing developments. Nevertheless, all of the techniques that are involved have emerged

F

from a finite set of basic ideas that guide sequence design, signal detection, data processing, and

O

hardware requirements. The goal of this article is to review the existing variety of short-T2 MRI methods from such a conceptual point of view, to work out commonalities and key differences

O

between techniques, and to place them in a shared framework. It strives to offer guidance for getting started in short-T2 method development and for choosing the right approach towards a given

PR

application.

The presentation starts with a brief recapitulation of the physics underlying signal decay in MRI. It

E-

then discusses the basic principles of short-T2 imaging, which serve as a framework for a subsequent overview of existing techniques. Particular attention is given to two classes of techniques that are

PR

currently most employed in vivo, the UTE and ZTE approaches. A separate section is dedicated to specific hardware requirements. The final section of this review surveys the realm of current and potential applications.

2.1

Preamble

AL

2 Relaxation in a nutshell

R N

The general topic of relaxation is a complex branch of NMR physics. For the purpose of this article, no attempt is made to survey it fully or at significant depth. The aim of this section is rather to set the scene for MRI with short-lived signals by a brief, simplified review of the chief causes of rapid signal

U

decay. For a thorough treatment of the topic we refer to NMR textbooks [9] as well as original work

JO

[10-12] and pertinent reviews [13-16]. 2.2

Overview

Relaxation is the process that drives a perturbed spin system back to thermal equilibrium with zero transverse

magnetization,

  ,

and

finite

longitudinal

magnetization,

  .

A

phenomenological, mono-exponential description of relaxation is given by the Bloch equations. In the absence of an RF field they reduce to

   



(1)

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI



    





(2)

with relaxation times T1 and T2, which vary greatly across different nuclei and their environments in given tissues and materials. The purpose of relaxation theory is to describe the mechanisms underlying relaxation and to relate T1 and T2 to more fundamental physical properties.

F

In most practical cases, relaxation of nuclear spins is based on two key ingredients: 1) some magnetic

O

interaction of the spins in question with their surroundings that depends on orientation or position of the molecules containing the spins and 2) quasi-random molecular tumbling. Tumbling motion

O

translates interaction into magnetic field fluctuations at the spin locations. These fluctuations cause continuous slight variations in the precession frequency of the spins and in their angle with respect to

PR

the static field B0, leading to irreversible loss of transverse magnetization and recovery of longitudinal magnetization, respectively. Interactions

E-

2.3

Relevant interactions occur within molecules as well as between spins and their wider environment.

PR

To contribute to relaxation, interactions must be spatially anisotropic or inhomogeneous. In particular, intra-molecular interactions must depend on the orientation of the molecule with respect to the B0 field whereas inter-molecular interactions must change with distance or relative orientation. The

AL

usually strongest and therefore most important interactions are: 1) Dipolar coupling: The dipolar field associated with a spin’s magnetic moment contributes to the local field at the sites of nearby spins. Dipolar coupling is a through-space effect that drops rapidly

R N

with distance and therefore acts mostly intra-molecularly. It will be discussed in more detail below. 2) Chemical shift anisotropy: In electron shells of molecules, the external B0 field induces currents, which in turn add to the total field experienced by nuclear spins. This effect is generally anisotropic

U

because the topology of a molecule favors certain electron paths. The strength and direction of the

JO

created field thus depend on the orientation of the molecule relative to B0. With respect to relaxation, this effect is most relevant in solids. In liquids and gases, it is typically canceled by motional averaging. 3) Electric quadrupole coupling: In nuclei with spin > 1/2, the distribution of charge is non-isotropic,

giving rise to an electric quadrupole moment that interacts with electric field gradients generated by surrounding electrons. This interaction depends on the relative orientation of the molecule and the nuclear moments and thus fluctuates upon tumbling in the B0 field. The effects of orientation-dependent interactions are schematically illustrated in Fig. 1 for the case of dipolar coupling between two equivalent nuclear spins in a molecule. Upon rotation of the molecule,

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

the spins remain aligned with B0 such that their dipolar interaction changes (Fig. 1a). The associated effect on NMR spectra depends on the configuration and mobility of molecules in a sample (Fig. 1b). For a rigid single crystal, all local lattices form the same angle with B0, leading to a single sharp peak. A second symmetric peak is created by opposite spin alignment. In a powder, all lattice orientations co-exist and give rise to a continuum of resonance frequencies forming the characteristic Pake pattern

F

of a powder spectrum [17, 18]. In liquids, tumbling of molecules partly cancels dipolar interaction by

O

temporal averaging. In the limit of rapid tumbling, the spectrum converges to a single sharp peak as would be obtained without any coupling. However, in real liquids, tumbling occurs on a finite time

O

scale, still averaging dipolar coupling and leading to a single spectral peak, yet with a finite linewidth (Fig. 1b), which is associated with relaxation as discussed in the following section.

PR

Free induction decay (FID) signals corresponding to these spectra are shown in Fig. 1c. The hypothetical sharp peak corresponds to a signal of constant amplitude. With relaxation, the signal decays exponentially. As this decay is driven by stochastic fields, it is irreversible and cannot be

E-

countered by spin-echo refocusing. The two sharp peaks of the single crystal correspond to oscillation with a constant envelope in the time domain while the broad, sharply delineated spectrum of the

PR

powder entails a rapid signal drop with subsequent ringing. In the ideal powder, the loss of coherence is not of stochastic nature. In principle, it could be fully recovered by refocusing and is thus not considered actual relaxation. Mechanism of relaxation

AL

2.4

The dependence of relaxation on random molecular motion is shown in the final plot of Fig. 2

R N

(“Relaxation”, bottom right): In liquids, high molecular mobility leads to long T2 and T1 whereas in solids restricted mobility causes short T2 associated with long T1. The whole Fig. 2 summarizes the mechanism of relaxation driven by orientation- or position-dependent interaction in combination with

U

tumbling motion. Quasi-random rotational or translational motion and the anisotropic or non-uniform field associated with interaction lead to a random field ΔB(t) at the location of a considered spin. The

JO

amplitude of the field fluctuations is given by the interaction whereas the temporal behavior is linked with the speed of the tumbling motion and characterized by the correlation time τc. The latter can be

determined by calculating the autocorrelation function of the random field which is usually well described by an exponential function. The Fourier transform (FT) of the autocorrelation function is the spectral density J(ω) of the tumbling motion, i.e., the energy available as a function of frequency. It has the shape of a Lorentz curve and depends on τc as illustrated by two examples. The spectral density is the key link between molecular motion and relaxation. Relaxation changes the state of the spin system. Longitudinal relaxation affecting Mz is associated with change of the populations of the different energy levels, as indicated for a two-spin system in

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

Fig. 2. The transition probabilities depend on the order of the state change and are proportional to the spectral density of the tumbling motion at multiples of the Larmor frequency ω0. 2.5

T1 and T2

The final plot in Fig. 2 shows the qualitative behavior of relaxation times as functions of the

F

correlation time of tumbling motion. At very short τc as found in small-molecule liquids, interactions are averaged effectively, resulting in very slow relaxation, i.e., long T1 and T2. As τc increases, both

O

types of relaxation initially speed up. T1 reaches a minimum at τc ≈ 1/ω0, reflecting dependence on

O

J(ω0) and J(2ω0). Hence, longitudinal relaxation is most effective when 1/τc approximately matches the Larmor frequency, similar to the resonance condition for RF excitation. It is instructive to view

PR

longitudinal relaxation as caused by the transverse components of the random field whose action on Mz starts to cancel itself once spins perform more than half a precession cycle during the correlation time [14]. In contrast, T2 has no minimum, reflecting the fact that for large τc it is dominated by J(0).

E-

In this regime, transverse relaxation is driven mostly by the longitudinal component of the random field. Longitudinal relaxation inherently entails a certain loss also of transverse coherence, limiting

2.6

PR

the transverse relaxation time to T2 ≤ 2T1. In most practical cases T2 ≤ T1 [19]. T2*

In real experiments, the decay of an FID signal obtained from a given volume is often faster than

AL

according to T2. The chief cause is spatial variation of the static magnetic within the volume. The associated frequency spread causes dephasing of the spins, i.e., loss of macroscopic transverse coherence and thus signal decay. As long as the field variations are static, the dephasing can, in

R N

principle, be undone by spin-echo refocusing. In contrast to pure random-field-driven relaxation, signal decay including macroscopic dephasing has been termed apparent transverse relaxation. The shape of such decay varies widely depending on the underlying field distribution. Nevertheless, it is

U

frequently approximated by an exponential with the reduced time constant T2*. In the following, three

JO

cases of apparent transverse relaxation are discussed more closely: 1) Dipolar coupling in a powder as described in Fig. 1. In this case, the shape of the FID deviates strongly from exponential behavior. Large frequency spreads on the order of typically tens to hundreds of kHz require large-bandwidth, high-power RF pulses for refocusing. Alternatively, the dipolar interaction may be averaged out by spinning the sample at the magic angle or eliminated by RF decoupling pulses [18]. Note that such approaches are usually impractical in vivo. 2) Microscopic structure: Biological tissues (e.g. trabecular bone) and other samples of interest (e.g. porous rock) often exhibit microscopic structure on a length scale below feasible MRI

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

resolution. Related microscopic variation of magnetic susceptibility causes corresponding field inhomogeneity and signal dephasing. In T2*-weighted MRI, the associated signal decay alters image intensity on a per-pixel basis and can be used to characterize the underlying structure and materials involved. 3) Macroscopic B0 inhomogeneity: Macroscopic compartments filled with materials of differing

F

bulk susceptibility such as soft tissue, bone, and air, give rise to long-range field variation

O

that also causes dephasing and signal decay. However, such dephasing is not characteristic of the tissue in which it occurs. For compartment sizes in the order of or larger than the spatial

O

resolution, macroscopic susceptibility effects rather hamper MRI by depiction errors, e.g., signal shifts or blurring, and unspecific signal dropout. These problems can be mitigated by

2.7

PR

increasing spatial resolution and the strength of gradient fields used for spatial encoding. Diffusion

E-

In the presence of static microscopic field inhomogeneity, signal decay is further influenced by diffusion of the molecules that bear the resonant spins. Random motion translates field inhomogeneity

PR

into random field fluctuation that is individual for each spin and thus affects the nature of dephasing among the spin ensemble. This mechanism differs from relaxation as discussed above in that the effective field fluctuations are smaller by many orders of magnitude and arise from fields on a greater length scale such that their correlation times are much longer. As a result, dephasing under diffusion

AL

is amenable, at least in part, to refocusing by radiofrequency pulses. In a linear field gradient with amplitude G and unbounded Gaussian diffusion, the signal amplitude of an FID decays according to



(3)

R N

 



   

with self-diffusion constant D and gyromagnetic ratio γ [20, 21]. In spin-echo experiments with echo

U

time TE, refocusing partly recovers the decay, resulting in an apparent relaxation rate of 







    



(4)

JO

 

Hence, MR experiments may be hampered by increased apparent relaxation due to diffusion but can also give access to the physical property D. 2.8

Examples

As discussed above, various mechanisms can cause loss of coherence and thus rapid signal decay in tissues and materials. Accordingly, short-T2 MRI has a large variety of potential imaging targets, including the following prominent examples that will be further discussed in Section 10:

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

1) Solids: In solids, molecular motion is strongly restricted, leading to long correlations times. As seen in the relaxation plot of Fig. 2, in such substances short T2 is associated with long T1. The latter has two key implications for sequences with rapidly repeated RF excitations as introduced later: 1) steady-state magnetization is generally low and 2) maximum steady-state signal is reached with small excitation flip angles. In the biomedical context, a solid of

F

particular interest is hydroxyapatite in the bone matrix [22].

O

2) Bound water: Water molecules are small, leading to short correlation times and slow relaxation in free water. However, when bound to macromolecules or solids, water molecules

O

exhibit reduced mobility and much shorter T2. An example is collagen-bound water in bone [23].

PR

3) Water in collagen fibers: Some tissues (e.g. tendons, cartilage) are composed of collagen fibers which are organized in directed, parallel bundles. Motion of water embedded in such highly ordered structures is strongly restricted, leading to short T2. In addition, there are

E-

preferred orientations along which bound water molecules line up. Thus the effect of dipolar coupling and hence T2 become a function of the angle between fibers and B0 [24-33]. In T2-

PR

weighted acquisitions, signal varies and peaks at the so-called magic angle of 55°. This needs to be considered for image interpretation but also holds potential for contrast manipulation. 4) Myelin: This membrane tissue, which plays a key role in the nervous system (see section 10), is largely composed of phospholipids arranged in liquid crystalline bilayers. Intra-molecular

AL

dipolar interactions cause rapid T2 decay of 1H signal, which varies with the orientation of the bilayer with respect to B0. As in typically probed volumes all orientations exist, the resulting

R N

signal is well described by a superposition of exponential decays with a wide range of T2 values, leading to so-called “Super-Lorentzian” lineshapes [34, 35].

5) Lung parenchyma: Signal decay of water in pulmonary tissue is an important example for the diffusion-driven mechanism described in section 2.7. Microscopic tissue-air interfaces at the

U

level of the alveoli create local field gradients that rapidly dephase spins in water diffusing in

JO

the tissue [36]. This results in fast signal decay in gradient-echo techniques, which can be partly reversed with spin-echo sequences, according to Eqs. 3 and 4.

3 Fundamentals of short-T2 MRI

In this Section, the specific requirements for depicting tissues with short T2 are expressed in terms of two guiding principles. Starting from these principles, concepts of sequence design are derived, leading to four basic templates of short-T2 techniques.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

3.1

Principles

Figure 3 illustrates the challenges of short-T2 imaging. Common echo sequences hardly capture any signal from short-T2 spins due to signal loss by a factor







until the echo is formed and acquired

around the echo time TE (Fig. 3a). For gradient echoes this holds equally for pure T2 decay and T2*

F

decay. In spin echoes, the latter is partly refocused, which can be beneficial when T2* is much shorter than T2. However, at very short T2*, RF refocusing requires commensurately large bandwidth and thus

O

power deposition, which limits the applicability of spin echoes for this purpose. Therefore this option

O

is not further considered here in favor of techniques that can equally handle short T2* and T2 . To contain the signal loss observed with echo techniques, dedicated short-T2 sequences must start

PR

gradient encoding and data acquisition quickly after signal excitation. Following this first principle renders short-T2 materials visible to begin with (Fig. 3b). However, signal decay during data acquisition still impairs spatial resolution. In mathematical terms, this is due to apodization in k-space,

E-

which is equivalent to convolution with a kernel of reciprocal width in the image domain [37]. To avoid such resolution loss, all image data must be acquired within a small time range after excitation.

PR

In Fig. 3c, this second principle is implemented by use of a stronger gradient, which shortens the readout and thus improves the sharpness of short-T2 features. More formally, the relationship between time and k-space in an MRI sequence is described by

AL

defining the time of acquisition  as the time after signal excitation at which a k-space point  is acquired. In echo-based MRI techniques, the value of this function at the k-space center,  , is the echo time TE. The same definition of TE has widely been adopted for short-T2 imaging even though

R N

most of the techniques work without actual echo formation (see Section Error! Reference source not found. for details). To keep with this convention, the same notion is used in the present work. The range of times across which data are collected shall be referred to as the acquisition range   

U

  !"  .

JO

Based on these definitions, the principles of short-T2 sequences can be narrowed down to 1) data acquisition must start early, i.e., !"  must be small, and 2) the acquisition range  must be

small. Both of these conditions must hold on the time scale of relevant T2. Table 1 summarizes these considerations along with characteristic sequence parameters and their impact on image quality. Note that !"  in short-T2 imaging is often, but not necessarily, equal to TE. To approach the encoded, nominal resolution, the acquisition range should approximately match the decay time constant, i.e., #  [38]. Shortening the acquisition further will yield diminishing returns in terms of resolution but comes at the expense of signal-to-noise ratio (SNR).

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

The principles have an immediate consequence on any short-T2 sequence, as illustrated in Fig. 4. Competing with fast signal decay by reducing TE and/or  requires rapid signal encoding. On the other hand, achieving high spatial resolution requires large time integrals of gradient amplitude. These demands can be fulfilled jointly only by using strong gradients and corresponding large spread of resonance frequency across the imaged object. As discussed in the following, such strong gradients

F

may be active during RF excitation and/or data acquisition, making large bandwidth a characteristic of

O

short-T2 MRI. This consideration also reveals the fundamental challenge that necessarily short observation times of rapidly decaying transverse magnetization are equivalent to high-bandwidth

Concepts

PR

3.2

O

acquisition. Hence, short-T2 imaging is generally more SNR-limited than conventional MRI.

At the level of sequence design, the principles discussed above lead to three chief concepts that recur across the variety of short-T2 techniques:

E-

A) No slice selection: Slice selection under a gradient of finite strength requires RF pulses of suitably small bandwidth and thus long duration, followed by slice rephasing. To avoid the associated delays,

PR

short-T2 imaging frequently works with non-selective volume excitation, omitting the slice selection gradient. All spatial encoding is then achieved by three-dimensional (3D) Fourier encoding. (An exception to this approach will be described later in Section 4.3.)

AL

B) Radial center-out encoding: To reach every targeted k-space point in the shortest possible time, kspace is traversed by straight radial trajectories that point outward from the center. This approach is



R N

used in two principally different ways:

Pure phase encoding: after each excitation, a single point in k-space is targeted and sampled once reached.

Pure frequency encoding: data acquisition is performed continuously along radial trajectories.

U



C) Excitation under encoding gradient: To avoid losing time ramping up the readout gradient after

JO

excitation, the gradient can be switched on first such that it is already present at when excitation is performed. Figure 5 visualizes the combinations of these design elements and indicates published implementations using the respective approaches. While all techniques share concepts A and B, they differ with respect to the options for concepts B and C. All current short-T2 techniques can be represented in this scheme and thus viewed in a shared conceptual framework. As the figure indicates, two options each with respect to concepts B and C give rise to four basic sequence templates, which will now be considered in more detail.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

3.3

Basic sequences

The four basic schemes of short-T2 sequences and their associated k-space sampling patterns are shown in Figs. 6 and 7, respectively. All techniques employ non-selective hard-pulse excitation followed by 3D radial center-out encoding in k-space. Pure phase encoding is used in constant time imaging (CTI) [39] (Fig. 6a) and SPI [5, 40, 41] (Fig.

F

6b). In both, a gradient of varying amplitude and direction is used to travel to one point in k-space per

O

excitation (Fig. 7a). Each k-space point is sampled at the same fixed  = TE, which must be chosen large enough to permit reaching the maximum k-space excursion that corresponds to the targeted

O

resolution. While k-space is traversed radially, the points sampled are typically chosen to lie on a

PR

Cartesian grid. By virtue of fixed  , all data points have identical T2 weighting. Therefore, signal decay prior to TE does not impair spatial resolution. However, it does affect the SNR in resulting images, most so for short-T2 material. In CTI, TE includes time for ramping the encoding gradient up

E-

and down such that it is off during excitation and acquisition. This has the advantage that the bandwidth of both can then be limited to the intrinsic frequency spread in the object. However, it

PR

exacerbates the SNR penalty of signal decay.

As a remedy, in SPI TE is reduced by performing RF excitation and data acquisition in the presence of the gradient (Fig. 6b). The bandwidth of both must then be increased by making them suitably short. As a result, the SNR yield of SPI is subject to competing influences of signal decay and acquisition

AL

bandwidth. It is optimal at    , requiring a maximum gradient amplitude of   $ %&

for given T2 and spatial resolution %& [40]. In addition to its short-T2 capabilities, the SPI approach

R N

offers relatively silent operation if gradients are not switched off between excitations but merely adjusted to a nearby value [42]. The key advantage of purely phase-encoded techniques is full nominal spatial resolution at short T2.

U

Their chief drawbacks are long scan times and low SNR efficiency associated with point-wise sampling, low acquisition duty cycle, and comparatively long TE. For these reasons, pure phase

JO

encoding is rarely used in vivo [43]. It is mainly applied for imaging materials [44, 45] and, in MR microscopy, to mitigate diffusion effects [46]. Recently, as discussed in Section 6, partial SPI acquisition has found additional use as a means of complementing other, more time-efficient techniques. The terminology of purely phase-encoded imaging is not entirely consistent, with somewhat variable use of the terms CTI and SPI in the literature. With pure frequency encoding, lines in k-space are sampled by continuous signal acquisition during radial encoding. In this fashion, shorter TE, faster overall scanning, and higher SNR efficiency are accomplished. For an image matrix size M, the gain in acquisition duty cycle relative to SPI leads to

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

an increase in SNR efficiency of  at equal maximum G and of  at optimized G (see 9.2). In both UTE [6, 47-54] (Fig. 6c) and ZTE [7, 55-58] (Fig. 6d) imaging, radial encoding is performed by playing out gradients of fixed amplitude but direction varying from repetition to repetition. Thus one spoke is obtained per excitation (Figs. 7b and 7c). In the UTE sequence, data acquisition is started concurrently with ramping up the gradient shortly after excitation. The minimum TE then amounts to

F

half the RF pulse duration plus the dead time required to change from transmit to receive operation.

O

Due to gradient ramping, the central k-space is traversed at lower, linearly increasing speed (Fig. 7b). This is avoided in the ZTE sequence where RF excitation is performed only after switching on the

O

gradient, thus achieving instant, full encoding speed. As a consequence, similar to SPI, RF excitation must be broadband and near-silent operation is possible by only adjusting and not switching off the

PR

gradient between excitations. In the ZTE case, spatial encoding starts immediately at the time of creating transverse magnetization; hence the time TE at which it leaves the k-space center is effectively zero. However, during the finite RF pulse duration and the RF dead time initial signal is

E-

lost, leading to a gap in central k-space (Fig. 7c), which needs to be addressed by additional measures (see 5.3 and 6.1).

PR

For in-vivo MRI, UTE and ZTE imaging and their relatives are currently the most suitable and most widely used techniques and will therefore be discussed in more detail in Sections 4 and 5. For all implementations discussed above, the concept of radial encoding is central to short-T2

AL

depiction performance. With common Cartesian encoding, the principles stated initially cannot be implemented as effectively. Notwithstanding this, Cartesian gradient echo imaging methods can also

R N

be optimized for short-T2 imaging [59, 60]. 4 Ultra-short echo time imaging The UTE strategy is the most established and most widely deployed short-T2 approach in use today.

U

While arguably simple in principle, UTE sequence design and image reconstruction deviate significantly from standard MRI. To obtain robust results, implementations must consider a number of

JO

subtleties and pitfalls as detailed below. 4.1

Gradient ramping

The basic idea of the UTE sequence as shown in Fig. 6c is to start data acquisition very soon after RF excitation while concurrently ramping up the radial readout gradient to its target amplitude. The gradient ramp plays a key role in UTE MRI because it affects image quality in several ways. Depending on the desired gradient strength and the available slew rate (see Section 9.2), the ramping part can engage a considerable portion of the total encoding time. For a given nominal spatial resolution, the encoding time must be extended by half the ramp duration relative to the ZTE case

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

with full gradient amplitude throughout. The effect of the gradient shape on depiction behavior in UTE imaging is illustrated in Fig. 8, using the point spread function (PSF) as a criterion. The time of acquisition of each k-space value (Fig. 8a) shows sections of increasing and constant k-space speed during gradient ramp and plateau, respectively. Longer ramps lead to larger acquisition range . In conjunction with T2 decay, the acquisition timing translates into the modulation transfer function

F

(MTF) (Fig. 8b), which is then Fourier-transformed to obtain the PSF (Fig. 8c). As ramp times

O

increases, apodization of the MTF causes progressive broadening of the main lobe of the PSF. The

4.2

O

normalized display (Fig. 8d) illustrates the associated loss of spatial resolution. Gradient fidelity

PR

Gradient ramping is also critical to UTE image quality in that the demanded gradient waveforms may not be produced with the necessary fidelity. Delays, finite gradient amplifier bandwidth, and eddy currents cause deviations from nominal gradient dynamics. These may include dynamic fields

E-

different from first spatial order, usually with a dominant zeroth-order component. The simulations in Fig. 9 illustrate effects of eddy currents on UTE image quality. The distorted gradient shape in Fig. 9c

PR

leads to severe alterations in signal intensity when image reconstruction is based on the nominal trajectory (Fig. 9d). So-called B0 eddy currents that generate zeroth-order field and thus spatially constant phase errors (Fig. 9e) result in edge artifacts (Fig. 9f). Both types of error are amenable to retrospective correction. Phase perturbation by B0 eddy currents can be removed from the acquired

AL

raw data while first-order field effects can be addressed by image reconstruction based on the effective, distorted k-space trajectory. In this way, a virtually artifact-free image is obtained,

R N

comparable to imaging with ideal gradient behavior (Fig. 9b). Models of system imperfection for retrospective correction may be set up at different levels of complexity ranging from mere gradient chain delays to general first-order field perturbation,

U

additional zeroth-order [61], and second- and higher-order field terms [62]. Various strategies have been described for determining the parameters of such models, including approaches based on raw

JO

UTE data itself [63, 64], calibration scans [65-67], impulse-response characterization of the gradient system [68], and measurement of field evolution using modified MRI sequences [69, 70] or dedicated field probes [71] (see Section 9.5). 4.3

RF excitation

In the 3D UTE sequence shown in Fig. 6c, spatially non-selective RF excitation is followed by 3D radial encoding of the whole object. Alternatively, unlike SPI and ZTE techniques, UTE permits slice excitation under a separate selection gradient, forming a two-dimensional variant. To avoid increasing TE by longer RF pulses and slice refocusing, it has been proposed to employ half-excitation pulses

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

with opposite slice gradients in two successive acquisitions [47, 72, 73]. Full slice selectivity is then obtained by summation of acquired data. On the downside, besides doubling the scan time this approach has been reported to increase sensitivity to eddy currents [74, 75] and motion. Relaxation during excitation must be considered when pulse durations are significant relative to T2 of interest. In-pulse relaxation reduces the actual flip angle and renders it T2-dependent [52]. The

F

presence of an RF field changes the physics of relaxation, which must be taken into account in

Echo time

O

4.4

O

sequence optimization and the interpretation of data [76, 77].

Two different definitions of TE are common in the UTE literature. The first, as shown in Fig. 6c and

PR

used throughout this article, assumes that relaxation starts at the magnetic center or focus point of the RF pulse [77, 78]. For symmetric waveforms, including the commonly employed hard pulses, this is the center of the pulse. The second definition assumes that relaxation starts only after the RF pulse

E-

[38]. Both pictures are incomplete in that they do not account for relaxation during excitation as discussed in the preceding section. Hence, the full characterization of a short-T2 sequence also

PR

requires RF shape and duration in addition to TE of either definition.

In short-T2 samples, the SNR that can be achieved with UTE imaging depends on the shortest possible TE, which is governed by sequence implementation and hardware constraints (Tab. 1). A lower bound

AL

for TE is given by the RF dead time, which is composed of the fraction of the RF pulse duration after its magnetic center, the time required to switch between RF transmission and reception, and the leadin time of bandpass filtering [79] (see 9.3). On current MR systems, typical TE values for 3D UTE

R N

MRI are in the order of several tens of microseconds. In the UTE scheme, TE can also be varied, which is a useful feature. It permits creating T2* contrast, emphasizing short-T2 tissues by image subtraction [80], and mapping T2* quantitatively [81].

U

Deliberate choice of somewhat longer than minimal TE can serve to avoid detection of signals that

JO

decay too fast for meaningful spatial encoding even with the UTE approach. Disturbing signal of such exceedingly short T2 most commonly originates from hardware parts and causes cloudy image background when unaddressed [82] (see 9.3). 4.5

Image reconstruction

Due to non-Cartesian sampling, image reconstruction from UTE data involves interpolation in kspace, which is usually performed with gridding procedures [83, 84]. In particular, this includes correction for non-uniform k-space density which can be accomplished by using a pre-calculated density correction function (DCF) [85-89]. One approach to obtain a DCF defines local k-space

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

density as the area surrounding a sampled point as calculated by Voronoi diagrams [86]. For radial sampling, this approach leads to the so-called rho filter [90] for which an implementation is provided in Section 12.1. Density correction is an approximation based on splitting up reconstruction into FT and a preceding filter, the DCF. More rigorously, reconstruction is viewed as a general inverse problem. It can be

F

solved efficiently using iterative algorithms [91] which inherently perform correction for non-uniform

O

k-space sampling. In this case, application of a DCF can serve for pre-conditioning, thus accelerating convergence. The choice of the DCF then affects only convergence but not the accuracy of the

O

reconstruction result.

PR

5 Zero echo time imaging

Image acquisition according to the ZTE scheme shown in Fig. 6d dates back several decades. Renewed recent attention has been driven by technical and methodological advances as well as

E-

increasing interest in short-T2 capability and silent MRI. To express its relation to the more established

5.1

Sequence properties

PR

UTE method, the approach has been dubbed ZTE imaging [7].

The basic ZTE sequence is an imaging analogue of the pulse-acquire NMR experiment. The only modification consists in spatial spreading of resonance frequencies by addition of a projection

AL

gradient that is present during RF excitation and FID acquisition. This approach was in fact used in the very first MRI experiment ever reported [92]. Only later it was identified and exploited as an

R N

approach particularly suitable for short-T2 imaging [55-58]. Relative to UTE and standard present-day MRI the distinguishing feature of ZTE imaging is RF excitation in the presence of the encoding gradient. While arguably a small difference in a sequence

U

diagram, it is of great consequence with regard to sequence considerations, hardware requirements, data acquisition, and image reconstruction. The chief objective of excitation under the encoding

JO

gradient is to remove the need for gradient ramping while short-T2 magnetization evolves. It saves

time, offers immediate full k-space speed, and minimizes gradient switching even after data acquisition. In this way, enhanced short-T2 sensitivity is obtained in combination with reduced

acoustic noise or even virtually silent operation and robustness against eddy currents. Removal of gradient ramps also reduces the time required per repetition and increases the acquisition duty cycle. In combination with large gradient strengths as required for short T2, very short repetition time (TR) can be achieved and scan times for ZTE imaging are comparatively short [93].

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

These capabilities come at the expense of shifted and partly increased hardware requirements. Continuous gradient operation for short-T2 imaging requires high gradient strength at full duty cycle. This demand is not commonly met by present-day gradient systems and has prompted developments in dedicated gradient hardware (see Section 9.2). On the other hand, the ZTE mode comes with reduced specifications in terms of gradient switching speed. In turn, it partly shifts hardware

Dead time

O

5.2

O

transmit-receive (T/R) switching are essential as will be detailed in the following.

F

challenges from the gradients to the RF subsystem. In particular, high-bandwidth excitation and rapid

In pulsed MR experiments, signal excitation and acquisition are separated in time to accommodate the

PR

large difference in the RF amplitudes involved. Changing between the respective modes of operation takes time, resulting in a dead time after RF excitation during which MR signal is not yet detected. In many MR experiments, especially those involving echo formation, acquiring data immediately after

E-

RF excitation is not of interest. However, in short-T2 imaging, encoding and acquiring fresh MR signal as soon as possible is one of the key principles (see 3.1). With immediate encoding at full k-

PR

space speed, the ZTE sequence is arguably the most rigorous implementation of this rule and thus particularly vulnerable to data loss during the dead time. In 3D k-space it causes a spherical gap at the center (Fig. 7c). To contain the gap size one must minimize the dead time, which has contributions from the RF pulse, T/R switching, and digital filtering (Fig. 10) [7]. The chief means of limiting the

AL

gap are thus the use of short RF pulses [93] (see Section 5.5), speeding up T/R switching (see Section 9.3), and the design and use of short digital filters (see Section 9.4). The residual k-space gap must be

R N

addressed either by means of image reconstruction (see Section 5.3), by acquisition of additional data (see Section 6.1), or a combination of the two. Ideally, RF dead time would be avoided altogether by performing signal acquisition also during

U

excitation. However, due to the large dynamic range between transmit and receive signal this

JO

approach is of limited practicality (see Section 9.3). 5.3

Algebraic Reconstruction

Straightforward image reconstruction from data with a gap in central k-space, using standard FT, leads to prohibitive artifacts of low spatial frequency. However, for small gap sizes the missing data can be recovered by taking two measures: 1) oversampling in the direct, radial dimension (Fig. 7c, Fig. 10) and 2) treating reconstruction as a general inverse problem set in the framework of linear algebra [94]. Standard sampling according to the Nyquist-Shannon theorem is associated with harmonic encoding functions that are mutually orthogonal within the field-of-view (FOV). In contrast, oversampling is associated with encoding functions that exhibit non-vanishing scalar products with

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

the Nyquist set. In the algebraic framework, these dependencies permit the implicit formation of missing low-frequency harmonics by linear combination of available encodings [95]. This approach is known as finite-support or band-limited extrapolation and was previously used for low-frequency restoration in MRI in a different context [96]. It is limited to a relatively small range of extrapolation, beyond which it becomes impractical due to ill-conditioning (see Section 5.4). Greater oversampling

F

improves conditioning somewhat but also increases the amount of resulting data. Therefore,

O

oversampling factors range between 2 and 4 in practice [95].

Formally, the algebraic approach models data acquisition by means of an encoding matrix [97]. Its

O

rows reflect the encoding functions actually implemented. Hence, the gap is represented by the absence of corresponding rows while oversampling elsewhere leads to an increased number of rows

PR

relative to standard Fourier encoding [7]. Finite support is implied by limiting the length of the rows to the FOV. Algebraic reconstruction is achieved by solving the encoding equation for the unknown image, eliminating the dead time gap and returning to regular sampling density in k-space. Figure 11

E-

illustrates this approach for the case of one-dimensional (1D) ZTE data. In principle, it can be straightforwardly expanded to 3D, working with one single encoding matrix. However, the latter will

PR

usually be too large for direct inversion and require iterative treatment. Alternatively, and more efficiently, algebraic processing can be restricted to 1D sub-problems, reconstructing full k-space lines through the origin from pairs of acquisitions with opposite gradient direction. The resulting full

3D image [7].

AL

set of radial k-space data is then merged by standard 3D gridding (see Section 4.5) to form the final

Data extrapolated to the gap in the described manner reflect the actual T2 decay during the dead time.

R N

In particular, the value at k = 0 is not affected by transverse relaxation. Hence, although these data points are not actually acquired, they correspond to signal obtained as if sampling was started at TE =

U

0, which has led to the “ZTE” notation. 5.4

Bandwidth limitations

JO

There are limitations to the gap size and hence the signal bandwidth that can be handled with the algebraic approach. These arise from deteriorating conditioning of the encoding matrix as gap size increases [94]. Ill-conditioning boosts error in the signal and the signal model, causing three kinds of image distortion: 1) enhancement and spatial correlation of noise, 2) aliased and amplified out-ofband signal, and 3) intra-voxel oscillation of the spatial response [95]. Usually, the first two are dominant in practice and will therefore be described in more detail below. As the reason for illconditioning is data missing in central k-space, all image distortions described are of low-frequency nature and may not be distinguished easily.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

1) Noise in common Cartesian k-space raw data exhibits constant variance which leads to the usual uncorrelated appearance of noise in MR images. In contrast, in reconstructed ZTE data, noise variance increases strongly towards the k-space center, thus causing noise correlation in image space. Therefore, at large gaps, noise in ZTE images can have an unusual artifact-like appearance (Fig. 12). Its random nature is only observed by varying shapes in the image resulting from different noise input

F

data. Importantly, ZTE image reconstruction benefits from intrinsic averaging in central k-space

O

associated with the sampling density of 3D radial data increasing strongly at small radii. Therefore, distortions appearing in 1D projections do not necessarily translate into deterioration of 3D image

O

quality.

2) Out-of-band signal violates the fundamental assumption of ZTE reconstruction in that it stems from

PR

locations outside the FOV. As a consequence, such signal will be assigned into the FOV, corresponding to aliasing. Furthermore, related to the bad conditioning it will appear strongly amplified in the image. This is a serious practical issue as ZTE acquisition not only captures the

E-

desired short-T2 signal of the object of interest but also unwanted contributions from nearby hardware parts, in particular RF coils (Fig. 13). The straightforward solution of targeting a larger FOV

PR

including all detected signal is of limited scope as it results in larger bandwidth, i.e., larger gap size, as well as longer scan times and larger data sets. Better ways of reducing such background signal are suppression techniques [7, 98], image subtraction [99], and designing RF coils with negligible content

AL

of the observed nucleus (see Section 9.3).

Within these boundaries and depending on setup and sequence parameters, practical gap sizes in ZTE experiments are usually limited to approximately 2 - 3.5 DW, where DW = 1/BW is the Nyquist dwell

R N

time for the bandwidth (BW) in the FOV. Using ZTE-optimized RF hardware, this can nevertheless translate into BW as high as 500 kHz [93], which is the key to high-resolution short-T2 MRI (see

U

Section 3.1). 5.5

RF excitation

JO

RF excitation in ZTE imaging plays a particular role as it occurs while the radial encoding gradient is already on (Fig. 6d). To ensure consistency between data of different radial directions, the spectrum of the excitation pulse must be largely uniform within the bandwidth spanned by the gradient across the object. In the basic ZTE sequence, a short block-shaped, so-called hard pulse is used. For greatest uniformity the pulse duration should not exceed approximately DW/3 [100]. However, in practice durations up to about DW can also provide useful results. For bandwidths of hundreds of kilohertz this limits hard pulse durations to a few microseconds, which is much shorter than RF pulses commonly applied on

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

clinical MR scanners. With the maximum B1 amplitudes usually available on such systems, flip angles of only a few degrees can be obtained (see Section 9.3). On the other hand, as little time is spent on RF pulses and gradient switching, the ZTE sequence can be operated with very short TR (see Section 7.1), even below 1 ms [93], thus enabling maximum steady-state magnetization even at very small flip angles.

F

If flip angles should nevertheless be increased due to low B1 or very high bandwidth, or to create T1

O

contrast, three options have been proposed, namely 1) longer hard pulses with intensity correction, 2) short modulated pulses, or 3) long modulated pulses with simultaneous excitation and acquisition.

O

1) By increasing the hard pulse duration, its bandwidth will be reduced. For durations above

PR

DW, the main lobe of the Sinc-shaped pulse spectrum exhibits relevant variation in the FOV. Hence, larger flip angles are primarily obtained in the central FOV, which may still be useful if this coincides with the location of interest. Inconsistencies in signal intensity between

E-

projections can be reduced by correcting for the pulse profile [101]. This approach is limited by singularities above durations of 2DW and by variations of local steady-state magnetization

PR

over time.

2) The hard pulse is replaced by an RF excitation with both modulated amplitude and frequency [102, 103]. If designed appropriately, such pulses enable high spectral uniformity at high flip angle efficiency and thus outperform hard pulses of equal duration (Fig. 14). Unlike hard

AL

pulses, their length is not limited by uniformity considerations but by the maximum acceptable dead time (see Section 5.2). 3) Further extending modulated pulses may lead to intolerable dead times. To reconcile long

R N

pulses and short dead time it has been suggested in the SWIFT (sweep imaging with Fourier transformation) technique to overlap RF excitation and acquisition [104, 105]. This is indeed possible by alternating the two RF operation modes rapidly (Fig. 15), thus offering enhanced

U

degrees of freedom in pulse definition. However, frequent switching comes at the expense of

JO

reduced acquisition duty cycle and thus SNR efficiency. Furthermore, feasible bandwidth is ultimately limited by the inverse of the dead time [79]. To image large bandwidths at reduced switching rate, the SWIFT approach has been expanded to a multi-band mode [106]. Alternatively, technically challenging approaches have been proposed where excitation and acquisition are performed truly simultaneously without gapping (see Section 9.3).

Apart from short hard pulses, all excitation options listed above generate some kind of spectral nonuniformity which needs to be considered during image reconstruction [101, 103, 104, 107]. A potential implementation obstacle may be the high temporal resolution required for creating highbandwidth pulses which is not always available on clinical scanners (see Section 9.4). A more

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

fundamental limitation is posed by the large power deposited by high-bandwidth excitations, leading to potentially excessive specific absorption rates (SAR). 5.6

Comparison with UTE

In summary, the ZTE option exhibits the following advantages (+) and disadvantages (-) with respect

F

to UTE imaging: Higher short-T2 sensitivity and resolution

+

Greater resilience against off-resonance effects

+

Greater robustness against eddy currents, no trajectory calibration required

+

Silent operation



High-bandwidth RF excitation: higher SAR and tighter limits on flip angles



Bandwidth limited by dead time



Only 3D



Limited intrinsic contrast



Higher sensitivity to background signal

O

PR

E-

PR

6 Intermediate techniques

O

+

Various modifications of the basic sequences discussed so far in Sections 3.3, 4, and 5 have been proposed to mitigate some of their limitations. This results in the intermediate techniques indicated in

AL

the classification scheme of Fig. 5. They implement transitions on the two axes corresponding to concept B (vertical axis) concerning phase and frequency encoding and concept C (horizontal axis)

R N

concerning the succession of RF excitation and gradient switching. Intermediate techniques are obtained by either combining data obtained with different techniques or by merging opposed concepts within a sequence.

Filling the ZTE dead-time gap

U

6.1

JO

Gap size limitations in algebraic ZTE have mainly been tackled by acquisition of additional data in the k-space center: 1) Sampling at k = 0 by additional acquisition in the absence of a gradient improves the conditioning [95, 108, 109]. However, this data point will have a finite TE equaling the dead time, thus leading to considerable signal decay for tissues and materials very short T2 values. This introduces an inconsistency with the bulk data which leads to image artifacts for badly conditioned algebraic reconstruction. 2) The previous approach is expanded in the PETRA (pointwise encoding time reduction with radial acquisition) technique [110] by filling the entire spherical void in k-space with

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

Cartesian SPI data (Fig. 16a). In this way, arbitrary gaps can be handled, though at the expense of scan time and SNR efficiency. The effect on these two parameters increases steeply with gap size. 3) Alternatively, in the technique known as “water- and fat-suppressed proton projection imaging” (WASPI) [111], the dead-time gap is filled with radial data acquired at a gradient

F

strength reduced by a factor 2 × dead time / DW (Fig. 16b), thus achieving Nyquist sampling

O

in the very center also. By adapting the angular density to the radius given by the gap size, the additional scan time is usually less than for PETRA.

O

The latter two options change the weighting obtained in k-space associated with T2 decay. The related effect of on image quality can be assessed via the PSF (Fig. 17). The PSFs show that compared to

PR

ZTE, the main lobe width is reduced for both options but the PETRA technique offers relatively benign side-lobe behavior whereas WASPI is prone to increased side lobes. It has been found that this tendency can lead to artifacts in the presence of components with very short T2 on the order of the gap

E-

size, stemming, e.g., from hardware parts (see Section 9.3) or solid tissues [109].

PR

4) In the hybrid filling (HYFI) technique [112], radial data from a WASPI acquisition with multiple gradient strengths and a Cartesian SPI acquisition in the very center similar to PETRA are combined. This approach enables improved PSF quality as compared to WASPI with less scan time increase than required for PETRA.

AL

In the techniques described above, data missing in k-space is complemented by performing additional gradient encoding. In a different approach following the principles of parallel imaging [97],

R N

supplementary information is obtained by detecting the MR signal with a receiver coil array. During image reconstruction, independent data from multiple coils facilitates deriving missing from existing neighboring data, thus improving the conditioning associated with large gap size [113]. A challenge of this approach is the need to solve the full 3D inverse problem for image reconstruction, differently

U

from the simplified procedure outlined in Section 5.3. Finally, the encoding potential of RF pulses

JO

may be exploited to reduce the effective gap size [114]. 6.2

Reducing the excitation bandwidth

Another class of intermediate techniques departs from single gradient ramps either before or after the RF pulse by performing excitation under a lower gradient which is then ramped to full strength for rapid encoding. Such techniques have been proposed as gradient-modulated PETRA [115] or RHE (ramped hybrid encoding) [116] as well as gradient-modulated SWIFT [117]. As an advantage, RF excitation can be performed at lower bandwidth, thus relaxing flip-angle and SAR constraints. On the

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

other hand, gradient ramping per se and its fidelity affect image quality in analogy to the UTE case (see Sections 4.1 and 4.2). 7 Further sequence aspects 7.1

Repetition time

F

The sequences that have been described above are usually operated with short TR and a steady state of

O

the magnetization. In this way, the many repetitions required for 3D acquisition can be performed within useful scan times. Assuming complete spoiling of transverse magnetization after each TR,

O

maximal transverse steady-state magnetization ' and thus optimal SNR are obtained by excitation

PR

at the Ernst angle ( [118]. For ) * , which holds for most short-T2 imaging scenarios, both ( and ' are approximately proportional to ) (see Section 12.2). ' + ) implies that, when operating in the Ernst-angle regime, changing TR has no effect on SNR efficiency [119]. For example,

E-

reducing TR by a factor of two while doubling the number of averages will yield the same SNR at equal scan time. Furthermore, ( + ) implies that changing TR does not affect the SAR generated

PR

by the excitation.

To achieve short TR, gradient spoiling must be implemented in an efficient way (see Section 0). Particularly fast repetition is possible with minimal gradient switching as permitted by SPI (see Section 3.3) and ZTE (see Section 5.1). For these techniques, as an additional benefit, reducing TR

AL

helps approach the Ernst angle which may otherwise not be reached (see Section 5.5). Hence, in highbandwidth ZTE scanning, TR may range down to only hundreds of microseconds whereas UTE

7.2

R N

techniques are often operated with TR values of several or even tens of milliseconds. Gradient spoiling

U

To obtain a strict steady state, the same net gradient area must be applied in all dimensions for each TR [120]. Hence, in radial sequences each readout gradient must be rewound and gradient spoilers

JO

must be applied along dimensions fixed in physical space to destroy unwanted coherences [38]. Therefore, despite sub-millisecond encoding times, TR often amounts to several milliseconds. To reduce TR, pseudo-steady-state operation can be employed by leaving out spoilers and instead randomizing gradient areas [65] and/or projection orders [7, 108]. The latter is common particularly in ZTE and SPI sequences where switching gradients is usually avoided to minimize TR (see Section 7.1) and to keep their near-silent nature. Spoiling is then performed by extending the readout gradient along with suitable projection ordering.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

7.3

Reducing scan time

3D radial scanning is relatively time-consuming, particularly at large matrix sizes. To fulfil the Nyquist criterion for center-out encoding, a factor of π more encoding steps are required compared with a conventional Cartesian sequence. Ways of reducing the scan time include generalized parallel imaging schemes [97, 121], angular undersampling [122] (see Section 7.4), and compressed sensing

F

[123]. Furthermore, the k-space area covered per TR may be increased by using other center-out

O

encoding schemes such as, e.g., spiral- [124], or cone-shaped [125-128] trajectories, within the limits given by gradient specifications and targeted T2. Another modification proposed for the UTE gradient

O

shape aims at SNR optimization by adapting the k-space speed to compensate for the non-uniform

7.4

PR

density of radial encoding [129]. Geometry

The specific way of performing spatial encoding with the described techniques has implications for

E-

scan geometry.

PR

1) For both pure phase and radial frequency encoding there is no particular readout dimension in which bandwidth-limited filtering can be applied to restrict the FOV. Although possible for individual 1D projections in the radial techniques, such filtering would lead to inconsistent 3D data similar to selective RF excitation (see Section 5.5). Hence, in all (usually three)

AL

spatial dimensions a FOV given by object shape and the sensitivity range of the RF coils needs to be fully encoded to avoid aliasing. Thus potentially more encoding effort is required than with conventional sequences.

R N

2) As opposed to the discrete nature of the PSF in Cartesian imaging, the PSF for radial encoding exhibits a dense pattern of low-amplitude streaks. Therefore, undersampling causes streak artifacts of intensities much lower than with Cartesian sampling, and reducing the

U

angular distance of the radial spokes by factors 2-4 leads to relatively benign and often

JO

tolerable artifacts.

3) For isotropic encoding, the FOV in 3D radial techniques is spherical. A simple adaptation to the object shape can be achieved by individual scaling of the gradients per dimension, resulting in an ellipsoidal FOV and anisotropic resolution. Disentangling these two parameters is possible with a generalized approach to radial encoding [130]. 4) Creating FOV off-centers with radial encoding is also more involved than with Cartesian techniques. For off-center excitation under a gradient, as in the SPI and ZTE sequences, the RF frequency must be adjusted to the respective gradient direction [7]. Similarly the receiver frequency needs to be set individually. In UTE encoding, due to gradient ramping, the receiver demodulation frequency must be changed also during the readout, which may be a

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

limitation with some spectrometers. Alternatively, after sufficiently broadband acquisition, off-centers are applied by data demodulation during image reconstruction. 7.5

Contrast

The utility of medical images hinges on contrast between tissues types or conditions. MR image

F

contrast is controlled by changing sequence parameters (e.g., excitation flip angle) or adding magnetization preparation (e.g., inversion), and further modified by combining multiple images (e.g.,

O

by subtraction). In principle, in short-T2 imaging, the same concepts can be applied to generate

O

contrast. However, a specific requirement is the separation of short-T2 signals from the often contrast generation in short-T2 techniques:

PR

dominant signals of long-T2 tissues. Furthermore, some specific considerations apply concerning

1) With negligible echo times, all tissues contribute according to their spin density and steadystate magnetization, often leading to images with relatively low intrinsic contrast.

E-

2) The effect of RF pulses on magnetization depends on their duration relative to T2 [52, 131] which needs to be taken into account for excitation, refocusing, and preparation. This may be

PR

limiting in some cases but can also be utilized to achieve T2 selectivity [76, 131-135]. 3) T2* weighting can be introduced in the UTE and SPI sequences by increasing TE, an option that obviously does not exist for the ZTE approach. A related effect may be achieved by

AL

reducing the encoding duration [7]. However, this approach emphasizes T2* decay at higher spatial frequencies.

4) T1 weighting using large flip angles is most feasible for the UTE technique but with long

R N

pulses, T2 decay during the pulse may need to be considered (see Section 4.3). For ZTE and SPI excitation, restrictions to flip angles exist associated with RF amplitude, dead time, and SAR (see Section 5.5).

U

5) Most modes of magnetization preparation require more time than the TR values typical of the steady-state sequences described above. Hence, efficient performance requires a segmented

JO

implementation as employed in the Cartesian technique called “magnetization prepared rapid gradient echo imaging” (MP-RAGE) [136]. Emphasizing contrast by specific acquisition ordering as in the Cartesian case is not possible for basic UTE or ZTE scanning, as in all TR intervals equal parts of k-space are acquired. However, SPI and some of the hybrid techniques (see Section 0) enable such contrast manipulation. For the potentially SAR-intensive ZTE and SPI sequences, adding preparation pulses may be limited.

6) Subtraction images based on multiple gradient echo acquisitions with different TE can provide excellent T2* contrast. However, on radial trajectories, refocusing is sensitive to eddycurrent-induced artifacts due to concentration of errors in central k-space and requires

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

correction measures [64, 66, 137-139]. For multi-echo ZTE scanning, eddy currents can be mitigated by employing the echo-shifting principle [140] to reduce gradient switching [141, 142]. 7) With zero or ultra-short TE, the image phase will be less affected by off-resonance than in conventional echo-based sequences. However, phase accrued during RF excitation and

F

gradient encoding still contains information about local susceptibility or chemical shifts that

O

can be used as a source of image phase contrast [143-145].

Within these constraints, contrast can be generated in short-T2 imaging. Several approaches have been

O

proposed to separate short- from long-T2 signal, in particular to generate images with short-T2 signal only [38, 80, 132, 133, 146-158]. This results in so-called positive contrast for short-T2 tissues as

PR

opposed to classical T2-weighted images where these tissues are identified by the absence of signal (Fig. 18). Moreover, contrast within short-T2 components has been obtained by weighting by or quantification of T2 and T2* [53, 81, 159-163], T1 [162, 164-169], T1ρ [170, 171], chemical shift [150,

E-

172, 173], or magnetization transfer [167, 174-178]. Importantly, the last option requires suitable measures to separate true magnetization transfer from direct saturation which is stronger for lines

PR

broadened by rapid transverse relaxation [175]. 8 Artifacts and post-processing

AL

In this section, further aspects of image quality are discussed in addition to the resolution limitations and sequence-specific artifacts described in Sections 3-5. Off-resonance

R N

8.1

In radial images, the off-resonance effects due to chemical shifts, global field variations, or local susceptibility gradients deviate from the related signal shifts observed in Cartesian gradient echo

U

techniques. With the center-out acquisition scheme often used for short-T2 MRI, off-resonance phase is accrued to form a cone-like shape in two-dimensional k-space or a corresponding structure for the

JO

3D case. This leads to blurring or ringing depending on the resonance offset. The transition between the two effects occurs for offsets between 1 and 2 PWB, where the pixel bandwidth of an image is defined as PBW = BW / N and N is the number of pixels in the FOV. To reduce off-resonance due to B0 non-uniformity, shimming is usually applied. The generic means of

mitigating off-resonance artifacts of any origin is using strong gradients, i.e., large PBW, as employed anyhow for encoding short-T2 tissues. Complementary, off-resonance correction can be applied to the data based on B0 mapping, using conjugate phase reconstruction [179-182] or iterative approaches [183, 184]. As in other sequences, a range of off-resonances values within a voxel leads to signal loss

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

until the start of data acquisition and signal decay during the readout, which can only be mitigated by shortening TE and increasing encoding gradients, respectively. A particular effect in vivo is created by the chemical shift between water and fat of approximately 3.5 ppm. This corresponds to hundreds of Hertz at usual field strengths which is on the order of commonly used PBWs. If both water and fat contribute to the MR signal, the demodulation frequency

F

chosen for data acquisition and/or image reconstruction can be on resonance for only one of them,

O

leading to image blurring for the respective other component. This has a strong impact on the appearance of tissue interfaces and thus image contrast as demonstrated in Fig. 19. It is important to

O

keep these effects in mind in practice because artifacts such as those seen in the figure might be misinterpreted. Besides using higher PBW, such effects may be avoided by suppressing either fat or

8.2

PR

water using magnetization preparation or separation techniques [144, 185-190]. Resolution enhancement

E-

As outlined in Section 3.1, resolution in MRI is limited by apodization in k-space caused by T2 decay affecting the acquisition range. However, in analogy to time-domain corrections in spectroscopy

PR

[191], such apodization can be countered by suitable correction factors in k-space, corresponding to a filter operation in the image domain. In particular, T2 filtering can be applied which compensates for exponential T2 decay (Fig. 20):

 ',!- &

AL

,!- &  

(5)

When setting  ',!- &   of a given target tissue, the spatial resolution of this tissue will not be

R N

limited by the encoding duration. However, this tempting approach has two serious limitations: 1) While resolution will be improved for all  .  ',!- &, high spatial frequencies will be overemphasized for any  /  ',!- &, leading to ringing artifacts. For samples with excessive

U

T2 range, the spatial variation of T2 can be included in the signal model [192]. However, this

JO

approach requires additional T2 mapping and more extensive computation.

2) Multiplication by the exponential in Eq. (5) boosts not only high-k signal but also noise, at the expense of image SNR.

Therefore, although moderate T2 filtering may be a useful means of optimizing image appearance, it cannot replace sufficiently rapid encoding in short-T2 MRI. 9 Hardware As discussed in preceding sections, imaging of tissues with short T2 differs from common MRI in terms of acquisition strategies and sequence timing. As a consequence, it alters requirements on the

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

instrumentation used. In several respects, short-T2 imaging calls for enhanced hardware performance while in others it demands lower or different specifications. For some short-T2 techniques, dedicated hardware is not essential. Nevertheless, for best results and for conquering ever shorter T2, tailored hardware is clearly beneficial if not indispensable. Magnet

F

9.1

Today, the strength of the static B0 field of most clinical MRI systems is 1.5 T or 3 T, whereas in

O

numerous research systems it is 7 T and reaches up to 10.5 T [193]. B0 affects mainly three aspects:

O

1) Baseline sensitivity: As usual in human MRI setups, the SNR increases approximately linearly with B0 [194]. SNR is particularly critical for tissue components with short T2 due to

PR

often low spin density and long T1 (see Section 2.5) as well as short, high-bandwidth acquisitions. Thus, with regard to baseline sensitivity, high B0 is favored. 2) Relaxation times: For field strengths and tissues of interest, T1 increases along with B0 while

E-

both T2 and T2* decrease [13], eating up part of the sensitivity gain. In turn, at high field, shorter T2 calls for potent short-T2 techniques and longer T1 entails lower optimal flip angles.

PR

3) B1 efficiency: The power efficiency of generating B1 declines as B0 increases, approximately according to 0 120 0 [195, 196] with 120 being the power absorbed in the exposed tissue. This dependency renders RF power demands and SAR at high field challenging

AL

despite slightly reduced flip angle requirements.

Optimal choice of the field strength accounts for all of these aspects and also depends on the tissue

R N

and application of interest.

Concerning B0 uniformity, short-T2 techniques pose reduced demands since they use strong gradients and high-bandwidth acquisition. Therefore, per-subject shimming is often not required. It still tends to

U

be relevant, however, when weaker or ramped gradients are employed during data acquisition, such as in WASPI, UTE imaging, or gradient-modulated variants (see Section 6.2). Furthermore,

JO

magnetization preparation may require more uniform B0 than data acquisition (see Section 7.5). 9.2

Gradient system

The gradient system [197] has two main types of components, amplifiers and coils, which together produce dynamic field gradients according to input waveforms. 1) Gradient strength: The maximum gradient strength G is a key parameter as it determines kspace speed, which competes with signal decay. For purely frequency-encoded techniques (e.g. UTE or ZTE imaging) the latter causes approximately exponential apodization (Figs. 8 and 17). To actually achieve a given nominal resolution %&, the readout gradient must

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

approximately obey  3 $ %& [37, 38]. For example, for T2 = 500 μs an actual resolution of 1 mm can be obtained with G = 23.5 mT/m, which is feasible on most clinical systems, providing G is in the range 20 - 80 mT/m. However, T2 = 50 μs requires the extreme value of G = 235 mT/m or will yield only 10 mm resolution at the lower G. For purely phaseencoded techniques (e.g. SPI), resolution is not limited by T2, but large G is still important to

F

limit TE for best SNR (see Section 3.3 and Fig. 17).

O

2) Slew rate: Techniques that set readout gradients before RF excitation (e.g., SPI and ZTE) may operate gradients quasi-continuously with minimized switching (see Section 5) and thus do

O

not depend on high slew rates. Still, faster gradient switching permits shorter TR, which is associated with smaller optimal flip angle, higher acquisition duty cycle, and less overall scan

PR

time (see Section 7.1). In contrast, when RF excitation is performed first (e.g. in UTE), the speed of gradient ramping is a key determinant of resolution and robustness against offresonance (see Section 4.1, Fig. 8). Slew rates in commercial scanners range up to 200

E-

mT/m/ms in combination with a gradient strength of, e.g., 40 mT/m. Importantly, to efficiently use an increase in gradient strength, the slew rate must grow quadratically [119].

PR

3) Duty cycle: Quasi-continuous operation of ZTE- and SPI-like sequences means that gradients are always on, often at constant net magnitude, but with changing direction. This can only be accomplished if full duty cycle is provided at the applied high strength, which is a rather

AL

uncommon design criterion for gradient systems. Although the same techniques can also be run non-continuously with gradient switching and reduced duty cycle, the continuous mode is clearly preferable as it permits shorter TR and strongly reduces eddy currents and acoustic

R N

noise. For methods with intrinsic gradient switching (e.g. UTE imaging), high duty cycle is less essential but still helps minimize TR.

4) Amplifier: Requirements for the gradient amplifier arise from the three specifications

U

discussed above. Gradient strength and slew rate are directly related to maximum current and voltage, respectively. Some amplifiers offer added flexibility by permitting multiple

JO

combinations of peak current and peak voltage [198]. Operation at full duty cycle requires that the current necessary for a desired gradient strength can be provided continuously along with the corresponding power for given ohmic losses in the gradient coil.

5) Fidelity: The spatiotemporal fidelity of gradient fields is generally hampered by delays, finite bandwidth, eddy currents, and amplifier imperfections. This affects image quality primarily in UTE techniques (see Section 4.2) and when additional echoes are generated for creating contrast (see Section 7.5). To a certain degree, these effects can be reduced by delay calibration and pre-emphasis. If necessary, residual field errors can be measured and accounted for at the level of image reconstruction (see Section 9.5).

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

6) Peripheral nerve stimulation: In vivo, gradient switching must avoid peripheral nerve stimulation (PNS) induced by rapid, large-amplitude field changes [199]. In this respect, short, localized gradient coils are advantageous over whole-body gradients [197]. Short-T2 techniques with minimized gradient switching are not a problem with regard to PNS. Sequences with switching need more attention but are typically less prone to PNS than, e.g.,

F

echo-planar imaging with long trains of shortly spaced rapid field changes.

O

In summary, fulfilling all of the listed criteria at the same time is not trivial but it is possible with

9.3

O

dedicated designs, thus also enabling extreme applications [200]. RF chains

PR

Short-T2 MRI relies on high-bandwidth, sub-millisecond acquisitions and thus also alters requirements on RF hardware. For UTE imaging, the factors to consider are rapid T/R switching and avoiding background signal. For ZTE and SPI techniques, which require excitation during readout

power is needed (see Section 5.5).

E-

gradients, these demands are further increased, and in addition, efficient RF transmission at high

PR

1) Amplifier: For ZTE and SPI, the bandwidth of RF amplifiers must be sufficiently large to enable high-bandwidth excitation using either short hard pulses of duration in the μs range or frequency-swept pulses with large sweep width (see Section 5.5). Moreover, rapid amplitude

AL

changes at the beginning and end of these pulses call for sub-microsecond rise times. In analogy to gradient amplifiers, attenuation in transition bands may be overcome by preemphasis as far as permitted by power margins. Moreover, accurate realization of modulated

R N

pulses requires good linearity which may also be improved by means of calibrated compensation. High-power transmission is particularly important to achieve adequate flip angles by hard-pulse excitation. As an example, when using an RF coil providing a typical B1

U

efficiency of 20 45 67, obtaining a flip angle of 5° with a hard pulse of 2 μs duration

JO

requires 2.9 kW. Furthermore, the available duty cycle must permit application of the desired pulses (few to hundreds of microseconds) at short TR (potentially sub-millisecond, see Section 7.1) and for long scan times (potentially tens of minutes). Amplifier ring-down after the RF pulse causes artifacts in high-bandwidth images when time constants are similar to T/R switching times [201]. To suppress this ring-down sufficiently, blanking after transmission must be fast and achieve a high degree of isolation in addition to the attenuation provided by T/R switching or de-tuning (see below).

2) Coils: MRI generally relies on either a single coil for both RF transmission and detection or separate devices for the two purposes. The latter setup usually consists of a receiver array

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

near the subject to maximize SNR and a surrounding transmit coil. However, those short-T2 techniques with high B1 requirements call for particularly high transmit efficiency, which can favor T/R coils, in particular T/R arrays. A further advantage of such single-coil setups are reduced space requirements along the radial direction. On the other hand, for safety reasons, the distance to the subject needs to be larger than for a receive-only coil. Furthermore, pre-

F

amplifier decoupling cannot be implemented. On the receive side, as usual, coil arrays offer

O

SNR benefit [202] and parallel imaging capability [97]. In ZTE techniques, complementary spatial encoding by array detection also holds potential to support filling the k-space gap (see

O

Section 6.1). Similar to RF amplifiers, the bandwidth of RF coils must be large enough to generate and detect, at sufficient uniformity, suitable excitation pulses and resonance signals,

PR

respectively. If not provided by conventional coils, dedicated designs are required. For residual imperfections, appropriate pulse compensation or data correction procedures may be employed. A further design criterion arises from the original aim to capture short-T2 signals

E-

from the subject. In addition to the latter, short-T2 detection is prone to capture MR signals also from hardware parts, which conventional MRI does not pick up due to their rapid decay.

PR

This issue is most prominent for materials close to coil conductors where B1 is high (see Section 5.4). Ideally, such contributions are avoided by choosing materials that do not contain the targeted nucleus. For 1H imaging this is achieved with, e.g., glass or fluorinated polymers such as polytetrafluoroethylene (PTFE) [82, 203-205].

AL

3) T/R switching: In pulsed MR, the receive chain must be sufficiently isolated from the transmit path in both operation modes. During transmission, destruction and preferably even saturation

R N

of receive pre-amplifiers must be prevented. During reception, residual ring-down of the power amplifier must be blocked. Depending on the coil configuration, either switching between the two operation modes by means of a T/R switch or de-tuning of the transmit coil and tuning of the receive coil are required. The particular challenge for short-T2 techniques is

U

to perform this rapidly, i.e., in the low microsecond range as opposed to several tens of

JO

microseconds in commercial systems, to minimize the RF dead time [206, 207] (see Section 5.2). To avoid image artifacts, switching transients must have decayed to a level well below that of the MR signal when data is acquired. Rapid switching is particularly challenging at high peak and average power. It is typically based on positive-intrinsic-negative (PIN) diodes [208], but has also been achieved with transistors [208, 209] or micro-electro-mechanical systems (MEMS) [210]. In a recent design based on PIN diodes, switching transients are cancelled by symmetric biasing [211].

4) Simultaneous operation: A tempting alternative to RF switching is simultaneous transmission and reception [212-214]. With this approach, dead time overheads are avoided and sequences

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

are simplified. Simultaneous transmission and reception was common practice before the advent of pulsed NMR, albeit in rigid setups and at much lower frequencies. For in-vivo imaging at present-day MR frequencies, sufficient isolation between transmit and receive paths is difficult to achieve without T/R switching or detuning. Isolation between the two RF paths can be achieved by geometric coil decoupling, additional decoupling by passive

F

electronics, and active cancelation of transmit leakage into the receive path. However, the

O

robustness of such setups is limited by residual variation in RF conditions, particularly by patient motion, which causes subtle changes in RF loading. Feedback control of cancelation

O

has been proposed to address this limitation [215, 216]. In an alternative approach, isolation is facilitated by sideband modulation, which permits performing transmission and reception in

9.4

PR

different frequency bands [217]. Console

E-

The console of an MRI system includes a sequencer generating RF and gradient waveforms, a receiver for data acquisition, and a reconstructor that calculates images from raw data. This

PR

architecture is the same for short-T2 imaging. However, for high performance in the short-T2 mode, the components should meet more stringent specifications. 1) Sequencer: The temporal resolution of the sequence description should be 1 μs or below. The

AL

definition and generation of RF pulses must be more flexible than in common clinical MRI systems, including a wider range of pulse durations (1 μs to tens of milliseconds) and TR (down to 100 μs, see Section 7.1). To perform off-center imaging with ZTE or SPI (see

R N

Section 7.4), per-shot adjustment of the excitation carrier frequency is necessary. The same is true for the receiver with any radial technique.

2) Receiver: The feasible acquisition bandwidth must cover the signal bandwidth of up to the

U

megahertz range. Even higher acquisition bandwidth is needed for oversampling to enable finite-support extrapolation (see Section 5.3). In conjunction with high spatial resolution and

JO

array detection, 3D acquisition with oversampling yields very large amounts of data that need to be handled and stored. High bandwidth and short TR render these data challenging also in terms of throughput. Short-T2 imaging benefits particularly from flexibility at the level of digital filtering, where short, customized, and multi-rate filters [218] permit minimizing RF dead times without bloating the raw data (see Section 5.2).

3) Reconstructor: 3D radial data sets tend to be very large and do not readily lend themselves to splitting into smaller portions for separate processing. Therefore, a computer used for image reconstruction should be equipped with ample memory to avoid time-consuming swapping on the hard disk. For example, the raw data for image matrix size 256 acquired with full Nyquist

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

angular sampling, 4-fold radial oversampling, and 8 channels amount to approximately 13 GB. For 3D gridding reconstruction, about 8 times that amount of memory needs to be allocated. Reconstruction time is longer than in standard Cartesian imaging and is determined by the numerical complexity of algebraic parts, density correction, gridding interpolation, and FT, and by the number of cycles in the case of iterative reconstruction (see Sections 4.5 and

F

5.3). For speed-up, reconstruction may be partly parallelized and distributed to different nodes

Calibration tools

O

9.5

O

on a single computer or a cluster of computers, or on graphical processing units (GPUs).

All magnetic fields – static, audio- and radio-frequency – experienced by the spins contribute to

PR

encoding of the MR signal. Hence, for proper image reconstruction, the spatial and temporal behavior of these fields must be accurately included in the underlying signal model. However, encoding fields are produced with finite accuracy by the respective hardware, i.e., the magnet, gradient subsystems,

E-

and RF chains. To determine actual fields and field evolutions, either the MR system itself or additional dedicated measurement hardware may be used to record field information at runtime or in a

PR

calibration step. Depending on the severity of deviations, the results can be used to improve image reconstruction or to pre-distort input waveforms for higher accuracy. 1) Gradients: Apart from the methods discussed in Section 4.2, spatial and temporal behavior of

AL

gradient fields per se as well as eddy current effects can be measured by a set of MR field probes [71]. A field probe is a small RF coil wrapped around an NMR sample and determines the variation over time of field strength at one position. Several such probes enable spatial

R N

expansion of the field dynamics in the imaging volume. Such information has been used to improve reconstruction of, e.g., UTE images [219, 220] and dual-echo ZTE data [139]. Instead of measuring field behavior for each protocol, by assuming linearity and time

U

invariance, the gradient system may also be characterized by its impulse response [221], enabling field calculations for arbitrary gradient shapes.

JO

2) RF: Effective RF waveforms are subject to variable degrees of distortion caused by RF amplifiers, lines, cables, and coils. These are most easily observed with a pick-up coil placed near the sample and connected to an MR spectrometer or a spectrum analyzer. In this way, effective RF pulses as experienced by the spins are observed, including all disturbing effects in the transmit chain. Alternatively, field probes may be used in an advanced setup to monitor RF and gradient fields simultaneously [222], including information on spatial variation. Knowledge of effective RF waveforms is particularly important for the use of amplitude- or phase-modulated pulses (see Section 5.5).

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

9.6

Summary

The hardware considerations compiled here give an indication of how dedicated systems for short-T2 MRI might differ from present-day instrumentation. For optimal UTE imaging, requirements are actually rather similar to those of conventional MRI. High B0 should be combined with high gradient performance in terms of amplitude, slew rate, and

F

fidelity, and careful consideration of PNS. With regard to RF, in a similar manner to regular MRI,

O

dense receiver arrays are desirable for optimal SNR. As one feature specific to the short-T2 scenario, a UTE system should offer rapid T/R switching to minimize TE. Integrated field monitoring equipment

O

would facilitate accurate gradient characterization.

PR

In contrast, a system designed for ZTE imaging would differ more fundamentally. Due to the more limiting role of RF power considerations, the optimal field strength will be lower for ZTE imaging than for UTE. B0 field homogeneity is less critical in the ZTE scenario and may be traded off against

E-

magnet length and cost. Gradients should offer high amplitude and, critically, unlimited duty cycle. However, neither high slew rate nor gradient fidelity are particularly important, which will alter and

PR

potentially simplify coil design and amplifier considerations. Likewise, PNS is of little concern. On the RF side, the ZTE mode is more demanding, calling not only for dense array detection but also for high-power transmit capability, efficient transmission, very rapid T/R switching, and 1H-free coil front-ends. A ZTE console should offer superior temporal definition, whereas integrated RF

AL

monitoring would improve waveform correction and thus expand the range of feasible pulse modulation schemes.

R N

10 Applications

Most applications of short-T2 MRI relate to its ability to depict solid tissues. This is in contrast to conventional MRI, which depicts such tissues mostly by the absence of signal, relying on so-called

U

negative contrast. By using dedicated short-T2 techniques instead, signals from solid tissues can be

JO

observed directly and thus with positive contrast. This approach has two major advantages: 1) Short-T2 tissues can be distinguished from actually void spaces.

2) Signal variation related to spin density or other sources of contrast (see Section 7.5) permits differentiation between different short-T2 tissues or between healthy and pathological tissues.

Table 2 lists several tissues and materials of interest along with T2 values as reported in the literature, indicating fields of application. With conventional MRI hardware, the short-T2 techniques described previously enable imaging of T2 values with several hundreds of microseconds at useful resolution. By employing dedicated RF and gradient hardware (see Section 9), T2 values below 100 μs also become realistic targets [201, 223].

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

1) Musculoskeletal (MSK) system: Short-T2 tissues are particularly prevalent in the MSK system, such as bone, tendons, ligaments, menisci and cartilage (Fig. 21a). The focus of many investigations is to image bone and quantify bone mineral density, either indirectly from bone water [49, 51, 53, 111, 151, 160, 168, 172, 176, 224-246] or directly from hydroxyapatite employing

31P-imaging

[6, 22, 244, 247-257], due to its potentially high diagnostic value,

F

e.g., in osteoporosis. This application has also been expanded to bone substitutes [258].

O

Furthermore, tendons [31, 38, 150, 170, 259-261], ligaments [49, 54], menisci [49, 51, 150, 170, 262], and cartilage [263-265] are of specific interest due to frequent injuries,

O

inflammations, or wear.

2) Positron emission tomography (PET): The ability to depict bone can be utilized in devices for

PR

combined PET and MRI. PET requires attenuation correction for bone and other tissues which may be based on segmentation of short-T2 MR images [266-273]. The same approach can be applied for attenuation caused by hardware parts of the MRI scanner, in particular RF

E-

coils.

3) Pulmonary MRI with short-T2 techniques [47, 99, 108, 162, 167, 274-301] (Fig. 21b) has

PR

large potential, e.g., to improve specificity in early diagnosis of certain diseases (e.g., fibrosis) as compared to computer tomography (CT). For this application, a particular advantage of the radial short-T2 techniques is high robustness with respect to motion. 4) Iron oxide nanoparticles (IONP) are used as contrast agents for, e.g., cell labeling or tracking

AL

drug delivery. Typically, in MR images IONPs are observed as signal voids due to reduced T2*. Alternatively, short-T2 techniques enable imaging them with positive contrast, thus

R N

utilizing reduced T1 instead of T2* [7, 149, 165, 302-306]. In this way, their identification and localization can be improved.

5) Dental MRI using short-T2 techniques allows depiction of teeth with high quality (Fig. 21c). Application in vivo holds promise for improved early diagnosis of, e.g., caries lesions and

U

benefits from dedicated hardware [307-315]. However, at present patient comfort, scan times,

JO

and cost cannot yet compete with standard X-ray techniques.

6) Myelin is a biomaterial that is essential in rapid electrical signal propagation along nerve fibers in white matter of the brain and the spinal cord. Therefore demyelination plays an important role in numerous neurological disorders. In a conventional approach, myelin is quantified indirectly via myelin water [316, 317]. Short-T2 techniques have the potential to directly visualize and quantify myelin but must rely on effective suppression of long-T2 tissue components [35, 147, 318-325].

7) Plaques in vessels due to calcification by arteriosclerosis reduce or obstruct blood flow in, e.g., peripheral, cardiac, or brain vasculature. Due to the short T2* of plaques, their assessment

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

is conventionally based on X-ray imaging. Short-T2 MRI offers a prospective alternative modality for plaque detection and specification [319, 326-329]. 8) Functional MRI (fMRI): MR techniques with ultra-short or zero TE are not intrinsically sensitive to blood-oxygenation-level dependent (BOLD) contrast. Nevertheless, they can be used for performing fMRI. Sensitivity to neuronal activation is either obtained by T2

F

preparation [330] or based on other mechanisms, e.g., perfusion [331, 332]. As a particular

O

advantage, the zero-TE techniques offer silent operation, avoiding brain activation related to acoustic noise generated by gradient switching.

O

9) Electric property tomography (EPT): Spatial mapping of electric permittivity and conductivity by EPT has potential for various applications [333]. EPT is based on evaluation

PR

of the image phase. Techniques with negligible TE may contribute in this field due a reduced sensitivity of their image phase to off-resonance and eddy currents [334]. 10) Sodium MRI: Sodium is of interest due to its important role in many biochemical processes

E-

and is a potential biomarker in the diagnosis of diseases such as stroke, cancer, and osteoarthritis. The NMR-active isotope 23Na has spin 3/2, leading to bi-exponential relaxation

PR

behavior (see Section 2.3). T2 of the shorter component is typically best addressed with dedicated short-T2 techniques [126, 129, 335, 336]. 11) Materials: Apart from in-vivo uses, short-T2 MRI has large potential for materials research [337-340] and inspection as well as for applications in the food industry [341-343]. The

AL

ability to depict short-T2 materials is also instrumental in the non-destructive examination of mummies [344].

R N

12) Inhomogeneous B0: The use of minimized TE and high bandwidth, short-T2 techniques also provides considerable resilience against signal loss and displacement from strong macroscopic B0 (see Section 2.6). This kind of robustness improves imaging in the presence

U

of tissue-air interfaces or metal implants [7, 345].

JO

11 Conclusions

In this article, we have reviewed the growing family of techniques for in-vivo MRI of tissues with rapid decay of transverse magnetization. Besides covering the variety of current approaches and implementations, this survey has taken a conceptual perspective, emphasizing fundamental challenges, derived principles, as well as commonalities and key differences among the various techniques. It has been shown that, despite their variety, short-T2 imaging techniques can be categorized according to a few fundamental design choices, leading to four basic sequence ideas, CTI, SPI, UTE, and ZTE. The Cartesian techniques CTI and SPI offer the attractive feature of enabling theoretically unlimited

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

spatial resolution even in the case of rapid signal decay. However, this option comes at the expense of poor SNR. Greater SNR can be achieved with the radial UTE and ZTE sequences, which are most widely used at this point. UTE imaging is feasible with conventional hardware but requires accurate gradient calibration. In comparison, ZTE approaches offer higher short-T2 sensitivity and resolution but are more demanding in terms of RF hardware and image reconstruction. The choice of technique

F

thus depends on the targeted application and technical resources. Offering guidance for these

O

considerations is one of the chief purposes of this article.

As an outlook, two roads can be foreseen in the emerging field of short-T2 MRI. First, clinical

O

usefulness needs to be investigated for applications that are readily accessible with existing commercial scanners. These relate to tissues with T2 above approximately 300 μs, such as in MSK

PR

and pulmonary imaging. The basis for this task will be relatively established implementations of both UTE and ZTE variants, whereas further effort will be required to verify reliable contrast options. Second, the feasibility of imaging in vivo at T2 values below a few hundreds of microseconds,

E-

possibly down to tens of microseconds (as, e.g., in myelin) still needs to be explored and requires additional methodological and technical efforts. With regard to imaging methods, it is expected that

PR

there is still considerable potential in fusing and optimizing current approaches. Furthermore, on the technical side, advances in high-bandwidth capabilities of MR scanner hardware with respect to gradient strength and duty cycles, RF switching speed, console specifications, and proton-free

12 Appendices

AL

materials will pave the way for further discoveries in the world of short-T2 tissues.

R N

12.1 Implementation of the rho filter

For sequences with radial encoding, a suitable DCF is efficiently obtained by calculation of the rho filter. For center-out encoding schemes as used in UTE and ZTE imaging, it may be implemented as

U

follows:

JO

89    

:

89!  !    !

(A1) :

!  ! 

89!  ! ! 

:

:

!  ! 

(A2) :

(A3)

Here, !  ;;! counts the radial k-space positions starting from the center, k-values are given in arbitrary units, and d is the dimensionality of the data. Note that the DCF calculated according to Eqs. (A1-A3) is valid also for  not being at    as well as non-equidistant k-space positions.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

12.2 Steady-state magnetization at short TR In a steady-state sequence with complete spoiling of transverse magnetization after each TR interval, the flip angle for maximized transverse steady-state magnetization ' is the Ernst angle [118] (  <=>>?@

(A4)

F

with )





(A5)

O

 

O

leading to [346]  

'  

PR

  

For ) *  the approximation [119] ) 

E-

 # 



( # &AABC 

) 

#

as well as

PR

leads to

) 



AL

' # 

 ) 

(A6)

(A7)

(A8)

(A9)

R N

which by neglecting the -1 further simplifies to ' # 

) 

(A10)

JO

U

Hence at very short TR, both the Ernst angle and the signal are approximately proportional to ).

Funding

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Declarations of interest MW None

F

Declarations of interest KP

O

Klaas Pruessmann holds a research agreement with and receives research support from Philips Healthcare. He is a shareholder of Gyrotools LLC and Skope

PR

O

Magnetic Resonance Technologies.

References

E-

[1] P. Mansfield, Pulsed NMR in solids, Prog Nucl Mag Res Sp, 8 (1971) 43-101. [2] E.R. Andrew, E. Szczesniak, A historical account of NMR in the solid state, Prog Nucl Mag Res Sp, 28 (1995) 11-36.

PR

[3] P. Jezzard, J.J. Attard, T.A. Carpenter, L.D. Hall, Nuclear magnetic resonance imaging in the solid state, Prog Nucl Mag Res Sp, 23 (1991) 1-41. [4] J.B. Miller, NMR imaging of materials, Prog Nucl Mag Res Sp, 33 (1998) 273-308. [5] S. Emid, J. Creyghton, High resolution NMR imaging in solids, Physica B+ C, 128 (1985) 81-83.

AL

[6] M.D. Robson, P.D. Gatehouse, M. Bydder, G.M. Bydder, Magnetic resonance: an introduction to ultrashort TE (UTE) imaging, J Comput Assist Tomogr, 27 (2003) 825-846.

R N

[7] M. Weiger, K.P. Pruessmann, MRI with zero echo time, in: eMagRes, John Wiley & Sons, Ltd, Chichester, United Kingdom, 2012, pp. 311-322. [8] G.M. Bydder, G.D. Fullerton, I.R. Young, MRI of tissues with short T2s or T2*s, John Wiley & Sons Ltd, Chichester, United Kingdom, 2012.

U

[9] M.H. Levitt, Spin dynamics: basics of nuclear magnetic resonance, John Wiley & Sons, Chichester, United Kingdom, 2001.

JO

[10] I. Solomon, Relaxation processes in a system of two spins, Phys Rev, 99 (1955) 559. [11] N. Bloembergen, E.M. Purcell, R.V. Pound, Relaxation effects in nuclear magnetic resonance absorption, Phys Rev, 73 (1948) 679. [12] A.G. Redfield, On the theory of relaxation processes, IBM J Res Dev, 1 (1957) 19-31. [13] D.D. Traficante, Relaxation: an introduction, in: eMagRes, John Wiley & Sons, Ltd, Chichester, United Kingdom, 2007. [14] J.C. Gore, A.W. Anderson, The physics of relaxation, in: eMagRes, John Wiley & Sons, Ltd, 2007. [15] A.G. Redfield, Relaxation theory: density matrix formulation, in: eMagRes, John Wiley & Sons, Ltd, 2007.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[16] P.A. Beckmann, Spectral densities and nuclear spin relaxation in solids, Physics reports, 171 (1988) 85-128. [17] G.E. Pake, Nuclear resonance absorption in hydrated crystals: fine structure of the proton line, J Chem Phys, 16 (1948) 327-336. [18] D.D. Laws, H.M.L. Bitter, A. Jerschow, SolidDState NMR spectroscopic methods in chemistry, Angew Chem Int Edit, 41 (2002) 3096-3129.

F

[19] D.D. Traficante, Relaxation. Can T2 be longer than T1?, Concept Magn Reson A, 3 (1991) 171177.

O

[20] H.C. Torrey, Bloch equations with diffusion terms, Phys Rev, 104 (1956) 563-565.

O

[21] P.T. Callaghan, Principles of nuclear magnetic resonance microscopy, Oxford University Press, 1993.

PR

[22] S. Anumula, J. Magland, S. Wehrli, H. Zhang, H. Ong, H. Song, F. Wehrli, Measurement of phosphorus content in normal and osteomalacic rabbit bone by solidDstate 3D radial imaging, Magn Reson Med, 56 (2006) 946-952. [23] F.W. Wehrli, M.A. Fernández-Seara, Nuclear magnetic resonance studies of bone water, Annals of biomedical engineering, 33 (2005) 79-86.

E-

[24] G. Fullerton, I. Cameron, V. Ord, Orientation of tendons in the magnetic field and its effect on T2 relaxation times, Radiology, 155 (1985) 433-435.

PR

[25] S. Peto, P. Gillis, V.P. Henri, Structure and dynamics of water in tendon from NMR relaxation measurements, Biophys J, 57 (1990) 71-84. [26] S.J. Erickson, R.W. Prost, M. Timins, The" magic angle" effect: background physics and clinical relevance, Radiology, 188 (1993) 23-25.

AL

[27] R.M. Henkelman, G.J. Stanisz, J.K. Kim, M.J. Bronskill, Anisotropy of NMR properties of tissues, Magn Reson Med, 32 (1994) 592-601. [28] W. Gründer, MRI assessment of cartilage ultrastructure, NMR Biomed, 19 (2006) 855-876.

R N

[29] G.D. Fullerton, A. Rahal, Collagen structure: the molecular source of the tendon magic angle effect, J Magn Reson Imag, 25 (2007) 345-361. [30] M. Bydder, A. Rahal, G.D. Fullerton, G.M. Bydder, The magic angle effect: a source of artifact, determinant of image contrast, and technique for imaging, J Magn Reson Imag, 25 (2007) 290-300.

U

[31] J. Du, B.C. Pak, R. Znamirowski, S. Statum, A. Takahashi, C.B. Chung, G.M. Bydder, Magic angle effect in magnetic resonance imaging of the Achilles tendon and enthesis, Magn Reson Imag, 27 (2009) 557-564.

JO

[32] H.J.C. Berendsen, Nuclear Magnetic Resonance Study of Collagen Hydration, J Chem Phys, 36 (1962) 3297-3305. [33] L.V. Krasnosselskaia, G.D. Fullerton, S.J. Dodd, I.L. Cameron, Water in tendon: Orientational analysis of the free induction decay, Magn Reson Med, 54 (2005) 280-288. [34] H. Wennerström, Proton nuclear magnetic resonance lineshapes in lamellar liquid crystals, Chem Phys Lett, 18 (1973) 41-44. [35] M.J. Wilhelm, H.H. Ong, S.L. Wehrli, C. Li, P.H. Tsai, D.B. Hackney, F.W. Wehrli, Direct magnetic resonance detection of myelin and prospects for quantitative imaging of myelin density, Proc Natl Acad Sci, 109 (2012) 9605-9610. [36] J.M. Wild, H. Marshall, M. Bock, L.R. Schad, P.M. Jakob, M. Puderbach, F. Molinari, E.J.R. Van Beek, J. Biederer, MRI of the lung (1/3): methods, Insights Imaging, 3 (2012) 345-353.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[37] P.T. Callaghan, C.D. Eccles, Sensitivity and resolution in NMR imaging, J Magn Reson, 71 (1987) 426-445. [38] J. Rahmer, P. Bornert, J. Groen, C. Bos, Three-dimensional radial ultrashort echo-time imaging with T2 adapted sampling, Magn Reson Med, 55 (2006) 1075-1082. [39] S. Choi, X.W. Tang, D. Cory, Constant time imaging approaches to NMR microscopy, Int J Imag Syst Tech, 8 (1997) 263-276.

F

[40] S. Gravina, D.G. Cory, Sensitivity and resolution of constant-time imaging, J Magn Reson Ser B, 104 (1994) 53-61.

O

[41] I.V. Mastikhin, B.J. Balcom, Centric SPRITE MRI of biomaterials with short T2, in: L. John Wiley & Sons (Ed.) eMagRes, Chichester, United Kingdom, 2012, pp. 783-788.

O

[42] B.J. Balcom, R.P. MacGregor, S.D. Beyea, D.P. Green, R.L. Armstrong, T.W. Bremner, Singlepoint ramped imaging with T1 enhancement (SPRITE), J Magn Reson Ser A, 123 (1996) 131-134.

PR

[43] T. Speidel, V. Rasche, Accelerated adaptive quasi-random single-point imaging for in vivo artefact reduction, In Proceedings of the 25th Annual Scientific Meeting of ISMRM, Honolulu, Hawaii, USA, 2017, p.4026

E-

[44] P.J. Prado, B.J. Balcom, S.D. Beyea, R.L. Armstrong, T.W. Bremner, Concrete thawing studied by single-point ramped imaging, Solid State Nucl Mag, 10 (1997) 1-8.

PR

[45] M. Bruschewski, H. Kolkmannn, K. John, S. Grundmann, Phase-contrast single-point imaging with synchronized encoding: a more reliable technique for in vitro flow quantification, Magn Reson Med, (2018). [46] M. Weiger, D. Schmidig, S. Denoth, C. Massin, F. Vincent, M. Schenkel, M. Fey, NMR microscopy with isotropic resolution of 3.0 μm using dedicated hardware and optimized methods, Concept Magn Reson B, 33B (2008) 84-93.

AL

[47] C.J. Bergin, J.M. Pauly, A. Macovski, Lung parenchyma - projection reconstruction MR imaging, Radiology, 179 (1991) 777-781. [48] G.H. Glover, J.M. Pauly, K.M. Bradshaw, B-11 Imaging with a 3-dimensional reconstruction method, J Magn Reson Imag, 2 (1992) 47-52.

R N

[49] P. Gatehouse, G. Bydder, Magnetic resonance imaging of short T2 components in tissue, Clinical radiology, 58 (2003) 1-19.

U

[50] J.E. Holmes, G.M. Bydder, MR imaging with ultrashort TE (UTE) pulse sequences: basic principles, Radiography, 11 (2005) 163-174.

JO

[51] M.D. Robson, G.M. Bydder, Clinical ultrashort echo time imaging of bone and other connective tissues, NMR Biomed, 19 (2006) 765-780. [52] D.J. Tyler, M.D. Robson, R.M. Henkelman, I.R. Young, G.M. Bydder, Magnetic resonance imaging with ultrashort TE (UTE) PULSE sequences: technical considerations, J Magn Reson Imag, 25 (2007) 279-289. [53] J. Du, G.M. Bydder, Qualitative and quantitative ultrashort-TE MRI of cortical bone, NMR Biomed, 26 (2013) 489-506. [54] E.Y. Chang, J. Du, C.B. Chung, UTE imaging in the musculoskeletal system, J Magn Reson Imag, 41 (2015) 870-883. [55] B.H. Suits, D. White, NMR imaging in solids, Solid State Commun, 50 (1984) 291-295. [56] S. Hafner, Fast imaging in liquids and solids with the back-projection low-angle shot (BLAST) technique, Magn Reson Imag, 12 (1994) 1047-1051.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[57] D.P. Madio, I.J. Lowe, Ultra-fast imaging using low flip angles and FIDs, Magn Reson Med, 34 (1995) 525-529. [58] D.O. Kuethe, A. Caprihan, E. Fukushima, R.A. Waggoner, Imaging lungs using inert fluorinated gases, Magn Reson Med, 39 (1998) 85-88. [59] F. Hennel, M. Weiger, M.B. Rösler, R. Luechinger, B. Wilm, K.P. Pruessmann, Cartesian UTE with echo time below 0.5 millisecond, In Proceedings of the 25th Annual Scientific Meeting of ISMRM, Honolulu, Hawaii, USA, 2017, p.1608

O

F

[60] X. Deligianni, P. Bär, K. Scheffler, S. Trattnig, O. Bieri, HighDresolution FourierDencoded subDmillisecond echo time musculoskeletal imaging at 3 Tesla and 7 Tesla, Magn Reson Med, 70 (2013) 1434-1439.

O

[61] I.C. Atkinson, A. Lu, K.R. Thulborn, Characterization and correction of system delays and eddy currents for MR imaging with ultrashort echoDtime and timeDvarying gradients, Magn Reson Med, 62 (2009) 532-537.

PR

[62] B.J. Wilm, C. Barmet, M. Pavan, K.P. Pruessmann, Higher order reconstruction for MRI in the presence of spatiotemporal field perturbations, Magn Reson Med, 65 (2011) 1690-1701.

E-

[63] P.M. Margosian, T. Takahashi, M. Takizawa, Practical Implementation of UTE Imaging, in: eMagRes, John Wiley & Sons, Ltd, Chichester, United Kingdom, 2007, pp. 297-309. [64] V. Rasche, D. Holz, R. Proksa, MR fluoroscopy using projection reconstruction multiDgradientDecho (prMGE) MRI, Magn Reson Med, 42 (1999) 324-334.

PR

[65] K.-H. Herrmann, M. Krämer, J.R. Reichenbach, Time efficient 3D radial UTE sampling with fully automatic delay compensation on a clinical 3T MR scanner, PloS one, 11 (2016) e0150371. [66] D.C. Peters, J.A. Derbyshire, E.R. McVeigh, Centering the projection reconstruction trajectory: reducing gradient delay errors, Magn Reson Med, 50 (2003) 1-6.

AL

[67] P. Latta, Z. Starčuk, M.L.H. Gruwel, M.H. Weber, B. Tomanek, K-space trajectory mapping and its application for ultrashort echo time imaging, Magn Reson Imag, 36 (2017) 68-76.

R N

[68] S.J. Vannesjo, N.N. Graedel, L. Kasper, S. Gross, J. Busch, M. Haeberlin, C. Barmet, K.P. Pruessmann, Image reconstruction using a gradient impulse response model for trajectory prediction, Magn Reson Med, 76 (2016) 45-58. [69] J.H. Duyn, Y.H. Yang, J.A. Frank, J.W. van der Veen, Simple correction method for k-space trajectory deviations in MRI, J Magn Reson, 132 (1998) 150-153.

U

[70] Y.T. Zhang, H.P. Hetherington, E.M. Stokely, G.F. Mason, D.B. Twieg, A novel k-space trajectory measurement technique, Magn Reson Med, 39 (1998) 999-1004.

JO

[71] C. Barmet, N. De Zanche, K.P. Pruessmann, Spatiotemporal magnetic field monitoring for MR, Magn Reson Med, 60 (2008) 187-197. [72] J. Pauly, D. Nishimura, A. Macovski, Slice-selective excitation for very short T2 species, In Proceedings of the 8th Annual Meeting of SMRM, Amsterdam, The Netherlands, 1989, p.28 [73] J.M. Pauly, Selective excitation for ultrashort echo time imaging, in: eMagRes, John Wiley & Sons, Ltd, Chichester, United Kingdom, 2007, pp. 381-387. [74] S. Josan, J.M. Pauly, B.L. Daniel, K.B. Pauly, Double half RF pulses for reduced sensitivity to eddy currents in UTE imaging, Magn Reson Med, 61 (2009) 1083-1089. [75] K.D. Harkins, M.D. Does, W.A. Grissom, Iterative Method for Predistortion of MRI Gradient Waveforms, IEEE Trans Med Imag, 33 (2014) 1641-1647.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[76] M. Carl, M. Bydder, J. Du, A. Takahashi, E. Han, Optimization of RF excitation to maximize signal and T2 contrast of tissues with rapid transverse relaxation, Magn Reson Med, 64 (2010) 481490. [77] F. Springer, G. Steidle, P. Martirosian, C.D. Claussen, F. Schick, Effects of in-pulse transverse relaxation in 3D ultrashort echo time sequences: analytical derivation, comparison to numerical simulation and experimental application at 3T, J Magn Reson, 206 (2010) 88-96. [78] F. Hennel, The effective phase of soft RF pulses, Concept Magn Reson A, 43 (2014) 127-137.

O

F

[79] M. Weiger, K.P. Pruessmann, F. Hennel, MRI with zero echo time: hard versus sweep pulse excitation, Magn Reson Med, 66 (2011) 379-389.

O

[80] J. Rahmer, U. Blume, P. Bornert, Selective 3D ultrashort TE imaging: comparison of "dual-echo" acquisition and magnetization preparation for improving short-T2 contrast, Magn Reson Mater Phy, 20 (2007) 83-92.

PR

[81] J. Du, Quantitative ultrashort TE (UTE) imaging of short T2 tissues, in: eMagRes, John Wiley & Sons, Ltd, Chichester, United Kingdom, 2007, pp. 801-810. [82] R.A. Horch, K. Wilkens, D.F. Gochberg, M.D. Does, RF coil considerations for short-T2 MRI, Magn Reson Med, 64 (2010) 1652-1657.

E-

[83] J.I. Jackson, C.H. Meyer, D.G. Nishimura, A. Macovski, Selection of a convolution function for Fourier inversion using gridding, IEEE Trans Med Imag, 10 (1991) 473-478.

PR

[84] P.J. Beatty, D.G. Nishimura, J.M. Pauly, Rapid gridding reconstruction with a minimal oversampling ratio, IEEE Trans Med Imag, 24 (2005) 799-808. [85] R.D. Hoge, R.K.S. Kwan, G.B. Pike, Density compensation functions for spiral MRI, Magn Reson Med, 38 (1997) 117-128.

AL

[86] V. Rasche, R. Proksa, Reconstruction of MR images from data sampled along arbitrary k-space trajectories, In Proceedings of the 6th Annual Meeting of ISMRM, Sydney, Australia, 1998, p.668 [87] J.G. Pipe, P. Menon, Sampling density compensation in MRI: rationale and an iterative numerical solution, Magn Reson Med, 41 (1999) 179-186.

R N

[88] K.O. Johnson, J.G. Pipe, Convolution kernel design and efficient algorithm for sampling density correction, Magn Reson Med, 61 (2009) 439-447. [89] N.R. Zwart, K.O. Johnson, J.G. Pipe, Efficient sample density estimation by combining gridding and an optimized kernel, Magn Reson Med, 67 (2012) 701-710.

U

[90] R.N. Bracewell, A.C. Riddle, Inversion of fan-beam scans in radio astronomy, Astrophys J, 150 (1967) 427-434.

JO

[91] K.P. Pruessmann, M. Weiger, P. Bornert, P. Boesiger, Advances in sensitivity encoding with arbitrary k-space trajectories, Magn Reson Med, 46 (2001) 638-651. [92] P.C. Lauterbur, Image formation by induced local interactions: examples employing nuclear magnetic resonance, Nature, 242 (1973) 190-191. [93] M. Weiger, D.O. Brunner, B.E. Dietrich, C.F. Muller, K.P. Pruessmann, ZTE imaging in humans, Magn Reson Med, 70 (2013) 328-332. [94] D.O. Kuethe, A. Caprihan, I.J. Lowe, D.P. Madio, H.M. Gach, Transforming NMR data despite missing points, J Magn Reson, 139 (1999) 18-25. [95] M. Weiger, D.O. Brunner, M. Tabbert, M. Pavan, T. Schmid, K.P. Pruessmann, Exploring the bandwidth limits of ZTE imaging: Spatial response, out-of-band signals, and noise propagation, Magn Reson Med, 74 (2015) 1236-1247.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[96] J. Jackson, A. Macovski, D. Nishimura, Low-frequency restoration, Magn Reson Med, 11 (1989) 248-257. [97] K.P. Pruessmann, Encoding and reconstruction in parallel MRI, NMR Biomed, 19 (2006) 288299. [98] W. Dreher, I. Bardenhagen, L. Huang, M. Bäumer, On the suppression of background signals originating from NMR hardware components. Application to zero echo time imaging and relaxation time analysis, Magn Reson Imag, 34 (2016) 264-270.

O

F

[99] M. Weiger, M. Wu, M.C. Wurnig, D. Kenkel, W. Jungraithmayr, A. Boss, K.P. Pruessmann, Rapid and robust pulmonary proton ZTE imaging in the mouse, NMR Biomed, 27 (2014) 1129-1134.

O

[100] D. Idiyatullin, C. Corum, S. Moeller, M. Garwood, Gapped pulses for frequency-swept MRI, J Magn Reson, 193 (2008) 267-273.

PR

[101] D.M. Grodzki, P.M. Jakob, B. Heismann, Correcting slice selectivity in hard pulse sequences, J Magn Reson, 214 (2012) 61-67. [102] C. Li, J.F. Magland, A.C. Seifert, F.W. Wehrli, Correction of excitation profile in zero echo time (ZTE) imaging using quadratic phase-modulated RF pulse excitation and iterative reconstruction, IEEE Trans Med Imag, 33 (2014) 961-969.

E-

[103] K. Schieban, M. Weiger, F. Hennel, A. Boss, K.P. Pruessmann, ZTE imaging with enhanced flip angle using modulated excitation, Magn Reson Med, 74 (2015) 684-693.

PR

[104] D. Idiyatullin, C. Corum, J.Y. Park, M. Garwood, Fast and quiet MRI using a swept radiofrequency, J Magn Reson, 181 (2006) 342-349. [105] M. Garwood, MRI of fast-relaxing spins, J Magn Reson, 229 (2013) 49-54. [106] D. Idiyatullin, C.A. Corum, M. Garwood, Multi-band-SWIFT, J Magn Reson, 251 (2015) 1925.

AL

[107] M. Weiger, F. Hennel, K.P. Pruessmann, Sweep MRI with algebraic reconstruction, Magn Reson Med, 64 (2010) 1685-1695.

R N

[108] D.O. Kuethe, N.L. Adolphi, E. Fukushima, Short data-acquisition times improve projection images of lung tissue, Magn Reson Med, 57 (2007) 1058-1064. [109] R. Froidevaux, M. Weiger, D.O. Brunner, B.E. Dietrich, B.J. Wilm, K.P. Pruessmann, Filling the dead-time gap in zero echo time MRI: Principles compared, Magn Reson Med, 79 (2018) 20362045.

U

[110] D.M. Grodzki, P.M. Jakob, B. Heismann, Ultrashort echo time imaging using pointwise encoding time reduction with radial acquisition (PETRA), Magn Reson Med, 67 (2012) 510-518.

JO

[111] Y.T. Wu, G.P. Dai, J.L. Ackerman, M.I. Hrovat, M.J. Glimcher, B.D. Snyder, A. Nazarian, D.A. Chesler, Water- and fat-suppressed proton projection MRI (WASPI) of rat femur bone, Magn Reson Med, 57 (2007) 554-567. [112] R.N. Froidevaux, M.B. Rösler, D.O. Brunner, M. Weiger, K.P. Pruessmann, HYFI: Hybrid filling of the dead-time gap for faster zero echo time imaging, In Proceedings of the 27th Annual Meeting of ISMRM, Montreal, Canada, 2019, p.943 [113] T. Oberhammer, M. Weiger, F. Hennel, K.P. Pruessmann, Prospects of parallel ZTE imaging, In Proceedings of the 18th Annual Meeting of ISMRM, Montreal, Canada, 2011, p.2890 [114] R. Froidevaux, M. Weiger, K.P. Pruessmann, Enabling long excitation pulses in algebraic ZTE imaging by dead-time reduction via dual acquisition with alternative RF modulations, In Proceedings of the 26th Annual Scientific Meeting of ISMRM, Paris, France, 2018, p.8741

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[115] N. Kobayashi, U. Goerke, L. Wang, J. Ellermann, G.J. Metzger, M. Garwood, Gradientmodulated PETRA MRI, Tomography, 1 (2015) 85. [116] H. Jang, C.N. Wiens, A.B. McMillan, Ramped hybrid encoding for improved ultrashort echo time imaging, Magn Reson Med, 76 (2016) 814-825. [117] J. Zhang, D. Idiyatullin, C.A. Corum, N. Kobayashi, M. Garwood, Gradient-modulated SWIFT, Magn Reson Med, 75 (2016) 537-546.

F

[118] R.R. Ernst, W.A. Anderson, Application of Fourier Transform Spectroscopy to Magnetic Resonance, Rev Sci Instrum, 37 (1966) 93-102.

O

[119] S.B. Reeder, E.R. McVeigh, The effect of high performance gradients on fast gradient echo imaging, Magn Reson Med, 32 (1994) 612-621.

O

[120] J. Frahm, W. Hanicke, K.D. Merboldt, Transverse coherence in rapid FLASH NMR imaging, J Magn Reson, 72 (1987) 307-314.

PR

[121] N. Seiberlich, F.A. Breuer, M. Blaimer, K. Barkauskas, P.M. Jakob, M.A. Griswold, NonDCartesian data reconstruction using GRAPPA operator gridding (GROG), Magn Reson Med, 58 (2007) 1257-1265.

E-

[122] V. Rasche, R.W.D. Boer, D. Holz, R. Proksa, Continuous radial data acquisition for dynamic MRI, Magn Reson Med, 34 (1995) 754-761. [123] M. Lustig, D. Donoho, J.M. Pauly, Sparse MRI: The application of compressed sensing for rapid MR imaging, Magn Reson Med, 58 (2007) 1182-1195.

PR

[124] J. Du, M. Bydder, A.M. Takahashi, C.B. Chung, Two-dimensional ultrashort echo time imaging using a spiral trajectory, Magn Reson Imag, 26 (2008) 304-312. [125] P.T. Gurney, B.A. Hargreaves, D.G. Nishimura, Design and analysis of a practical 3D cones trajectory, Magn Reson Med, 55 (2006) 575-582.

AL

[126] F.E. Boada, J.S. Gillen, G.X. Shen, S.Y. Chang, K.R. Thulborn, Fast three dimensional sodium imaging, Magn Reson Med, 37 (1997) 706-715.

R N

[127] L. Wan, W. Zhao, Y. Ma, S. Jerban, A.C. Searleman, M. Carl, E.Y. Chang, G. Tang, J. Du, Fast quantitative 3D ultrashort echo time MRI of cortical bone using extended cones sampling, Magn Reson Med, 82 (2019) 225-236. [128] R.K. Robison, A.G. Anderson III, J.G. Pipe, Three-dimensional ultrashort echo-time imaging using a FLORET trajectory, Magn Reson Med, 78 (2017) 1038-1049.

U

[129] S. Konstandin, A.M. Nagel, P.M. Heiler, L.R. Schad, TwoDdimensional radial acquisition technique with density adaption in sodium MRI, Magn Reson Med, 65 (2011) 1090-1096.

JO

[130] P.E.Z. Larson, P.T. Gurney, D.G. Nishimura, Anisotropic field-of-views in radial imaging, IEEE Trans Med Imag, 27 (2008) 47-57. [131] M.S. Sussman, J.M. Pauly, G.A. Wright, Design of practical T2-selective RF excitation (TELEX) pulses, Magn Reson Med, 40 (1998) 890-899. [132] P.E.Z. Larson, P.T. Gurney, K. Nayak, G.E. Gold, J.M. Pauly, D.G. Nishimura, Designing longT2 suppression pulses for ultrashort echo time imaging, Magn Reson Med, 56 (2006) 94-103. [133] P.E.Z. Larson, S.M. Conolly, J.M. Pauly, D.G. Nishimura, Using adiabatic inversion pulses for long-T2 suppression in ultrashort echo time (UTE) imaging, Magn Reson Med, 58 (2007) 952-961. [134] J. Du, A.M. Takahashi, M. Bydder, C.B. Chung, G.M. Bydder, Ultrashort TE imaging with offresonance saturation contrast (UTE-OSC), Magn Reson Med, 62 (2009) 527-531.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[135] M. Carl, M. Bydder, J. Du, E. Han, Radiofrequency pulses for simultaneous short T2 excitation and long T2 suppression, Magn Reson Med, 65 (2011) 531-537. [136] J.P. Mugler, J.R. Brookeman, ThreeDdimensional magnetizationDprepared rapid gradientDecho imaging (3D MP RAGE), Magn Reson Med, 15 (1990) 152-157. [137] E.K. Brodsky, J.L. Klaers, A.A. Samsonov, R. Kijowski, W.F. Block, Rapid measurement and correction of phase errors from B0 eddy currents: Impact on image quality for non-cartesian imaging, Magn Reson Med, 69 (2013) 509-515.

O

F

[138] H. Lee, M. Weiger, K.P. Pruessmann, Exploring various radial trajectories for optimal relative phase correction in dual-echo subtraction ZTE MRI, In Proceedings of the 25th Annual Scientific Meeting of ISMRM, Honolulu, Hawaii, USA, 2017, p.1046

O

[139] H. Lee, M. Weiger, K.P. Pruessmann, High resolution dual-echo subtraction ZTE imaging at 7T with TE estimation based on multiple scaled trajectories, In Proceedings of the 25th Annual Scientific Meeting of ISMRM, Honolulu, Hawaii, USA, 2017, p.4023

PR

[140] C.T. Moonen, G. Liu, P.V. Gelderen, G. Sobering, A fast gradientDrecalled MRI technique with increased sensitivity to dynamic susceptibility effects, Magn Reson Med, 26 (1992) 184-189. [141] F. Wiesinger, A. Menini, A.B. Solana, Looping Star, Magn Reson Med, 81 (2019) 57-68.

E-

[142] R.F. Schulte, G. Buonincontri, M. Costagli, A. Menini, F. Wiesinger, A.B. Solana, Silent T2* and T2 encoding using ZTE combined with BURST, Magn Reson Med, 81 (2019) 2277-2287.

PR

[143] M. Carl, J.-T.A. Chiang, Investigations of the origin of phase differences seen with ultrashort TE imaging of short T2 meniscal tissue, Magn Reson Med, 67 (2012) 991-1003. [144] R.N. Froidevaux, M. Weiger, P.-J. Lu, K.P. Pruessmann, Signal-model-based water-fat separation in zero echo time (ZTE) MRI, In Proceedings of the 25th Annual Scientific Meeting of ISMRM, Honolulu, Hawaii, USA, 2017, p.772

AL

[145] L.J. Lehto, M. Garwood, O. Gröhn, C.A. Corum, Phase imaging in brain using SWIFT, J Magn Reson, 252 (2015) 20-28. [146] N.M. Szeverenyi, M. Carl, Contrast manipulation in MR imaging of short T2 and T2* tissues, in: eMagRes, John Wiley & Sons, Ltd, Chichester, United Kingdom, 2012, pp. 581-586.

R N

[147] J. Du, G. Ma, S. Li, M. Carl, N.M. Szeverenyi, S. VandenBerg, J. Corey-Bloom, G.M. Bydder, Ultrashort echo time (UTE) magnetic resonance imaging of the short T2 components in white matter of the brain using a clinical 3T scanner, NeuroImage, 87 (2014) 32-41.

U

[148] Q. Peng, Long-T2 suppression with low flip-angle b-SSFP ultra-short TE (LA-bUTE), In Proceedings of the 22nd Annual Meeting of ISMRM, Milan, Italy, 2014, p.819

JO

[149] L.N. Wang, W. Tang, Z.P. Zhen, H.M. Chen, J. Xie, Q. Zhao, Improving detection specificity of iron oxide nanoparticles (IONPs) using the SWIFT sequence with long T2 suppression, Magn Reson Imag, 32 (2014) 671-678. [150] G.E. Gold, J.M. Pauly, A. Macovski, R.J. Herfkens, MR spectroscopic imaging of collagen tendons and knee menisci, Magn Reson Med, 34 (1995) 647-654. [151] C. Li, J.F. Magland, X. Zhao, A.C. Seifert, F.W. Wehrli, Selective in vivo bone imaging with long-T2 suppressed PETRA MRI, Magn Reson Med, 77 (2017) 989-997. [152] J. Du, M. Bydder, A.M. Takahashi, M. Carl, C.B. Chung, G.M. Bydder, Short T2 contrast with three-dimensional ultrashort echo time imaging, Magn Reson Imag, 29 (2011) 470-482. [153] M. Weiger, M. Wu, M.C. Wurnig, D. Kenkel, A. Boss, G. Andreisek, K.P. Pruessmann, ZTE imaging with long-T2 suppression, NMR Biomed, 28 (2015) 247-254.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[154] M. Carl, G.M. Bydder, J. Du, UTE imaging with simultaneous water and fat signal suppression using a time-efficient multispoke inversion recovery pulse sequence, Magn Reson Med, 76 (2016) 577-582. [155] E.M. Johnson, U. Vyas, P. Ghanouni, K.B. Pauly, J.M. Pauly, Improved cortical bone specificity in UTE MR Imaging, Magn Reson Med, 77 (2017) 684-695.

F

[156] Y.-J. Ma, Y. Zhu, X. Lu, M. Carl, E.Y. Chang, J. Du, Short T2 imaging using a 3D double adiabatic inversion recovery prepared ultrashort echo time cones (3D DIR-UTE-Cones) sequence, Magn Reson Med, 79 (2018) 2555-2563.

O

[157] L. Soustelle, J. Lamy, F. Rousseau, J.P. Armspach, P. Loureiro de Sousa, A diffusion-based method for long-T2 suppression in steady state sequences: Validation and application for 3D-UTE imaging, Magn Reson Med, 80 (2018) 548-559.

O

[158] C. Li, J.F. Magland, H.S. Rad, H.K. Song, F.W. Wehrli, Comparison of optimized soft-tissue suppression schemes for ultrashort echo time MRI, Magn Reson Med, 68 (2012) 680-689.

PR

[159] S. Kirsch, L.R. Schad, Single-slice mapping of ultrashort T2, J Magn Reson, 210 (2011) 133136.

E-

[160] M.C. Wurnig, M. Calcagni, D. Kenkel, M. Vich, M. Weiger, G. Andreisek, F.W. Wehrli, A. Boss, Characterization of trabecular bone density with ultra-short echo-time MRI at 1.5, 3.0 and 7.0 T – comparison with micro-computed tomography, NMR Biomed, 27 (2014) 1159-1166.

PR

[161] A.S. Chaudhari, B. Sveinsson, C.J. Moran, E.J. McWalter, E.M. Johnson, T. Zhang, G.E. Gold, B.A. Hargreaves, Imaging and T2 relaxometry of short-T2 connective tissues in the knee using ultrashort echo-time double-echo steady-state (UTEDESS), Magn Reson Med, 78 (2017) 2136-2148. [162] J. Guo, X. Cao, Z.I. Cleveland, J.C. Woods, Murine pulmonary imaging at 7T: T2* and T1 with anisotropic UTE, Magn Reson Med, 79 (2018) 2254-2264.

AL

[163] M. Marcon, D. Keller, M.C. Wurnig, C. Eberhardt, M. Weiger, D. Eberli, A. Boss, Separation of collagenDbound and porous bone water transverse relaxation in mice: proposal of a multiDstep approach, NMR Biomed, 29 (2016) 866-872.

R N

[164] M. Weiger, D.O. Brunner, M. Wyss, B.E. Dietrich, B.J. Wilm, K.P. Pruessmann, ZTE imaging with T1 contrast, In Proceedings of the 21st Annual Meeting of ISMRM, Milan, Italy, 2014, p.4262 [165] J. Zhang, R. Chamberlain, M. Etheridge, D. Idiyatullin, C. Corum, J. Bischof, M. Garwood, Quantifying iron-oxide nanoparticles at high concentration based on longitudinal relaxation using a three-dimensional SWIFT look-locker sequence, Magn Reson Med, 71 (2014) 1982-1988.

JO

U

[166] J. Chen, E.Y. Chang, M. Carl, Y. Ma, H. Shao, B. Chen, Z. Wu, J. Du, Measurement of bound and pore water T1 relaxation times in cortical bone using three-dimensional ultrashort echo time cones sequences, Magn Reson Med, 77 (2017) 2136-2145. [167] M.C. Wurnig, M. Weiger, M. Wu, D. Kenkel, W. Jungraithmayr, K.P. Pruessmann, A. Boss, In vivo magnetization transfer imaging of the lung using a zero echo time sequence at 4.7 Tesla in mice: Initial experience, Magn Reson Med, 76 (2016) 156-162. [168] M. Marcon, D. Keller, M.C. Wurnig, M. Weiger, D. Kenkel, C. Eberhardt, D. Eberli, A. Boss, Separation of collagen-bound and porous bone-water longitudinal relaxation in mice using a segmented inversion recovery zero-echo-time sequence, Magn Reson Med, 77 (2017) 1909-1915. [169] Y.J. Ma, W. Zhao, L. Wan, T. Guo, A. Searleman, H. Jang, E.Y. Chang, J. Du, Whole knee joint T1 values measured in vivo at 3T by combined 3D ultrashort echo time cones actual flip angle and variable flip angle methods, Magn Reson Med, 81 (2019) 1634-1644.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[170] J. Du, M. Carl, E. Diaz, A. Takahashi, E. Han, N.M. Szeverenyi, C.B. Chung, G.M. Bydder, Ultrashort TE T1rho (UTE T1rho) imaging of the Achilles tendon and meniscus, Magn Reson Med, 64 (2010) 834-842. [171] Y.J. Ma, M. Carl, A. Searleman, X. Lu, E.Y. Chang, J. Du, 3D adiabatic T1rho prepared ultrashort echo time cones sequence for whole knee imaging, Magn Reson Med, 80 (2018) 14291439.

F

[172] J. Du, G. Hamilton, A. Takahashi, M. Bydder, C.B. Chung, Ultrashort echo time spectroscopic imaging (UTESI) of cortical bone, Magn Reson Med, 58 (2007) 1001-1009.

O

[173] J. Du, A.M. Takahashi, C.B. Chung, Ultrashort TE spectroscopic imaging (UTESI): application to the imaging of short T2 relaxation tissues in the musculoskeletal system, J Magn Reson Imag, 29 (2009) 412-421.

O

[174] F. Springer, P. Martirosian, F. Schick, Magnetization transfer – ultrashort echo time (MT-UTE) imaging, in: eMagRes, John Wiley & Sons, Ltd, Chichester, United Kingdom, 2007, pp. 597-604.

PR

[175] F. Springer, P. Martirosian, J. Machann, N.F. Schwenzer, C.D. Claussen, F. Schick, Magnetization transfer contrast imaging in bovine and human cortical bone applying an ultrashort echo time sequence at 3 Tesla, Magn Reson Med, 61 (2009) 1040-1048.

E-

[176] Y.-J. Ma, E.Y. Chang, M. Carl, J. Du, Quantitative magnetization transfer ultrashort echo time imaging using a time-efficient 3D multispoke Cones sequence, Magn Reson Med, 79 (2017) 692-700.

PR

[177] Y.J. Ma, A. Tadros, J. Du, E.Y. Chang, Quantitative two-dimensional ultrashort echo time magnetization transfer (2D UTE-MT) imaging of cortical bone, Magn Reson Med, 79 (2017) 19411949. [178] M. Marcon, M. Weiger, D. Keller, M.C. Wurnig, C. Eberhardt, D. Eberli, A. Boss, Magnetization transfer imaging of cortical bone in vivo using a zero echo time sequence in mice at 4.7 T: a feasibility study, Magn Reson Mater Phy, 29 (2016) 853-862.

AL

[179] I. Shenberg, A. Macovski, Inhomogeneity and multiple dimension considerations in magnetic resonance imaging with time-varying gradients, IEEE Trans Med Imag, 4 (1985) 165-174.

R N

[180] A. Maeda, K. Sano, T. Yokoyama, Reconstruction by weighted correlation for MRI with timevarying gradients, IEEE Trans Med Imag, 7 (1988) 26-31. [181] D.C. Noll, C.H. Meyer, J.M. Pauly, D.G. Nishimura, A. Macovski, A homogeneity correction method for magnetic resonance imaging with time-varying gradients, IEEE Trans Med Imag, 10 (1991) 629-637.

U

[182] L.C. Man, J.M. Pauly, A. Macovski, Multifrequency interpolation for fast offDresonance correction, Magn Reson Med, 37 (1997) 785-792.

JO

[183] B.P. Sutton, D.C. Noll, J.A. Fessler, Fast, iterative image reconstruction for MRI in the presence of field inhomogeneities, IEEE Trans Med Imag, 22 (2003) 178-188. [184] C. Barmet, J. Tsao, K.P. Pruessmann, Sensitivity encoding and B0 inhomogeneity - A simultaneous reconstruction approach, In Proceedings of the 13th Annual Meeting of ISMRM, Miami Beach, Florida, USA, 2005, p.682 [185] E. de Kerviler, A. Leroy-Willig, O. Clément, J. Frija, Fat suppression techniques in MRI: an update, Biomed Pharmacother, 52 (1998) 69-75. [186] E.M. Delfaut, J. Beltran, G. Johnson, J. Rousseau, X. Marchandise, A. Cotten, Fat suppression in MR imaging: techniques and pitfalls, RadioGraphics, 19 (1999) 373-382. [187] J. Ma, Dixon techniques for water and fat imaging, J Magn Reson Imag, 28 (2008) 543-558.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[188] M. Engström, C. Cozzini, M. Carl, G.C. McKinnon, F. Wiesinger, Perfect in-phase zero TE for musculoskeletal imaging, In Proceedings of the 27th Annual Meeting of ISMRM, Montreal, Canada, 2019, p.1287 [189] H. Eggers, K. Nehrke, P. Börnert, J. van den Brink, Quiet dixon imaging with looping star sequence, In Proceedings of the 27th Annual Meeting of ISMRM, Montreal, Canada, 2019, p.4638 [190] H. Jang, M. Carl, Y. Ma, S. Jerban, T. Guo, W. Zhao, E.Y. Chang, J. Du, Fat suppression for ultrashort echo time imaging using a single-point Dixon method, NMR Biomed, 32 (2019) e4069.

O

F

[191] R.R. Ernst, G. Bodenhausen, A. Wokaun, Principles of nuclear magnetic resonance in one and two dimensions, in, Clarendon Press Oxford, 1987, pp. 106-109.

O

[192] D. Mentrup, H. Eggers, Signal decay correction in 2D ultra-short echo time imaging, Magn Reson Mater Phy, 19 (2006) 62-70. [193] A. Nowogrodzki, The strongest scanners, Nature, 563 (2018) 24-26.

PR

[194] W.A. Edelstein, G.H. Glover, C.J. Hardy, R.W. Redington, The intrinsic signal-to-noise ratio in NMR imaging, Magn Reson Med, 3 (1986) 604-618. [195] P.A. Bottomley, R.W. Redington, W.A. Edelstein, J.F. Schenck, Estimating radiofrequency power deposition in body NMR imaging, Magn Reson Med, 2 (1985) 336-349.

E-

[196] P. Röschmann, Radiofrequency penetration and absorption in the human body: Limitations to high-field whole-body nuclear magnetic resonance imaging, Medical physics, 14 (1987) 922-931.

PR

[197] R. Kimmlingen, Magnetic field gradients, in: Magnetic resonance technology: hardware and system component design, The Royal Society of Chemistry, 2016, pp. 208-263. [198] P.R. Harvey, E. Katznelson, Modular gradient coil: A new concept in highDperformance wholeDbody gradient coil design, Magn Reson Med, 42 (1999) 561-570. [199] W. Irnich, F. Schmitt, Magnetostimulation in MRI, Magn Reson Med, 33 (1995) 619-623.

R N

AL

[200] M. Weiger, J. Overweg, M.B. Rösler, R. Froidevaux, F. Hennel, B.J. Wilm, A. Penn, U. Sturzenegger, W. Schuth, M. Mathlener, M. Borgo, P. Börnert, C. Leussler, R. Luechinger, B.E. Dietrich, J. Reber, D.O. Brunner, T. Schmid, L. Vionnet, K.P. Pruessmann, A high-performance gradient insert for rapid and short-T2 imaging at full duty cycle, Magn Reson Med, 79 (2018) 32563266.

U

[201] R. Froidevaux, M. Weiger, M.B. Rösler, D.O. Brunner, B. Wilm, B.E. Dietrich, J. Reber, K.P. Pruessmann, Pushing the limits of short-T2 MRI using 200 mT/m gradient strength and 2 MHz bandwidth, In Proceedings of the 26th Annual Scientific Meeting of ISMRM, Paris, France, 2018, p.3947

JO

[202] P.B. Roemer, W.A. Edelstein, C.E. Hayes, S.P. Souza, O.M. Mueller, The NMR phased array, Magn Reson Med, 16 (1990) 192-225. [203] T. Eichhorn, U. Ludwig, E. Fischer, J. Grobner, M. Gopper, A.K. Eisenbeiss, T. Flugge, J. Hennig, D. von Elverfeldt, J.B. Hovener, Modular coils with low hydrogen content especially for MRI of dry solids, PloS one, 10 (2015). [204] M.B. Rösler, M. Weiger, T. Schmid, D.O. Brunner, R. Froidevaux, K.P. Pruessmann, Ultrasonic soldering on glass for the construction of MRI coils with minimized background signal in short-T2 images, In Proceedings of the 33rd Annual Scientific Meeting of ESMRMB, Vienna, Austria, 2016, p.87 [205] M. Weiger, D.O. Brunner, T. Schmid, R. Froidevaux, M.B. Rösler, S. Gross, K.P. Pruessmann, A virtually 1H-free birdcage coil for zero echo time MRI without background signal, Magn Reson Med, 78 (2017) 399-407.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[206] M. Alecci, J. Brivati, G. Placidi, L. Testa, D. Lurie, A. Sotgiu, A submicrosecond resonator and receiver system for pulsed magnetic resonance with large samples, J Magn Reson, 132 (1998) 162166. [207] J. Brittain, A. Shankaranarayanan, V. Ramanan, A. Shimakawa, C.H. Cunningham, S. Hinks, R. Francis, R. Turner, J.W. Johnson, K.S. Nayak, S. Tan, J.M. Pauly, G.M. Bydder, Ultra-short TE imaging with single-digit (8μs) TE, In Proceedings of the 11th Annual Meeting of ISMRM, Kyoto, Japan, 2004, p.629

F

[208] Agilent Technologies, Understanding RF/microwave solid state switches and their applications, Application Note, 5989-7618EN (2010).

O

[209] M. Twieg, M.A.d. Rooij, M.A. Griswold, Active Detuning of MRI Receive Coils with GaN FETs, IEEE T Microw Theory, 63 (2015) 4169-4177.

O

[210] D. Spence, M. Aimi, Custom MEMS switch for MR surface coil decoupling, In Proceedings of the 23rd Annual Scientific Meeting of ISMRM, Toronto, Canada, 2015, p.704

PR

[211] D.O. Brunner, L. Furrer, M. Weiger, W. Baumberger, T. Schmid, J. Reber, B.E. Dietrich, B.J. Wilm, R. Froidevaux, K.P. Pruessmann, Symmetrically biased T/R switches for NMR and MRI with microsecond dead time, J Magn Reson, 263 (2016) 147-155.

E-

[212] G.R. Davies, D.J. Lurie, J.M.S. Hutchison, S.J. McCallum, I. Nicholson, Continuous-wave magnetic resonance imaging of short T2 materials, J Magn Reson, 148 (2001) 289-297.

PR

[213] D. Idiyatullin, S. Suddarth, C.A. Corum, G. Adriany, M. Garwood, Continuous SWIFT, J Magn Reson, 220 (2012) 26-31. [214] A.C. Özen, M. Bock, E. Atalar, Active decoupling of RF coils using a transmit array system, Magn Reson Mater Phy, 28 (2015) 565-576. [215] S.M. Sohn, J.T. Vaughan, R.L. Lagore, M. Garwood, D. Idiyatullin, In vivo MR imaging with simultaneous RF transmission and reception, Magn Reson Med, 76 (2016) 1932-1938.

AL

[216] A.C. Ozen, E. Atalar, J.G. Korvink, M. Bock, In vivo MRI with concurrent excitation and acquisition using automated active analog cancellation, Scientific reports, 8 (2018) 10631.

R N

[217] D.O. Brunner, B.E. Dietrich, M. Pavan, K.P. Pruessmann, MRI with sideband excitation: Application to continuous SWIFT, In Proceedings of the 20th Annual Meeting of ISMRM, Melbourne, Australia, 2012, p.150

U

[218] J. Marjanovic, M. Weiger, J. Reber, D.O. Brunner, B.E. Dietrich, B.J. Wilm, R. Froidevaux, K.P. Pruessmann, Multi-rate acquisition for dead time reduction in magnetic resonance receivers: application to imaging with zero echo time, IEEE transactions on medical imaging, 37 (2018) 408416.

JO

[219] D.O. Brunner, B.E. Dietrich, C. Barmet, S. Gross, B.J. Wilm, K.P. Pruessmann, Automatic reconstruction of gradient echo imaging sequences by concurrent and continuous monitoring of gradient and RF waveforms, In Proceedings of the 21st Annual Meeting of ISMRM, Salt Lake City, Utah, USA, 2013, p.553 [220] A.P. Aitken, D. Giese, C. Tsoumpas, P. Schleyer, S. Kozerke, C. Prieto, T. Schaeffter, Improved UTE-based attenuation correction for cranial PET-MR using dynamic magnetic field monitoring, Medical physics, 41 (2014) 012302. [221] S.J. Vannesjo, M. Haeberlin, L. Kasper, M. Pavan, B.J. Wilm, C. Barmet, K.P. Pruessmann, Gradient system characterization by impulse response measurements with a dynamic field camera, Magn Reson Med, 69 (2013) 583-593.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[222] D.O. Brunner, B.E. Dietrich, M. Çavuşoğlu, B.J. Wilm, T. Schmid, S. Gross, C. Barmet, K.P. Pruessmann, Concurrent recording of RF pulses and gradient fields – comprehensive field monitoring for MRI, NMR Biomed, 29 (2016) 1162-1172. [223] R.N. Froidevaux, M. Weiger, M.B. Rösler, B. Wilm, F. Hennel, R. Luechinger, B.E. Dietrich, J. Reber, K.P. Pruessmann, Ultra-high-bandwidth, high-resolution MRI of fast relaxing spins, In Proceedings of the 25th Annual Scientific Meeting of ISMRM, Honolulu, Hawaii, USA, 2017, p.4037

F

[224] W.C. Bae, P.C. Chen, C.B. Chung, K. Masuda, D. D'Lima, J. Du, Quantitative ultrashort echo time (UTE) MRI of human cortical bone: correlation with porosity and biomechanical properties, J Bone Miner Res, 27 (2012) 848-857.

O

O

[225] H.H. Cao, J. Ackerman, K. Ecklund, S. Steigman, M. Glimcher, Y.T. Wu, Water and fat suppressed proton projection MRI (WASPI)/micro CT measurements of bone mineralization, Bone, 43 (2008) S95.

PR

[226] H.H. Cao, J.L. Ackerman, M.I. Hrovat, L. Graham, M.J. Glimcher, Y.T. Wu, Quantitative bone matrix density measurement by water- and fat-suppressed proton projection MRI (WASPI) with polymer calibration phantoms, Magn Reson Med, 60 (2008) 1433-1443.

E-

[227] E. Diaz, C.B. Chung, W.C. Bae, S. Statum, R. Znamirowski, G.M. Bydder, J. Du, Ultrashort echo time spectroscopic imaging (UTESI): an efficient method for quantifying bound and free water, NMR Biomed, 25 (2012) 161-168.

PR

[228] M.A. FernándezDSeara, S.L. Wehrli, M. Takahashi, F.W. Wehrli, Water content measured by protonDdeuteron exchange NMR predicts bone mineral density and mechanical properties, J Bone Miner Res, 19 (2004) 289-296. [229] C. Li, A.C. Seifert, H.S. Rad, Y.A. Bhagat, C.S. Rajapakse, W. Sun, S.C.B. Lam, F.W. Wehrli, Cortical bone water concentration: dependence of MR imaging measures on age and pore volume fraction, Radiology, 272 (2014) 796-806.

AL

[230] M.K. Manhard, K.D. Harkins, D.F. Gochberg, J.S. Nyman, M.D. Does, 30-Second bound and pore water concentration mapping of cortical bone using 2D UTE with optimized half-pulses, Magn Reson Med, 77 (2017) 945-950.

R N

[231] H.S. Rad, S.C.B. Lam, J.F. Magland, H. Ong, C. Li, H.K. Song, J. Love, F.W. Wehrli, Quantifying cortical bone water in vivo by three-dimensional ultra-short echo-time MRI, NMR Biomed, 24 (2011) 855-864.

U

[232] I.L. Reichert, M.D. Robson, P.D. Gatehouse, T. He, K.E. Chappell, J. Holmes, S. Girgis, G.M. Bydder, Magnetic resonance imaging of cortical bone with ultrashort TE pulse sequences, Magn Reson Imag, 23 (2005) 611-618.

JO

[233] A. Techawiboonwong, H.K. Song, M.B. Leonard, F.W. Wehrli, Cortical bone water: in vivo quantification with ultrashort echo-time MR imaging, Radiology, 248 (2008) 824-833. [234] F.W. Wehrli, Magnetic resonance of calcified tissues, J Magn Reson, 229 (2013) 35-48. [235] F.W. Wehrli, H.K. Song, P.K. Saha, A.C. Wright, Quantitative MRI for the assessment of bone structure and function, NMR Biomed, 19 (2006) 731-764. [236] M. Weiger, M. Stampanoni, K.P. Pruessmann, Direct depiction of bone microstructure using MRI with zero echo time, Bone, 54 (2013) 44-47. [237] Y. Wu, J.L. Ackerman, H. Cao, T.G. Reese, M.I. Hrovat, M.J. Glimcher, In vivo bone matrix density measurements by water and fat suppressed projection MRI (WASPI), J Bone Miner Res, 23 (2008) S154-S154.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[238] Y.T. Wu, J.L. Ackerman, D.A. Chesler, L. Graham, Y. Wang, M.J. Glimcher, Density of organic matrix of native mineralized bone measured by water- and fat-suppressed proton projection MRI, Magn Reson Med, 50 (2003) 59-68. [239] Y.T. Wu, M.I. Hrovat, J.L. Ackerman, T.G. Reese, H.H. Cao, K. Ecklund, M.J. Glimcher, Bone matrix imaged in vivo by water- and fat-suppressed proton projection MRI (WASPI) of animal and human subjects, J Magn Reson Imag, 31 (2010) 954-963.

F

[240] E.C. Argentieri, M.F. Koff, R.E. Breighner, Y. Endo, P.H. Shah, D.B. Sneag, Diagnostic accuracy of zero-echo time MRI for the evaluation of cervical neural foraminal stenosis, Spine, 43 (2018) 928-933.

O

O

[241] R.E. Breighner, Y. Endo, G.P. Konin, L.V. Gulotta, M.F. Koff, H.G. Potter, Zero echo time imaging of the shoulder: Enhanced osseous detail by using MR imaging, Radiology, 286 (2018) 960966.

PR

[242] A. Nazaran, M. Carl, Y. Ma, S. Jerban, Y. Zhu, X. Lu, J. Du, E.Y. Chang, Three-dimensional adiabatic inversion recovery prepared ultrashort echo time cones (3D IR-UTE-Cones) imaging of cortical bone in the hip, Magn Reson Imag, 44 (2017) 60-64. [243] A.C. Seifert, C. Li, S.L. Wehrli, F.W. Wehrli, A surrogate measure of cortical bone matrix density by long T2-suppressed MRI, J Bone Miner Res, 30 (2015) 2229-2238.

E-

[244] A.C. Seifert, F.W. Wehrli, Solid-state quantitative 1H and 31P MRI of cortical bone in humans, Current osteoporosis reports, 14 (2016) 77-86.

PR

[245] J. Chen, S.P. Grogan, H. Shao, D. D'Lima, G.M. Bydder, Z. Wu, J. Du, Evaluation of bound and pore water in cortical bone using ultrashort-TE MRI, NMR Biomed, 28 (2015) 1754-1762. [246] X. Lu, S. Jerban, L. Wan, Y. Ma, H. Jang, N. Le, W. Yang, E.Y. Chang, J. Du, Threedimensional ultrashort echo time imaging with tricomponent analysis for human cortical bone, Magn Reson Med, 82 (2019) 348-355.

AL

[247] Y. Wu, J.L. Ackerman, D.A. Chesler, J. Li, R.M. Neer, J. Wang, M.J. Glimcher, Evaluation of bone mineral density using three-dimensional solid state phosphorus-31 NMR projection imaging, Calcif Tissue Int, 62 (1998) 512-518.

R N

[248] D. Gillies, B. Newling, E. Randall, Phosporus-31 solid-state NMR in high-field gradients: prospects for imaging bone using the long echo-train summation technique (LETS), J Magn Reson, 151 (2001) 235-241.

U

[249] Y.T. Wu, T.G. Reese, H.H. Cao, M.I. Hrovat, S.P. Toddes, R.A. Lemdiasov, J.L. Ackerman, Bone mineral imaged in vivo by P-31 solid state MRI of human wrists, J Magn Reson Imag, 34 (2011) 623-633.

JO

[250] A.C. Seifert, A.C. Wright, S.L. Wehrli, H.H. Ong, C. Li, F.W. Wehrli, 31P NMR relaxation of cortical bone mineral at multiple magnetic field strengths and levels of demineralization, NMR Biomed, 26 (2013) 1158-1166. [251] M.D. Robson, Ultrashort echo time imaging of phosphorus in man, in: eMagRes, John Wiley & Sons, Ltd, Chichester, United Kingdom, 2012, pp. 757-761. [252] M.D. Robson, P.D. Gatehouse, G.M. Bydder, S. Neubauer, Human imaging of phosphorus in cortical and trabecular bone in vivo, Magn Reson Med, 51 (2004) 888-892. [253] Y. Wu, J.L. Ackerman, D.A. Chesler, J. Wang, M.J. Glimcher, In vivo solid state 31P MRI of human tibia at 1.5 T, In Proceedings of the 7th Annual Meeting of ISMRM, Philadelphia, USA, 1999, p.313

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[254] Y.T. Wu, D.A. Chesler, M.J. Glimcher, L. Garrido, J.X. Wang, H.J. Jiang, J.L. Ackerman, Multinuclear solid-state three-dimensional MRI of bone and synthetic calcium phosphates, Proc Natl Acad Sci, 96 (1999) 1574-1578. [255] X. Zhao, H.K. Song, A.C. Seifert, C. Li, F.W. Wehrli, Feasibility of assessing bone matrix and mineral properties in vivo by combined solid-state 1H and 31P MRI, PloS one, 12 (2017) e0173995. [256] X. Zhao, H.K. Song, F.W. Wehrli, In vivo bone 31P relaxation times and their implications on mineral quantification, Magn Reson Med, 80 (2018) 2514-2524.

O

F

[257] A.C. Seifert, C. Li, C.S. Rajapakse, M. Bashoor-Zadeh, Y.A. Bhagat, A.C. Wright, B.S. Zemel, A. Zavaliangos, F.W. Wehrli, Bone mineral 31P and matrix-bound water densities measured by solidstate 31P and 1H MRI, NMR Biomed, 27 (2014) 739-748.

O

[258] Y. Sun, M. Ventura, E. Oosterwijk, J.A. Jansen, X.F. Walboomers, A. Heerschap, Zero echo time magnetic resonance imaging of contrast-agent-enhanced calcium phosphate bone defect fillers, Tissue Eng Part C Methods, 19 (2013) 281-287.

PR

[259] G.E. Gold, T. Wren, K. Nayak, D. Nishimura, G. Beaupre, In vivo short echo time imaging of Achilles tendon, In Proceedings of the 9th Annual Meeting of ISMRM, Glasgow, Scotland, UK, 2001, p.244

E-

[260] M. Robson, M. Benjamin, P. Gishen, G. Bydder, Magnetic resonance imaging of the Achilles tendon using ultrashort TE (UTE) pulse sequences, Clinical radiology, 59 (2004) 727-735.

PR

[261] J. Du, A.J.-T. Chiang, C.B. Chung, S. Statum, R. Znamirowski, A. Takahashi, G.M. Bydder, Orientational analysis of the Achilles tendon and enthesis using an ultrashort echo time spectroscopic imaging sequence, Magn Reson Imag, 28 (2010) 178-184. [262] Y.H. Lee, J.S. Suh, D. Grodzki, Ultrashort echo (UTE) versus pointwise encoding time reduction with radial acquisition (PETRA) sequences at 3 Tesla for knee meniscus: A comparative study, Magn Reson Imag, 34 (2016) 75-80.

AL

[263] A. Williams, Y. Qian, C.R. Chu, UTE-T2 mapping of human articular cartilage in vivo: a repeatability assessment, Osteoarthr Cartil, 19 (2011) 84-88.

R N

[264] C. Pauli, W.C. Bae, M. Lee, M. Lotz, G.M. Bydder, D.L. D’Lima, C.B. Chung, J. Du, Ultrashort–echo time MR imaging of the patella with bicomponent analysis: correlation with histopathologic and polarized light microscopic findings, Radiology, 264 (2012) 484-493. [265] G. Andreisek, M. Weiger, T2* mapping of articular cartilage: current status of research and first clinical applications, Invest Radiol, 49 (2014) 57-62.

U

[266] I. Bezrukov, F. Mantlik, H. Schmidt, B. Schölkopf, B.J. Pichler, MR-based PET attenuation correction for PET/MR imaging, Semin Nucl Med, 43 (2013) 45-59.

JO

[267] V. Keereman, P. Mollet, Y. Berker, V. Schulz, S. Vandenberghe, Challenges and current methods for attenuation correction in PET/MR, Magn Reson Mater Phy, 26 (2013) 81-98. [268] G. Wagenknecht, H.-J. Kaiser, F.M. Mottaghy, H. Herzog, MRI for attenuation correction in PET: methods and challenges, Magn Reson Mater Phy, 26 (2013) 99-113. [269] A. Boss, M. Weiger, F. Wiesinger, Future image acquisition trends for PET/MRI, Semin Nucl Med, 45 (2015) 201-211. [270] F. Wiesinger, L.I. Sacolick, A. Menini, S.S. Kaushik, S. Ahn, P. Veit-Haibach, G. Delso, D.D. Shanbhag, Zero TEMR bone imaging in the head, Magn Reson Med, 75 (2016) 107-114. [271] H. Jang, F. Liu, T. Bradshaw, A.B. McMillan, Rapid dual-echo ramped hybrid encoding MRbased attenuation correction (dRHE-MRAC) for PET/MR, Magn Reson Med, 79 (2018) 2912-2922.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[272] F. Wiesinger, M. Bylund, J. Yang, S. Kaushik, D. Shanbhag, S. Ahn, J.H. Jonsson, J.A. Lundman, T. Hope, T. Nyholm, P. Larson, C. Cozzini, Zero TE-based pseudo-CT image conversion in the head and its application in PET/MR attenuation correction and MR-guided radiation therapy planning, Magn Reson Med, 80 (2018) 1440-1451. [273] H. Lee, X. Zhao, H.K. Song, R. Zhang, S.P. Bartlett, F.W. Wehrli, Rapid dual-RF, dual-echo, 3D ultrashort echo time craniofacial imaging: A feasibility study, Magn Reson Med, (2018).

F

[274] G. Bauman, K.M. Johnson, L.C. Bell, J.V. Velikina, A.A. Samsonov, S.K. Nagle, S.B. Fain, Three-dimensional pulmonary perfusion MRI with radial ultrashort echo time and spatial–temporal constrained reconstruction, Magn Reson Med, 73 (2015) 555-564.

O

[275] C. Bergin, G. Glover, J. Pauly, Lung parenchyma: magnetic susceptibility in MR imaging, Radiology, 180 (1991) 845-848.

PR

O

[276] A. Bianchi, F. Lux, O. Tillement, Y. Cremillieux, Contrast enhanced lung MRI in mice using ultra-short echo time radial imaging and intratracheally administrated Gd-DOTA-based nanoparticles, Magn Reson Med, 70 (2013) 1419-1426. [277] G. Dournes, D. Grodzki, J. Macey, P.-O. Girodet, M. Fayon, J.-F. Chateil, M. Montaudon, P. Berger, F. Laurent, Quiet submillimeter MR imaging of the lung is feasible with a PETRA sequence at 1.5 T, Radiology, 276 (2015) 258-265.

E-

[278] S.L. Gewalt, G.H. Glover, L.W. Hedlund, G.P. Cofer, J.R. Macfall, G.A. Johnson, MR microscopy of the rat lung using projection reconstruction, Magn Reson Med, 29 (1993) 99-106.

PR

[279] F. Gibiino, L. Sacolick, A. Menini, L. Landini, F. Wiesinger, Free-breathing, zero-TE MR lung imaging, Magn Reson Mater Phy, (2014) 1-9. [280] K.M. Johnson, S.B. Fain, M.L. Schiebler, S. Nagle, Optimized 3D ultrashort echo time pulmonary MRI, Magn Reson Med, 70 (2013) 1241-1250.

AL

[281] N. Kobayashi, D. Idiyatullin, C. Corum, J. Weber, M. Garwood, D. Sachdev, SWIFT MRI enhances detection of breast cancer metastasis to the lung, Magn Reson Med, 73 (2015) 1812-1819. [282] S. Köhler, M. Weiger, U. Molkentin, F. Hennel, Towards routine lung MRI in small animals, In Proceedings of the 17th Scientific Meeting of ISMRM, Honolulu, Hawai’i, USA, 2009, p.2006

R N

[283] S.J. Kruger, S.B. Fain, K.M. Johnson, R.V. Cadman, S.K. Nagle, OxygenDenhanced 3D radial ultrashort echo time magnetic resonance imaging in the healthy human lung, NMR Biomed, 27 (2014) 1535-1541.

U

[284] Y. Ohno, H. Koyama, T. Yoshikawa, K. Matsumoto, M. Takahashi, M. Van Cauteren, K. Sugimura, T2* measurements of 3-T MRI with ultrashort TEs: capabilities of pulmonary function assessment and clinical stage classification in smokers, Am J Roentgenol, 197 (2011) W279-W285.

JO

[285] L.E. Olsson, M. Lindahl, P.O. Onnervik, L.B. Johansson, M. Palmer, M.K. Reimer, L. Hultin, P.D. Hockings, Measurement of MR signal and T2*in lung to characterize a tight skin mouse model of emphysema using single-point imaging, J Magn Reson Imag, 25 (2007) 488-494. [286] M.D. Shattuck, S.L. Gewalt, G.H. Glover, L.W. Hedlund, G.A. Johnson, MR microimaging of the lung using volume projection encoding, Magn Reson Med, 38 (1997) 938-942. [287] K.W. Stock, Q. Chen, H. Hatabu, R.R. Edelman, Magnetic resonance T2 measurements of the normal human lung in vivo with ultra-short echo times, Magn Reson Imag, 17 (1999) 997-1000. [288] M. Takahashi, O. Togao, M. Obara, M. van Cauteren, Y. Ohno, S. Doi, M. Kuro-o, C. Malloy, C.C. Hsia, I. Dimitrov, Ultra-short echo time (UTE) MR imaging of the lung: comparison between normal and emphysematous lungs in mutant mice, J Magn Reson Imag, 32 (2010) 326-333.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[289] M. Tibiletti, A. Bianchi, A. Kjorstad, S. Wundrak, D. Stiller, V. Rasche, Respiratory self-gated 3D UTE for lung imaging in small animal MRI, Magn Reson Med, 78 (2017) 739-745. [290] M. Tibiletti, J. Paul, A. Bianchi, S. Wundrak, W. Rottbauer, D. Stiller, V. Rasche, Multistage three-dimensional UTE lung imaging by image-based self-gating, Magn Reson Med, 75 (2016) 13241332.

F

[291] O. Togao, Y. Ohno, I. Dimitrov, C.C. Hsia, M. Takahashi, Ventilation/perfusion imaging of the lung using ultra-short echo time (UTE) MRI in an animal model of pulmonary embolism, J Magn Reson Imag, 34 (2011) 539-546.

O

[292] O. Togao, R. Tsuji, Y. Ohno, I. Dimitrov, M. Takahashi, Ultrashort echo time (UTE) MRI of the lung: assessment of tissue density in the lung parenchyma, Magn Reson Med, 64 (2010) 14911498.

O

[293] J. Yu, Y. Xue, H.K. Song, Comparison of lung T2* during freeDbreathing at 1.5 T and 3.0 T with ultrashort echo time imaging, Magn Reson Med, 66 (2011) 248-254.

PR

[294] D. Kenkel, Y. Yamada, M. Weiger, W. Jungraithmayr, M.C. Wurnig, A. Boss, Magnetization transfer as a potential tool for the early detection of acute graft rejection after lung transplantation in mice, J Magn Reson Imag, 44 (2016) 1091-1098.

E-

[295] J. Delacoste, J. Chaptinel, C. Beigelman-Aubry, D. Piccini, A. Sauty, M. Stuber, A double echo ultra short echo time (UTE) acquisition for respiratory motion-suppressed high resolution imaging of the lung, Magn Reson Med, 79 (2018) 2297-2305.

PR

[296] W. Jiang, F. Ong, K.M. Johnson, S.K. Nagle, T.A. Hope, M. Lustig, P.E.Z. Larson, Motion robust high resolution 3D free-breathing pulmonary MRI using dynamic 3D image self-navigator, Magn Reson Med, 79 (2018) 2954-2967. [297] D.O. Kuethe, J.M. Hix, L.E. Fredenburgh, T1 , T1 contrast, and Ernst-angle images of four ratlung pathologies, Magn Reson Med, 81 (2019) 2489-2500.

AL

[298] L. Mendes Pereira, T. Wech, A.M. Weng, C. Kestler, S. Veldhoen, T.A. Bley, H. Kostler, UTESENCEFUL: first results for 3D high-resolution lung ventilation imaging, Magn Reson Med, 81 (2019) 2464-2473.

R N

[299] A. Bianchi, M. Tibiletti, Å. Kjørstad, G. Birk, L.R. Schad, B. Stierstorfer, V. Rasche, D. Stiller, Three-dimensional accurate detection of lung emphysema in rats using ultra-short and zero echo time MRI, NMR Biomed, 28 (2015) 1471-1479.

U

[300] J. Delacoste, H. Feliciano, J. Yerly, V. Dunet, C. Beigelman-Aubry, G. Ginami, R.B. van Heeswijk, D. Piccini, M. Stuber, A. Sauty, A black-blood ultra-short echo time (UTE) sequence for 3D isotropic resolution imaging of the lungs, Magn Reson Med, 81 (2019) 3808-3818.

JO

[301] M.M. Willmering, R.K. Robison, H. Wang, J.G. Pipe, J.C. Woods, Implementation of the FLORET UTE sequence for lung imaging, Magn Reson Med, 0. [302] H.L. Ring, J. Zhang, N.D. Klein, L.E. Eberly, C.L. Haynes, M. Garwood, Establishing the overlap of IONP quantification with echo and echoless MR relaxation mapping, Magn Reson Med, 79 (2018) 1420-1428. [303] S. Magnitsky, J. Zhang, D. Idiyatullin, G. Mohan, M. Garwood, N.E. Lane, S. Majumdar, Positive contrast from cells labeled with iron oxide nanoparticles: Quantitation of imaging data, Magn Reson Med, 78 (2017) 1900-1910. [304] X. Lu, Y. Ma, E.Y. Chang, Q. He, A. Searleman, A. von Drygalski, J. Du, Simultaneous quantitative susceptibility mapping (QSM) and R2* for high iron concentration quantification with 3D ultrashort echo time sequences: An echo dependence study, Magn Reson Med, 79 (2018) 23152322.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[305] W. Dou, S. Mastrogiacomo, O. Koshkina, A. Veltien, M. Srinivas, X.F. Walboomers, A. Heerschap, 19F ZTE MR imaging of 19F labeled calcium phosphate cement at 11.7T, In Proceedings of the 24th Annual Meeting of ISMRM, Singapore, 2016, p.3877 [306] F. Schmid, C. Höltke, D. Parker, C. Faber, Boosting 19F MRI—SNR efficient detection of paramagnetic contrast agents using ultrafast sequences, Magn Reson Med, 69 (2013) 1056-1062.

F

[307] A.K. Bracher, C. Hofmann, A. Bornstedt, S. Boujraf, E. Hell, J. Ulrici, A. Spahr, B. Haller, V. Rasche, Feasibility of ultraDshort echo time (UTE) magnetic resonance imaging for identification of carious lesions, Magn Reson Med, 66 (2011) 538-545.

O

[308] D. Idiyatullin, C. Corum, S. Moeller, H.S. Prasad, M. Garwood, D.R. Nixdorf, Dental magnetic resonance imaging: making the invisible visible, J Endod, 37 (2011) 745-752.

O

[309] J.-B. Hövener, S. Zwick, J. Leupold, A.-K. Eisenbeiβ, C. Scheifele, F. Schellenberger, J. Hennig, D.v. Elverfeldt, U. Ludwig, Dental MRI: imaging of soft and solid components without ionizing radiation, J Magn Reson Imag, 36 (2012) 841-846.

PR

[310] M. Weiger, K.P. Pruessmann, A.K. Bracher, S. Kohler, V. Lehmann, U. Wolfram, F. Hennel, V. Rasche, High-resolution ZTE imaging of human teeth, NMR Biomed, 25 (2012) 1144-1151.

E-

[311] A. Bracher, C. Hofmann, A. Bornstedt, E. Hell, F. Janke, J. Ulrici, B. Haller, M. Geibel, V. Rasche, Ultrashort echo time (UTE) MRI for the assessment of caries lesions, Dentomaxillofac Radiol, 42 (2013) 20120321.

PR

[312] Y. Sun, O. Brauckmann, D.R. Nixdorf, A. Kentgens, M. Garwood, D. Idiyatullin, A. Heerschap, Imaging human teeth by phosphorus magnetic resonance with nuclear Overhauser enhancement, 6 (2016) 30756. [313] Y. Kang, B. Wu, S. Fu, N. Hong, Can zero TE imaging be a viable alternative to micro CT in dentistry imaging : application in tooth implanting and extraction?, In Proceedings of the 24th Annual Meeting of ISMRM, Singapore, 2016, p.704

AL

[314] W. Dou, S. Mastrogiacomo, A. Veltien, H.S. Alghamdi, X.F. Walboomers, A. Heerschap, Visualization of calcium phosphate cement in teeth by zero echo time 1H MRI at high field, NMR in biomedicine, 31 (2018) e3859.

R N

[315] K.M. Rychert, G. Zhu, M.M. Kmiec, V.K. Nemani, B.B. Williams, A.B. Flood, H.M. Swartz, B. Gimi, Imaging tooth enamel using zero echo time (ZTE) magnetic resonance imaging, SPIE, 2015. [316] E. Alonso-Ortiz, I.R. Levesque, G.B. Pike, MRI-based myelin water imaging: a technical review, Magn Reson Med, 73 (2015) 70-81.

U

[317] C. Laule, I.M. Vavasour, S.H. Kolind, D.K.B. Li, T.L. Traboulsee, G.R.W. Moore, A.L. MacKay, Magnetic resonance imaging of myelin, Neurotherapeutics, 4 (2007) 460-484.

JO

[318] R.A. Horch, J.C. Gore, M.D. Does, Origins of the ultrashortDT2 1H NMR signals in myelinated nerve: a direct measure of myelin content?, Magn Reson Med, 66 (2011) 24-31. [319] A. Waldman, J. Rees, C. Brock, M. Robson, P. Gatehouse, G. Bydder, MRI of the brain with ultra-short echo-time pulse sequences, Neuroradiology, 45 (2003) 887-892. [320] V. Sheth, H. Shao, J. Chen, S. Vandenberg, J. Corey-Bloom, G.M. Bydder, J. Du, Magnetic resonance imaging of myelin using ultrashort echo time (UTE) pulse sequences: phantom, specimen, volunteer and multiple sclerosis patient studies, NeuroImage, 136 (2016) 37-44. [321] M.J. Wilhelm, H.H. Ong, S.L. Wehrli, P.H. Tsai, D.B. Hackney, F.W. Wehrli, Prospects for quantitative imaging of myelin with dual-echo short inversion time 3D UTE MRI, In Proceedings of the 19th Annual Meeting of ISMRM, Montréal, Québec, Canada, 2011, p.2460

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[322] A.C. Seifert, C. Li, M.J. Wilhelm, S.L. Wehrli, F.W. Wehrli, Towards quantification of myelin by solid-state MRI of the lipid matrix protons, NeuroImage, 163 (2017) 358-367. [323] J. Du, V. Sheth, Q. He, M. Carl, J. Chen, J. Corey-Bloom, G.M. Bydder, Measurement of T1 of the ultrashort T2* components in white matter of the brain at 3T, PloS one, 9 (2014) e103296. [324] T. Boucneau, P. Cao, S. Tang, M. Han, D. Xu, R.G. Henry, P.E.Z. Larson, In vivo characterization of brain ultrashort-T2 components, Magn Reson Med, 80 (2018) 726-735.

O

F

[325] M. Weiger, R.N. Froidevaux, D.O. Brunner, M.B. Rösler, K.P. Pruessmann, Advances in direct myelin imaging, In Proceedings of the 27th Annual Meeting of ISMRM, Montreal, Canada, 2019, p.1118

O

[326] G.E. Gold, J.M. Pauly, G.H. Glover, J.C. Moretto, A. Macovski, R.J. Herfkens, Characterization of atherosclerosis with a 1.5DT imaging system, J Magn Reson Imag, 3 (1993) 399407.

PR

[327] R.R. Edelman, O. Flanagan, D. Grodzki, S. Giri, N. Gupta, I. Koktzoglou, Projection MR imaging of peripheral arterial calcifications, Magn Reson Med, 73 (2015) 1939-1945.

E-

[328] S. Sharma, S. Boujraf, A. Bornstedt, V. Hombach, A. Ignatius, A. Oberhuber, V. Rasche, Quantification of calcifications in endarterectomy samples by means of high-resolution ultra-short echo time imaging, Invest Radiol, 45 (2010) 109-113.

PR

[329] C.F. Chan, N.G. Keenan, S. Nielles-Vallespin, P. Gatehouse, M.N. Sheppard, J.J. Boyle, D.J. Pennell, D.N. Firmin, Ultra-short echo time cardiovascular magnetic resonance of atherosclerotic carotid plaque, J Cardiov Magn Reson, 12 (2010) 17. [330] A.B. Solana, A. Menini, L.I. Sacolick, N. Hehn, F. Wiesinger, Quiet and distortion-free, whole brain BOLD fMRI using T2-prepared RUFIS, Magn Reson Med, 75 (2016) 1402-1412.

AL

[331] S. Mangia, R. Chamberlain, F. De Martino, S. Moeller, C. Corum, T. Kim, C. Kalavagunta, S. Michaeli, M. Garwood, Functional MRI with SWIFT, In Proceedings of the 20th Annual Meeting of ISMRM, Melbourne, Australia, 2012, p.326

R N

[332] L.J. Lehto, D. Idiyatullin, J. Zhang, L. Utecht, G. Adriany, M. Garwood, O. Gröhn, S. Michaeli, S. Mangia, MB-SWIFT functional MRI during deep brain stimulation in rats, NeuroImage, 159 (2017) 443-448. [333] U. Katscher, T. Voigt, C. Findeklee, P. Vernickel, K. Nehrke, O. Doessel, Determination of electric conductivity and local SAR via B1 mapping, IEEE Trans Med Imag, 28 (2009) 1365-1374.

U

[334] S.K. Lee, S. Bulumulla, F. Wiesinger, L. Sacolick, W. Sun, I. Hancu, Tissue electrical property mapping from zero echo-time magnetic resonance imaging, IEEE Trans Med Imag, 34 (2015) 541550.

JO

[335] G. Madelin, J.-S. Lee, R.R. Regatte, A. Jerschow, Sodium MRI: methods and applications, Prog Nucl Mag Res Sp, 79 (2014) 14-47. [336] Y. Blunck, B.A. Moffat, S.C. Kolbe, R.J. Ordidge, J.O. Cleary, L.A. Johnston, Zero-gradientexcitation ramped hybrid encoding (zGRF-RHE) sodium MRI, Magn Reson Med, 81 (2019) 11721180. [337] K. Borzutzki, G. Brunklaus, Chapter Three - Magnetic Resonance Imaging Studies of the Spatial Distribution of Charge Carriers, in: G.A. Webb (Ed.) Annual Reports on NMR Spectroscopy, Academic Press, 2017, pp. 115-141. [338] W.P. Węglarz, A. Krzyżak, G. Machowski, M. Stefaniuk, ZTE MRI in high magnetic field as a time effective 3D imaging technique for monitoring water ingress in porous rocks at sub-millimetre resolution, Magn Reson Imag, 47 (2018) 54-59.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[339] W.P. Węglarz, A. Krzyżak, M. Stefaniuk, ZTE imaging of tight sandstone rocks at 9.4T — Comparison with standard NMR analysis at 0.05T, Magn Reson Imag, 34 (2016) 492-495. [340] C.B. Kennedy, B.J. Balcom, I.V. Mastikhin, Three-dimensional magnetic resonance imaging of rigid polymeric materials using single-point ramped imaging with T1 enhancement (SPRITE), Can J Chem, 76 (1998) 1753-1765. [341] J.-P. Renou, P.S. Belton, G.A. Webb, Magnetic resonance in food science: an exciting future, Royal Society of Chemistry, 2011.

O

F

[342] H.T. Fabich, A.J. Sederman, D.J. Holland, Development of ultrafast UTE imaging for granular systems, J Magn Reson, 273 (2016) 113-123.

O

[343] D. Groß, K. Zick, G. Guthausen, Chapter four - Recent MRI and diffusion studies of food structures, in: G.A. Webb (Ed.) Annual Reports on NMR Spectroscopy, Academic Press, 2017, pp. 145-197.

PR

[344] A.C. Ozen, U. Ludwig, L.M. Ohrstrom, F.J. Ruhli, M. Bock, Comparison of ultrashort echo time sequences for MRI of an ancient mummified human hand, Magn Reson Med, 75 (2016) 701708.

E-

[345] C.N. Wiens, N.S. Artz, H. Jang, A.B. McMillan, K.M. Koch, S.B. Reeder, Fully phase-encoded MRI near metallic implants using ultrashort echo times and broadband excitation, Magn Reson Med, 79 (2018) 2156-2163.

PR

[346] K. Sekihara, Steady-state magnetizations in rapid NMR imaging using small flip angles and short repetition intervals, IEEE Trans Med Imag, 6 (1987) 157-164. [347] P. Fantazzini, R.J.S. Brown, G.C. Borgia, Bone tissue and porous media: common features and differences studied by NMR relaxation, Magn Reson Imag, 21 (2003) 227-234. [348] J.S. Nyman, Q. Ni, D.P. Nicolella, X. Wang, Measurements of mobile and bound water by nuclear magnetic resonance correlate with mechanical properties of bone, Bone, 42 (2008) 193-199.

AL

[349] A.C. Seifert, S.L. Wehrli, F.W. Wehrli, BiDcomponent T2* analysis of bound and pore bone water fractions fails at high field strengths, NMR Biomed, 28 (2015) 861-872.

R N

[350] J. Du, E. Diaz, M. Carl, W. Bae, C.B. Chung, G.M. Bydder, Ultrashort echo time imaging with bicomponent analysis, Magn Reson Med, 67 (2012) 645-649. [351] R. Biswas, W. Bae, E. Diaz, K. Masuda, C.B. Chung, G.M. Bydder, J. Du, Ultrashort echo time (UTE) imaging with bi-component analysis: bound and free water evaluation of bovine cortical bone subject to sequential drying, Bone, 50 (2012) 749-755.

JO

U

[352] N. Kokabi, W. Bae, E. Diaz, C.B. Chung, G.M. Bydder, J. Du, Ultrashort TE MR imaging of bovine cortical bone: The effect of water loss on the T1 and T2* relaxation times, Magn Reson Med, 66 (2011) 476-482. [353] R.A. Horch, J.S. Nyman, D.F. Gochberg, R.D. Dortch, M.D. Does, Characterization of 1H NMR signal in human cortical bone for magnetic resonance imaging, Magn Reson Med, 64 (2010) 680-687. [354] R.A. Horch, D.F. Gochberg, J.S. Nyman, M.D. Does, Clinically compatible MRI strategies for discriminating bound and pore water in cortical bone, Magn Reson Med, 68 (2012) 1774-1784. [355] M. Robson, P. Gatehouse, T. He, D. Firmin, G. Bydder, S. Neubauer, Human in vivo imaging of phosphorus in cortical bone using ultrashort TE pulse sequences, In Proceedings of the 11th Annual Meeting of ISMRM, Toronto, Ontario, Canada, 2003, p.791 [356] L. Li, R.A. Kruger, P-31 NMR imaging of solid bone with solid echoes combined with refocused gradients, Phys Med Biol, 35 (1990) 1153-1158.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[357] U. Grosse, F. Springer, T. Hein, G. Grozinger, C. Schabel, P. Martirosian, F. Schick, R. Syha, Influence of physical activity on T1 and T2* relaxation times of healthy Achilles tendons at 3T, J Magn Reson Imag, 41 (2015) 193-201. [358] V. Juras, S. Apprich, P. Szomolanyi, O. Bieri, X. Deligianni, S. Trattnig, Bi-exponential T2* analysis of healthy and diseased Achilles tendons: an in vivo preliminary magnetic resonance study and correlation with clinical score, Eur Radiol, 23 (2013) 2814-2822.

F

[359] G.H. Filho, J. Du, B.C. Pak, S. Statum, R. Znamorowski, P. Haghighi, G. Bydder, C.B. Chung, Quantitative characterization of the Achilles tendon in cadaveric specimens: T1 and T2* measurements using ultrashort-TE MRI at 3 T, Am J Roentgenol, 192 (2009) W117-W124.

O

[360] E.Y. Chang, W.C. Bae, S. Statum, J. Du, C.B. Chung, Effects of repetitive freeze–thawing cycles on T2 and T2* of the Achilles tendon, Eur J Radiol, 83 (2014) 349-353.

PR

O

[361] V. Juras, S. Zbyn, C. Pressl, L. Valkovic, P. Szomolanyi, I. Frollo, S. Trattnig, Regional variations of T2* in healthy and pathologic achilles tendon in vivo at 7 Tesla: preliminary results, Magn Reson Med, 68 (2012) 1607-1613. [362] D.A. Reiter, P.C. Lin, K.W. Fishbein, R.G. Spencer, Multicomponent T2 relaxation analysis in cartilage, Magn Reson Med, 61 (2009) 803-809.

E-

[363] P.J. Lattanzio, K.W. Marshall, A.Z. Damyanovich, H. Peemoeller, Macromolecule and water magnetization exchange modeling in articular cartilage, Magn Reson Med, 44 (2000) 840-851.

PR

[364] L. Schreiner, I. Cameron, N. Funduk, L. Miljković, M. Pintar, D. Kydon, Proton NMR spin grouping and exchange in dentin, Biophys J, 59 (1991) 629-639. [365] N. Funduk, D. Kydon, L. Schreiner, H. Peemoeller, L. Miljković, M. Pintar, Composition and relaxation of the proton magnetization of human enamel and its contribution to the tooth NMR image, Magn Reson Med, 1 (1984) 66-75.

AL

[366] K.S. Nayak, J.M. Pauly, G.E. Gold, D.G. Nishimura, Imaging ultra-short T2 species in the brain, In Proceedings of the 8th Annual Meeting of ISMRM, Denver, USA, 2000, p.509

R N

[367] E. Ercan, P. Boernert, A. Webb, I. Ronen, Whole-brain tissue-based assessment of the ultrashort T2 component using 3D UTE MRI relaxometry, In Proceedings of the 20th Annual Meeting of ISMRM, Melbourne, Australia, 2013, p.4279 [368] Q. He, Y. Ma, S. Fan, H. Shao, V. Sheth, G.M. Bydder, J. Du, Direct magnitude and phase imaging of myelin using ultrashort echo time (UTE) pulse sequences: A feasibility study, Magn Reson Imag, 39 (2017) 194-199.

JO

U

[369] S.J. Fan, Y. Ma, E.Y. Chang, G.M. Bydder, J. Du, Inversion recovery ultrashort echo time imaging of ultrashort T2 tissue components in ovine brain at 3 T: a sequential D2O exchange study, NMR in biomedicine, 30 (2017) e3767. [370] S.J. Fan, Y. Ma, Y. Zhu, A. Searleman, N.M. Szeverenyi, G.M. Bydder, J. Du, Yet more evidence that myelin protons can be directly imaged with UTE sequences on a clinical 3T scanner: Bicomponent T2* analysis of native and deuterated ovine brain specimens, Magn Reson Med, 80 (2018) 538-547. [371] A. Ramani, A. Aliev, G. Barker, P. Tofts, Another approach to protons with constricted mobility in white matter: pilot studies using wideline and high-resolution NMR spectroscopy, Magn Reson Imag, 21 (2003) 1039-1043. [372] M.A. Schmidt, G.Z. Yang, P.D. Gatehouse, D.N. Firmin, FID-based lung MRI at 0.5 T: theoretical considerations and practical implications, Magn Reson Med, 39 (1998) 666-672. [373] M. Kveder, I. Zupančič, G. Lahajnar, R. Blinc, D. Šuput, D. Ailion, K. Ganesan, C. Goodrich, Water proton NMR relaxation mechanisms in lung tissue, Magn Reson Med, 7 (1988) 432-441.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

[374] S.M. Triphan, F.A. Breuer, D. Gensler, H.U. Kauczor, P.M. Jakob, Oxygen enhanced lung MRI by simultaneous measurement of T1 and T2* during free breathing using ultrashort TE, J Magn Reson Imag, 41 (2015) 1708-1714. [375] E.D. Pracht, J.F.T. Arnold, T.T. Wang, P.M. Jakob, Oxygen-enhanced proton imaging of the human lung using T2*, Magn Reson Med, 53 (2005) 1193-1196.

F

[376] H. Hatabu, D.C. Alsop, J. Listerud, M. Bonnet, W.B. Gefter, T2* and proton density measurement of normal human lung parenchyma using submillisecond echo time gradient echo magnetic resonance imaging, Eur J Radiol, 29 (1999) 245-252.

O

[377] R.J. Theilmann, T.J. Arai, A. Samiee, D.J. Dubowitz, S.R. Hopkins, R.B. Buxton, G.K. Prisk, Quantitative MRI measurement of lung density must account for the change in T2* with lung inflation, J Magn Reson Imag, 30 (2009) 527-534.

O

[378] N. Beckmann, B. Tigani, L. Mazzoni, J.R. Fozard, MRI of lung parenchyma in rats and mice using a gradient-echo sequence, NMR in biomedicine, 14 (2001) 297-306.

PR

[379] A. Zhu, D. Hernando, K.M. Johnson, S.B. Reeder, Quantification of short-T2* signal components in the liver using radial 3D UTE chemical shift-encoded MRI, In Proceedings of the 25th Annual Meeting of ISMRM, Honolulu, Hawaii, USA, 2017, p.120

E-

[380] E.K. Doyle, K. Toy, B. Valdez, J.M. Chia, T. Coates, J.C. Wood, Ultra-short echo time images quantify high liver iron, Magn Reson Med, 79 (2018) 1579-1585.

PR

[381] S. Codd, M. Mallett, M. Halse, J.H. Strange, W. Vennart, T. Van Doorn, A three-dimensional NMR imaging scheme utilizing doubly resonant gradient coils, J Magn Reson B, 113 (1996) 214-221.

JO

U

R N

AL

[382] M. Fernandez-Seara, S. Wehrli, F. Wehrli, Multipoint mapping for imaging of semi-solid materials, J Magn Reson, 160 (2003) 144-150.

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

Glossary 1D: one-dimensional 3D: three-dimensional AC: autocorrelation

F

AR: algebraic reconstruction

O

B0: static magnetic field

O

B1: radio-frequency magnetic field

PR

BLAST: back-projection low angle shot BOLD: blood-oxygenation-level dependent BW: bandwidth

E-

CT: computer tomography

PR

CTI: constant time imaging DCF: density correction function

EPT: electrical property tomography

AL

FID: free induction decay

fMRI: functional magnetic resonance imaging

R N

FOV: field-of-view

FT: Fourier transform

GM: gradient-modulated

U

GPU: graphical processing unit

JO

HYFI: hybrid filling IONP: iron oxide nanoparticle MEMS: micro-electro-mechanical system MP mapping: multipoint mapping MP RAGE: magnetization prepared rapid gradient echo imaging MR: magnetic resonance MRI: magnetic resonance imaging

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

MSK: musculoskeletal MTF: modulation transfer function NMR: nuclear magnetic resonance PBW: pixel bandwidth

F

PET: positron emission tomography

O

PETRA: pointwise encoding time reduction with radial acquisition

O

PIN: positive-instrinsic-negative

PR

PNS: peripheral nerve stimulation PSF: point spread function PTFE: polytetrafluoroethylene

E-

RF: radio-frequency

PR

RHE: ramped hybrid encoding RUFIS: rotating ultra-fast imaging sequence SAR: specific absorption rate

AL

SPI: single-point imaging

SPRITE: single-point ramped imaging with T1 enhancement

R N

SNR: signal-to-noise ratio

SWIFT: sweep imaging with Fourier transformation T1: longitudinal relaxation time

U

T2: transverse relaxation time

JO

T2*: apparent transverse relaxation time TE: echo time : acquisition range, or range of times across which data are collected

TR: repetition time T/R: transmit-receive UTE: ultra-short echo time WASPI: water- and fat-suppressed proton projection imaging

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

ZTE: zero echo time

Fig. 1: Dipolar coupling as an example of an orientation-dependent interaction. a) The dipolar field of one spin (red arrow) affects the total field at the location of another spin (blue arrow). The local field contribution %0E >?@ E  depends on the orientation of the molecule as the spins keep their

F

alignment with B0 upon molecular rotation, thus changing the position of the affected spin within the

O

dipolar field. The relevant angle θ is the angle between the B0 field vector and the axis connecting the two spins (dashed line). Molecular tumbling changes the orientation of the molecule over time and the

O

local field as seen by the spin fluctuates. b) Effects of dipolar coupling on NMR spectra. In a single crystal, all local lattices form the same angle with B0 (θ = 20° in this case), leading to a distinct peak.

PR

The second, symmetric peak arises from spin pairs in the same sample with anti-parallel instead of parallel mutual alignment (one of the spins in a) flipped). In a powder, all lattice orientations co-exist; thus the spectrum reflects all corresponding frequency shifts. The latter are weighted with @FG E

E-

according to their occurrence for an equal distribution of orientations, leading to the peculiar shape peaking at θ = 90°, the so-called Pake pattern. Again, both ways of alignment contribute. Note that

PR

the field shift vanishes at the so-called magic angle of θ ≈ 55°. In a hypothetical liquid with arbitrarily fast molecular rotation, the dipolar field is averaged to zero over time, leading to a single narrow line. In a real liquid, averaging also occurs, but in addition, the temporal fluctuations cause T1 and T2 relaxation, the latter leading to line broadening. Note that the spectra are scaled for presentation

AL

purposes. c) Free induction decay signals obtained for the spectra shown in b). The hypothetical liquid gives a constant on-resonant signal whereas the liquid with relaxation shows exponential decay. The

R N

off-resonant peaks of the single crystal lead to oscillations with constant envelope. The signal associated with the powder spectrum drops rapidly due to superposition of a wide range of frequencies. In this case, the observed loss of transverse magnetization is a non-random, reversible

U

process and therefore not considered to be actual relaxation. Fig. 2: Relaxation in a nutshell. In the presence of an orientation- or position-dependent interaction

JO

(e.g. dipolar coupling, see Fig. 1), molecular tumbling leads to randomly fluctuating fields ΔB(t) at the spins’ locations. The temporal behavior of the random field is determined by the speed of the tumbling motion and characterized by the correlation time τc, whereas its amplitude depends on the strength of the interaction. Two typical random fields are shown with rapid (τc = 0.05 ns, red) and

slow fluctuations (τc = 0.4 ns, blue). The autocorrelation (AC) function calculated from the average of such a random field decays exponentially with the constant τc. Consequently, the Fourier transform (FT) of the AC function is a Lorentz function, also known as the spectral density J(ω) of the random field, i.e., the energy available at different frequencies, which will be required below. Relaxation of longitudinal magnetization involves changes of the populations of different spin states. This is

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

illustrated in the energy diagram drawn for an isolated pair of equivalent spins, where transitions between combined states occur with probabilities W0, W1, and W2 where the subscripts indicate the corresponding “step size” in the diagram. The transition probabilities scale with the spectral density taken at multiples of the Larmor frequency ω0. How they contribute to T1 is given by the particular interaction causing the relaxation. T1 relaxation is most efficient for τc ≈ 1/ω0. Transverse relaxation

F

involves losses in coherences. Therefore, J(0) contributes and dominates T2 for large τc, leading to a

O

continuous T2 reduction for slower motion. The dashed vertical red and blue lines indicate the relaxation times obtained for the example values of τc.

O

Fig. 3: Principles of short-T2 MRI, illustrated by three different acquisitions performed on a phantom containing mineral oil with long T2 surrounded by pieces of rubber with short T2. All scans share the

PR

same nominal isotropic resolution, repetition time, and flip angle. (a) The gradient echo image is dominated by the oil signal while rubber is hardly visible. This is due to signal losses occurring before the echo is formed and acquired around the echo time TE. (b) In contrast, with zero echo time (ZTE)

E-

imaging, frequency encoding and acquisition start immediately after signal excitation, so that TE = 0. Hence, rubber appears in the image at full intensity. However, signal decay during the acquisition

PR

range Tk leads to image blurring, i.e., reduced actual spatial resolution for rubber. (c) Speeding up encoding by using a stronger gradient, thus creating a higher bandwidth BW, reduces the effects of signal decay and sharpens the edges of the rubber in the image. Hence, the two basic principles for

AL

encoding and acquisition in short-T2 MRI are 1) to start as early as possible and 2) to minimize the range Tk to avoid detrimental effects of signal decay. Fig. 4: High bandwidth as a requirement for short-T2 MRI. Imaging short-T2 samples at high spatial

R N

resolution requires both rapid encoding and large gradient (G) time-integrals. This can be achieved simultaneously only by using strong gradients, thus resulting in high resonance bandwidth. Fig. 5: Concepts and techniques for short-T2 MRI. The concepts follow the two principles of (1) an

U

early start and (2) a small range of data acquisition, as illustrated in Fig. 3 and indicated by the red

JO

labels. (1) is achieved by A) using 3D imaging without slice selection, B) employing radial center-out gradient encoding, and C) performing signal excitation in the presence of the encoding gradient. For concept B, using either pure Cartesian phase or pure radial frequency encoding patterns favors one or the other principle. Combinations of the different excitation and encoding schemes lead to the four basic techniques for short-T2 MRI, namely constant time, single point, ultra-short TE, and zero TE

imaging, along with their intermediate relatives. CTI = constant time imaging, SPI = single point imaging, SPRITE = single-point ramped imaging with T1 enhancement, MP mapping = multipoint mapping [382], UTE = ultra-short echo time, BLAST = back-projection low angle shot, RUFIS = rotating ultra-fast imaging sequence, ZTE = zero echo time, WASPI = water- and fat-suppressed

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

proton projection imaging, SWIFT = sweep imaging with Fourier transformation, PETRA = pointwise encoding time reduction with radial acquisition, HYFI = hybrid filling, GM = gradient-modulated, RHE = ramped hybrid encoding. Fig. 6: Sequence diagrams of the four basic techniques for short-T2 MRI. They all employ 3D centerout k-space encoding from which one TR interval is depicted. (a) In constant time imaging (CTI), pure

F

phase encoding is performed using a gradient of constant length and variable amplitude. One data

O

point is acquired per excitation after identical TE, including k = 0. Thus, spatial resolution is not impaired by T2 decay but signal is reduced related to TE. (b) Shorter TE is achieved with single point

O

imaging (SPI) where gradient ramps are excluded from encoding by performing both excitation and acquisition of high bandwidth while the gradient is at full strength. As indicated in the scheme by the

PR

grey bars, gradients need not be switched off between successive TR intervals, thus enabling nearsilent operation. (c) In contrast to the previous techniques, ultra-short echo time (UTE) imaging employs pure frequency encoding. Thus higher SNR efficiency is achieved but T2 decay may affect

E-

spatial resolution. After non-selective excitation, data acquisition is started concurrently with gradient ramping, the minimum TE being determined by the RF dead time. (d) In zero echo time (ZTE)

PR

imaging, high-bandwidth excitation is performed after ramping up the gradient. TE is actually zero because k = 0 coincides with the center of the RF pulse. Data acquisition starts only after the RF dead time, thus leading to a gap in central k-space. Equivalently to SPI, ZTE imaging can be performed

AL

with strongly reduced acoustic noise.

Fig. 7: Patterns in k-space of the four basic techniques for short-T2 MRI, schematically showing a central plane of 3D k-space for each case. The distance between neighboring k-space points according

R N

to the Nyquist criterion is given as the dwell time DW = 1/BW for the bandwidth BW in the FOV. (a) Both CTI and SPI are usually performed with Cartesian sampling. One k-space point is obtained per excitation after radial encoding of varying speed and direction. Thus, all points are equally T2-

U

weighted but with a relatively long TE. (b) In UTE imaging, data are acquired along the radial

JO

trajectory. In central k-space (yellow), encoding speed increases while ramping up the gradient, leading to a non-uniform radial sampling density. (c) In ZTE imaging, full, constant encoding speed is achieved but data are missing in central k-space (red) due to finite RF dead time. Additional data points are acquired by radial oversampling. Fig. 8: Effect of gradient ramping on depiction behavior in UTE imaging. A 1D simulation is shown using different ramp durations associated with different gradient slew rates. All times and durations are normalized with respect to ' which is the acquisition range at constant gradient without ramping, i.e., as in the ZTE case shown in Fig. 4d. Using the parameters FOV 240 mm, spatial resolution 1 mm, and gradient strength 25 mT/m results in ' = 0.47 ms. A T2 value shorter than

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

' of 0.20 ms is assumed to emphasize the effects of signal decay. As a reference, the ideal case is shown where all data are obtained without any delay after signal generation, i.e., unaffected by T2. (a) Time of acquisition  of the k-space positions of one radial readout. During the ramp, as indicated by the dashed lines, k-space speed increases until the target gradient strength is achieved. Compared to ', the UTE acquisition range  is increased by half the ramping time. (b) Modulation transfer

F

function (MTF) describing the apodization in k-space due to exponential T2 decay according to the

O

delayed acquisition times shown in a). (c) Point-spread-function (PSF) obtained by symmetrically complementing, zero-filling, and Fourier-transforming the MTF in b). The displayed PSFs describe

O

the distribution of signal from a pixel center (position 0) to other reconstructed pixels in the image and are depicted for the first 12 neighboring pixels. In the ideal case, the PSF is 1 at the location of origin

PR

and 0 at all other pixel positions. The Sinc shape reflects finite sampling up to kmax. For UTE encoding, the PSFs show a broadening of the main lobe and damped Sinc side lobes. (d) Normalizing

slew rates further emphasize this effect.

E-

the UTE PSFs to their maximum value reveals the associated loss of spatial resolution. Low gradient

Fig. 9: Impact of gradient fidelity on UTE image quality. (a) Ideal gradient shape with ramp. (b)

PR

Simulated image of a sphere obtained from data acquired with the ideal shape in a). (c) Distortion of gradient shape obtained for an exponential first-order, so-called gradient eddy current with time constant 0.05 '. (d) Image obtained when reconstructing data acquired with the distorted gradient

AL

in c) using the nominal, ideal trajectory. (e) Phase evolution from a zeroth-order field generated by an exponential, so-called B0 eddy current with time constant 0.5 '. (f) Image obtained from data with B0 eddy current in e) induced by the gradient coil encoding the horizontal direction. Correcting

R N

data for the B0-eddy-current-induced phase and reconstructing based on the k-space trajectory distorted by the gradient eddy current results in an image virtually identical to b). Fig. 10: Contributions to the initial RF dead time in pulsed MR. For symmetric RF pulses of constant

U

frequency, the focal point is located in the pulse center, so that only half of its duration contributes to

JO

the dead time. During T/R switching, the MR signal is corrupted by spikes and transients. After these RF events, fully filtered and uncorrupted data are only available after the lead-in time of the digital filter, corresponding to half of its length in the case of a symmetric filter with finite impulse response. The dwell time DW corresponds to the bandwidth in the target FOV according to Nyquist (see Fig. 7), whereas DWFOV is smaller and associated with radial oversampling. Fig. 11: Illustration of 1D algebraic reconstruction (AR) of a ZTE image. (a) Simulated ZTE k-space data obtained for a box-shaped object, exhibiting a gap of 3 Nyquist dwell times on the radius and two-fold oversampling. (b) Image obtained by FT of the data after filling the gap with zeroes, leading to low-frequency artifacts. (c) Image obtained by AR which is free of artifacts. (d) Data in k-space

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

obtained as FT of the reconstructed ZTE image, demonstrating recovery of data in the gap along with elimination of the oversampling. Fig. 12: Noise behavior in ZTE imaging. (a) 1D image reconstructed from pure noise time-domain data. Bad conditioning at a gap size of 3 Nyquist dwell times DW = 1/BW leads to noise enhancement and correlation in the image domain, which manifests as a low-frequency image distortion. Another

F

realization of the noise input data would lead to a different shape in the image. (b) Using as input the

O

same noise plus the signal from Fig. 11a results in a superposition of the two components also in the image. (c) Reconstruction of a 3D image from noise projections with gap = 3 DW as a). Due to high

O

inherent oversampling of 3D radial data in central k-space, noise enhancement is averaged out. (d) Only at larger gap = 4, further enhanced and correlated noise patterns are observed also in the 3D

PR

image.

Fig. 13: Background signal in ZTE imaging. Images of a human head were acquired at 7 T using a

E-

T/R birdcage RF coil. All kinds of signals were detected, including signals from very short-T2 materials of the RF coil and cables, the head support, and the patient table. (a) With a small dead-time

PR

gap and a FOV adapted to the head size, ordinary aliasing of radial data is manifested at the image center and borders. (b) A FOV twice as large contains most of the background signal and aliasing is largely eliminated. However, scan time is four-fold, and bandwidth BW and the gap size given in Nyquist dwell times DW = 1/BW are doubled. (c) For a slightly larger gap, worse conditioning

AL

amplifies the aliased out-of-band signal. (d) Choosing a larger FOV does not help in this case as the increased gap of 3 DW also leads to amplification of residual aliased signal.

R N

Fig. 14: RF excitation for ZTE imaging for a given signal bandwidth BW in the FOV. (a) Time course of short and long block-shaped pulses and a modulated pulse. Durations are given in Nyquist dwell times DW = 1/BW. The frequency sweep of the modulated pulse approximately covers BW. (b) Spectral amplitudes exhibiting high uniformity for the short block (= hard) pulse and strong variations

U

for the longer block pulse. The modulated pulse achieves about three times the amplitude of the hard

JO

pulse for a comparable uniformity. (c) ZTE image of a sphere using the short block pulse, showing uniform excitation. (d) A longer block pulse of duration 2.5 DW results in strong variations in image intensity. (e) A pulse of equal duration but with amplitude and frequency modulation provides uniformity comparable to the short block pulse. Fig. 15: RF excitation using a long modulated pulse as done in the SWIFT technique. To avoid missing signals generated and encoded during earlier parts of the pulse, excitation and acquisition are performed quasi-simultaneously. This is achieved by periodically alternating between the two operation modes. Thus both excitation and acquisition are run at a reduced duty cycle which is further

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

decreased by switching dead times. To also fully encode signals created at the end of the pulse, the encoding gradient and acquisition need to be continued after excitation [107]. Fig. 16: Patterns in k-space of two techniques filling the dead-time gap in ZTE imaging. (a) In PETRA, the gap area is covered by additional Cartesian SPI acquisitions. (b) In WASPI, additional radial projections are acquired at reduced gradient strength to limit the gap size to half a Nyquist

F

distance. To limit extra scan time, the angular density is chosen to meet the Nyquist criterion only at

O

the interface between the two parts of the data.

Fig. 17: Image fidelity of ZTE imaging and related techniques, illustrated in analogy to Fig. 8. (a) In

O

SPI, the time of acquisition is constant as in the ideal case but delayed according to the largest

PR

encoding. ZTE encoding starts at 0, increases linearly, and ends at the SPI delay. Note that with the present gap size of 50 μs = 13 DW, the ZTE case is not actually possible due to ill-conditioning and is only shown to illustrate the purely exponential decay. With PETRA, the additional SPI data in the gap

E-

are acquired after the dead time, leading to a reduced acquisition range Tk. In WASPI, a discontinuity is created by the different gradient strengths. (b) T2 decay illustrated by the MTF shows no effect for

PR

SPI, exponential weighting for ZTE acquisition, a flattened initial part for PETRA, and a discontinuous behavior for WASPI. (c) The PSF shows that at the original pixel location the strongest signal reduction is experienced by SPI, whereas ZTE, WASPI, and PETRA exhibit a comparable decrease. (d) The normalized PSF shows ideal resolution for SPI and an increase of the main lobe

AL

width for the three ZTE relatives. For both ZTE and PETRA, a resolution loss is observed from the distribution of signals to neighboring pixels. With WASPI, this effect is stronger and reaches further

R N

out, thus leading to side lobes.

Fig. 18: Generating contrast in ZTE imaging by magnetization preparation demonstrated on an excised lamb joint. (a) Non-prepared ZTE image showing primarily proton density and some T1 contrast. (b) In gradient echo imaging, signals from short-T2 tissues are reduced, in particular from

U

cortical bone. (c) In ZTE imaging with fat suppression, both short- and long-T2 water signals remain.

JO

(d) Selectively suppressing long-T2 water in the ZTE sequence removes signals from cartilage and

muscle, and fat signal dominates the image. (e) Applying suppression of both fat and long-T2 water results in positive contrast for short-T2 tissues, showing primarily trabecular and cortical bone. Remaining bright fat signal (right) indicates imperfect suppression due to local off-resonance. (f) In the CT reference image, strong similarities with the MRI data in e) are observed. Reproduced with permission from Ref. [153]. Fig. 19: Off-resonance effects in radial imaging, demonstrated for the specific case of chemical shifts in the presence of both water and fat. (a) Simulated images with a pixel bandwidth (PBW) corresponding to the frequency difference between water and fat. Setting the demodulation frequency

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

f to the water resonance results in sharp edges of the water compartment and blurred depiction of the fat parts. At interfaces of the two compartments, destructive and constructive interference lead to edge enhancement (arrow). Demodulation at the fat frequency changes the order of low and high signal intensities at interfaces whereas at the intermediate frequency, edges are not enhanced but blurred. (b) In-vivo ZTE data of a human knee acquired at 7 T with PBW = 870 Hz, corresponding to 0.87 times

F

the water-to-fat-shift. At different demodulation frequencies, the equivalent effects to the simulations

O

in a) are observed, as indicated by arrows.

Fig. 20: Resolution enhancement by exponential T2 filtering in simulated radial MRI with center-out

O

encoding. (a) Two ellipse-shaped objects without T2 filtering. For a relatively long T2 equal to the acquisition range Tk, the object is sharply depicted. For a significantly shorter T2, blurring is observed.

PR

(b) Moderate T2 filtering enhances resolution of the short-T2 object but leads to edge artifacts for the object with longer T2. (c) With full compensation of T2 decay, the short-T2 object is sharply depicted.

E-

However, strong artifacts are observed for the long-T2 object and noise is enhanced.

PR

Fig. 21: Selected short-T2 MRI applications. All images were obtained using the 3D ZTE technique. (a) MSK: Human knee acquired at 7 T with isotropic spatial resolution of 0.75 mm in 4 m 17 s. Reproduced with permission from Ref. [205]. (b) Pulmonary MRI: Mouse lung acquired at 4.7 T with a resolution of 0.31 mm in 4 m 30 s using respiratory gating. Reproduced with permission from Ref.

AL

[99]. (c) Dental MRI: Excised human tooth acquired at 11.7 T with a resolution 0.13 × 0.13 × 0.18

R N

mm3 in 21 m 18 s. Reproduced with permission from Ref. [310].

U

Table 1: Principles for short-T2 MRI.  = time of acquisition, TE = echo time, SNR = signal-tonoise ratio,  = range of times of acquisition.

JO

Principle for data acquisition 1. Early start

2. Small range

Characterized by

Affects

!"  , TE

SNR



Spatial resolution

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

Table 2: Examples of short-T2 tissues and materials. Relaxation times were derived and selected from single- or multi-component analysis. An asterisk * indicates a T2* value as opposed to T2. Unless stated otherwise, the values associated with tissues refer to 1H signal in water molecules. B0 [T]

Origin

T2 or T2* [μs]

References

Cortical bone

0.47

Porcine, bovine, specimen

400

[347]

(bone water)

0.64

Human, specimen

400

[348]

1.5

Human

420 – 500 *

Human, specimen

401 *

Human

260 – 786 *

Human, specimen

237 – 576 *

[224, 233, 349]

Bovine, specimen

290 – 368 *

[175, 227, 350-352]

4.7

Human, specimen Human, specimen

7.0

O

O

290 *

[354]

431 *

[168]

290 *

[236]

Human, specimen

302 *

[349]

Human, specimen

540, 400 *

[243]

Porcine, specimen

597 *

[258]

Human

176 – 207 *

[252, 355]

Bovine, specimen

190 *

[356]

Human

160 *

[249]

Human

178 *

[256]

Lamb, specimen

189 *

[250]

Synthesized

448

[248]

7.0

Chicken, specimen

675

[248]

9.4

Lamb, specimen

119 *

[250]

Rabbit, specimen

90 *

[22]

0.47

Porcine, specimen

650 – 3800

[25]

1.5

Human

1200 – 2000 *

[150, 259]

3.0

Human

5200 – 6000

[161]

Human

600 – 2320 *

[31, 54, 173, 227, 260, 357, 358]

Human, specimen

1180 – 2180 *

Bovine, specimen

1280 *

Human

340

AL

1.5

PR

Bovine, specimen

3.0

R N U JO

[353, 354]

[349]

5.6

Tendons

400 – 416

368 *

11.7

(bone matrix, hydroxyapatite 31P)

[53, 166, 172, 231]

Human, specimen

9.4

Cortical bone

[349]

E-

Mouse

[232]

PR

3.0

F

Tissue/material

7.0

[359, 360] [350] [361]

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

Collagen protein

Dentine

[51]

3.0

Human

5800 – 9300

[161]

Human

3260 *

[173]

Caprine, specimen

1530 *

[350]

1.5

Human

3000 *

[150]

3.0

Human

9700 – 11400

[161]

Human

4190 *

Human, specimen

1800 – 2100 *

1.5

Human

5000 – 10000

3.0

Human

1340 *

Human, specimen

480 – 990 *

9.4

Bovine, specimen

0.47

O [51]

O

[173]

[264, 350]

Porcine, bovine, specimen

12 *

[347]

0.7

Bovine, specimen

20 *

[363]

1.4

Rat tendon, specimen

90

[32]

4.7

Human, specimen

57, 12 *

[353]

9.4

Rabbit, specimen

20 – 40 *

[228]

Human, specimen

38000

[364]

Human, specimen

12, 65, 200, 1000 *

[364]

Human

324

[307]

Human, specimen

12, 65, 230 *

[365]

0.9

AL

0.9 1.5

Human

150 – 350 *

[366]

3.0

Human

358 *

[320]

Human

420 *

[147]

Human

500 – 700 *

[324]

Human

260 *

[367]

Human

336 *

[368]

Human, MS specimen

216 *

[320]

Ovine, specimen

200 – 300 *

[369]

Ovine, specimen

150 – 400 *

[370]

Bovine, specimen

114 *

[320]

4.7

Rat, specimen

50 – 1000 *

[318]

7.0

Human

200 – 300 *

[324]

9.4

Human, specimen

38 – 57

[371]

Rat, bovine, specimen

8 – 1000 *

[35]

U JO

[227, 350]

[362]

R N

Myelin lipids

[173]

2300

3.0 Dental enamel

F

4000 – 10000

PR

Cartilage

Human

E-

Meniscus

1.5

PR

Ligaments

JOURNAL PRE-PROOF Weiger et al., Short-T2 MRI

Human

4000 *

[372]

0.68

Rat, specimen

2340 *

[373]

1.5

Human

1100 – 2110 *

[287, 293, 374-377]

2.0

Rat, specimen

830 *

[373]

3.0

Human

460 – 820 *

[284, 293]

Mouse

910 – 1000 *

[288, 292]

4.7

Mouse

460 – 500 *

7.0

Mouse

395 *

7.9

Mouse

770 *

Liver

3.0

Human

300 – 450 *

Rubber

2.5

Polyvinylchloride (PVC)

2.4

Acrylonitrile butadiene styrene (ABS)

2.5

R N U JO

O

[162] [167]

O

[379, 380] [381]

180 – 1350 *

[147, 172, 223]

38 *

[340]

98 *

[340]

E-

2.4

PR

Polyethylene (PE)

AL

2.4

[285, 378]

300 – 1300 *

3.0 Polymethylmethacrylate (PMMA, acrylic glass)

F

0.5

PR

Lung

442 *

[340]

760 *

[381]

JOURNAL PRE-PROOF

a

Molecular tumbling or crystal rotation Dipolar field

B0 θ

F

θ

B0 + ΔB(θ)

Spectrum

E-

Liquid

PR

b

O

O

Molecule

Single crystal

θ = 20°

PR

Liquid with relaxation

θ = 90°

Powder

AL

θ = 55°

f

R N

0

θ = 0°

Free induction decay

JO

U

c

0

Fig. 1

t

Fig. 2

ΔB

Amplitude

Interaction (e.g. dipolar coupling)

τc

t [s]

τc = 0.05 ns

τc = 0.4 ns

W0

W2

AL W1

W 1 ~ J(ω0)

W 2 ~ J(2ω0)

Transition probabilities

W 0 ~ J(0)

FT

[s]

Liquids

F

O

τc [s]

τc = 1/ω0

Relaxation

ω [rad/s]

T2

T1

Spectral density

0 ω0 2ω0

J

O

PR

E-

Weights

PR τ [s]

Autocorrelation function

Energy diagram spin pair

AC

R N

U

JO

Random field

Solids

JOURNAL PRE-PROOF

JO

U

Gradient echo

t

TE

AQ

Signal loss

Rubber T2 ≈ 0.4 ms

Oil T2 ≈ 30 ms

b ZTE

AQ

1 Early start

TE = 0

Tk

c ZTE

F

2 Small range

O

AQ

Tk

Tk = 0.32 ms, BW = 200 kHz

O

PR Signal decay

E-

PR

Tk = 1.28 ms, BW = 50 kHz

AL

R N

TE = 1.32 ms, BW = 100 kHz

Fig. 3

G

RF

Sig

a

JOURNAL PRE-PROOF

Fig. 4

+

High bandwidth

F

O

O

PR

Large G-integral

High resolution

E-

PR

Strong G

AL

Rapid encoding

Short T2

R N

U

JO JOURNAL PRE-PROOF

Fig. 5

1

Radial center-out encoding

B

1

Frequency encoding Radial

Intermediate

2

Phase encoding Cartesian

C

A

AL

Ultra-short TE UTE

Constant Time Imaging CTI

GM-SWIFT

MP mapping WASPI PETRA HYFI

Single Point Imaging SPI, SPRITE

RF during G

3D

2

1

1

Zero TE BLAST, RUFIS, ZTE SWIFT

F

O

O

PR

RHE, GM-PETRA

E-

Intermediate

No slice selection

PR

RF before G

R N

U

Concept

JO JOURNAL PRE-PROOF

TE

AQ

Tk

PR

TE = 0

F

AQ

Tk

TE

O

O

Dead time

d ZTE

b SPI

E-

PR

k=0

AQ

Tk = 0

AL

R N

TE

U

JO

Ramp

c UTE

Fig. 6

G

RF

G

RF

a CTI

AQ

Tk = 0

JOURNAL PRE-PROOF

gradient encoding

b

k-space points:

DW

Nyquist

k0

F

O

missed

Dead time gap

ZTE

oversampled

O

c

PR acquired

E-

PR

UTE

Gradient ramp

AL

R N

U

JO

CTI / SPI

Equal, long TE

Fig. 7

a

JOURNAL PRE-PROOF

Fig. 8

b

ramp

AL

R N

U

JO

a

Tk

d

F

O

O

PR

E-

PR

plateau

c

JOURNAL PRE-PROOF

Fig. 9

AL

R N

U

JO

F

O

O

PR

E-

PR

JOURNAL PRE-PROOF

Fig. 10

RF pulse

F

DW

O

O

PR

E-

PR

T/R switching Filter lead-in

AL

Dead time

R N

U

JO

DWov

JOURNAL PRE-PROOF

Fig. 11

d

a

Dead time gap

IFT

FT

c

b

AL

R N

U

JO

k-Space

F

O

O

PR

E-

PR

Image

JOURNAL PRE-PROOF

c

Fig. 12

3D

1D

a

gap = 3

gap = 3

Noise

d

gap = 4

gap = 3

AL

Signal + noise

R N

b

U

JO

F

O

O

PR

E-

PR

JOURNAL PRE-PROOF

Fig. 13

FOV 530 mm BW 500 kHz DW 2 μs

FOV 265 mm BW 250 kHz DW 4 μs

Gap 2

b

Gap 1

a

AL

R N

Dead time 4 μs

U

JO

F

O

O

PR Gap 3

E-

PR d

Gap 1.5

c

Dead time 6 μs

JOURNAL PRE-PROOF

Fig. 14

b

a

-1

-0.5

0

BW = FOV

0

0.5

f/BW

1 t/DW

d

c

2.5 DW

0.4 DW

e 2.5 DW

F

O

O

PR

E-

PR

BW

AL

Block short Block long Mod. amp. Mod. freq.

R N

U

JO JOURNAL PRE-PROOF

Freq

Fig. 15

G

AQ

RF

Amp

U

RF AQ Dead interval interval time

F

O

O

PR

E-

PR

Encoding time

AL

R N

Pulse duration

JO JOURNAL PRE-PROOF

Fig. 16

a

SPI

b

F

O

O

PR

E-

Lower gradient

PR

WASPI

AL

R N

U

JO

PETRA

JOURNAL PRE-PROOF

Fig. 17

b

gap

AL

R N

U

JO

a

d

F

O

O

PR

E-

PR

Tk

c

JOURNAL PRE-PROOF

fat

Gradient echo

d

c

Water sup

f

CT

F

O

O

PR

E-

e Fat + water sup

PR

cartilage

AL

Fat sup

R N

U

trab. cort. bone bone

ZTE

muscle

Fig. 18

b

a

JO JOURNAL PRE-PROOF

fat water

F

O

f = ffat

O

PR

E-

PR

f = (fwater+ffat)/2

AL

R N

U

JO

PBW = fwater - ffat

Fig. 19

b

a

f = fwater

JOURNAL PRE-PROOF

b

c

F

T2,filter = T2,short

O

O

PR

E-

PR

T2,filter = 2 T2,short

AL

R N

U

T2,short = 0.2 Tk

No T2 filtering

T2,long = Tk

Fig. 20

a

JO JOURNAL PRE-PROOF

Fig. 21

a

patellar tendon

b

AL

R N

U

JO

F

O

O

pulp

PR

E-

PR

caries lesion

c

bone

dentine

enamel

JOURNAL PRE-PROOF

JO

U

R N

AL

PR

E-

PR

O

O

F

JOURNAL PRE-PROOF

JOURNAL PRE-PROOF Short-T2 MRI: Principles and recent advances Highlights Overarching concepts structure the variety of short-T2 MRI techniques Relaxation in a nutshell

F

Detailed description of ZTE imaging and its variants

O

Dedicated hardware requirements for all scanner components

JO

U

R N

AL

PR

E-

PR

O

Comprehensive literature overview including T2 and T2* values of potential imaging targets