$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
PET imaging of oxygen and glucose metabolism using inhaled [15O]CO and [15O]O2 gases, intravenous injection of [15O]H2O, and intravenous injection of [18F]FDG have significant historical priors based on imaging accumulated from older generations of PET scanners14,15,16,17,26,27. These data commonly relied on transmission scanning with germanium/gallium sources for attenuation mapping, luminescence-limited bismuth germanate (BGO) scintillator arrays, high-activity dosing, and 2D acquisitions with ringed septa to accommodate limited scintillator luminescence and filtered back-projection reconstructions that are computationally inexpensive. Newer scanners, including that described in this protocol, use high-resolution CT for attenuation correction, highly luminescent lutetium-yttrium oxyorthosilicate (LYSO) scintillator arrays arranged in spatially dense scintillator configurations, 3D acquisitions that multiplex the efficiency of photon detection, TOF detection electronics that improve spatial localization of annihilation events, and computationally advanced reconstruction methods with 3D ordinary-Poisson ordered subsets of expectation maximization (OP-OSEM). New scanners provide superior image quality at a lower radiation dose29,30. Comparable exemplars from three generations of PET scanners are detailed in Figure 6.
This protocol takes advantage of several additional modifications to advance the quality of information from PET scanning from newer scanners: feet-first positioning, meticulously timed dose administrations, and automated arterial input function measurements25. These permit collection of arterial input functions and brain emissions data for four distinct tracers of oxygen and glucose metabolism within a 2-3 h scanning session. These data can then be submitted to traditional and novel kinetic modeling methods, as well as partial volume correction algorithms, to produce high-resolution, high-accuracy quantitative estimates of brain metabolism. Importantly, this protocol not only permits quantitative measurements of brain oxygen and glucose metabolism but also the calculation of brain AG. Notably, superior image quality from state-of-the-art scanners enables more realistic image-derived arterial input functions as well as novel analytical approaches that were previously limited by noise or sensitivity to timing differences between arterial sampling and brain emissions19.
Notably, imaging quality is more challenging to improve with 15O radiotracers than [18F]FDG. The ephemeral 122 s half-life of 15O causes severe loss of information. However, the use of gaseous radiotracers may also be pertinent, as gases may accumulate in the ambient space between the participant's head and the detection rings of the PET-CT. There, emissions from radiotracer gases can form significant sources of randoms and scattering, thereby confounding vendor-implemented scatter correction methods that use relative scaling of scattering sources compared to regions identified to be air. Even intravenously injected [15O]H2O may pose difficulties if significant amounts are exhaled as [15O]H2O vapor. Methods for scavenging gaseous radiotracers have been reported24 but not implemented in this work. Conventional scatter correction models can be inadequate for imaging 15O radiotracers, and the best performance for these tracers has depended on scattering models that account for the absolute sources of annihilation photons in the field of view, without rescaling, paired with accurate maps of attenuation. Even minor omissions of hardware, such as mirrors attached to head-coils of PET-MR scanners, upon which 15O gases commonly form condensates, can produce substantial artifacts in emissions reconstructions. Related artifacts arising from mismatched maps of attenuation following participant movement have also been reported31. As the interactions of scattering models with point-spread characteristics of gaseous radiotracers are yet poorly characterized, this work has omitted point-spread modeling methods which otherwise are essential for high-resolution PET. Optimizing reconstruction parameters for listmode for each tracer is best practice, but this work provides exemplars using simply console reconstructions that were identically applied to all tracers, as possible, as described in protocol section 9.1. Exemplars necessarily demonstrate frame timing schedules that have been adapted to the variability of tracer-kinetic timescales. The exemplars emphasize simplicity, comparability of counting statistics, and effective spatial resolutions available for each tracer, and minimizing bias that could arise while optimizing reconstruction parameters for each tracer. For comparability with existing studies of 15O, the representative results of this work use 8 iterations of 5 subsets, previously reported to be appropriate for modern PET-CT with time-of-flight acquisitions30,32. Time-of-flight may yield improved imaging resolution with 4 iterations of 5 subsets, and these optimizations are under active investigation for the tracers described in this work.
The longer scanner bore of modern PET-CT scanners encumbers investigators providing radiotracer gases. Confined gantry geometries increase the difficulties of reaching the participant's mouth when administering gases, communicating instructions to participants, and ensuring that gas administration devices do not leak radiotracers into the ambient space. Thereby, feet-first positioning of participants demonstrates advantages. Most scanner hardware makes feet-first positioning incompatible with the placement of neuroshielding devices, which typically place lead shielding that surrounds the participant's neck and partitions the participant's head from the body. Neuroshielding will reduce randoms entering the field-of-view from the body. The relative benefits of neuroshielding compared to avoidance of ambient gas radiotracers are presently poorly understood.
Quantifying brain metabolism
This paper demonstrates the results of analyses using traditional models for hemodynamics described by Raichle et al. (CBF)16 and Martin et al. (CBV)17, and for metabolism described by Mintun et al. (CMRO2)15 and Huang et al. (CMRGlc)14. These analyses utilize the arterial input functions and emissions data measured using the protocol described in this paper. Alternative means to measure the arterial input function, such as an image-derived method, are being actively pursued and tested against the invasive method. Traditional models also make full use of dynamic imaging to estimate kinetic model parameters. We have preferred the use of polynomial representations of CBF and CMRO226,27. We have favored the methods of Huang et al. for estimating CMRGlc but employing Bayesian parameter estimation28,33. AG was more recently conceived as a measure of glucose that enters glycolytic pathways but contributes to cellular functions other than oxidative phosphorylation4. Calculation of brain AG then simply requires conversion of CMRGlc and CMRO2 to molar equivalents and subtraction of one-sixth of the latter from the former.
Partial volume correction
We commonly perform regional partial volume correction on SUVR images using the symmetric geometric transfer matrix algorithm (sGTM)34. This relies on the high-resolution data obtained by MRI, which can include anatomic sequences alone, though the current protocol also obtains Human Connectome Project (HCP)-like sequences for more precise regional parcellation. We have used both anatomic (e.g., Desikan-Killiany35) and functional (e.g., Schaeffer36 or Glasser37) atlases. The former atlases often have larger regions of interest that are more appropriate for PET images obtained at a lower resolution and have been more widely used in prior PET research. However, the newer atlases are functionally and cytoarchitecturally more well-defined and may now be more applicable with higher-resolution PET imaging, such as that obtained with the Siemens Vision scanner. SUVR images for 15O scans are calculated as have been previously described and referenced to whole brain values4. Whole-brain CMRGlc and CMRO2 are then calculated using the kinetic modeling as per above to convert the SUVR images to absolute quantitative values. From this, partial volume corrected regional AG can then be calculated. While our exemplars for dynamic emissions imaging do not demonstrate partial volume correction for simplicity of presentation, partial volume correction methods that account for dynamic imaging are now being actively investigated.
Generalizability and adaptability of the protocol
This protocol describes methodological details pertinent to acquiring the highest achievable accuracy, precision, reproducibility, and interpretability from measurements of oxygen and glucose metabolism by PET in humans. The highest known accuracy and precision of measurements belong to analyses based upon principles and practices of tracer kinetics, which include direct measurements of arterial input functions, determination of regions of interest by co-registration with high-resolution structural MRI, and highly demanding use of computational resources for statistical models. In particular, invasive radial-artery cannulation, and its maintenance over serially repeated emission scans, is costly and bears finite risks of serious injury that demand the involvement of clinically experienced teams. The use of tracers based on 15O is yet another challenge for cyclotrons and their support teams. These provide correspondence with historically reported metabolic measurements in humans and ensure that interpretations of quantitative findings are built upon prior accumulated information. However, these detailed and invasive measures also provide means for validation of less invasive, less demanding methods that may be sufficient for many questions of human metabolism. Following appropriate validation, radial artery cannulation may be suitably replaced with image-derived sources of input functions38 for many specialized studies. Avoiding radial artery cannulation also precludes direct measurements of the primary metabolic product of [15O]O2, namely [15O]H2O, which can be estimated from samples of radial artery that have been centrifuged to estimate [15O]H2O of metabolism from plasma and [15O]O2 from centrifuged red blood cells. However, computationally intensive methods can be used to parametrically solve for the emergence of metabolic [15O]H2O in the arterial circulation, a linear model of which has been proposed by Mintun et al15. Many methodologies using SUVR may also be amenable to validation following comparisons to invasively obtained data. Currently, active research programs seek validation of oxygenation estimates based on the blood oxygenation level-dependent effect, and contemporaneous measurements with 15O can provide further validation39. Additionally, while only mentioned in this protocol without exemplars, detailed optimizations of image reconstructions from listmode are likely to improve the quality of measurements of metabolism. Further, many specialized methods such as motion correction of listmode data and partial volume corrections, may merit effort when study questions demand them, for example, when occasional participants have excessive head movements or studies of cortical atrophy are confounded by partial volume averaging with cerebral spinal fluid.
Cautions for PET scanning
This protocol exposes participants to ionizing gamma radiation. The current estimates of radiation in participants who have undergone the total PET session described here (including four to six 15O scans and one FDG scan), result in an average total effective dose, inclusive of all radiotracer administrations, of approximately 7 mSv (0.7 rem). This is comparable to the radiation dose of other radiotracers used in research and clinical practice and represents 14% of the total radiation dose allowable to a radiation worker in one year. All participants (and/or their legally acceptable representative) are carefully informed both in writing and verbally of the potential radiation risks prior to obtaining their consent for study enrollment.
Cautions for arterial lines
Complications of arterial cannulation include bruising (common), pain, swelling, hematoma (uncommon), bleeding at the site of insertion, vasospasm, thrombosis (rare), and reaction to heparinized saline when used (very rare). We routinely employ an interventional radiology service to place arterial lines and manage any complications if they should arise.
Cautions for MRI
The US Food and Drug Administration requires all medical devices that enter the magnetic field of an MRI scanner to undergo standardized safety testing. If a participant's device is labeled MRI Safe, then we can proceed using standard FDA-approved scanning conditions without additional risks to the participant. If a participant's device is MRI Conditional, we review the manufacturer, make, and model of the device, and then adjust the MRI scanning procedures accordingly. Imaging also has the potential risk of revealing incidental findings, which can incur additional risks such as anxiety, financial loss from additional work-up, and, rarely, complications arising from subsequent, clinically indicated invasive procedures.
Summary of limitations
Unavoidable limitations of measuring oxygen and glucose metabolism are specific to the radiotracer themselves. 15O-tracers have short half-lives that stipulate many of the most demanding requirements of this protocol, including the presence of highly specialized and accessible cyclotron facilities. Gaseous 15O-tracers introduce the complexity of gas management such as scavenging of expired gases24 and specialized methods for scatter corrections. [18F]FDG exchanges the considerable challenges of tracking metabolites in [11C]-glucose tracers for confounding by lumped-constant estimations. This protocol also has limitations arising from the confining geometry of integrated PET-CT scanners. While feet-first positioning of the participant ameliorates confinement problems, it may be incompatible with the use of neuroshielding devices. Finally, the use of invasive radial artery cannulation is itself a significant limitation. The radial artery is a proxy for the primary arterial supply of the brain. However, the radial artery may have discrepant flow characteristics from arterial supplies to the brain. Furthermore, invasive and painful procedures alter the cognitive state of study participants, limiting the scope and generalizability of neuroimaging studies that can be pursued by this protocol.