Ultra-High Contrast MRI of the Brain and Spinal Cord using Directly Acquired and Synthetic Bipolar Filter (BLAIR) Images

Daniel M. Cornfeld; Paul Condron; Mark Bydder; Tracy R Melzer; Emanuele Pravatà; John D Port; Samantha J. Holdsworth; Graeme Bydder

doi:10.20944/preprints202505.2380.v3

Submitted:

30 May 2025

Posted:

02 June 2025

Read the latest preprint version here

Abstract

Ultra-high contrast (UHC) MRI is a term used to describe MR imaging that shows abnormalities with high contrast when little or no abnormality is seen on conventional MR images. It is achieved without the use of increased static or gradient magnetic fields. UHC can be accomplished by using tissue properties such as T1 twice or more in the same sequence to synergistically contribute to overall contrast as with directly acquired divided subtracted inversion recovery (dSIR) T1-bipolar sequences. These sequences also provide high contrast and high spatial resolution signal boundaries. UHC MRI can also be achieved by applying synthetic bipolar filters (BLAIRs) to high quality tissue property maps. Illustrative cases of mild traumatic brain injury, multiple sclerosis, methamphetamine substance use disorder, Parkinson's disease and white matter changes associated with a cerebral tumour are shown. Patients showed widespread abnormalities with directly acquired and synthetic dSIR sequences in areas where little or no abnormality was seen with conventional T2-wSE and/or T2-FLAIR sequences. New signs seen with dSIR sequences include the whiteout sign and grayout signs as well as the bilaminar cortex sign and bubble signs. dSIR and other bipolar sequences create high contrast from small changes in T1 as well as other tissue properties, and are complementary to conventional clinical sequences which create high contrast from larger changes in T1 and/or T2. Comparison is made with other sequences which use single or combined inversion recovery images. dSIR and other tissue property bipolar filter (BLAIR) sequences utilise sequences that already exist on most MR machines and are easy to implement.

Keywords:

Ultra-high contrast

;

divided subtracted inversion recovery (dSIR) sequence

;

logarithmic subtracted inversion recovery (lSIR) sequence

;

bipolar filter (BLAIR)

;

synthetic MRI

;

tissue boundaries

;

multiple sclerosis

;

whiteout sign

;

grayout sign

Subject:

Medicine and Pharmacology - Neuroscience and Neurology

Introduction

Ultra-high contrast MRI is a term used to describe MR imaging that shows abnormalities with high contrast when little or no abnormality is seen using conventional state-of-the-art MRI sequences. It is achieved without the use of increased static or gradient magnetic fields. For soft tissues, Xray CT can be regarded as high contrast, conventional state-of-the-art MRI as very high contrast and MRI which shows abnormalities with obvious contrast when they are not seen with conventional MRI as ultra-high contrast (UHC) [1].

A common mechanism for achieving UHC MRI is synergistic contrast in which changes in a single tissue property such as T₁ are used twice or more in the same sequence to increase contrast rather than just once which is the case with most conventional sequences [2]. An example is the divided subtracted inversion recovery (dSIR) sequence in which the signals from two inversion recovery (IR) sequences with different inversion times (TIs) are subtracted to increase the contrast produced by small changes in T₁. This contrast is increased further by division of the subtraction by the sum of the signals from the two sequences [3].

Another way of achieving UHC MRI is to use a high quality single tissue property map such as a T₁ or T₂ map and retrospectively apply a shaped filter to the map. The filter can increase contrast in one or more domains of tissue property values and maintain anatomical detail in the rest of the image. This is unlike narrow windowing of images which, when used to increase contrast, is accompanied by saturation of the display gray scale at the upper and lower signal thresholds. Narrow windowing leads to loss of anatomical detail and limits the degree to which contrast can be increased by windowing.

The principal clinical focus of UHC MRI to date has been on small changes in tissue properties from normal that are insufficient to generate useful contrast with conventional state-of-the-art sequences. When abnormalities are already seen with high contrast using conventional sequences, there is no particular clinical gain in demonstrating them with even greater contrast. As a result, attention has been directed at normal or near normal appearing tissues as seen with conventional state-of-the-art imaging.

Small changes in T₁ may produce obvious abnormalities with the dSIR sequence as illustrated in a case of mild traumatic brain injury (mTBI) in an 18-year-old male patient shown in Figure 1. No abnormality is seen on the T₂-FLAIR image in the patient (Figure 1A). The dSIR sequence in a normal age matched control shows normal white matter in the cerebral hemispheres as low signal (dark) (arrows) in Figure 1B. The same dSIR sequence in the patient shows abnormal high signal in his entire white matter (arrows) (Figure 1C). High signal (light) well defined boundaries are seen between normal white and gray matter on the dSIR image in Figure 1B. These are less obvious in Figure 1C because of the high signal in the abnormal white matter.

High contrast abnormalities in white matter have been seen using dSIR T₁-bipolar filter (BLAIR) images of the type illustrated in Figure 1 in areas where little or no abnormality has been seen with T₂-FLAIR or T₂-wSE sequences in cases of mTBI [4], multiple sclerosis (MS) [5], methamphetamine substance use disorder [6] and Grinker’s myelinopathy [7].

The purpose of this educational review is to describe the basic physics underlying UHC MRI using BLAIR images and illustrate use of these images in patients in order to facilitate wider use of them in clinical practice.

Basic Physics

Tissue Property Filters (TP-filters) and the Inversion Recovery (IR) Sequence

The first part of this paper is a condensed and updated version of work published previously [1,2,3]. It is included to make the paper self-contained and not require concurrent reference to other papers to understand the contrast of clinical images included in the paper.

Instead of illustrating image contrast in the usual way by plotting longitudinal and then transverse magnetisation (M_Z and then M_XY) against time using the Bloch equations as in Figure 2, it is possible to use plots of signals produced by the exponential recovery of longitudinal magnetisation and exponential decay of transverse magnetisation against T₁ and T₂ respectively (rather than against time). These plots are tissue property filters (TP-filters); the tissue properties may be mobile proton density, T₁, T₂, T₂*, D*, etc [8,9]. In this paper the term tissue is used to include fluids unless otherwise specified.

For TR >> T₁, the T₁ dependent part of the IR sequence has the signal equation:

S_T1 = 1-2e^-TI/T₁

(1)

where S_T1 is the signal from the T₁-filter, TI in the inversion time and T₁ is the longitudinal relaxation time. The T₁-filter of the IR sequence (with a constant value of TI) is shown in phase-sensitive or signed form in Figure 3A where it is a monotonic low pass filter, and in magnitude form in Figure 3B where it is a negative unipolar filter. The maximum size of the slopes of these two filters are shown with red lines in Figs. 3A and 3B. The magnitude forms of the IR filter S_T1 = (1-2e^-TI/T₁) for T₁ < and S_T1 = (2e^-TI/T₁ -1) for T₁ > are used in the subsequent text unless otherwise specified.

When TI is increased, the magnitude IR T₁-unipolar filter shifts to the right as shown in Figure 4. Figure 4A (left) shows the IR T₁-filter with a short TI_s (e.g., the Short TI IR or STIR sequence) for brain where gray matter (G) has a higher signal than white matter (W). The slope of the T₁-filter between W and G is moderately positive (red line).

When the TI is increased to an intermediate TI, TI_i, as in Figure 4B (centre), the T₁-filter is shifted to the right. W and G are fixed in the same position on the X axis, and W now has a higher signal than G. The slope of the T₁-filter between W and G is strongly negative (red line).

When the TI is increased further to a long TI, TI_l, the T₁-filter is displaced further to the right, as in Figure 4C (right) where W has a slightly higher signal than G, and the slope of the T₁-filter between them is mildly negative (red line).

Using the small change approximation of differential calculus, the contrast (difference in signal) produced by each T₁-filter from the increase in T₁ from W to G (horizontal green arrows in Figs. 4A, 4B and 4C) is the size of this increase multiplied by the slopes of the respective T₁-filters (i.e. by their first derivatives, the red lines). This contrast is shown by the vertical blue arrows in Figure 4 and is moderately positive in Figure 4A, highly negative in Figure 4B and mildly negative in Figure 4C. The slopes of the T₁-filters show the contribution of the sequence to contrast in each example.

Subtracted IR (SIR) and divided Subtracted IR (dSIR) Bipolar Filters

Two conventional magnitude IR T₁-unipolar filters with different TIs, namely TI_short = TI_s and TI_intermediate = TI_i are shown in Figure 5A. The subtraction: first TI_s T₁-filter minus the second TI_i T₁-filter produces the subtracted IR (SIR) T₁-bipolar filter shown in Figure 5B. The vertical dashed lines at the null points of the two IR T₁-filters shown in Figure 5A divide the X axis in Figure 5B into the lowest Domain (lD), the middle Domain (mD) and the highest Domain (hD). In the mD in Figure 5B (between the two dashed lines) the size of the slope of the SIR T₁-bipolar filter is about double that of the IR T₁-filters shown in Figure 5A. This is because the slope of the SIR filter in its mD in Figure 5B is the positive slope of the TI_s filter in the mD shown in Figure 5A minus the negative slope of the TI_i filter in the mD also shown in Figure 5A.

The two IR T₁-filters shown in Figure 5A can also be added to give the Added IR (AIR) T₁-filter shown in Figure 5C. In its mD, which is bounded by the vertical dashed lines, the signal is reduced to about 0.20 compared with its value of 2 at T₁ = 0 (i.e. about one tenth).

Figure 6A shows a divided SIR (dSIR) T₁-bipolar filter in which the SIR T₁-bipolar filter in Figure 5B is divided by the AIR T₁-filter in Figure 5C i.e. . The signal (S_dSIR) of the dSIR filter is given by:

S_{d S I R} = \frac{S_{T I s} - S_{T I i}}{S_{T I S} + S_{T I i}}

(2)

where S_TIs is the signal from the shorter TI filter and S_TIi is the signal from the intermediate TI_i filter.

The resulting dSIR T₁-bipolar filter (Figure 6A) shows a negative pole and a positive pole. Its mD shows a highly positive nearly linear slope which is about ten times greater than the slopes of the IR T₁-filters shown in Figure 5A. The dSIR filter has maximum and minimum values of ±1. (The contributions from mobile proton density and T₂ to conventional IR sequences cancel out with dSIR sequences which are therefore T₁ maps; the whole sequence can be accurately represented by its T₁-bipolar filters.)

Figure 6B compares the contrast produced by a conventional S_TIs IR T₁-unipolar filter (pink) to that from the SIR T₁-bipolar filter (blue) shown in Figure 5B from the same increase in T₁ (horizontal positive green arrow, DT₁). DT₁ is multiplied by the slopes of the respective S_TI IR and SIR filters (red lines) to produce the differences in signal DS, i.e. the contrast from each of them. These are shown by the vertical pink and blue arrows on the right. The SIR T₁-bipolar filter generates about twice the contrast (blue arrow) of the S_TIs IR T₁-unipolar filter (pink arrow) from the same increase in T₁, DT₁.

Figure 6C compares the contrast produced by the S_TIs IR T₁-unipolar filter (pink) to that produced by the dSIR T₁-bipolar filter shown in Figure 6A (blue). The increase in T₁ (horizontal positive green arrow, DT₁) is multiplied by the slopes of the respective S_TIs IR and dSIR filters (red lines) to produce the differences in signal DS, or contrast generated by the two filters. These are shown by the vertical pink and blue arrows on the right. For the same increase in T₁ (DT₁), the dSIR T₁-bipolar filter produces about ten times greater contrast than the S_TIs IR T₁-unipolar filter. The S_TIs IR T₁-filter is that of a conventional TI IR sequence such as MP-RAGE (magnetisation prepared-rapid acquisition gradient echo).

To produce the large increase in contrast shown in Figure 6C, the dSIR sequence is typically targeted at the small increase in the T₁ of normal tissue produced by disease (shown by the horizontal green arrow) and this is positioned within the steeply sloping mD of the dSIR T₁-bipolar filter. This is done by choosing appropriate values of TI_s and TI_i. The target is often small increases in the T₁ of white matter. High contrast can also be produced by difference or change in T₁ within the lower part of the hD and the upper part of the lD where the dSIR T₁-bipolar filter has relatively steep slopes. More detail is available in reference [3].

Contrast at Tissue Boundaries

MRI tissue boundaries take different forms such as a gradual change in signal from one tissue to another as well as sharply defined high and low signal boundaries between tissues. In the boundary region between two pure tissues (such as between white and gray matter) the T₁s of voxels which contain mixtures of the two tissues typically span the range of T₁ values between those of the two pure tissues. If a narrow mD dSIR T₁-bipolar filter (e.g., with TI_s nulling normal white matter, and TI_i longer than TI_s but less than that needed to null gray matter) is used there are mixtures of the two pure tissues in voxels in the boundary region with intermediate T₁ values (T_1W,G) which correspond with the peak of the T₁-bipolar filter as shown in Figure 7. This produces a high signal boundary between white matter and gray matter. The dSIR T₁-bipolar filter shows much higher contrast than that produced between white matter and gray matter by a conventional white matter nulled IR sequence such as MP-RAGE (Figure 7). The mechanism shown in Figure 7 produces the high signal S_W,G at the boundary between normal white and normal gray matter shown in Figure 1B. High spatial resolution boundaries are also seen at the junction between CSF and white matter at ventricular boundaries, and can be seen at the boundary between gray matter and CSF with wide mD sequences [3]. The high and low signal boundaries on dSIR images are iso-T₁ contours.

T₁ maps and Qualitative - Quantitative MRI

To better understand the dSIR T₁-bipolar filter, a linear equation of the form y = mx + c can be used to approximate the filter in its mD. The equation is produced by fitting a straight line between the first and last points of the mD (i.e. first point x = and y = -1, and last point x = and y = +1). In the mD, the signal of the dSIR sequence S_dSIR, is given by:

S_{dSIR} \approx \frac{\ln 4}{∆ T I} T_{1} - \frac{Σ T I}{∆ T I}

(3)

where ΔTI = TI_i - TI_s (i.e., second TI minus first TI) which is positive, and ΣTI = TI_s + TI_i which is also positive. Note that because DTI is positive, the slope

\frac{\ln 4}{∆ T I}

is positive. The offset is negative.

The expression in Eq. (2) illustrates four key features of the dSIR T₁-bipolar filter, firstly, the near linear change in signal (i.e. S_dSIR) with T₁ in the mD, secondly, the filter has a slope in the mD equal to , thirdly the filter shows high sensitivity to small changes in T₁ when the size of ΔTI is small. When DT₁ is small, the size of DTI can be decreased to match it and so scale up the sensitivity of the filter. The reduction in DTI increases the steepness of the T₁-filter in the mD and thus the amplification of contrast produced from DT₁. This compensates for the small value of DT₁ until the sequence becomes SNR and/or artefact limited. This is not the case with conventional IR sequences where, if DT₁ is decreased, contrast decreases.

Fourthly, Eq. (3) can be used to calculate T₁ values in the mD so that:

T_{1} \approx \frac{∆ T I}{l n 4} S_{dSIR} + \frac{Σ T I}{l n 4}

(4)

The linear approximation is only valid in the mD. Also, it is assumed that TR is long compared to tissue T₁ values. An equivalent expression that allows for incomplete recovery of longitudinal magnetisation during TR can be formulated.

Outside of the mD in the lD and hD, the magnitude of the dSIR T₁-bipolar filter signal decreases towards zero at minimum and maximum values of T₁. T₁s in the mD occupy the full range of the display gray scale from black to white. T₁s from the shortest and longest T₁s in the lD and hD respectively only occupy parts of the display gray scale.

Log Then Subtracted Inversion Recovery (lSIR) Sequences

Another bipolar filter is the natural logarithm (ln) then subtracted inversion recovery (lSIR) filter [10]. The signal of this filter (S_lSIR) is half of the subtraction: ln short S_TIs T₁-filter minus ln intermediate S_TIi T₁-filter. The ln S_TIs and ln S_TIi unipolar negative T₁-filters are shown in Figure 8 and have sharper negative poles than those of unipolar S_TIs and S_TIi T₁-filters shown in Figure 5A.

Figure 9 shows a dSIR T₁-bipolar filter in blue and the corresponding lSIR T₁-bipolar filter with the same TIs in orange. The two filters appear similar in the lowest region of the lD and in the highest region of the hD as well as in the central region of the mD. The slope of the dSIR filter is essentially constant and equal to . The magnitude of the slope of the lSIR filter is the same that of the dSIR filter at the centre of the mD but greater around the lower and higher nulling TI values. The lSIR filter asymptotically approaches negative and positive infinity at the two corresponding T₁s. The lSIR filter increases its contrast amplification as the change in T₁ from normal becomes smaller. Like the dSIR filter, the lSIR filter essentially eliminates the mobile proton density and T₂ dependence of the full IR sequence and so can be used to model the full sequence behaviour.

The lSIR filter is an inverse hyperbolic tangent of the dSIR filter in the mD (S_lSIR = atanh S_dSIR), and an inverse hyperbolic cotangent of the dSIR filter in the lD and hD.

From a practical point of view, the lSIR T₁-filter has 2-3 times the slope of the dSIR T₁-filter for very small differences or changes in T₁ around the null points and thus has 20-30 times greater contrast than conventional IR sequences for the same very small differences or changes in T₁. The maximum and minimum values of the lSIR filter are usually shown as ±2-3 (rather than ±infinity). This compares with values of ±1 for the dSIR filter. The low and high signal boundaries of the lSIR filter are narrower than those of the corresponding dSIR filter.

The lSIR filter shows a selective increase in slope, or sharpening, at iso-T₁ contours close to the nullpoints. It is not generalised edge enhancement of all boundaries on an image.

Composite (c) bipolar filters (T₁ as well as T₂, T₂*, D*, c and/or MT)

It is possible to introduce an attenuation filter S_OF (other filter) and apply it to one of the two IR filters used for dSIR T₁-bipolar filter imaging in Figs. 5 and 6. S_OF may equal e^–DTE/T2, e^–DTE/T2* or e^–DbD* which introduce T₂, T₂* and D* dependence respectively. This provides additional contrast to the T₁ contrast of the dSIR bipolar filter. When DTE equals zero, or Db equals zero, S_OF = 1. If S_OF is set to 0.5 and 0.25, the curves shown in Figs. 10A and 10B result when using a narrow mD dSIR sequence (TIs = 350 ms and 500 ms). The nulling times are the same for all values of S_OF.

In Figure 10A attenuation of the first IR filter (with TI_s) is shown with S_OF set at 1 (blue), 0.5 (orange) and 0.25 (yellow). The yellow and orange curves are wider around the outside of the first (negative) pole at the first TI_s and are narrower and inside the second (positive) pole at TI_i. This results in a broadening and loss of contrast and lower signal white matter around the first negative pole as well as increased contrast and narrower boundaries between white and gray matter around the positive pole. S_OF can be set negative values to reverse the sign of the relevant contrast and make it synergistic with the T₁ contrast.

In (A), as S_OF is decreased to 0.5 and 0.25 the filters around the first negative pole become wider and are positioned outside the blue filter. The attenuated filters around the second positive pole show higher contrast and become narrower. They are located inside the blue filter. (B) shows the conventional T₁-bipolar filter (blue) with S_OF = 1, and attenuated second TI_i filter S_OF = 0.5 (red) and a more attenuated second TI_s filter (S_OF = 0.25) (yellow). As S_OF is decreased, the filters around the first negative pole become narrower and show higher contrast. They are positioned within the blue filter. The attenuated filters around the second positive pole become wider and show lower contrast. They are positioned outside of the blue filter. For T₂, a short T₂ increases attenuation relative to a longer T₂ and this corresponds to the narrower and sharper peak with greater negative contrast and spatial resolution of the contrast (A).

In Figure 10B, attenuation of the second IR filter (with TI_i) is shown. As S_OF is decreased from 1.0 to 0.5 and 0.25, the filter narrows around the first (negative) pole and widens around the second (positive) pole. The negative contrast produced for the same difference in T₁ is increased around the first negative pole. Composite (c) forms of the bipolar filters can be designated as cSIR, cdSIR, clSIR etc.

Susceptibility and magnetisation transfer can also be used to attenuate signals and produce supplementary contrast. The combination of T₁ and T₂* contrast may be of particular value in demonstrating the central vein sign and paramagnetic rim sign in multiple sclerosis (MS) as well as late sequelae of haemorrhage and fMRI in combination with perfusion.

Synthetic dSIR and lSIR T₁-bipolar Filter Images

Synthesizing narrower mD images from wider mD dSIR and lSIR T₁-bipolar filter images

It is possible to synthetically create narrower mD dSIR images from wider mD dSIR images using Eq. 1 and the relationships illustrated in Figure 6A. The TIs of the narrower mD images are within those of the wider mD when using magnitude forms of the IR sequences.

Synthesizing dSIR and lSIR T₁-bipolar filter images from T₁ maps

dSIR and lSIR T₁-bipolar filter images with particular values of TI can also be synthesised from T₁ maps using Eq.1 and the relationships shown in Figs. 6A, 8 and 9. T₁ and other tissue property maps (T₂, T₂*, D*) have linear TP-filters with maximum signal values at the high end of the tissue property X axis. As a result, CSF has high signals on the maps and this can dominate their appearance. The problem can be avoided when appropriate T₁-bipolar filters are applied to T₁ maps to create T₁-bipolar images in which CSF signal is reduced. The equivalent mapping using MP2RAGE acquisitions has also been described previously. While in principle the same contrast is obtained, there are wide discrepancies in the T₁s measured by different techniques suggesting T₁ is not independent of the method of acquisition. The unaccounted variance may mask small white matter abnormalities observed using dSIR images that were not previously seen with general purpose T₁ mapping.

Synthesizing T₂-, T₂*-, D*andc-bipolar filter images using purpose designed bipolar filters together with T₂, T₂*, D* andc maps

It is also possible to create TP-bipolar filter images from tissue properties other than T₁ using purpose designed synthetic TP-bipolar filters. These filters are typically based on the appearance of the dSIR T₁-bipolar filter shown in Figs. 6A and 6C. They can have a linear function in their mD and a scaled reciprocal function (e.g., y = ) or similar function in their lDs and hDs. These functions join the linear function in the mD at signal values of -1 and +1, respectively. Synthetic TP-bipolar filters can then be used with T₂, T₂*, D* and c maps to create T₂-, T₂*-, D*-and c-bipolar filter images.

The lSIR inverse hyperbolic tangent and cotangent functions can also be used in the same way to create synthetic lSIR TP-bipolar filter images from tissue property maps.

Directly acquired T₁-bipolar filter images and synthetic T₁-, T₂-, T₂*-, D*-, c- and MT-bipolar images with different TPs can be multiplied together to produce synergistic multi TP-bipolar filter images. Under log then subtraction (lSIR) the multiplicative contrast become additive. The signs of the slopes of the TP-bipolar filters (negative or positive) can be matched with the signs of the changes in tissue properties (negative or positive) to ensure the contrast produced by the different tissue properties is synergistic [2].

Window Levels and Widths of Displayed Images

Conventional changes in window level and width of MR images

When an image is displayed, the window level sets the brightness (signal level) and the window width sets the contrast (signal difference). Changes in the window level and width of images change their signals on the gray scale display from S to S_D and can be described using a high pass filter as seen in Figure 11A. This plots the display signal S_D against the image signal S. This is a signal (S) filter not a TP-filter. Contrast (the difference in signal between the S_Ds of P and Q) may be increased by narrowing the window width as shown in Figure 11B so that P and Q are now respectively shown at the bottom and top ends of the display scale. With narrow windowing, signals with values at the upper end of the X axis from Q and beyond are all shown with the same high value of S_D. As a result, there is no contrast between these voxels and they coalesce into blocks. This leads to a loss of anatomical coherence. The same also happens at the lower end of the S scale where values of S less than P have the same signal in Figure 11B and coalesce. This saturation of signals establishes a practical limit to how tightly images can usefully be windowed (i.e., how narrow the window width can usefully be made) and thus how much the contrast of the displayed image can be increased before the loss of anatomical coherence leads to a lack of credibility of the images. This is in addition to limits imposed by SNR considerations and image artefact levels.

dSIR images

dSIR images can decrease DTI and so increase contrast. Maximum and minimum values are shown as high and low signal boundaries as in Figure 11C, T₁ values lower than the lower nulling T₁ value and greater than the greater nulling T₁ value do not appear on the same plateaux as with conventional narrow windowing in Figure 11B, but vary in signal because of the sloping regions in the lDs and hDs of the dSIR T₁-bipolar filter (Figure 11C). As a consequence, contrast is maintained and anatomical detail is preserved in all three Domains when DTI is decreased. As a result, contrast can be increased by decreasing DTI until it becomes SNR and/or artefact limited.

Methods

With approval from the New Zealand Health and Disability Ethics Committee 20/CEN/246 (26 June 2020), 21/CEN/246 (1 May 2021); 20/NTB/14 (21 May 2020); UAHPEC 018466 (20 Dec 2019); and the Northern B Health and Disability Ethics Committee 2022 EXP 13106 (18 Nov 2022) as well as informed consent from each subject, MR scans were performed on four adult normal controls, seven adult patients with mild traumatic brain injury (mTBI), MS, methamphetamine substance use disorder, Grinker’s myelinopathy, Parkinson’s disease and white matter in a patient with a glioma. The research was conducted in accordance with the Declaration of Helsinki. 3T scanners (General Electric Healthcare Premier and Siemens Healthineers Skyra) were used. 2D IR FSE sequences were performed with a TI_s chosen to null the shortest T₁ normal white matter (or slightly shorter) and a longer TI_i chosen to produce narrow mD dSIR images targeted at small increases in the T₁ of white matter from normal as illustrated in Figure 6C. 3D wide mD dSIR images with a short TI_s and a long TI_l were also obtained [5]. Longer TE 2D IR images were acquired to produce cdSIR images. Positionally matched 2D and 3D T₂-FLAIR images were obtained as well as susceptibility weighted filtered images as described in Table 1.

Illustrative Cases

Mild Traumatic Brain Injury (mTBI)

Conventional MRI images usually show little or no change in mTBI but narrow mD dSIR images frequently show extensive high signal areas in white matter (the whiteout sign) (Figure 1). This may resolve within as little as two days or persist over years. It is generally attributed to neuroinflammation but oedema, demyelination and degeneration may also have a role.

In addition to the whiteout sign, there may be loss of contrast in the gray matter of the thalamus due to an increase in T₁. This is manifest as a low contrast appearance between the medial and lateral thalamus (Figure 12A) (arrows) (a grayout sign). After remission normal high gray white matter contrast is restored (Figure 12B) (arrows). Contrast between the lateral thalamus (arrows) and the posterior limb of the internal capsule (PLIC) is low in (Figure 12A) but very high in (Figure 12B). The low contrast between the thalamus and the PLIC in (A) results from a reduction in signal in the gray matter of the lateral thalamus due to an increase in its T₁ in the hD as well as an increase in signal in the white matter of the PLIC due to an increase in its T₁ in the mD. This is reversed in (Figure 12B) where the normal high contrast boundary is restored.

Figure 13A shows a normal control with obvious contrast between the heads of the caudate nuclei and the adjacent CSF as well as between the cortex and CSF. Figure 13B shows a patient with mTBI who has a whiteout sign. There are grayout signs in the thalamus and putamen. In addition, contrast is lost between the heads of the caudate nuclei and CSF as well as between cortex and CSF which appear isointense. These are grayout signs. Thus, the patient has a combination of whiteout and grayout signs.

There may also be gray matter changes elsewhere in mTBI. The normal red nucleus has a relatively short T₁ resulting in its normal signal being mid-gray in the mD (arrows) (Figure 14A). It shows a higher signal due to an increase in T₁ in the case of mTBI (arrows) (Figure 14B).

Multiple Sclerosis (MS)

Focal lesions, as well as patchy white matter changes, may be seen in areas that appear normal on T₂-wFSE images in a patient in remission (Figure 15) (compare with normal low signal appearance of normal peripheral white matter in Figure 1B). In another case during a relapse, one lesion is seen on the T₂-FLAIR image in Figure 16A (long arrow). This lesion is also seen in Figure 16B (long arrow). There are an additional six lesions (small arrows) in Figure 16B. Some of the MS lesions in Figure 16B have high signal boundaries. This can be due to a large increase in T₁ in the lesion or a decrease in the abnormal T₁ of a lesion due to the presence of paramagnetic free iron. Some of the lesions with high signal boundaries shown in Figure 16B have paramagnetic rims on the filtered susceptibility weighted images (arrows, Figure 16C).

In addition, there are widespread abnormal areas in white matter which are only seen on the dSIR images (Figure 16). These changes are typically bilateral, symmetrical and have an increased signal. They have the features of a whiteout sign (Figure 16B) and to date have only been seen in MS patients during a relapse.

Composite T₁ and T₂ cdSIR images may show higher contrast than corresponding dSIR images (Figure 17A,B).

In the spinal cord in another patient with MS, the T₂-FLAIR image only shows an ill-defined smudge (arrow) (Figure 18A). The corresponding wide mD dSIR image (Figure 18B) shows an extensive, well defined, high contrast abnormality (arrows). Another lesion is seen in the medulla on the dSIR image (arrow) but not on the T₂-FLAIR image.

Methamphetamine Substance Use Disorder

Patients with methamphetamine substance use disorder may show a heterogenous pattern of white matter abnormalities with dSIR sequences but they may also show clearly defined whiteout signs and these can remit with continued abstinence (Figure 19) (compare with normal appearance in Figure 1B). The situation may be confounded in the case shown by a previous mTBI.

Delayed Post-Hypoxic Leukoencephalopathy (Grinker’s Myelinopathy)

This is thought to be a rare syndrome in which classically patients make an initial recovery after a hypoxic and/or ischemic episode, and then deteriorate with severe neurological signs such as Parkinsonism and akinetic mutism. Conventional MRI shows widespread abnormalities in the white matter of the cerebral hemispheres. In this type of clinical situation, it is possible to see extensive abnormalities in white matter with dSIR sequences where no abnormalities are seen with conventional sequences [7]. The patients typically have cognitive symptoms of less severity than those included in classical descriptions. This less severe form of Grinker’s myelinopathy may be much more common than the classical form but not be recognised using conventional MRI sequences.

Parkinson’s Disease

Parkinson’s disease may show more obvious bilaminar signs in the cortex than in age matched controls (i.e., an increased signal in the outer layer of the cortex) with narrow mD dSIR sequences (white arrows) (Figure 20). The inner higher signal is from the boundary between white matter and gray matter and the next layer outwards is lower signal in the inner cortex. The next outer layer beyond this is the higher signal layer in the outer cortex.

Bubble signs in which circular areas of similar size are seen in the gray matter in the basal ganglia and thalamus can also be seen (e.g., black arrows) (Figure 20).

White Matter Associated with Cerebral Tumours

In cases of cerebral tumours, extensive white matter changes have been seen in areas that appeared normal on T₂-FLAIR images after chemotherapy (Figure 21A and 21B). When typical radiation or chemotherapy changes have been seen on T₂-FLAIR images, additional changes have been observed in the form of whiteout signs.

Normal Control and MS Patient with dSIR and lSIR Images

dSIR and lSIR images (Figure 22A and 22B) are compared in a 46-year-old male normal control. White matter gray matter boundaries and the bilaminar cortex sign (arrows) are seen with higher contrast and higher spatial resolution on the lSIR image. Bubble signs are also better seen in the thalamus and putamen on the lSIR image. A very small decrease in T₁ through the sharp peak of the lSIR T₁-bipolar filter results in finer boundaries and narrower margin bubble signs with the lSIR T₁-bipolar filter compared to the dSIR T₁-bipolar filter.

In a 41-year-old female patient with MS, the dSIR image (Figure 23A) shows a blurred leukocortical lesion in the right medial frontal region with no evidence of a white matter/gray matter boundary within it. The matched lSIR image (Figure 23B) shows a disrupted high signal boundary between white matter and gray matter within the lesion. Profiles (plots of signal against position in mm) through the lesion at the level shown by the arrows for the dSIR (blue) and lSIR (orange) images shows a higher signal and steeper slope for the lSIR image (Figure 23C). The spatial resolution of the contrast i.e., change in signal with change in position is generally higher on the lSIR image.

Discussion

UHC MRI using Bipolar Filters

UHC MRI is a term used to describe MRI techniques which show much greater contrast than conventional state-of-the-art MRI images while having similar spatial resolutions and comparable acquisition times. This is achieved with no increase in static or gradient magnetic field. UHC MRI using bipolar filters employs sequences existing on most MRI machines with minimal image processing. It can be implemented in a very short time at very little cost.

UHC MRI is distinguished from ultra-high field MRI which uses static magnetic field strengths of 7T or higher to improve image signal-to-noise ratio. This improvement may be used to increase contrast, spatial resolution and/or the speed of MRI scans [11].

UHC MRI is also distinguished from ultra-high spatial resolution MRI where typically increased gradient strength is used to improve spatial resolution [12].

The three main types of directly acquired bipolar filter (BLAIR) used to achieve UHC MRI in this paper are summarised in Table 2 with some of their variants. Their signal equations and other relevant functions are included in Table 3. This table also provides information for synthetic imaging.

The dSIR sequence is the BLAIR sequence most used to date. It has been validated in quantitative phantom and human studies [13,14].

While the emphasis in this paper has been on white matter in the mD and grey matter in the hD, very high contrast can also be seen in the central grey matter of the brain when this is in the mD rather than in the hD. In addition, iron containing tissues have only been studied with dSIR using the single tissue property T₁, not the composite sequence cdSIR using the two tissue properties T₁ and T₂*.

Understanding the contrast produced by TP-bipolar filter sequences is greatly helped by use of their TP-filters. These filters relate differences or changes in tissue properties to contrast using their slopes and explain signal and contrast patterns that can otherwise appear inexplicable using conventional approaches [8,9].

Artefacts from too high a value of TI_i, misregistration, partial volume effects and other causes are described in reference [5]. Means of avoiding these and/or reducing their effects are included.

Tissue and Fluid Boundaries

dSIR and lSIR images are characterised by sharply defined high contrast, higher spatial resolution boundaries between white matter and gray matter as well as between white matter and CSF, and between gray matter and CSF. The site of the boundaries as well as their width can be changed by varying the TIs of the sequence [5]. In addition, there are well defined low signal (dark) boundaries between white matter, gray matter and CSF as modelled by TP-bipolar filters.

High signal in-plane boundaries on dSIR images usually have a uniform and consistent appearance although these are subject to partial volume effects and may simulate lesions particularly with 2D images acquired with relatively thick slices (e.g., 2-4 mm). Through-plane partial volume effects also produce high signal from boundaries and these are more prominent with thick slices.

A systematic treatment of boundaries including variation in signal with position using TP-filters, the partial derivative change in T₁ with tissue fraction and the partial derivative change of tissue fraction with position is described in reference [5]. Using this approach, differences in signal (i.e., contrast) can be related to differences in position within the boundary region. The dSIR and lSIR filters produce higher contrast and higher spatial resolution of the contrast at boundaries between tissues and fluids. The boundaries are specific and located at iso-T₁ contours. They are not part of a generalised pattern of edge enhancement.

Hybrid Quantitative and Qualitative Imaging

dSIR sequences produce single images that can be used both quantitatively and qualitatively. The images show high contrast and this can be used for conventional qualitative image interpretation. The images are also T₁ maps so they can be used to directly measure T₁ values in regions of interest placed on them for quantitative studies. In comparison with the conventional approaches to T₁ quantitation which use separate images for qualitative interpretation and T₁ mapping, dSIR images: (i) do not require working with images of two different types; (ii) allow precise placement of ROIs in relation to lesions on the images used for radiological interpretation, and (iii) require no extra time for additional acquisitions.

T₁ values of dSIR sequences have shown comparable accuracy to conventional T₁ mapping in the mD but, as expected, differences in the lD and hD [13,14].

Magnetisation Transfer (MT)

With 2D multislice acquisitions, there may be incidental magnetisation transfer (MT) effects from off-resonance 180° or similar pulses used during FSE acquisitions. Using two pool modelling, MT decreases the observed mobile proton density in the free pool but also decreases the observed T₁ of the free pool in proportion to the reduction in observed mobile proton density [15]. This reduction may be quite substantial and decrease normal observed T₁s (T_1obs) by as much as 50%. In general, in disease MT effects are decreased so the reduction in T_1obs in diseased tissue due to MT is generally less than that in normal tissue. This is manifest as an apparent increase in T_1obs in diseased tissue relative to normal tissue in addition to any increases in T₁ in diseased tissue for other reasons e.g. pathological change such as oedema.

2D FSE images with higher echo train lengths (ETLs) require more 180° or similar pulses and thus produce a corresponding increase in MT effects during multislice acquisitions. These produce a greater reduction in values of T_1obs and so need shorter nulling TIs. MT effects are less with gradient echo acquisitions as usually used with 3D acquisitions. As a result, T_1obs may be longer with 3D gradient echo acquisitions than with 2D FSE acquisitions and so their nulling TIs may need to be correspondingly longer. The contrast between normal and abnormal tissues may also be lower with gradient echo acquisitions.

MT following IR can also be used as an attenuating filter S_OF in composite bipolar filters.

Synthetic T₁-bipolar Filters

Synthetic bipolar filters can be of two types as shown in Table 2. The first uses Eq. 1 and the relationships shown in Figs. 6, 8 and 9 to synthesise dSIR and lSIR images with different TIs from two directly acquired images. The second type uses a synthetic TP-bipolar filter which may be the same or similar to the dSIR and lSIR filters and applies these to tissue property maps. This approach has two aspects: (i) acquiring the tissue property map which should be of high quality in the tissue property Domains or parts of Domains of interest, and (ii) using a visualization function such as a bipolar filter to increase contrast without saturating voxel values. These two aspects are incorporated in directly acquired dSIR and lSIR imaging.

Signs Produced with dSIR Sequences

The Whiteout Sign

This is a generalised bilateral symmetrical increase in signal in white matter of the cerebral and cerebellar hemispheres as well as the brainstem [1]. The increase in signal follows the pattern of normal white matter signals and can involve all, or almost all, of the white matter in a particular slice. It has specific features such as sparing of the anterior and posterior central corpus callosum (genu and splenium) as well as peripheral white matter in the cerebral hemispheres (Figs. 1, 12, 13, 14, 19, 21). The white matter signal can increase to that of the high signal boundary between white and gray matter and can be divided into five grades [1]. It is due to widespread, small increases in T₁ and is thought to be particularly associated with neuroinflammation but may also be due to oedema, demyelination, degeneration and other causes. It has also been seen in mTBI, hypoxic/ischaemic injury to the brain, multiple sclerosis and white matter associated with tumours with/without treatment.

Grayout Signs

With grayout signs there is a loss of contrast in gray matter such as is seen in the thalamus and basal ganglia after mTBI (Figs. 12, 13), or a loss of contrast between gray matter and white matter, or between gray matter and CSF (Figure 13B). Grayout signs are generally more subtle than the whiteout sign. They are associated with an increase in T₁ in gray matter.

In young subjects, the normal lateral thalamus gray matter (TG_lat) tends to have a shorter T₁ than the normal medial thalamus gray matter (TG_med) (upper horizontal negative solid green arrow). This leads to higher signal in the lateral thalamus with a narrow mD dSIR sequence (first vertical blue arrow) (Figure 24). The T₁ of the lateral thalamus may increase in disease (middle horizontal positive green striped arrow in Figure 24) which produces a loss of contrast with the medial thalamus (blue circle). With reversion back to normal on recovery, the T₁ may decrease back to normal (lowest horizontal negative green arrow) restoring contrast (second vertical blue arrow).

The whiteout sign may be accompanied by grayout signs in conditions such as mTBI as described in this paper. The term post-insult leukoencephalopathy syndromes (PILS) used in reference [1] only refers to whiteout signs after mTBI, Grinker’s myelinopathy etc. A term that includes both the whiteout sign and grayout signs as well as the possibility of ongoing effects (rather than just post-insult effects) is reactive encephalopathy (RE). Depending on the observed features (presence/absence of grayout signs), the pattern may be described as PILS/RE. Experimental studies to determine histological origins of the signs are in progress.

The Bilaminar Cortex Sign

This is an increase in signal in the outer layer of the cortex which is often barely seen in young adults but becomes more obvious with increasing age as well as in Parkinson’s disease. It is seen with narrow mD dSIR images in which there is a high signal at the white matter gray matter junction, a lower signal layer in the inner cortex and a higher signal layer in the outer cortex (Figs. 20, 22). The higher signal generally corresponds to layers 1-3 and the lower signal to layers 4-6 of the cortex. The increase in signal is usually most obvious at the crest of gyri and to a lesser extent in the banks of sulci. It may be obscured by partial volume effects particularly on 2D images of 4 mm thickness compared with the estimated 2-4 mm thickness of the cortex. The sign may be due to increased free iron content with a decrease in T₁ of the normal cortex with increase in age, with an exaggeration of this in Parkinson’s disease [16]. It is distinguished from the double cortex sign [17] which is due to iron deposited in a subcortical location. A bilaminar cortical appearance has been seen with MP2RAGE sequences at 7T after detailed image processing and has been attributed to differences in myelin [18]. The bilaminar cortex sign differs from the high signal that can be seen at the cortical margin on T₂-FLAIR images. This signal is attributed to increased T₂ relative to brain in cortical veins/pia arachnoid (as opposed to a decreased T₁ in the bilaminar cortex sign). It is usually more evident on 2D T₂-FLAIR images with longer TIs and TRs than 3D T₂-FLAIR images though the latter may have very long effective TEs. It can be present in young subjects who do not show a bilaminar cortex sign.

Bubble Signs

With the narrow mD dSIR sequence nulling white matter, bubble signs are seen with an increase in T₁ that crosses the high signal boundary (“overshoots”) and enters the hD (Figure 16B) [5]. There is an outer high signal margin and a lower signal central region [3]. Bubble signs may also be seen in the gray matter of the thalamus and basal ganglia with narrow mD dSIR sequences when the T₁ in the gray matter decreases and crosses the high T₁-bipolar filter peak. Bubble signs become more obvious with age and are seen in Parkinson’s disease. The decrease in T₁ in gray matter sufficient to cross the high signal boundary of narrow mD dSIR sequences is illustrated in Figure 20 and probably reflects heterogeneous distribution of free iron. Bubble signs may also be produced in MS by release of free iron at the rim of lesions decreasing the T₁ in abnormally increased T₁ lesions.

Figure 25 illustrates the origin of the bubble sign in the thalamus. The T₁ of lateral gray matter is decreased from normal (TG_latnorm) to abnormal (TG_latab). In doing so, it passes through the peak signal where there is a partial volume effect between TG_latnorm and TG_latab. This results in a high signal on the periphery and a lower signal centrally in the bubble as shown by the vertical blue arrows in Figure 25.

Primary Cerebral Hemisphere Cortices

The primary motor, somato-sensory and visual cortices have relatively greater myelin than elsewhere in the cortex, and this results in a shorter T₁ and so higher signal in the hD relative to cortical signal elsewhere when using narrow mD dSIR sequences.

Clinical Use of T₂-FLAIR and Tissue Property BipoLAr fIlteR (TP-BLAIR) Sequences

T₂-FLAIR and T₁-BLAIR sequences are complementary. The T₂-FLAIR images show contrast due to larger changes in T₂ and T₁-BLAIR sequences such as dSIR show abnormalities due to smaller changes in T₁ in white and/or gray matter in areas that appear normal or near normal with conventional T₂-FLAIR sequences (e.g., Figs. 16A and 16B). In many diseases there are concurrent increases in both T₁ and T₂ so that the T₁-BLAIR sequence to match the T₂-FLAIR sequence is sensitive to changes in T₁ and does not have to be specifically sensitive to changes in T₂.

Other bipolar sequences can be used besides directly acquired dSIR sequences. These include synthetic dSIRs as well as T₂ and other tissue property synthetic TP- bipolar filters. Directly acquired T₁-bipolar filter sequences may be made sensitive to T₂, T₂* and D* as well as T₁ as composite bipolar filters e.g. T₁,T₂-BLAIR. There is an advantage in T₁-BLAIR imaging with dSIR sequences since they can be directly acquired and do not require acquisition of a tissue property map.

Other Sequences

The Double Inversion Recovery (DIR) Sequence

The DIR sequence multiples two IR sequences together and is one of the multiplied added subtracted and/or divided IR (MASDIR) group of sequences [5]. It has a unipolar T₁-filter as well as T₂ and mobile proton density filters. Initially, the two IR sequences were used to null fat and CSF [19] but in later versions either white or gray matter was nulled and so was CSF [20]. Lesions which increase T₁ and T₂ contrast produce synergistic contrast. When white matter is nulled with the conventional DIR sequences, this can be of particular value for imaging lesions of the cerebral cortex.

In a study using a gray matter and CSF nulled DIR sequence in patients with MS so that only white matter signal remained, Tillema et al found a dark rim around lesions in white matter with longer T₁s and this helped in the diagnosis of MS [21]. The rim was due to nulling of T₁s in the boundary region between normal white matter and lesions which was evident if the T₁ of the MS lesions was increased beyond the T₁ of nulled gray matter. This was seen more often in MS than other diseases. The nulling effect is dependent on the T₁ of the mixture of tissues matching the first nulling TI of the DIR sequence. It corresponds to the high signal bubble sign seen in white matter (Figure 16B). Release of paramagnetic free iron could also have a role.

Magnetisation Prepared 2 Rapid Acquisition Gradient Echo (MP2RAGE) [22]

This sequence multiplies the signals of two IR sequences together and divides this by the sum of the squares of the signals from the two sequences. Complex multiplication is employed. Signal values vary between ± ½. It is a negative unipolar T₁-filter for brain and has a moderate negative slope between white matter and gray matter with a shallower positive slope between gray matter and CSF (Figure 26, blue curve labelled FLAWS-uni which is identical with MP2RAGE). The sequence places white matter at the upper end of the display scale, gray matter at the lower end of the display scale and CSF at an intermediate level which is similar to the divided reverse SIR (drSIR) sequence but the reverse of the dSIR sequence.

Fluid and White Matter Suppression high contrast (FLAWS-hc), FLAWS- high contrast opposite (FLAWS-hco) and FLAWS-minimum (FLAWS-min) [23,24]

The FLAWS-hc and FLAWS-hco sequences have the same general mathematical form as the dSIR and drSIR sequences i.e. subtraction of two IR sequences divided by their sum i.e. , but differ in important ways: firstly, the TIs used are widely separated as opposed to dSIR and drSIR sequences where the TIs are usually narrowly separated in order to specifically target small differences or changes in T₁ in tissues. Secondly, complex multiplication is used and the images are in phase-sensitive form not magnitude form. In Figure 26 the FLAWS-hc filter (green) has a positive slope over a very wide T₁ domain and the FLAWS-hco filter (orange) has a negative slope over the same very wide Domain. The magnitude of the slopes of the two filters are generally greater than that of the FLAWS-min/MP2RAGE filter. The monotonic form of the FLAWS-hc and FLAWS-hco T₁-filters is strikingly different from the bipolar form of dSIR and lSIR T₁-filters which have well defined, sharp discontinuities at their null points due to the use of magnitude reconstruction. They also have steep positive slopes within narrow mDs rather than the broad lower magnitude slopes of the FLAWS-hc and FLAWS-hco filters over wide T₁ Domains (Figure 26).

The FLAWS-min filter is the lower signal from two IR T₁-filters, one nulling white matter and the other nulling CSF. This is divided by the sum of the signals from the two filters and has a maximum value of a ½ and a minimum value of zero for the brain. It follows the shorter TI nulling IR filter in the lower part and the longer TI nulling IR filter in the higher part of the T₁ Domain. It is a positive unipolar filter with the lowest signal nulling white matter, a peak at a T₁ just greater than the T₁ of gray matter, and CSF as zero signal (Figure 26). It has similarities to the narrow mD dSIR sequence but has a lower display range and a lower slope. Müller et al described a halo sign when imaging MS lesions in the cortex with the FLAWS-min sequence (Figure 5 in reference [25]). This is due to an increase in cortex T₁ in lesions passing through the peak of the FLAWS-min filter into longer T₁ regions. The slope of the T₁-filter of the FLAWS-min sequence around its maximum value is relatively low and this may account for the greater width of the halo sign compared with the width of the bubble sign found with the narrow mD dSIR sequence when there is an increase in T₁ through the peak of that sequence (e.g. Figure 16B). There are also high signal boundaries between ventricular white matter and CSF as well as between white matter and a periventricular lesion with the FLAWS-min sequence (in Figure 5, reference [25]) as expected from its T₁-filter.

A similar pattern is seen in the leukocortical lesion shown in Figure 2B of reference [26] using T₁w-FLAWS which is the same as FLAWS-min. However, the cortical lesion shown in the same figure imaged with T₁w-FLAWS shows low contrast with the surrounding normal gray matter. Increased signal is seen at the boundary between gray matter and CSF. These features are consistent with Figure 26 (blue curve for FLAWS-min). FLAWS-min images do not show high signal between white matter and gray matter because the peak signal on the FLAWS-min T₁-filter is not between these two tissues. If the second nulling TI (of CSF) is reduced, the FLAWS-min sequence becomes more like a narrow mD dSIR filter but with lower magnitude slopes.

In summary, the MP2RAGE and FLAWS T₁-filters are unipolar or monotonic, not bipolar. They have fixed TIs and wide T₁ Domains with lower slopes than narrow mD dSIR T₁-bipolar filters in their mDs. The FLAWS-uni/MP2RAGE T₁-filter has a moderate slope, reduced CSF signal and shows white and gray matter at the upper and lower ends of the display scale. FLAWS-hc has a higher slope than FLAWS-uni/MP2RAGE, with CSF very high signal and a wide display range from WM to CSF. FLAWS-hco also has a higher magnitude negative slope than FLAWS-uni/MP2RAGE with CSF very low signal and a wide display range from CSF to white matter. FLAWS-min has a positive T₁-unipolar filter with a moderate slope and peak signal at a T₁ just greater than that of gray matter. The T₁-filters shown in Figure 26 provide a succinct and accessible way of comparing the signal levels and contrast performance of these sequences.

White Matter Nulled MP-RAGE [27] and Fast Gray Matter IR (FGATIR) [28]

These are single IR sequences that null white matter so that most of the remaining signal present in the brain is from gray matter. They have negative unipolar T₁-filters. The slopes of their T₁-filters are shown in Figure 6B and 6C when TI is chosen to null white matter. Small differences or changes in white matter from normal show much less contrast with white matter nulled MP-RAGE than with white matter nulled narrow mD dSIR sequences (Figure 6C). Boundaries between white matter and gray matter also show higher contrast with narrow mD sequences (Figure 7).

3D-Edge Enhancing Gradient Echo (3D-EDGE) [29] and EDGE-MP2RAGE [30] Sequences

The EDGE sequence has a single IR gradient echo with a negative unipolar T₁-filter and a TI chosen to produce zero signal at the boundary between white and gray matter with equal signals in each tissue. It uses magnitude reconstruction. The boundary contains voxels with a mixture of the two tissues and a T₁ between those of white and gray matter. It is of value in localising boundaries and diagnosing focal cortical dysplasia. EDGE-MP2RAGE acquires two separate images, the first of which is nulled in the same way as for the 3D-EDGE sequence. FGATIR, 3D-EDGE and EDGE-MP2RAGE images can be synthetically created from MP2RAGE derived T₁ maps [31].

Ethical Approval

All procedures performed in the studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards.

Informed consent

Written informed consent was obtained from all participants in this study.

Statement and Declarations

Conflict of Interest: GMB is a consultant to Magnetica, Brisbane, Australia. No other authors have any conflict of interest to disclose.

Funding

We have received support from the Fred Lewis Enterprise Foundation, the Hugh Green Foundation, Manaaki Moves, NZ P Pull, the JN & HB Williams Foundation, the Mangatawa Beale Williams Memorial Trust, Neurological Foundation, Trust Tairāwhiti, Friends of Mātai Blue Sky Fund, GE Healthcare, Mātai Ngā Māngai Maori, and Kānoa— Regional Economic Development & Investment Unit of New Zealand.

Abbreviations

BLAIR = BipoLAr fIlteR, bipolar filter

cdSIR = composite divided Subtracted Inversion Recovery

clSIR = composite logarithmic then Subtracted Inversion Recovery

dSIR = divided Subtracted Inversion Recovery

hD = highest Domain

IR = Inversion Recovery

lSIR = logarithmic then Subtracted Inversion Recovery

mD = middle Domain

MP2RAGE = Magnetisation Prepared 2 Rapid Acquisition Gradient Echo

MT = Magnetisation Transfer

SIR = Subtracted Inversion Recovery

S_OF = Signal from Other Filter

T₁D = T₁ Derivative

TP-bipolar filter = Tissue Property-bipolar filter

TP-filter = Tissue Property-filter

UHC = Ultra-High Contrast

References

Condron, P.; Cornfeld, D.M.; Scadeng, M.; Melzer, T.R.; Newburn, G.; Bydder, M.; Kwon, E.E.; McGeown, J.P.; Handsfield, G.G.; Emsden, T.; et al. Ultra-High Contrast MRI: The Whiteout Sign Shown with Divided Subtracted Inversion Recovery (dSIR) Sequences in Post-Insult Leukoencephalopathy Syndromes (PILS). Tomography 2024, 10, 983–1013. [Google Scholar] [CrossRef] [PubMed]
Ma, Y.-J.; Shao, H.; Fan, S.; Lu, X.; Du, J.; Young, I.R.; Bydder, G.M. New options for increasing the sensitivity, specificity and scope of synergistic contrast magnetic resonance imaging (scMRI) using Multiplied, Added, Subtracted and/or FiTted (MASTIR) pulse sequences. Quant. Imaging Med. Surg. 2020, 10, 2030–2065. [Google Scholar] [CrossRef]
Ma, Y.-J.; Moazamian, D.; Port, J.D.; Edjlali, M.; Pruvo, J.-P.; Hacein-Bey, L.; Hoggard, N.; Paley, M.N.J.; Menon, D.K.; Bonekamp, D.; et al. Targeted magnetic resonance imaging (tMRI) of small changes in the T1 and spatial properties of normal or near normal appearing white and gray matter in disease of the brain using divided subtracted inversion recovery (dSIR) and divided reverse subtracted inversion recovery (drSIR) sequences. Quant. Imaging Med. Surg. 2023, 13, 7304–7337. [Google Scholar] [CrossRef]
Newburn, G.; McGeown, J.P.; Kwon, E.; Tayebi, M.; Condron, P.; Emsden, T.; Holdsworth, S.J.; Cornfeld, D.M.; Bydder, G.M. Targeted MRI (tMRI) of Small Increases in the T1 of Normal Appearing White Matter in Mild Traumatic Brain Injury (mTBI) Using a Divided Subtracted Inversion Recovery (dSIR) Sequence. OBM Neurobiol. 2023, 07, 1–27. [Google Scholar] [CrossRef]
Ma Y-J, Moazamian D, Cornfeld DM, et al. Improving the understanding and performance of clinical MRI using tissue property filters and the central contrast theorem, MASDIR pulse sequences and synergistic contrast MRI. Quant Imaging Med Surg 2022, 12, 4658–4690. [Google Scholar] [CrossRef]
Cornfeld D, Condron P, Newburn G, et al Ultra-high contrast MRI: Using divided subtracted inversion recovery (dSIR) and divided echo subtraction (dES) sequences to study the brain and musculoskeletal system. Bioengineering 2024, 11, 441.
Newburn G, Condron P, Kwon EE, et al Diagnosis of delayed post-hypoxic leukoencephalopathy (Grinker’s myelinopathy) with MRI using divided subtracted inversion recovery (dSIR) sequences: time for reappraisal of the syndrome? Diagnostics 2024, 14, 418. [CrossRef]
Yokoo T, Bae WC, Hamilton G, et al A quantitative approach to sequence and image weighting. J Comput Assist Tomogr 2010, 34, 317–331. [CrossRef] [PubMed]
Young IR, Szeverenyi NM, Du J, Bydder GM Pulse sequences as tissue property filters (TP-filters): a way of understanding the signal, contrast and weighting of magnetic resonance images. Quant Imaging Med Surg 2020, 10, 1080–1120. [CrossRef]
Bydder, M. , Cornfeld, D. M., Melzer, T.R. et al Log subtracted inversion recovery. Magn Reson Imaging 2025, 117, 110328. [Google Scholar]
Ugurbil K Magnetic resonance imaging at ultrahigh fields. IEE Trans Biomed Eng 2014, 61, 1364–1379. [CrossRef] [PubMed]
Feinberg DA, Beckett AJS, Vu AT, et al. Next-generation MRI scanner designed for ultra-high-resolution human brain imaging at 7 Tesla. Nat Methods 2023, 20, 2048–2057. [Google Scholar] [CrossRef]
Bydder M, Condron P, Cornfeld DM, et al Validation of an ultra-high contrast divided subtracted inversion recovery technique using a standard T₁ phantom. NMR Biomed 2024, 2024, e5269.
Losa L, Peruzzo D, Galbiati S, Locatelli F, Agarwal N Enhancing T₁ signal of normal-appearing white matter with divided subtracted inversion recovery: a pilot study in mild traumatic brain injury. NMR Biomed 2024, 37, e5175.
Hajnal JV, Baudouin CJ, Oatridge A, Young IR, Bydder GM Design and implementation of magnetisation transfer pulse sequences for clinical use. J Comput Assist Tomogr 1992, 16, 7–18. [CrossRef]
Nürnberger L, Gracien R-M, Hok P, et al Longitudinal changes of cortical microstructure in Parkinson’s disease assessed with T₁ relaxometry. NeuroImage Clinical 2017, 13, 405–414. [CrossRef]
Roeben B, Zeltner L, Hagberg GE, Scheffler K, Schols L, Bender B Susceptibility-weighted imaging reveals subcortical iron deposition in PLA2G6- associated neurodegeneration: the “Double Cortex Sign”. Mov Disord 2023, 38, 904–906. [CrossRef]
Mueller SG 7T MP2RAGE for cortical myelin segmentation: impact of aging. PLoS One 2024, 19, e0299670.
Bydder GM, Young IR MR imaging: clinical use of the inversion recovery sequence. J Comput Assist Tomogr 1985, 9, 659–675. [CrossRef]
Redpath TW, Smith FW Technical note: use of a double inversion recovery pulse sequence to image selectively grey or white brain matter. Br J Radiol 1994, 67, 1258–1263. [CrossRef]
Tillema J-M, Weigand SD, Dayan M, et al Dark rims: novel sequence enhances diagnostic specificity in multiple sclerosis. AJNR Am J Neuroradiol 2018, 39, 1052–1058. [CrossRef]
Marques JP, Kober T, Krueger G, van der Zwaag W, Van de Moortele P-F, Gruetter R MP2RAGE, a self bias-field corrected sequence for improved segmentation and T₁-mapping at high field. Neuroimage 2010, 49, 1271–1281. [CrossRef]
Beaumont J, Saint-Jalmes H, Acosta O, et al Multi T₁-weighted contrast MRI with fluid and white matter suppression at 1. 5T. Magn Reson Imaging 2019, 63, 217–225. [Google Scholar]
Beaumont J, Gambarota G, Saint-Jalmes H, et al. High-resolution multi-T₁-weighted contrast and T₁ mapping with low B₁ sensitivity using the fluid and white matter suppression (FLAWS) sequence at 7T. Magn Reson Med 2021, 85, 1364–1378. [Google Scholar] [CrossRef]
Müller J, La Rosa F, Beaumont J, et al. Fluid and white matter suppression. New sensitive 3T magnetic resonance imaging contrasts for cortical lesion detection in multiple sclerosis. Invest Radiol 2022, 57, 592–600. [Google Scholar] [CrossRef] [PubMed]
Durozard P, Maarouf A, Zaaraoui W, et al. Cortical lesions as an early hallmark of Multiple Sclerosis: visualization by 7T MRI. Invest Radiol 2024, 59, 747–753. [Google Scholar] [CrossRef]
Saranathan M, Tourdias T, Bayram E, Ghanouni P, Rutt BK Optimization of white-matter-nulled magnetization prepared rapid gradient echo (MP-RAGE) imaging. Magn Reson Med 2015, 73, 1786–1794. [CrossRef] [PubMed]
Grewal SS, Middlebrooks EH, Kaufmann TJ, et al Fast gray matter acquisition T₁ inversion recovery MRI to delineate the mammillothalamic tract for preoperative direct targeting of the anterior nucleus of the thalamus for deep brain stimulation in epilepsy. Neurosurg Focus 2018, 45, E6. [CrossRef]
Middlebrooks EH, Lin C, Westerhold E, et al. Improved detection of focal cortical dysplasia using a novel 3D imaging sequence: Edge-Enhancing Gradient Echo (3D-EDGE) MRI. NeuroImage 2020, 28, 102449. [Google Scholar] [CrossRef] [PubMed]
Tao S, Zhou E, Greco V, et al. Edge-Enhancing Gradient-Echo MP2RAGE for clinical epilepsy imaging at 7T. AJNR Am J Neuroradiol 2023, 44, 268–270. [Google Scholar] [CrossRef]
Middlebrooks EH, Tao S, Zhou X, et al Synthetic inversion image generation using MP2RAGE T₁ mapping for surgical targeting in deep brain stimulation and lesioning. Stereotact Funct Neurosurg 2023, 101, 326–331. [CrossRef] [PubMed]

Figure 1. Positionally matched images of the brain in a 24-year-old male patient with mTBI (A and C) and a normal control (B). (A) is a T₂-FLAIR image in the patient which shows no abnormality. (B) is a narrow mD dSIR image of the brain in the normal control. The white matter in the central region of this image has a normal low signal (dark) appearance (arrows). (C) is a narrow mD dSIR image performed in the patient with the same sequence as in the control. This image shows all the patient’s white matter with abnormal high signal (light) (arrows) rather than the normal low signal (dark) appearance in the control (arrows) (B). There is a night and day difference in signal between normal and abnormal white matter in (B) and (C). High signal boundaries are seen between normal white matter and normal gray matter in (B). These are less obvious in (C) because of the high signal in the abnormal white matter.

Figure 2. Plot of M_Z/M_XY against time for a T₂-weighted version of the SE sequence for two tissues P (with a shorter T₁ and T₂) and Q (with a longer T₁ and T₂). Overall T₁ and T₂ dependent contrast between P and Q is shown with the positive blue arrow on the right.

Figure 3. IR T₁-filters with plots of signal S_T1 against T₁ in phase-sensitive (A) and magnitude forms (B). (A) is a low pass filter and (B) is a negative unipolar filter. (A) shows both positive and negative values for S_T1 whereas (B) shows negative values “reflected” across the X axis so they become positive. The maximum size of the slopes of the two T₁-filters which are tangential to the slope are shown as red lines. The slopes are negative in both cases. .

Figure 4. The long TR IR sequence. Negative T₁-unipolar filters for short TI_s (A, left), intermediate TI_i (B, centre) and long TI_l (C, right) including white (W) and gray (G) matter. Signal S_T1 is plotted against T₁ in each case (black curves). The positions of W and G are the same at each of the three TIs. TI is increased from TI_s (left) to TI_i (centre) and increased further to TI_l (right). In (A), (B) and (C), the increase in T₁ from W to G (horizontal positive green arrows) is multiplied by the relevant slopes of the filters (red lines) to produce contrast. This is moderately positive, strongly negative, and mildly negative contrast in (A), (B) and (C) respectively (vertical blue arrows).

Figure 5. Unipolar IR (A), subtracted IR (SIR) (B) and Added IR (AIR) (C) T₁-filters. T₁ is shown along the X axes in ms. (A) shows a TI_s T₁-unipolar filter (pink) and a TI_i T₁-unipolar filter (blue), (B) shows the subtraction (S_TIs - S_TIi) IR or SIR T₁-bipolar filter, and (C) shows the addition (S_TIs + S_TIi) IR or AIR T₁- filter. The vertical lines divide the X axes into lowest, middle and highest Domains (lD, mD and hD). The mD contains T₁ values between the vertical dashed lines. In (B) the slope of the curve of the SIR T₁-filter in the mD is about double that of the S_TIs T₁-filter (pink in [A]). In (C) the signal at T₁=0 is doubled to 2.0, and the signal in the mD is reduced to about 0.20 in the nearly linear, slightly downward sloping central part of the AIR T₁-filter (i.e., in the mD).

Figure 6. Figure 6A shows division of the SIR T₁-bipolar filter in Figure 5B by the addition T₁-filter AIR in Figure 5C to give the dSIR T₁-bipolar filter. T₁ is shown along the X axis in ms. Figure 6B shows comparison of the conventional IR S_TIs T₁-unipolar filter (pink) with the SIR T₁-bipolar filter (blue) for a small increase in T₁ (horizontal positive green arrow, DT₁). Figure 6C is a comparison of the S_TIs T₁-unipolar filter (pink) with the dSIR T₁-bipolar filter (blue) for the same small increase in T_1. In 6B the contrast produced by the SIR T₁-bipolar filter is about twice that produced by the IR T₁-unipolar filter (blue and pink arrows). In 6C the contrast produced by the dSIR T₁-bipolar filter is about ten times greater than that produced by the IR T₁-unipolar filter (blue and pink arrows). The display gray scale has no negative values so white matter appears on it as zero signal (black) at the lowest part of the display scale when using both the T₁-filters shown in each of 6B and 6C.

Figure 7. Boundaries. This image shows a narrow mD dSIR T₁-bipolar filter with a mD extending from white matter (W) to a T_1W,G between the TIs of W and gray matter (G) (blue curve) as well as a white matter nulled conventional IR T₁-unipolar filter e.g., MP-RAGE (pink curve). The X axis is shown in ms. The contrast produced from the difference in signal between W and G by the T₁-unipolar filter is (S_G minus S_W, pink curve) and is shown by the vertical pink arrow. With the T₁-bipolar filter (blue curve) there is a partial volume effect between W and G producing a T_1W,G between the T1s of W and G which results in the high signal, SW,G shown on the T₁-bipolar filter (blue curve). The contrast between this high signal and white matter is the difference (S_W,G minus S_W, in blue) and is shown by the vertical blue arrow. The contrast produced by the T₁-bipolar filter at the boundary between W and G is far greater than that produced at the boundary between W and G by the T₁-unipolar filter.

Figure 8. Plots of ln S_TIs (TI_s = 350 ms) vs T₁ (pink), and ln S_TIi (TI_i = 500 ms) vs T₁ (blue). The plots have a negative unipolar form with steeply sloping signal gradients around their negative poles. The pink and blue ln S_TIs and ln S_TIi curves are steeper than the corresponding magnitude curves shown in Figure 5A. The lSIR filter S_lSIR = ½(ln S_TIs - ln S_TIi) reverses the sign of the S_TIi filter so it becomes positive as seen in the following figure (Figure 9) (orange curve).

Figure 9. Plot of signal vs T₁ in ms for a dSIR T₁-bipolar filter (blue curve) and the corresponding lSIR T₁-bipolar filter with the same null points (orange curve) (TI_s = 350 ms and TI_i = 500 ms, DTI = 150 ms). The dSIR filter shows the usual pattern with maximum values of ±1. It has an essentially constant slope in the mD. The lSIR filter has similar values to the dSIR filter in the regions of lowest and highest values of T₁ as well as in the centre of the mD. However around the lower and upper T₁ nulling values it has much steeper slopes and proceeds asymptotically to minus and plus infinity respectively (values of S along the Y axis in the range of ± 2 are shown). From a practical point of view, the contrast for very small differences or changes in T₁ is 2-3 times greater with the lSIR filter compared with the corresponding dSIR filter. For the same difference in T₁ in the regions close to the null points the lSIR filter is steeper than the dSIR filter. This results in greater contrast and spatial resolution.

Figure 10. Composite dSIR filters (cdSIRs with TI_s = 350 ms and TI_i = 500 ms). Attenuated first IR T₁-filter (TI_s) in (A) and attenuated second IR T₁-filter (TI_i) in (B). (A) shows the conventional T₁-bipolar filter (blue) with S_OF = 1, an attenuated first TI_s filter (S_OF = 0.5) (red) and a further attenuated first TI_s filter (S_OF = 0.25) (yellow).

Figure 11. Comparison of conventional windowing (A and B) and a dSIR T₁-bipolar filter (C). (A) plots the display signal S_D against image signal S for two tissues P and Q over the full image gray scale range shown along the Y axis. There is a linear relationship between S_D and S over the full gray scale range. In (B) the image is narrowly windowed to place the signals from P and Q at the lower and upper ends of the S_D range shown on the Y axis. The contrast between P and Q (difference in signal level shown on the Y axis) is increased in (B) compared with (A). In (B) signal values less than P on the X axis are all shown at the lowest signal level on the Y axis, and signal values along the X axis greater than Q are all shown at the highest signal level on the Y axis (horizontal lines). In these regions, voxels have the same signals and show no contrast between them. They form isointense blocks and anatomical detail is lost within the blocks. (C) plots signal against T₁ using a dSIR T₁-bipolar filter. Contrast between P and Q is high in the mD as in (B). Values of T₁ less than P and greater than Q do not form plateaux as in (B) but follow the T₁-bipolar filter in the lD and hD. There are therefore differences in signal between voxels so contrast is maintained and anatomical detail is preserved outside of the mD bounded by P and Q.

Figure 12. 18-year-old patient with mTBI 21h (A) and 64h (B) post injury imaged with the same narrow mD dSIR sequence. In (A) the patient shows a whiteout sign (grade 4 out of 5) with high signal in most of the white matter in the cerebral hemispheres except for the anterior and posterior central corpus callosum. The posterior limb of the internal capsule is high signal. The thalami show low internal contrast from medial to lateral (arrows on lateral margins of the thalami). There is also low external contrast between the lateral margins of the thalami and the adjacent posterior limbs of the internal capsule. On the follow up examination at 64h (B) the whiteout sign has resolved and there is low signal in the white matter including the posterior limbs of the internal capsule (except for the corticospinal tracts). The thalami now show high internal contrast from medial to lateral (arrows on lateral margins of the thalamus) which is the normal appearance at this age. There is now also very high external contrast between the lateral margin of the thalamus and the adjacent posterior limb of the internal capsule. Image (A) shows the grayout sign which is a reduction in the high contrast between the medial and lateral gray matter of the thalamus. The high contrast is restored in (B). No abnormality was seen on the corresponding T₂-FLAIR images obtained at both time points.

Figure 13. Normal control (A) and patient with mTBI (B) showing a whiteout sign and grayout signs (narrow mD dSIR images). The normal control shows the heads of the caudate nuclei with higher signal than the adjacent CSF. Contrast is also seen between the cortex and CSF. In (B), the patient shows a whiteout sign. There are grayout signs in the thalami and putamina. In addition, contrast between the heads of the caudate nuclei and CSF is lost and there is little or no contrast between cortex and CSF which are isointense. These are also grayout signs. No abnormality was seen on the T₂-FLAIR images in the normal control or patient.

Figure 14. Red nuclei in an 18-year-old normal control (A) and an 18-year-old male patient with mTBI (B). 2D narrow mD dSIR images. In (A) the normal control shows low signal in the white matter of the cerebral hemisphere, cortical spinal tracts and the ascending sensory tracts. The red nuclei (arrows) have an intermediate mid-gray signal. In (B) the patient shows high signal in the cerebral white matter, the corticospinal tracts and the ascending sensory tracts (whiteout sign, grade 4 out of 5). In addition, the red nuclei are higher signal than in (A) (arrows). No abnormality was seen on the T₂-FLAIR images in the normal control or patient.

Figure 15. 41-year-old female patient with MS in remission. 2D T₂-wSE (A) and narrow mD dSIR (B) images at the same level. No abnormality is seen in (A). A focal lesion is seen in (B) (long narrow white arrow) and the corticospinal tracts show a high signal (short narrow white arrows). In addition there is patchy increased signal in white matter (short thick white arrows) with only a small region showing a normal or near normal low signal appearance (long black arrow). High contrast and high spatial resolution contrast are seen at the boundaries between normal white matter and normal gray matter in (B). These features are less obvious in areas where the white matter has abnormal high signal.

Figure 16. 32-year-old female with MS during a relapse. T₂-FLAIR (A), synthetic narrow mD dSIR (B) and filtered gradient echo (C) images. On the T₂-FLAIR image (A), one lesion is seen (long arrow). The surrounding white matter appears normal. On the dSIR image (B), the lesion shown on the T₂-FLAIR image is seen (long arrow) as well as six other lesions (short arrows). High signal boundaries around lesions are also seen in (B). Some of these lesions show paramagnetic rims on the filtered gradient echo image (arrows) in (C). In addition, most of the white matter in (B) is high signal corresponding to a high grade 4/5 whiteout sign [1]. The whiteout sign is not seen on the T₂-FLAIR image (A).

Figure 17. 24-year-old female patient with MS in remission (TI_s = 350 ms, TI_i = 500 ms). Narrow mD dSIR (A) and composite (T₁ and T₂) cdSIR (TI_s/TE_s = 350/7 and 80 ms; TI_i/TE_i = 500/7 ms) (B) images. A leukocortical lesion is seen on the dSIR image (A) and with higher contrast on the cdSIR image in (B) (arrows). White matter and gray matter boundaries are higher contrast and narrower on the composite filter cdSIR image (B). Also white matter is more uniformly low signal in (B). These features are consistent with the attenuation of the first IR TI_s filter signal increasing contrast and narrowing boundaries at the positive pole and decreasing these at the negative pole as shown in Figure 10A.

Figure 18. 76-year-old female patient in remission with a diagnosis of MS. Sagittal 3D T₂-FLAIR (A) and 3D wide mD dSIR (B) images. The T₂-FLAIR image shows a poorly defined area of increased signal in the cervical cord (arrow). The dSIR image shows a high contrast lesion with sharply defined (“punched out”) boundaries in the cervical cord (lower three arrows). This is much more extensive than in (A). An additional lesion is seen in the medulla on the dSIR image in the region of the area postrema (highest arrow) (B) but not on the T₂-FLAIR image (A). The extended lesion in the cervical cord and the lesion in the medulla raise the possibility of neuromyelitis spectrum disorder. Other conventional sequences such as MP-RAGE, T₂-wSE or STIR may perform better than T₂-FLAIR in the cervical cord.

Figure 19. 51-year-old male patient with methamphetamine substance use disorder after one month’s abstinence (A) and after nine months’ abstinence (B). Matched narrow mD dSIR images. In (A) there is extensive high signal in the white matter of the cerebral hemispheres with only small areas of normal or near normal white matter (dark) at the periphery of the white matter (whiteout sign, grade 4 out of 5). After nine months’ abstinence (B), the high signal areas in (A) have markedly regressed. There is some intermediate signal in the more central white matter and lower signal in the peripheral white matter (whiteout sign grade 2 out of 5, where grade 1 is normal). No abnormality was seen in either examination on the corresponding T₂-FLAIR images. A previous mTBI in the patient may have been a confounder.

Figure 20. 72-year-old male patient with Parkinson’s disease. 2D narrow mD dSIR image (A) with insets of anterior left cortex (B) and thalami (C). There is high signal in the superficial layer of the cortex at multiple sites (bilaminar cortex signs) which are more obvious on the inset (white arrows) (B). There are circular appearances in the thalami, putamen and heads of the caudate nuclei (bubble signs) which are more obvious on the inset (e.g., black arrows) (C). The bubble sign is due to focal reductions in T₁ probably from the presence of free iron. The thalamic signal is slightly higher medially compared with laterally unlike the appearance of the normal thalamus in the 18-year-old patient shown in Figure 12B.

Figure 21. 65-year-old female with a glioma (not shown) post chemotherapy. Parasagittal T₂-FLAIR (A) and synthetic narrow mD dSIR (B) images. On the T₂-FLAIR image no abnormality is seen. On the dSIR image (B), the white matter has a high signal corresponding to a high grade whiteout sign. This includes the cerebellar white matter which is similar in intensity to the cerebral white matter. There is some sparing of the peripheral white matter in (B).

Figure 22. 46-year-old normal control. dSIR (A) and lSIR (B) images with insets of anterior left cortex dSIR in (C) and anterior left cortex lSIR in (D). The boundaries between white matter and gray matter are seen with higher contrast and the contrast is seen with higher spatial resolution in (B) and (D) compared with (A) and (C). Bilaminar cortex signs are also seen with higher contrast and spatial resolution in (B) and (D) (white arrows). There are also small bubble signs in the putamen and medial thalamus which are better seen on the lSIR image (B) compared with the dSIR image (A).

Figure 23. 41-year-old female patient with MS. A leukocortical lesion in the right medial frontal region is shown on the narrow mD dSIR image (A) and a matching lSIR image (B) (arrows). There are also left to right profiles with signal plotted against position (in mm) for the dSIR (blue) and lSIR (orange) images (C) at the level of the horizontal arrows shown in (A) and (B). No boundary between white matter and gray matter is seen within the lesion in (A). A disrupted high signal boundary between white matter and gray matter is seen in the lesion in (B). The lSIR profile (orange) has higher signal and steeper slopes than the dSIR profile (blue) in (C). The difference in signal (or contrast) achieved for the same change in position is generally greater with the lSIR filter i.e., the contrast shown on the lSIR image generally has a higher spatial resolution.

Figure 24. Thalamus. dSIR narrow mD (TI_s = 350 ms and TI_i = 500 ms). T₁-bipolar filter for normal appearance, abnormal increase in T₁ in lateral thalamic gray matter within the hD with loss of contrast with medial thalamic gray matter (striped green arrow), and subsequent decrease in T₁ to normal with restoration of contrast with medial thalamic gray matter (solid green arrow). The normal thalamus gray matter shows a decrease in T₁ from medial (TG_med) to lateral (TG_lat) (highest horizontal negative solid green arrow) and corresponding positive contrast (first vertical positive blue arrow). With disease there is an increase in T₁ in the lateral thalamus (middle horizontal positive striped green arrow) with a loss of contrast with the medial thalamus (blue circle). With later regression of this and return of the T₁ difference between the lateral and medial thalamus to normal (lowest negative solid green arrow), normal positive contrast is seen between the medial and lateral thalamus (second vertical positive blue arrow).

Figure 25. Thalamus. Narrow mD (TI_s = 350 ms and TI_i = 500 ms). dSIR T₁-bipolar filter for a larger abnormal decrease in T₁ in lateral thalamic gray matter from normal (TG_latnorm) to abnormal (TG_latab) extending from the hD through the signal peak to the mD (horizontal negative green arrow). Signal in the lateral thalamus is highest where there is partial volume effect between abnormal and normal gray matter of the thalamus (TG_latab/TG_latnorm) (longer vertical blue arrow). Signal from abnormal lateral thalamus in the middle Domain (TG_latab) (shorter blue arrow) is less than the signal in the lateral thalamus at the boundary between normal and abnormal thalamus (longer blue arrow). The higher and lower signals (vertical blue arrows) result in a bubble sign with higher signal in the periphery where there are partial volume effects and lower signal centrally where there is only abnormal gray matter. The same argument applies to the medial thalamus.

Figure 26. Plots of signal vs T₁ (T₁-filters) for the FLAWS-uni/MP2RAGE (blue), FLAWS-min (yellow), FLAWS-hc (green) and FLAWS-hco (orange) T₁-filters from Beaumont et al Figure 1B in reference [24]. The FLAWS-uni/MP2RAGE T₁-filter (blue) is a negative unipolar filter for brain and CSF with a moderate negative slope between white matter (WM) and gray matter (GM), and a moderate positive slope between GM and CSF at 4000 ms. The CSF and WM signals are similar. The display range is from GM to WM. (GP is the globus pallidus).The FLAWS-min T₁-filter (yellow) is a positive unipolar filter with a moderate positive slope from WM to beyond GM and a moderative negative slope from this point to CSF with CSF low signal. The display range is from WM and CSF to just higher than GM.The FLAWS-hc T₁-filter (green) is monotonic for brain and CSF. It has a higher positive slope than the magnitude of the FLAWS-min/MP2RAGE slope from WM to GM and CSF has a very high signal. WM is at the lowest level on the display range and CSF is at the highest level. The FLAWS-hco T₁-filter is monotonic with a higher size of slope than FLAWS-uni/MP2RAGE between WM and GM and a very low signal for CSF. The display range is from CSF to white matter.(Reproduced with permission: from Beaumont et al [24]).

Table 1. Pulse sequences and pulse sequence parameters used at 3T. Only 2D scans were performed on the Siemens Healthineers Skyra system.

#	Sequence	TR (ms)	TI (ms)	TE (ms)	Matrix size Voxel size (mm)	Number of slices	Slice thickness (mm)
1*	2D FSE IR (for white nulling matter)	9,192	350	7, 80	256 x 224 0.9 x 0.1 Z512 0.4 x 0.4	26	4
2*	2D FSE IR (used with #1 for narrow mD dSIR)	5,796	500	7, 80	256 x 224 0.9 x 0.1 Z512 0.4 x 0.4	26	4
3*	2D FSE IR (for longer T₁ nulling)	5,796	600	7	256 x 224 0.9 x 0.1 Z512 0.4 x 0.4	26	4
4*	2D FSE IR (used with #1 for wide mD dSIR)	5,796	800	7	256 x 224 0.9 x 0.1 Z512 0.4 x 0.4	26	4
5	3D BRAVO (for white matter nulling)	2,000	400	6	256 x 256 0.8 x 0.8 Z512	220	0.8
6	3D BRA VO (used with #5 for wide mD dSIR)	2,000	800	6	256 x 256 0.8 x 0.8 Z512	220	0.8
7*	2D T₂-FLAIR	6,300	1,851	102	320 x 240 0.7 x 0.7 Z512 0.4 x 0.4	26	4
8	3D T₂-FLAIR	6,300	1,850	102	256 x 256 0.8 x 0.8 Z512 0.6 x 0.6	252	0.8
9	3D susceptibility weighted	40	-	32	300 x 300 0.8 x 0.8 Z512	110	2

Z = zipped * = positionally matched for slice position, slice thickness, slice number, field of view and interpolated matrix size.

Table 2. Directly acquired tissue property bipolar filters.

Bipolar filter	Reverse bipolar filter	Tissue properties single (upper) composite (lower)
SIR	rSIR	T₁
dSIR cdSIR	drSIR cdrSIR	T₁ T₁; T₂, T₂, D*, c and/or MT
lSIR clSIR	lrSIR clrSIR	T₁ T₁; T₂, T₂, D, c and/or MT

Table 3. Signal equations for bipolar filters and other functions.

#	Filter, other functions	Signal Equation	Fig. #
1	IR, TI_s	S_TIs = 1-2e^-TIs/T₁	3,4,5
2	IR, TI_i	S_TIi = 1-2e^-TIi/T₁	3,4,5
3	SIR	S_SIR = S_TIs - S_TIi	5
4	dSIR		6,9
5	dSIR	(in lD and hD) ΔTI = TI_i – TI_s (in mD)	6,9
6	cdSIR		10
7	cdSIR, S_OF	S_OF = ±e^-DTE/T2, ± e^-DTE/T2, ±e^-DbD, etc	10
8	dSIR, S_dSIR	$S_{dSIR} \approx \frac{\ln 4}{∆ T I} T_{1} - \frac{Σ T I}{∆ T I}$ (in mD)	6
9	dSIR, T₁	$T_{1} \approx \frac{∆ T I}{l n 4} S_{dSIR} + \frac{Σ T I}{l n 4}$ (in mD)	6
10	lSIR	S_lSIR = ½(ln S_TIs - ln S_TIi)	8,9
11	clSIR
12	clSIR, lSIR	, ± , ±DbD*, etc
13	lSIR, dSIR†	S_lSIR = atanh S_dSIR (in mD)	9
14	dSIR, T₁D††	(in mD)	6,9
15	lSIR, T₁D †††	(in mD)	9

† S_lSIR is the inverse hyperbolic tangent of S_dSIR in the mD. †† is the T₁ derivative (T₁D) of SdSIR and gives the slope of dSIR filter in the mD which is essentially constant and equal to (Figure 9). ††† is the T₁ derivative (T₁D) of the S_lSIR filter in the mD which varies. It has values of minus or plus infinity when and respectively (Figure 9).

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.