A comparative analysis of the sensitivity and BOLD contamination of the VASO response at 3 Tesla: ME-DEPICTING vs. ME-EPI readouts

Ratnamanjuri Devi; Jöran Lepsien; Toralf Mildner; Harald E. Möller

PMC · DOI:10.1162/imag_a_00333·October 28, 2024

A comparative analysis of the sensitivity and BOLD contamination of the VASO response at 3 Tesla: ME-DEPICTING vs. ME-EPI readouts

Ratnamanjuri Devi, Jöran Lepsien, Toralf Mildner, Harald E. Möller

PDF

Open Access

TL;DR

The study compares two MRI techniques for measuring blood volume changes in the brain, finding that signal corrections are crucial to avoid overestimation.

Contribution

The paper introduces a novel comparison of ME-DEPICTING and ME-EPI for VASO signal correction and sensitivity at 3 Tesla.

Findings

01

Uncorrected VASO data from the first echo showed significant differences in sensitivity between the two readouts.

02

ME extrapolation corrected VASO data showed similar sensitivity for detecting cerebral blood volume changes.

03

Dynamic division correction led to a near-perfect linear dependence on echo time due to intravascular BOLD contributions.

Abstract

‘Non-BOLD fMRI’ data acquired at non-zero echo time (TE) suffer from contamination by the Blood Oxygenation Level Dependent (BOLD) signal due to the unavoidable signal decay caused by transverse relaxation. This contamination further reduces their already low inherent functional sensitivities and makes their correction essential. The Slice-Saturation Slab-Inversion Vascular Space Occupancy (SS-SI–VASO), for instance, cancels out BOLD contributions from VASO data, reflecting cerebral blood volume (CBV) changes, via a dynamic division approach. Alternatively, multi-echo (ME) data provide the possibility of extrapolating toTE=0. Acquisitions at very shortTEwould minimize the need for such corrections. The center-out EPI variant (‘DEPICTING’) is one such readout which allows for shortTE. The ME 2D DEPICTING was compared here against a traditional ME 2D EPI for its sensitivity to functional…

Linked entities

Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.

Chemicals1

Blood Oxygenation

Figures11

Click any figure to enlarge with its caption.

Fig. 1 — SS-SI-VASO implementation with ME-DEPICTING and ME-EPI as imaging readouts. For the slice package centered at the calcarine sulcus, the blood nulling condition at 3 T was fulfilled for the 6th slice (in orange) of the VASO images. The interleaved acquisition of VASO and Control images was accomplished with both ME-DEPICTING and ME-EPI readouts forn= 3 echoes. The sequence diagrams of the two readouts are provided. The double-shot nature of DEPICTING results in a shorter TE when compared to ME-EPI.

Fig. 2 — Simultaneously acquired VASO and BOLD responses in participant P1. (A) VASO activation maps obtained from the BOLD-uncorrectedSbn(TE1), BOLD-correctedSbn(TE→0), andSdd(TE1). Note the larger activation areas obtained with ME-DEPICTING. (B) BOLD activation maps obtained fromSctr(sum)and fittedR2,ctr*data. TheδsVASO(%),δsBOLD(%), andΔR2*(s–1) in the activation maps, shown here for 3 slices, have been overlaid over the corresponding slices of the pre-processedSbn(TE→0)of the two readouts.

Fig. 3 — Functional sensitivity comparison of the two readouts (group-averaged data). Results obtained with the ME-DEPICTING readout are presented in orange and those with the ME-EPI readout in blue. Differences between the two readouts are given in % (orange and blue text boxes) and their significance is marked with asterisks (*, , and * forp<0.05,0.01, and0.001, respectively). Red and blue arrows indicate sensitivity gains for ME-DEPICTING and ME-EPI, respectively. Error bars denote standard deviations across participants. (A) Sensitivity for VASO signal detection: At minimumTE,DEPICTING yielded 72% more voxels than EPI. Correction for BOLD contamination withSbn(TE→0)led to 30% and approx. 100% increases in the activated area for ME-DEPICTING and ME-EPI, respectively. Dynamic division yielded unexpectedly greater boosts of 78% and 300%, respectively. (B) VASO CNR in commonly activated regions of both readouts:Sbn(TE1)revealed a 27% higher CNR for DEPICTING compared to EPI (1.4 ± 0.3 vs. 1.1 ± 0.2). The CNR values obtained from the dynamically divided data were once again higher than those of the extrapolated data, yielding gains of 35% and 100% for ME-DEPICTING and ME-EPI, respectively. (C) BOLD CNRs based onSctr(sum)andR2,ctr*: The CNR of ME-EPI was only a 12% higher than that of ME-DEPICTING. (D)Sctr(TEi) showing the expected increase of the BOLD sensitivity with TE for both readouts.

Fig. 4 — VASO activation maps obtained from the nulled slice in two participants (P1 and P5). The maps based on theSbn(TE1),Sbn(TE→0)andSdd(TEi) data show BOLD-uncorrectedδsVASO,bn(TE1)as well as BOLD-corrected results obtained by extrapolation to zeroTE,δsVASO,bn(TE→0), and by dynamic division,δsVASO,dd(TEi)for all echo times. The activation maps have been overlaid over the 6th slice of the respective pre-processed VASO images at the specified TE. For both BOLD-correction strategies, the increase of the size and strength of activation compared to uncorrected data is evident. An increase in the intensity of the VASO signal change,δsVASO,ddwith TE is indicated by the growing blue/violet area for theSdd(TE2) andSdd(TE3) data of both readouts.

Fig. 5 — Cycle-averaged time courses of percent VASO signal changes. (A) ME-DEPICTING results, (B) ME-EPI results; both obtained in participant P1 (nulled slice; see top panel inFig. 4). Red squares: BOLD response from the control data,δsBOLD,ctr(TEi); open green squares: uncorrected VASO response,δsVASO,bn(TEi); filled green squares: VASO response corrected by dynamic division,δsVASO,dd(TEi). Yellow triangles represent the difference [δsVASO,bn(TEi)−δsBOLD,ctr(TEi)]. The yellow-shaded region denotes the stimulus duration.

Fig. 6 — Echo-time dependence of dynamic-division corrected VASO signals. MeanδsVASO,dd(A, D), calculated over the common ROI of VASO activation in the nulled slices of ME-DEPICTING (A) and ME-EPI (D), and the correspondingδcbv obtained withEq. 2(B, E) are plotted for all individual subjects against the correspondingTEi.δsBOLD,ctr(C, F) serves as references for the TE dependence. The subject-averaged values are shown in red (also seeSupplementary Tables S5 and S6).

Fig. 7 — Linear relationships of subject-averaged VASO and BOLD responses withTE. Subject-averagedδsVASO,dd(A, C) andδsBOLD, ctr(B, D) were obtained fromSdd(TEi)andSctr(TEi)recorded with ME-DEPICTING (A, B) and ME-EPI (C, D). The linear relations and their corresponding goodness of fit,r2, are given within the rectangular boxes. For comparison purposes, the subject-averaged δsVASO,bn(TE→0)values are indicated in (A) and (C) by black diamonds on they-axes. Residual BOLD signals extrapolated to zeroTEwere 0.02±0.14% and 0.07±0.03% for EPI and DEPICTING, respectively. These residual BOLD-like time-dependent signals and the correspondingδsVASO,bn(TE→0)signals (for reference) for each individual participant are shown inSupplementary Figure S2.

Fig. 8 — SS-SI-VASO acquisition with ME-EPI with 5 echoes. (A) Significant activation observed in the 6th slice (thresholded atp< 10–6) of the uncorrected VASO response fromSbn(TE1), the BOLD-corrected VASO responses from the extrapolatedSbn(TE→0), and the dynamically divided echoesSdd(TEi). The maps have been overlaid over the respective pre-processed VASO images at the specified TE. (B) Corresponding uncorrected (δsVASO,bn(TE1))and corrected cycle-averaged time courses ofδsVASO,bn(TE→0)andδsVASO,dd(TEi)for 1≤i≤5. The shaded yellow region denotes stimulus duration.

Fig. 9 — ME-DEPICTING SS-SI-VASO acquisition at higher resolution. (A) Significant VASO activation (p< 10–4) as seen on the nulled slice (8th slice) of the uncorrectedSbn(TE1)(2nd column, within the black square box), the ME-extrapolatedSbn(TE→0)(1st column), and the dynamically dividedSdd(TEi)(columns 3 to 5), overlaid over the corresponding unprocessed images. (B) Comparison of cycle-averaged time courses of VASO signal changes at different spatial resolutions:δsVASO,bn(TE1)(in green),δsVASO,bn(TE→0)(in blue), andδsVASO,dd(TE1)(gray broken line). The signals were extracted from the region defined by activation in theSbn(TE→0)data. Shaded yellow region denotes stimulus duration.

Fig. 10 — T2*,rest,ΔR2*andfEVmaps for participant P16 (6th slice, common ROI-1). Voxel-wise results for (A)T2,EV+IV*,restandT2,EV*,rest, (B)ΔR2,EV+IV*andΔR2,EV*, and (C) the extravascularΔR2*fraction,fEVobtained with ME-DEPICTING (top) and ME-EPI (bottom). The maps have been overlaid over the corresponding structural slice.

Fig. 11 — Dependence ofΔR2,IV*onT2,IV*,rest. (A) Results obtained with subject-averaged input values for ME-DEPICTING (orange) and ME-EPI (blue). (B) High-resolution versus low-resolution data obtained in a single participant (P1) with ME-DEPICTING. (C) Single-subject results (P1) obtained with ME-EPI at different TE. The shaded green area represents the venous range withT2,IV*,rest≈ 15–30 ms (Donahue et al., 2011).

Tables3

Table 1.. Values ofT2*,restandΔR2*in the EV+IV and EV compartments (obtained from time series ofR2,ctr*andR2,bn*, respectively) and the corresponding extravascularΔR2*fraction. Results obtained with both readouts are shown separately for ROI-1 and ROI-2 located in the blood-nulled slice.

Participant	ROI size	Parenchyma				GM
		$T_{2, EV+IV}^{*, rest}$ (ms)		$Δ R_{2, EV+IV}^{*}$ (s ^–1 )		$T_{2, EV}^{*, rest}$ (ms)		$Δ R_{2, EV}^{*}$ (s ^–1 )		Extravascular $Δ R_{2}^{*}$ fraction $f_{EV}$ (%)
		EPI	DEPICTING	EPI	DEPICTING	EPI	DEPICTING	EPI	DEPICTING	EPI	DEPICTING
ROI-1: Slice 6 of functional parenchymal visual cortex mask
P1	270	41.9	36.3	–0.69	–0.83	40.9	33.9	–0.35	–0.45	50.7	54.3
P3	195	41.2	40.1	–0.69	–0.96	39.7	38.6	–0.45	–0.53	71.2	56.1
P4	98	40.5	38.1	–0.81	–0.91	39.7	36.3	–0.50	–0.56	61.0	61.3
P5	141	48.5	42.4	–0.53	–0.74	46.4	39.8	–0.18	–0.34	33.0	44.6
P9	157	47.5	42.6	–0.46	–0.67	45.8	39.7	–0.36	–0.46	78.3	68.7
P10	130	42.4	36.2	–0.40	–0.78	42.3	34.9	–0.21	–0.44	53.6	53.4
P11	55	47.5	40.1	–0.52	–0.77	49.0	36.7	–0.41	–0.45	78.4	57.1
P12	163	52.1	46.5	–0.68	–0.58	49.2	42.0	–0.28	–0.15	42.6	25.9
P13	158	48.9	46.6	–0.50	–0.52	47.5	43.2	–0.25	–0.24	51.0	42.0
P14	75	47.2	45.0	–0.47	–0.57	45.5	42.7	–0.17	–0.26	38.5	44.6
P16	164	49.7	43.6	–0.45	–0.51	48.5	40.9	–0.23	–0.25	51.2	45.6
Mean ± SD	146 ± 59	46.1 ± 3.9	41.6 ± 3.7	–0.56 ± 0.13	–0.71 ± 0.16	45 ± 3.7	39 ± 3.2	–0.31 ± 0.11	–0.38 ± 0.13	55 ± 15	50 ± 12
ROI-2: Slice 6 of functional extravascular visual cortex mask
P1	71	41.5	33.4	–0.86	–1.04	41.4	31.4	–0.49	–0.62	59.4	60.7
P3	97	41.5	42.3	–0.87	–1.16	39.3	39.9	–0.58	–0.73	70.1	63.5
P4	74	41.7	40.6	–0.88	–0.98	41.2	39.2	–0.58	–0.67	68.7	60.4
P5	17	50.1	48.1	–0.71	–0.86	46.5	45.8	–0.35	–0.45	49.5	54.0
P9	73	47.2	40.4	–0.57	–0.76	45.6	38.2	–0.46	–0.61	83.5	83.2
P10	30	36.6	35.1	–0.61	–1.14	37.6	35.2	–0.37	–0.74	60.2	62.6
P11	39	45.6	39.9	–0.59	–0.91	48.2	37.0	–0.47	–0.56	78.4	57.1
P12	8	48.8	37.1	–1.39	–1.20	48.2	36.0	–0.50	–0.49	35.1	40.6
P13	45	46.7	42.9	–0.89	–0.81	45.9	39.4	–0.45	–0.43	51.0	52.2
P14	–	–	–	–	–	–	–	–	–	–	–
P16	43	51.7	46.9	–0.55	–0.66	49.9	44.0	–0.36	–0.40	65.5	61.8
Mean ± SD	45 ± 31	45.1 ± 4.7	40.7 ± 4.7	–0.79 ± 0.25	-0.95 ± 0.18	44 ± 4.2	38.6 ± 4.2	–0.46 ± 0.08	–0.57 ± 0.12	62 ± 14	60 ± 11

Table 2.. Summary of results from the EV+IV and EV compartments forT2*,restandΔR2*, and the corresponding extravascularΔR2*fraction,fEV, reported in literature for the visual cortex and from the current work.

Publication	Methodology	TE (ms)	Sample size	Voxel size (mm ³ )	ROI	EV+IV		EV
Publication	Methodology	TE (ms)	Sample size	Voxel size (mm ³ )	ROI	$T_{2}^{*, rest}$ (ms)	$Δ R_{2}^{*}$ (s ^–1 )	$T_{2}^{*, rest}$ (ms)	$Δ R_{2}^{*}$ (s ^–1 )	Extravasc. $Δ R_{2}^{*}$ fraction $f_{EV}$ (%)
Lu and van Zijl (2005)	ME-VASO	14, 34.5, 55, 75.9	4	2 × 2 × 5	Single slice	45.5 ± 3.5	–0.58 ± 0.18	47.4 ± 3.1	–0.38 ± 0.10	67 ± 12
Donahue et al. (2011)	ME GRE-BOLD EPI with/without bipolar crushers	32.7, 44.6, 57.6, 70.7	7	3.5 × 3.5 × 3.5	Occipitallobe		–0.74 ± 0.13		–0.52 ± 0.19	70 ± 29
This work	SS-SI-VA	1.7, 10.7, 19.7	11	3 × 3 × 4	ROI-1	41.6 ± 3.7	–0.71 ± 0.16	39.0 ± 3.2	–0.38 ± 0.13	50 ± 12
	SO ME-DEPICTING				ROI-2	40.7 ± 4.7	–0.95 ± 0.18	38.6 ± 4.2	0.57 ± 0.12	60 ± 11
			1	1.7 × 1.7 × 1.7	ROI-1	47.5	–1.03	47.2	-0.80	80
					ROI-2	49.3	–1.13	47.2	-1.0	95
	SS-SI-VA	7.5, 20.7, 33.9	11	3 × 3 × 4	ROI-1	46.1 ± 3.9	–0.56 ± 0.13	45.0 ± 3.7	–0.31 ± 0.11	55 ± 15
	SO ME-EPI				ROI-2	45.1 ± 4.7	–0.79 ± 0.25	44.0 ± 4.2	–0.46 ± 0.08	62 ± 14
		10, 25.3, 40.6, 55.9, 71.2	1		ROI-1	43.1	–0.54	42.6	–0.31	59
					ROI-2	43.2	–0.59	42.8	–0.37	67

Table 3.. List of experimental values used to estimate the intravascularΔR2,IV*for the two readouts based on the results in ROI-2.δsVASOvalues were extracted from the dynamically divided data at the respective echo time.

	TE _mid (ms)	${CBV}^{act}$ (ml/ml)	$δ s_{VASO,dd} ({TE}_{mid})$ (%)	$T_{2, EV}^{*, rest}$ (ms)	$T_{2, EV}^{*, act}$ (ms)
ME-DEPICTING
11 subjects, averaged	10.7	0.0553	–0.9	38.6	39.5
P1 (3 × 3 × 4 mm ³ )	10.7	0.0573	–1.1	31.4	32.1
P1 (1.7 × 1.7 × 1.7 mm ³ )	15.1	0.0635	–1.0	49.0	51.5
ME-EPI
11 subjects, averaged	20.7	0.0567	–1.5	44.4	45.3
P1 (3 echoes)	20.7	0.0579	–1.7	41.4	42.3
P1 (5 echoes)	40.6	0.0553	–1.3	42.8	43.5

Funding1

—Max Planck Society

Keywords

multi-echo EPIcenter-out EPIVASOΔ R 2 *bloodextravascular fractionintravascular BOLD response

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Taxonomy

TopicsAdvanced MRI Techniques and Applications · Advanced Neuroimaging Techniques and Applications · Cardiac Imaging and Diagnostics

Full text

Introduction

1

VAscular Space Occupancy (VASO) functional magnetic resonance imaging (fMRI) provides indirect measures of changes in cerebral blood volume (CBV; in units of ml blood per ml of tissue) by nulling all blood signal at the time of acquisition (Lu et al., 2003), such that the VASO signal is related to the CBV according to (Huber, Ivanov, et al., 2014):

[eqn]

Since its inception in 2003, the technique has evolved into many variants, including Multiple Acquisitions with Global Inversion Cycling (MAGIC) VASO, which allowed for multi-slice and whole brain imaging (Lu, Van Zijl, et al., 2004;Scouten & Constable, 2007); inflow-based VASO (iVASO) for arterial and arteriolar quantitative CBV measurements (Hua et al., 2011); or Slab-selective Inversion (SI) VASO (Jin & Kim, 2008) and Slice-Saturation Slab-Inversion (SS-SI) VASO (Huber, Ivanov, et al., 2014) for improved sensitivity. Of these, SS-SI-VASO has caught the most traction in the past few years due to its applicability in high-resolution high-field (≥7 T) layer-fMRI studies (Huber et al., 2015). The higher sensitivity and temporal resolution of SS-SI-VASO can also be exploited at lower field strengths. Corresponding implementations, however, were delayed due to the lower demand for layer-fMRI studies at 3 T. There are at present only three published SS-SI-VASO studies at 3 T (Guidi et al., 2023;Huber, Kronbichler, et al., 2023;Knudsen et al., 2023).

Alternatives to blood oxygenation level-dependent (BOLD) fMRI techniques, such as VASO fMRI, also serve as proxies for neuronal activation and help supplement a more holistic understanding of BOLD fMRI. The irony, however, is in the fact that at non-zero echo time (TE), these non-BOLD methods themselves suffer from contamination by the BOLD response (Hetzer et al., 2011;Huber et al., 2019). The functional sensitivity and accuracy of these measurements could then benefit greatly from short-TE readouts. One such readout is the Double-shot Echo Planar Imaging with Center-out Trajectories and Intrinsic NaviGation or DEPICTING (Hetzer et al., 2011). The multi-echo (ME) version of DEPICTING with very short first TE ( $[eqn]$ ) and inter-echo time was recently found to substantially reduce BOLD contamination in pseudo-Continuous Arterial Spin Labeling (pCASL) (Alsop et al., 2015;Dai et al., 2008;Lorenz et al., 2018) measurements of cerebral blood flow (CBF) changes, while providing reliable measures of the simultaneous BOLD response (Devi et al., 2022). SS-SI-VASO generally takes care of the inherent BOLD contamination at non-zero TE by employing a dynamic division strategy, wherein the blood-nulled image is divided by its consecutive non-nulled (BOLD-weighted) control image. A complete correction for BOLD contamination, however, relies on a number of factors, most importantly being the assumption of a similar BOLD contribution in both thenullingandcontrolcondition, that is, extravascular (EV) BOLD contributions being equivalent to those of the combined extravascular plus intravascular (IV) BOLD response in the non-nulled parenchyma (Huber, Ivanov, et al., 2014).

In the present study, the feasibility of ME-DEPICTING as a potential readout for SS-SI-VASO at 3 T was investigated. Owing to its shorter TEs, a higher functional sensitivity is expected, as compared to the ME variant of the commonly used Echo Planar Imaging (EPI). All assessments of ME-DEPICTING-based SS-SI-VASO were, hence, contrasted to those of ME-EPI-based SS-SI-VASO. Apart from the provision for simultaneous BOLD and CBV measurements, the ME readouts also allowed for a comparison of the BOLD-correction strategies. In particular, a correction provided by ME extrapolation to $[eqn]$ was tested against the correction provided by the original dynamic-division strategy.

Materials & Methods

2

Participants

2.1

Sixteen healthy volunteers (30 ± 5 years, 9 female) gave written informed consent before undergoing the experiments that had been approved by the Ethics Committee at the Medical Faculty of Leipzig University. All participants were right-handed and had normal or corrected-to-normal vision.

Functional paradigm

2.2

A full-field 8-Hz flickering black-and-white radial checkerboard was used for visual stimulation. Each functional cycle lasted 20 repetitions, starting with a rest block, which consisted of a blank gray screen of 12 repetitions (24 s), followed by the task block of 8 repetitions (16 s). The paradigm was programmed using Presentation (v17.2, Neurobehavioral Systems, Berkeley, CA, USA). A central, colored fixation point was present throughout the experiment. Subjects were instructed to focus on this dot and press a button whenever it changed color. Their attention was monitored by visually tracking their responses.

Magnetic resonance acquisitions

2.3

SS-SI-VASO was implemented on a 3-T MAGNETOM Skyra^fit^scanner (Siemens Healthineers, Erlangen, Germany) equipped with a 32-channel receive head coil. Data were obtained from 10 slices of 4 mm thickness (no slice gap, nominal in-plane resolution 3 mm × 3 mm, field of view 192 mm, matrix 64 × 64, bandwidth 2232 Hz/Px) located along the calcarine sulcus and acquired in descending slice order. Selective inversion for the VASO scans (nulling condition) was achieved over a 34-cm slab centered at the middle of the slice package through a hyperbolic secant inversion radiofrequency (RF) pulse of 8 ms duration, bandwidth-time product of 10, and RF field peak amplitude of $[eqn]$ = 14 μT. The blood-nulling condition, corresponding to an inversion time, TI = 1153 ms, for an assumed longitudinal relaxation time of blood of $[eqn]$ (Lu, Clingman, et al., 2004), was fulfilled for the 6th slice. A schematic of the SS-SI-VASO implementation is given inFigure 1along with the sequence diagrams for ME-DEPICTING and ME-EPI. The specifications for the two two-dimensional (2D) ME readout modules were as follows:

SS-SI-VASO implementation with ME-DEPICTING and ME-EPI as imaging readouts. For the slice package centered at the calcarine sulcus, the blood nulling condition at 3 T was fulfilled for the 6th slice (in orange) of the VASO images. The interleaved acquisition of VASO and Control images was accomplished with both ME-DEPICTING and ME-EPI readouts forn= 3 echoes. The sequence diagrams of the two readouts are provided. The double-shot nature of DEPICTING results in a shorter TE when compared to ME-EPI.

ME-DEPICTING: $[eqn]$ / $[eqn]$ / $[eqn]$ = 1.7 ms/10.7 ms/19.7 ms; repetition time, $[eqn]$ = 2 s; GRAPPA (GeneRalized Autocalibrating Partial Parallel Acquisition) factor 2.ME-EPI: $[eqn]$ / $[eqn]$ / $[eqn]$ = 7.5 ms/20.7 ms/33.9 ms; $[eqn]$ = 2 s; GRAPPA factor 2; partial-Fourier factor 6/8.

Functional runs consisted of ten functional cycles (200 repetitions, i.e., 100 nulling/control condition pairs; effective TR, $[eqn]$ = 4 s) and were of approximately 7 min duration. ME-EPI and ME-DEPICTING acquisitions were recorded during the same session. Their order was shuffled across participants to avoid primacy bias.

In addition to these two functional runs, two auxiliary functional runs of similar duration were recorded in one of the subjects, with ME-EPI and ME-DEPICTING specifications different from those stated above. They were:

(i)ME-EPI with 5 echoes: $[eqn]$ / $[eqn]$ / $[eqn]$ / $[eqn]$ / $[eqn]$ = 10 ms/25.3 ms/40.6 ms/55.9 ms/71.2 ms; $[eqn]$ = 2 s; GRAPPA factor 2; partial-Fourier factor 7/8.(ii)High-resolution ME-DEPICTING: 15 slices, 1.7 mm isotropic nominal resolution; $[eqn]$ / $[eqn]$ / $[eqn]$ = 2.8 ms/15.1 ms/27.4 ms; $[eqn]$ = 2 s; GRAPPA factor 2.

For co-registration purposes, a 2D spoiled gradient-recalled echo (GRE) scan (1.5 mm nominal in-plane resolution; $[eqn]$ = 3.1 ms; $[eqn]$ = 1300 ms; flip angle 90°) was obtained at the start of each session, with slice geometry identical to the functional scans.

Data preprocessing and analysis

2.4

Data were preprocessed and analyzed using Statistical Parametric Mapping (SPM12; Wellcome Trust Centre for Neuroimaging, UCL, London, UK) implemented in MATLAB 2022b (The MathWorks, Natick, MA, USA), FMRIB Software Library (FSL, (Jenkinson et al., 2012)), and additional scripts written in Interactive Data Language (IDL 8.1, Exelis Visual Information Solutions, Boulder, CO, USA).

Preprocessing

2.4.1

All functional data were preprocessed in an identical manner. The data were split into VASO and BOLD datasets using ‘*fslsplit’*and ‘fslmerge’. The series acquired at $[eqn]$ of each dataset was realigned, and the resulting realignment parameters applied to the corresponding remaining echoes. To allow for a meaningful comparison of functional sensitivity between the two readouts, only participants who exhibited minimal motion, defined as translation <1 mm and rotation <0.015 radians in both scans, were included in the analysis. The mean correlations between the realignment parameters of the VASO and BOLD scans in these subjects were also evaluated. The realigned time series were temporally high-pass filtered at a cut-off frequency equal to two effective functional cycles [1/(20 $[eqn]$ )] and 3D-Gaussian filtered at a full width at half maximum (FWHM) equal to the nominal voxel size (3 mm or 1.7 mm). Extrapolated ME signals at TE = 0, $[eqn]$ (please refer^1^) and $[eqn]$ , as well as the effective transverse relaxation rates, $[eqn]$ and $[eqn]$ , of the VASO and BOLD datasets, respectively, were then extracted from the preprocessed ME image volumes acquired at $[eqn]$ with voxel intensity $[eqn]$ via linear regression (LINFIT) of the expression $[eqn]$ .

Parameters

2.4.2

Percent VASO signal changes were then evaluated from: (i)uncorrectedfirst echo of the blood-nulled dataset, $[eqn]$ , (ii) blood-nulled data of all echoes corrected by dynamic division, such that $[eqn]$ $[eqn]$ , and (iii) the ME-extrapolated intercept, $[eqn]$ . The simultaneously acquired BOLD signal was evaluated (i) in terms of $[eqn]$ from the $[eqn]$ and $[eqn]$ regression results, and (ii) in terms of percent BOLD signal changes from the weighted summation of all echoes, $[eqn]$ (Poser et al., 2006;Posse et al., 1999). The weights were computed from $[eqn]$ and the underlying BOLD model according to the fitted $[eqn]$ of each voxel (seeEq. 6ofPoser et al., 2006).

Analysis

2.4.3

A general linear model (GLM) with the canonical hemodynamic response function (double gamma function) was implemented in IDL for statistical analyses and applied to all data.

Quantification

2.4.4

Activation-related percent signal changes were quantified based on the resulting $[eqn]$ coefficients, such that $[eqn]$ . Percent VASO signal changes were converted into relative changes in cerebral blood volume by assuming a resting value of $[eqn]$ = 0.05 ml/ml and negligible contribution from cerebrospinal fluid (CSF) (Scouten & Constable, 2007). The imperfect blood nulling in the remaining slices was then compensated based onEq. 2btherein, rewritten here as:

[eqn]

with $[eqn]$ for non-nulled slices; $[eqn]$ for the nulled slice; and water densities of blood, $[eqn]$ 0.86 ml water/ml blood, and parenchyma, $[eqn]$ 0.89 ml water/ml parenchyma. $[eqn]$ and $[eqn]$ are the longitudinal magnetizations per unit of water in blood and parenchyma, respectively, at the time of acquisition; and $[eqn]$ is the absolute $[eqn]$ change in ml blood/ml of parenchyma. For calculation of the quantitative rate changes $[eqn]$ in units of s^−1^, $[eqn]$ was taken directly from the estimated $[eqn]$ parameter of the respective $[eqn]$ time series.

Sensitivity of VASO and BOLD signal changes

2.5

Significant regions of VASO and BOLD activation were identified based on a voxel-based significance threshold of $[eqn]$ , except for the high-resolution data, which owing to the lower signal to noise ratio (SNR) was thresholded at $[eqn]$ . Slices cropped during the realignment process due to inter-session motion were disregarded from all scans of the particular session.

The sensitivity of the VASO response was evaluated from $[eqn]$ , $[eqn]$ and $[eqn]$ , while that of the BOLD response was evaluated from $[eqn]$ and $[eqn]$ . The comparison of $[eqn]$ and $[eqn]$ VASO responses, inadvertently, allowed a comparison of the two BOLD correction strategies. The functional sensitivities were assessed in terms of (i) number of suprathreshold voxels, and (ii) the more reliable metric of temporal CNR (Geissler et al., 2007). The latter was evaluated from the common region of interest (ROI) between the ‘VASO activation’ obtained from the ME-extrapolated $[eqn]$ and the ‘BOLD activation’ taken from the $[eqn]$ images of both readouts. Comparisons of relative VASO or BOLD signal changes across TEs were also based on this ROI.

Extravascular and intravascular BOLD contributions

2.6

The parenchymal and extravascular BOLD contributions were estimated to evaluate the BOLD correction by dynamic division, which relies on the assumption $[eqn]$ $[eqn]$ $[eqn]$ . The intravascular BOLD contribution was then estimated in terms of $[eqn]$ , written as $[eqn]$ . These analyses were based on two ROIs in the nulled slice (i.e., the 6th/8th slice for standard-/high-resolution data):

(i)ROI-1 was based on common BOLD activation in the $[eqn]$ data of the two readouts. $[eqn]$ evaluated from $[eqn]$ was expected to represent parenchymal $[eqn]$ and is henceforth written as Δ $[eqn]$ .(ii)ROI-2 was based on common BOLD activation in the $[eqn]$ data of the two readouts. Δ $[eqn]$ evaluated from $[eqn]$ was expected to represent extravascular $[eqn]$ and is henceforth written as Δ $[eqn]$ .

The related effective transverse relaxation times during the resting condition, $[eqn]$ and $[eqn]$ , were also extracted. Prior to the extraction of the blood-nulled slice from the visual cortex masks (-fslslice), the multi-slice mask was multiplied by a brain mask that had been extracted from the structural image (-betwith*-f 0.2*) (Smith, 2002) and linearly registered with the $[eqn]$ images of each readout (-flirtwith*-applyxfm*,-usesqform& -noresampblur) (Jenkinson et al., 2002;Jenkinson & Smith, 2001). The relative extravascular rate change was evaluated from ROI-1 and ROI-2, such that

[eqn]

The effective intravascular BOLD contribution of both readouts was estimated solely from ROI-2 according to:

[eqn]

An expression for $[eqn]$ is derived in theAppendix A1(Eq. A6). Experimental values of $[eqn]$ and $[eqn]$ were taken from mean values in ROI-2. Similarly, $[eqn]$ was calculated from $[eqn]$ and the VASO signal from the dynamically divided middle echo $[eqn]$ with $[eqn]$ = 10.7 ms and 20.7 ms for ME-DEPICTING and ME-EPI, respectively, and $[eqn]$ = 40.6 ms for the ME-EPI with 5 echoes. The other relevant parameters had the following assumed values: $[eqn]$ = 1330 ms (Wansapura et al., 1999); relative proton densities, $[eqn]$ = 0.87 ml water/ml blood; $[eqn]$ = 0.89 ml water/ml parenchyma (Donahue et al., 2009;Herscovitch & Raichle, 1985;Lu et al., 2002). $[eqn]$ was estimated for $[eqn]$ varying from 15 ms to 40 ms in steps of 5 ms, that is from a more venous to a more arterial regime (Zhao et al., 2007). The influence of resolution and longer TE on $[eqn]$ and its dependence on $[eqn]$ was also examined for a single subject with standard versus high-resolution ME-DEPICTING data and 3-echoes versus 5-echoes ME-EPI data.

Results

3

Four participants were excluded from the evaluation based on the motion criterion and another one due to minimal activation.Supplementary Figure S1shows examples of plots of the translational and rotational displacements of an excluded and an included subject. The mean correlation between the realignment parameters of VASO and control scans over the remaining 11 participants averaged at 0.86 ± 0.07 and 0.86 ± 0.09 for EPI and DEPICTING scans, respectively (Supplementary Table S1).

Sensitivity comparison of the readouts

3.1

Significant VASO and BOLD activations were identified in all remaining participants. Examples are presented inFigure 2. The higher sensitivity of $[eqn]$ of DEPICTING at the subject level (Fig. 2A) was confirmed for the group-averaged data (Fig. 3A). The number of significant voxels showing VASO, BOLD and their common activation (VASO $[eqn]$ BOLD, from $[eqn]$ and $[eqn]$ respectively) for each participant can be assessed fromSupplementary Table S2. The temporal VASO CNR values are shown inFigure 3BandSupplementary Table S3. Interestingly, after extrapolation to zero TE, the activated areas of both readouts became rather similar and the CNR was almost equal.

Simultaneously acquired VASO and BOLD responses in participant P1. (A) VASO activation maps obtained from the BOLD-uncorrectedSbn(TE1), BOLD-correctedSbn(TE→0), andSdd(TE1). Note the larger activation areas obtained with ME-DEPICTING. (B) BOLD activation maps obtained fromSctr(sum)and fittedR2,ctrdata. TheδsVASO(%),δsBOLD(%), andΔR2*(s–1) in the activation maps, shown here for 3 slices, have been overlaid over the corresponding slices of the pre-processedSbn(TE→0)of the two readouts.*

Functional sensitivity comparison of the two readouts (group-averaged data). Results obtained with the ME-DEPICTING readout are presented in orange and those with the ME-EPI readout in blue. Differences between the two readouts are given in % (orange and blue text boxes) and their significance is marked with asterisks (, , and *** forp<0.05,0.01, and0.001, respectively). Red and blue arrows indicate sensitivity gains for ME-DEPICTING and ME-EPI, respectively. Error bars denote standard deviations across participants. (A) Sensitivity for VASO signal detection: At minimumTE,DEPICTING yielded 72% more voxels than EPI. Correction for BOLD contamination withSbn(TE→0)led to 30% and approx. 100% increases in the activated area for ME-DEPICTING and ME-EPI, respectively. Dynamic division yielded unexpectedly greater boosts of 78% and 300%, respectively. (B) VASO CNR in commonly activated regions of both readouts:Sbn(TE1)revealed a 27% higher CNR for DEPICTING compared to EPI (1.4 ± 0.3 vs. 1.1 ± 0.2). The CNR values obtained from the dynamically divided data were once again higher than those of the extrapolated data, yielding gains of 35% and 100% for ME-DEPICTING and ME-EPI, respectively. (C) BOLD CNRs based onSctr(sum)andR2,ctr: The CNR of ME-EPI was only a 12% higher than that of ME-DEPICTING. (D)Sctr(TEi) showing the expected increase of the BOLD sensitivity with TE for both readouts.

It is to be noted that even though the ROI selection relied on BOLD activation from $[eqn]$ (Supplementary Table S2), the sizes of the ROI were similar even when the significant BOLD area was based on $[eqn]$ (not shown). The sensitivity difference for the detection of BOLD signal changes was found to be negligible (~7%) between the $[eqn]$ and $[eqn]$ for both readouts (Fig. 3CandSupplementary Table S4). The CNR of BOLD data obtained from each $[eqn]$ for the two readouts from the ROI is plotted inFigure 3D.

TE dependence of the BOLD correction by dynamic division

3.2

Investigation into the startling gain in VASO sensitivity of the dynamic division-corrected data revealed further increases in both the activation area (Supplementary Table S2) and VASO CNR (Supplementary Table S3) for $[eqn]$ ) and $[eqn]$ ) for both readouts. Significantly activated areas in two arbitrary participants are presented inFigure 4along with corresponding maps obtained from $[eqn]$ and $[eqn]$ . However, in case of the dynamic-division method, along with the increase in activation area with TE, the intensity of the VASO signal change, $[eqn]$ was also found to increase with TE. The increase of $[eqn]$ with $[eqn]$ is also evident in their cycle-averaged time courses (green line with filled squares inFig. 5).

VASO activation maps obtained from the nulled slice in two participants (P1 and P5). The maps based on theSbn(TE1),Sbn(TE→0)andSdd(TEi) data show BOLD-uncorrectedδsVASO,bn(TE1)as well as BOLD-corrected results obtained by extrapolation to zeroTE,δsVASO,bn(TE→0), and by dynamic division,δsVASO,dd(TEi)for all echo times. The activation maps have been overlaid over the 6th slice of the respective pre-processed VASO images at the specified TE. For both BOLD-correction strategies, the increase of the size and strength of activation compared to uncorrected data is evident. An increase in the intensity of the VASO signal change,δsVASO,ddwith TE is indicated by the growing blue/violet area for theSdd(TE2) andSdd(TE3) data of both readouts.

Cycle-averaged time courses of percent VASO signal changes. (A) ME-DEPICTING results, (B) ME-EPI results; both obtained in participant P1 (nulled slice; see top panel inFig. 4). Red squares: BOLD response from the control data,δsBOLD,ctr(TEi); open green squares: uncorrected VASO response,δsVASO,bn(TEi); filled green squares: VASO response corrected by dynamic division,δsVASO,dd(TEi). Yellow triangles represent the difference [δsVASO,bn(TEi)−δsBOLD,ctr(TEi)]. The yellow-shaded region denotes the stimulus duration.

The time courses inFigure 5were extracted separately for each readout from a common ROI of VASO activation in the nulled slices of ME-DEPICTING and ME-EPI data. As expected, $[eqn]$ switches from a negative to a positive response with increasing $[eqn]$ for both readouts, while $[eqn]$ grows with $[eqn]$ , as evidenced by the amplitudes of the primary signals and post-stimulus undershoots. Interestingly*,* $[eqn]$ can be seen mirroring the corresponding BOLD signals. This applies to the evolution of both the signal intensities and the shapes of the time courses as a function of $[eqn]$ for both readouts. Another finding is that $[eqn]$ agreed rather well with the difference between the uncorrected VASO response during the nulling condition and the BOLD response during the control condition, $[eqn]$ .

Figures 6and7further demonstrate the TE dependence of the dynamic division-corrected VASO data for both readouts. With a near-perfect linear increase ( $[eqn]$ ), the relation of $[eqn]$ bears a striking resemblance to the TE dependence of $[eqn]$ . The slopes of the linear fits of $[eqn]$ and $[eqn]$ averaged across subjects were found to be comparable between the two readouts.

Echo-time dependence of dynamic-division corrected VASO signals. MeanδsVASO,dd(A, D), calculated over the common ROI of VASO activation in the nulled slices of ME-DEPICTING (A) and ME-EPI (D), and the correspondingδcbv obtained withEq. 2(B, E) are plotted for all individual subjects against the correspondingTEi.δsBOLD,ctr(C, F) serves as references for the TE dependence. The subject-averaged values are shown in red (also seeSupplementary Tables S5 and S6).

Linear relationships of subject-averaged VASO and BOLD responses withTE. Subject-averagedδsVASO,dd(A, C) andδsBOLD, ctr(B, D) were obtained fromSdd(TEi)andSctr(TEi)recorded with ME-DEPICTING (A, B) and ME-EPI (C, D). The linear relations and their corresponding goodness of fit,r2, are given within the rectangular boxes. For comparison purposes, the subject-averaged δsVASO,bn(TE→0)values are indicated in (A) and (C) by black diamonds on they-axes. Residual BOLD signals extrapolated to zeroTEwere 0.02±0.14% and 0.07±0.03% for EPI and DEPICTING, respectively. These residual BOLD-like time-dependent signals and the correspondingδsVASO,bn(TE→0)signals (for reference) for each individual participant are shown inSupplementary Figure S2.

Furthermore, the $[eqn]$ dependence of $[eqn]$ was also identified for the auxiliary single-subject data with 5 echoes (Fig. 8). A reduction in the extent of activation can be seen with later echoes (Fig. 8A). The amplitude of the primary signal of VASO responses grows rapidly at smaller TEs before appearing to stabilize at $[eqn]$ ≥ 25.3 ms. A shift in the shape of the time courses can be seen for $[eqn]$ ; specifically, in a faster return to baseline compared to $[eqn]$ and $[eqn]$ . An increase in the amplitude of the post-stimulus overshoot is also evident (Fig. 8B).Figure 9demonstrates a similar behavior of the dynamic division-corrected data (increase of $[eqn]$ with $[eqn]$ ) for higher-resolution data. The shape of the signals from $[eqn]$ deviates again from that of uncorrected $[eqn]$ and ME-extrapolated $[eqn]$ . The larger amplitude of VASO signal changes at higher resolution is also evident.

SS-SI-VASO acquisition with ME-EPI with 5 echoes. (A) Significant activation observed in the 6th slice (thresholded atp< 10–6) of the uncorrected VASO response fromSbn(TE1), the BOLD-corrected VASO responses from the extrapolatedSbn(TE→0), and the dynamically divided echoesSdd(TEi). The maps have been overlaid over the respective pre-processed VASO images at the specified TE. (B) Corresponding uncorrected (δsVASO,bn(TE1))and corrected cycle-averaged time courses ofδsVASO,bn(TE→0)andδsVASO,dd(TEi)for 1≤i≤5. The shaded yellow region denotes stimulus duration.

ME-DEPICTING SS-SI-VASO acquisition at higher resolution. (A) Significant VASO activation (p< 10–4) as seen on the nulled slice (8th slice) of the uncorrectedSbn(TE1)(2nd column, within the black square box), the ME-extrapolatedSbn(TE→0)(1st column), and the dynamically dividedSdd(TEi)(columns 3 to 5), overlaid over the corresponding unprocessed images. (B) Comparison of cycle-averaged time courses of VASO signal changes at different spatial resolutions:δsVASO,bn(TE1)(in green),δsVASO,bn(TE→0)(in blue), andδsVASO,dd(TE1)(gray broken line). The signals were extracted from the region defined by activation in theSbn(TE→0)data. Shaded yellow region denotes stimulus duration.

Extra- and intravascular BOLD contributions

3.3

Extravascular

ΔR2*

fraction

3.3.1

Averaged functional $[eqn]$ and $[eqn]$ values resulted in estimates of the extravascular $[eqn]$ fraction, over the parenchymal ROI-1 (146 ± 59 voxels), of $[eqn]$ = 55±15% and 50 ± 12% for ME-EPI and ME-DEPICTING, respectively. A substantial increase in this fraction was obtained in the smaller extravascular ROI-2 (45 ± 31 voxels): $[eqn]$ = 62 ± 14% and 60 ± 11%, respectively. The sizes of ROI-1 and ROI-2 for each participant and the corresponding $[eqn]$ , $[eqn]$ and $[eqn]$ are provided inTable 1for both readouts along with the $[eqn]$ and $[eqn]$ values. $[eqn]$ was found to be slightly shorter than their $[eqn]$ counterparts with both EPI and DEPICTING in both ROIs. An example of the voxel-wise distribution of $[eqn]$ , $[eqn]$ , $[eqn]$ , $[eqn]$ and the extravascular $[eqn]$ fraction over ROI-1 (participant P16: average $[eqn]$ of 51% and 46% for EPI and DEPICTING, respectively) is shown inFigure 10.

T2,rest,ΔR2andfEVmaps for participant P16 (6th slice, common ROI-1). Voxel-wise results for (A)T2,EV+IV,restandT2,EV*,rest, (B)ΔR2,EV+IVandΔR2,EV, and (C) the extravascularΔR2fraction,fEVobtained with ME-DEPICTING (top) and ME-EPI (bottom). The maps have been overlaid over the corresponding structural slice.

Increased estimates of the extravascular $[eqn]$ fraction were obtained from the single-subject high-resolution DEPICTING data with $[eqn]$ = 80% (ROI-1, 458 voxels) and 95% (ROI-2, 237 voxels). A somewhat higher extravascular $[eqn]$ fraction was also obtained with the 5-echo EPI data. These results along with those reported in recent literature for the visual cortex are summarized inTable 2.

Estimates of

ΔR2,IV*

3.3.2

Table 3lists the experimentally obtained data employed in the estimation of intravascular $[eqn]$ withEq. 4. Subject-averaged values have been provided for the multi-subject main study with 3-echoes and standard resolution. Values, corresponding to the participant (P1) from this study have also been provided along with their values for the single-subject study with 5 echoes (ME-EPI) and higher resolution (ME-DEPICTING). All values were based on ROI-2. The $[eqn]$ ( $[eqn]$ ) values were obtained by assuming $[eqn]$ = 0.05 ml/ml and are based on the $[eqn]$ data. Subject-averaged $[eqn]$ values of 0.0053 ml/ml and 0.0067 ml/ml corresponded to percent CBV changes of $[eqn]$ = 13 ± 3 and 11 ± 4% for EPI and DEPICTING, respectively. The difference between the results obtained with both readouts was insignificant (paired two-tailedt-test, $[eqn]$ = 0.08). Similarly, $[eqn]$ = 11% was obtained for the 5-echo ME-EPI while the high-resolution DEPICTING data yielded the highest $[eqn]$ change, $[eqn]$ 0.0135 ml/ml or $[eqn]$ = 27%.

The results for the intravascular $[eqn]$ estimations are plotted inFigure 11for assumed $[eqn]$ values ranging from 15 to 50 ms. A nonlinear relationship between $[eqn]$ and $[eqn]$ is evidenced in all the plots. A steeper decline of $[eqn]$ is observed in the region of venous $[eqn]$ . The differing $[eqn]$ contributions between the two readouts can also be seen converging at higher resting $[eqn]$ values (Fig. 11A). The within-sequence comparisons demonstrate a lower intravascular contribution (i.e., smaller $[eqn]$ ) with higher spatial resolution (Fig. 11B) and longer $[eqn]$ (Fig. 11C).

Dependence ofΔR2,IVonT2,IV*,rest. (A) Results obtained with subject-averaged input values for ME-DEPICTING (orange) and ME-EPI (blue). (B) High-resolution versus low-resolution data obtained in a single participant (P1) with ME-DEPICTING. (C) Single-subject results (P1) obtained with ME-EPI at different TE. The shaded green area represents the venous range withT2,IV*,rest≈ 15–30 ms (Donahue et al., 2011).*

Discussion

4

Readouts with short TE come with the promise of better sensitivity for non-BOLD contrasts (Hetzer et al., 2011;Huber et al., 2019) as BOLD contamination at $[eqn]$ <2 ms is expected to be minimal. Our recent pCASL study at 3 T had proven this to be the case for functional CBF changes in the visual cortex (Devi et al., 2022). The area of activation obtained with pCASL-prepared ME-DEPICTING at $[eqn]$ = 1.7 ms exceeded that obtained with ME-EPI at $[eqn]$ = 8 ms by 40%. The present study confirms the same for the VASO contrast. The sensitivity for detecting uncorrected VASO signal changes at $[eqn]$ of the blood-nulled dataset, $[eqn]$ was found to be substantially higher than that of ME-EPI, $[eqn]$ . However, comparing the functional sensitivity of data BOLD-corrected by ME extrapolation to $[eqn]$ (Fig. 3B), an almost equivalent VASO CNR was obtained for both readouts. This differs from the results of the pCASL study, wherein a higher CBF-CNR was obtained for the ME-DEPICTING acquisitions extrapolated to zero $[eqn]$ [Fig. 6C and 6E inDevi et al. (2022)] with approximately the same TEs as employed in the current work. Also contrary to the VASO results was a higher CNR for CBF changes obtained with $[eqn]$ data for both readouts compared to the corresponding $[eqn]$ data. An explanation for the improvement of the quality of the fitted data at $[eqn]$ in SS-SI-VASO could be the fact that the data were first split into blood-nulling and control images prior to extrapolation, whereas the pCASL data were fit as a whole, with the fluctuations between the control and label data, possibly degrading the quality of the resulting $[eqn]$ estimate. Nonetheless, as with the pCASL study, the sensitivity for the BOLD response obtained from the $[eqn]$ and $[eqn]$ data was found to be comparable for both readouts (Fig. 3CandSupplementary Table S4). For instance, $[eqn]$ yields 1500 ± 448 suprathreshold voxels with ME-DEPICTING vs. 1420 ± 439 voxels with ME-EPI ( $[eqn]$ = 0.60; CNR 3.2 $[eqn]$ 0.9 vs. 3.7 $[eqn]$ 1.1, $[eqn]$ = 0.05). This finding, along with the prediction of a shorter $[eqn]$ providing a better approximation of CBV (Genois et al., 2021), makes ME-DEPICTING a promising substitute of the traditional EPI readout for the simultaneous measurement of CBV changes accompanying the BOLD response.

BOLD correction by dynamic division

4.1

A much higher gain in VASO sensitivity was deduced for the $[eqn]$ data corrected for BOLD contributions by dynamic division compared to the correction by extrapolation to zero TE for both readouts. Further investigation into this remarkable improvement revealed a linear TE dependence of the BOLD-corrected $[eqn]$ data (Figs. 6A,6D,7Aand7C). This TE dependence bears a close resemblance to that of the BOLD data, the difference being that the latter’s intercept is close to zero (Figs. 6C,6F,7Band7D). This suggests a failure of the dynamic division strategy and, hence, the presence of residual BOLD contamination in the $[eqn]$ data, as also hinted by the similar but mirrored shape of their time courses. The VASO measures obtained from extrapolated and even uncorrected data followed temporal dynamics expected for a functional CBV response: a slower return to baseline and less prominent post-stimulus transients (Lu et al., 2003;Mandeville et al., 1998). Dynamic division-corrected data, on the other hand, showed an earlier return to baseline and post-stimulus overshoots that increased in peak amplitude with TE (Figs. 5and8B). This finding mimics that of the BOLD response in ME acquisitions (Havlicek et al., 2017) and was even more evident in the ME-EPI data with extended echo train length ( $[eqn]$ = 71.2 ms) (Fig. 8B).

Intravascular BOLD contribution

4.1.1

Relating VASO signal changes to the CBV response in SS-SI-VASO (Eq. 1) relies on the assumption of equivalent BOLD contributions in the blood-nulling and control condition (Huber, Ivanov, et al., 2014). This assumption holds in situations where the BOLD signal is mostly of extravascular origin. The extent of intravascular BOLD contributions to the BOLD signal, however, relies on a number of factors, including field strength ( $[eqn]$ ), TE, readout sequence, and diffusion weighting (Duong et al., 2003;Silvennoinen et al., 2003;Uludaǧ et al., 2009). The intravascular $[eqn]$ fraction is expected to decrease with $[eqn]$ with predicted intravascular BOLD contributions from the microvasculature of approximately 57%, 36%, 11%, and 5% at the respective field strengths of 1.5 T, 3 T, 4 T, and 4.7 T and $[eqn]$ = $[eqn]$ (Uludaǧ et al., 2009). Experimentally obtained extravascular $[eqn]$ fractions are approximately 70% at 3 T (Donahue et al., 2011;Lu et al., 2003) and approximately 90% at 7 T (Cheng et al., 2015;Donahue et al., 2011) in the human visual cortex. The results from the current study are well within the range of what has been reported (Table 2). Our estimated $[eqn]$ values (Fig. 11) agree also with previousin-vitromeasurements in bovine blood of varying hematocrit (Hct) levels ( $[eqn]$ = –8.16 s^–1^,—14.3 s^–1^, and—16.6 s^–1^for Hct = 0.21, 0.44, and 0.57, respectively), with the assumption of changes in blood oxygen saturation fraction,Y, from 0.61 to 0.73 (Zhao et al., 2007).

The observed linear dependence on $[eqn]$ of the VASO data corrected by dynamic division align with simulations based on a vascular anatomical model (Genois et al., 2021), wherein the intravascular BOLD contributions were taken into account. The linear dependence of $[eqn]$ on $[eqn]$ ranging from 0 to 30 ms as depicted inFigure 5AofGenois et al. (2021)for 3-T data with an extravascular fraction of 72% resemble our plots inFigure 6Aand6D.

Our results, however, differ from those reported in the original SS-SI-VASO paper (Huber, Ivanov, et al., 2014). The 3-T results therein, with data from two participants, found $[eqn]$ to be independent of TE (14 ms and 30 ms) in the visual cortex. This is remarkable as at a nominal resolution of 3 × 3 × 4 mm^3^, partial voluming effects are expected to be at play and introduce significant intravascular BOLD contributions. For instance, with the assumption of a 5% increase in oxygen metabolism, 57.5% increase in CBF, TE = 50 ms, and TR = 1s,Zhao et al. (2007)estimated an intravascular BOLD contribution of 13% from pure parenchyma, which increased to a substantial 42% for parenchyma contaminated by 2% veins. Consistently, a reduction in intravascular BOLD contributions from larger veins at higher resolution is evident in our single-subject data (Fig. 11B). Interestingly, despite an estimated 80–90% extravascular $[eqn]$ fraction for this data, the intravascular fraction was sufficient to bring about a TE dependence of $[eqn]$ .

A TE independence of 3-T ME-EPI data (12 ms and 48 ms) corrected by dynamic division has been reported in a recent layer-fMRI VASO study at sub-millimeter resolution (Huber, Kronbichler, et al., 2023), where contamination by larger veins is expected to be largely reduced. It should be noted that even in the absence of intravascular BOLD signal, any CBV change leads to a TE-dependent SS-SI-VASO signal induced by the TE dependence of the resting signals of corresponding (arterial or venous) blood compartments relative to the tissue signal. Consequently, an SS-SI-VASO signal increasing with $[eqn]$ was also found at 7 T in a study employing a resolution of 1.5 mm (Fig. 5AinHuber, Goense, et al., 2014). Quantification of CBV changes during activation at $[eqn]$ will always be governed by such biases. Under certain circumstances, however, these biases might cancel out for $[eqn]$ , because such arterial and venous contributions tend to increase and decrease the SS-SI-VASO signal, respectively.

Relationship between SS-SI-VASO signals

4.1.2

The echo-wise evaluation indicated a remarkable similarity of amplitude and transients of $[eqn]$ and the difference signal, defined as $[eqn]$ (Fig. 5). Differences in the extravascular $[eqn]$ fraction and, subsequently, in intravascular BOLD contaminations do not explain this finding. As shown inAppendix A2, the similarity of both quantities is easily explained by the relation $[eqn]$ (seeEq. A9). Because $[eqn]$ holds (positive BOLD response) and given the small stimulus-related signal change (~1–2%), $[eqn]$ is expected to be slightly smaller than $[eqn]$ . This relationship, however, could explain the stabilization of $[eqn]$ at longer $[eqn]$ in our 5-echo example. Since $[eqn]$ contains an intravascular signal whereas $[eqn]$ does not, their difference is expected to vary faster at shorter $[eqn]$ , because at 3 T, the intravascular signal increases significantly more than the extravascular (Triantafyllou et al., 2011).

Imaging sequence-related effects

4.1.3

Differences between the EPI readouts employed in previous studies and the DEPICTING readout, such as different $[eqn]$ -sensitivities of the specific k-space trajectories (Devi et al., 2022;Patzig et al., 2021), may contribute to the failure of the dynamic-division strategy in the present work. Although this aspect is likely to be of only secondary order because a TE dependence was introduced in our $[eqn]$ data and corresponding $[eqn]$ acquired with EPI as well (Fig. 6A,6B,6Dand6E), it cannot be ignored. As shown inAppendix A3, additional $[eqn]$ -weighting due to the k-space trajectory would result in non-zero $[eqn]$ and may be captured by an offset term $[eqn]$ leading to an effective echo time, $[eqn]$ . With the intercepts obtained inFigure 7Band7Dand the results fromTable 2,Eq. A11yields estimates of $[eqn]$ ≈ 0.96 ms and 0.35 ms for DEPICTING and EPI, respectively. The larger offset $[eqn]$ for DEPICTING is due to the different $[eqn]$ -sensitivity of its k-space trajectory (Devi et al., 2022;Patzig et al., 2021). A significant deviation from the nominal TE results only at $[eqn]$ of DEPICTING ( $[eqn]$ ≈ 0.56 $[eqn]$ ), but no relevant bias at longer TEs and for EPI. This effect leads to the slightly lower $[eqn]$ value estimated from $[eqn]$ of DEPICTING (11 ± 4%) compared to that of EPI (13 ± 3%) and also explains somewhat larger BOLD-like fluctuations in its extrapolated $[eqn]$ signal (cf., participants P11 and P12 inSupplementary Fig. S2). The different zero-crossing TE of the $[eqn]$ of the two readouts (Supplementary Fig. S3: approximately 11 ms and 17 ms for DEPICTING and EPI, respectively) could also be attributed to this difference in BOLD weighting at equivalent TEs.

A distinction lies in the TEs employed by earlier studies, which were comparatively longer and, hence, less likely to be significantly influenced by intravascular BOLD contributions. A consistent tendency could be the near-plateauing of $[eqn]$ values at longer TE of the ME-EPI data, while the ME-DEPICTING data with shorter TEs appear to follow an increasing trend over the entire range (Fig. 6Aand6D). This is also evident in time courses of later echoes extracted from the nulled slice of ME-EPI inFigures 5and9B. Further ME studies might provide additional insight into this matter. Such studies are, however, presently rare, even with SS-SI-VASO being the current workhorse of CBV-based layer-fMRI studies. Apart from the work mentioned above, all other studies used single-echo acquisitions at 7 T with 15 ms ≤ $[eqn]$ ≤ 28 ms (Beckett et al., 2020;de Oliveira et al., 2023;Faes et al., 2023;Guidi et al., 2016,2020;Huber et al., 2015,2021;Huber, Kassavetis, et al., 2023;Liu et al., 2022;Oliveira et al., 2021,2022;Pizzuti et al., 2023) and with 27 ms at 3 T (Knudsen et al., 2023). We note that we do not expect relevant differences between the 2D readouts used here and those of the 2D Simultaneous Multi-Slice (SMS) EPI and 3D EPI variants (Huber et al., 2018), given that evaluations were exclusively based on the nulled slice.

Short TEs can also be achieved with spiral acquisitions. With recent developments in their implementation, artifact correction, and reconstruction, these can be expected to be a promising alternative to EPI sequences with Cartesian trajectories for non-BOLD fMRI (Glover, 2012;Kasper et al., 2022). A preliminary report of a two-fold improvement in temporal SNR (tSNR) compared to 3D-EPI was recently presented for VASO fMRI with spiral readouts at 7 T (Monreal-Madrigal et al., 2023).

Limitations

4.2

SS-SI-VASO was implemented in the present study for TI/TR = 1153 ms/2000 ms. This fulfills the requirements of the original paper [see Fig. 1a ofHuber, Ivanov, et al. (2014)] by choosing a TI shorter than the arterial transit time of blood and a ‘period III’ (i.e., $[eqn]$ ) that prevents contamination with blood that was inverted more than once. Our timing, however, does not null the CSF contribution to the signal. A higher extravascular signal intensity was preferred over CSF nulling by our choice of a shorter TR compared to that required for simultaneous nulling of blood and CSF contributions (TR = 2.75 s). Pilot experiments had revealed a more focal activation with CSF nulling. However, $[eqn]$ values within a common ROI were found to be very similar. Nonetheless, partial volume effects due to CSF contribution cannot be disregarded in our results. The TI of 1153 ms was based on a blood $[eqn]$ value of 1664 ms measured from bovine blood at 3 T (Lu, Clingman, et al., 2004), which is lower than recently published $[eqn]$ values (~1800 ms) of human blood (Li et al., 2016,2017). With our approximate slice acquisition time of 71 ms, our choice of a TI might has nulled the 7th rather than the 6th slice. This is expected to result in only minimal blood contributions and, hence in significant over- or underestimations of $[eqn]$ changes.

The ROI definition in the current study could also have been attributed to deviations from the original report. The voxel-based thresholding employed in the present work could bias the selection to voxels that showed the maximum percent signal changes and, hence, towards voxels with more partial voluming with blood vessels. On the other hand, the more widely employed cluster-based thresholding is sensitive to weaker and more diffuse signal changes, while suffering from lower spatial specificity (Woo et al., 2014).

The very short TE achieved with ME-DEPICTING was of primary interest for the current work and as such, 2D readouts were compared. While the error due to incompletely nulled slices was corrected for evaluating $[eqn]$ , such corrections are not easily obtained for $[eqn]$ values. Consequently, estimations of $[eqn]$ , and hence also $[eqn]$ , were limited to the single blood-nulled slice. Future work with 2D SMS-EPI (Huber et al., 2016) and SMS-DEPICTING or corresponding 3D readouts could solve this issue. At 7 T, 2D SMS-EPI was found to outperform 3D EPI (Poser et al., 2010) at lower resolutions due to lesser contributions from physiological noise (Huber et al., 2018).

Only a single pilot acquisition at higher resolution, which is likely to be less biased by intravascular BOLD contributions, and another one with more echoes have been provided here. Increasing the number of participants in each could further solidify these findings. Assessing the TE dependence may benefit from a further increase of the number of echoes. Extending the evaluation beyond the visual system, such as the motor cortex, is also warranted. The impact of intravascular BOLD contributions might also be investigated by combining 3-T experiments with those at 7 T, where the extravascular BOLD effect is expected to be dominant. Further experimentation is, indeed, needed to comment on suitable ranges of $[eqn]$ and spatial resolutions, within which the dynamic division method would be feasible at the respective field strengths. This would help identify the impact of the dynamic-division correction in SS-SI-VASO-based layer fMRI studies.

Presently, in the absence of a ground truth, $[eqn]$ and, subsequently, $[eqn]$ values obtained by extrapolation to zero TE appear to deliver the most reliable estimations, especially, when quantification is the goal. However, the assumed mono-exponential $[eqn]$ decay with TE is likely to be too simplistic and may introduce some bias.

Conclusion

5

SS-SI-VASO was implemented at 3 T with two 2D readouts. The feasibility of the ME-DEPICTING sequence as a potentially advantageous implementation compared to ME-EPI was investigated, motivated by the shorter $[eqn]$ and inter-echo time capabilities of this double-shot EPI readout. Unwanted BOLD contamination was, in fact, drastically reduced at the very short $[eqn]$ of 1.7 ms. The extrapolation to zero TE, however, resulted in equivalent VASO sensitivity for the two readouts. Interestingly, the dynamic division approach was found to bring about a higher gain in VASO sensitivity for ME-EPI than for ME-DEPICTING. VASO signal changes corrected for BOLD with this strategy were, however, found to exhibit a distinct $[eqn]$ - dependence. This is probably due to the influence of intravascular BOLD contributions at 3 T. The correction for BOLD contamination by extrapolation to zero TE using multiple echoes does not rely on a negligible intravascular BOLD signal and is, therefore, recommended, if an overestimation of functional CBV changes is to be avoided. Consequently, in the absence of ME data, the short $[eqn]$ of the DEPICTING readout still provides an alternative to conventional EPI if better sensitivity and accuracy of the VASO data are desired.

Supplementary Material

Bibliography62

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1Alsop , D. C. , Detre , J. A. , Golay , X. , Günther , M. , Hendrikse , J. , Hernandez-Garcia , L. , Lu , H. , Macintosh , B. J. , Parkes , L. M. , Smits , M. , Van Osch , M. J. P. , Wang , D. J. J. , Wong , E. C. , & Zaharchuk , G. ( 2015 ). Recommended implementation of arterial spin-labeled Perfusion mri for clinical applications: A consensus of the ISMRM Perfusion Study group and the European consortium for ASL in dementia . Magn Reson Med , 73 , 102 – 116 . 10.10 · doi ↗ · pubmed ↗
2Beckett , A. J. S. , Dadakova , T. , Townsend , J. , Huber , L. , Park , S. , & Feinberg , D. A. ( 2020 ). Comparison of BOLD and CBV using 3D EPI and 3D GRASE for cortical layer functional MRI at 7 T . Magn Reson Med , 84 , 3128 – 3145 . 10.1002/mrm.28347 32557752 PMC 7721980 · doi ↗ · pubmed ↗
3Cheng , Y. , van Zijl , P. C. M. , & Hua , J. ( 2015 ). Measurement of parenchymal extravascular R 2* and tissue oxygen extraction fraction using multi-echo vascular space occupancy MRI at 7 T . NMR Biomed , 28 , 264 – 271 . 10.1002/nbm.3250 25521948 PMC 4297270 · doi ↗ · pubmed ↗
4Dai , W. , Garcia , D. , De Bazelaire , C. , & Alsop , D. C. ( 2008 ). Continuous flow-driven inversion for arterial spin labeling using pulsed radio frequency and gradient fields . Magn Reson Med , 60 , 1488 – 1497 . 10.1002/mrm.21790 19025913 PMC 2750002 · doi ↗ · pubmed ↗
5de Oliveira , Í. A. F. , Siero , J. C. W. , Dumoulin , S. O. , & van der Zwaag , W. ( 2023 ). Improved selectivity in 7 T digit mapping using VASO-CBV . Brain Topogr , 36 , 23 – 31 . 10.1007/s 10548-022-00932-x 36517699 PMC 9834127 · doi ↗ · pubmed ↗
6Devi , R. , Lepsien , J. , Lorenz , K. , Schlumm , T. , Mildner , T. , & Möller , H. E. ( 2022 ). Multi-echo investigations of positive and negative CBF and concomitant BOLD changes . Neuroimage , 263 , 119661 . 10.1016/j.neuroimage.2022.119661 36198353 · doi ↗ · pubmed ↗
7Donahue , M. J. , Blicher , J. U. , Østergaard , L. , Feinberg , D. A. , Mac Intosh , B. J. , Miller , K. L. , Günther , M. , & Jezzard , P. ( 2009 ). Cerebral blood flow, blood volume, and oxygen metabolism dynamics in human visual and motor cortex as measured by whole-brain multi-modal magnetic resonance imaging . J Cereb Blood Flow Metab , 29 , 1856 – 1866 . 10.1038/jcbfm.2009.107 19654592 · doi ↗ · pubmed ↗
8Donahue , M. J. , Hoogduin , H. , Van Zijl , P. C. M. , Jezzard , P. , Luijten , P. R. , & Hendrikse , J. ( 2011 ). Blood oxygenation level-dependent (BOLD) total and extravascular signal changes and ΔR 2* in human visual cortex at 1.5, 3.0 and 7.0 T . NMR Biomed 24 , 25 – 34 . 10.1002/nbm.1552 21259367 · doi ↗ · pubmed ↗