PET-BIDS, an extension to the brain imaging data structure for positron emission tomography.

Background

Positron Emission Tomography (PET) was developed in the late 1950s with the ultimate goal of measuring and visualizing physiological processes in vivo such as metabolism, blood flow and the concentration of proteins in various receptor systems[1-3]. Since then, PET has been used extensively in pre-clinical and clinical settings, mostly for oncological purposes[4], but increasingly also in areas of cardiology and neurology (for investigation of, e.g., synaptic plasticity, neuroinflammation, and neurodegeneration[5]). For brain imaging specifically, PET has largely been applied together with high affinity radiolabeled molecules to quantify the brain’s distribution, concentration and drug occupancy of proteins, using pharmacokinetic models[6]. The outcomes of this work have provided significant insights into the complex neurobiology of receptor systems in the healthy and diseased brain, as well as advancements in our understanding of pharmacological treatments[6]. However, compared to imaging modalities such as structural Magnetic Resonance Imaging (MRI), experiments using PET often involve data from multiple sources including, e.g., chemical characteristics of radiolabelled molecules and blood and metabolite data acquired during the imaging procedure (Fig. 1). The acquisition and availability of different types of data related to the PET experiment depends on the biological target of interest, and complicates the standardization of the acquired data for storage, analysis and sharing. Furthermore, due to differences in PET scanner vendors and blood sampling devices across PET centers, data are often stored in different file formats, and with varying nomenclature describing the same type of data. Historically, the field of PET brain imaging attained maturity with the adoption of a consensus nomenclature for pharmacokinetic modeling of PET data in 2007[7]. Recently in 2020, similar to the Organization for Human Brain Mapping recommendations from the Committee on Best Practices in Data Analysis and Sharing (COBIDAS) for MRI and MEEG, consensus guidelines for describing the content and format of PET brain data in publications and archives has emerged[8].

Fig. 1

Overview of a common PET experiment. This example includes blood measurements and is defined on a common time scale. Note, “time zero” can either be defined as time of injection or scan start, and all the PET and blood data should be decay-corrected to this time point.

The Brain Imaging Data Structure (BIDS) standard was originally developed for MRI as a community standard for organizing and sharing brain imaging study data within and between laboratories[9]. While the focus was originally targeting MRI, and especially structural and functional MRI (fMRI), BIDS has since rapidly expanded to include many different imaging modalities, including magnetoencephalography (MEG)[10], electroencephalography (EEG)[11], intracranial electroencephalography (iEEG)[12], and linking brain imaging to genetics[13]. BIDS mainly addresses the heterogeneity of data organization by following the FAIR principles (findability, accessibility, interoperability, and reusability)[14]. Findability and reusability are addressed in BIDS by providing rich metadata in dedicated sidecar files. Interoperability is addressed by using the existing Neuroimaging Informatics Technology Initiative (NIfTI) format for storing brain imaging data, text files arranged according to the JavaScript Object Notation (JSON) format for complementary metadata, and Tab Separated Value (TSV) format for tabular data. While accessibility is not directly addressed within the BIDS standard itself, the existence of such a standard facilitates the development of public data repositories. The largest of these repositories, OpenNeuro (https://openneuro.org), is already built around the BIDS standard, and a new repository fully dedicated to PET scans (OpenNeuroPET) is currently under development (https://openneuro.org/pet). BIDS also fosters interoperability and reuse of already acquired data by defining how to structure and store data using naming conventions and dedicated metadata files. Because BIDS data follow a common data structure and description, the proliferation of BIDS datasets incentivizes the creation of analysis pipelines that target this structure, the adoption of which promotes verifiable and reproducible research practices.

In this work, we present the main features of the extension of BIDS to include PET data (PET-BIDS). This extension largely builds upon the original BIDS specification[9], and the guidelines for the content and format of brain PET data in publications and archives[8]. The full documentation of the PET-BIDS extension can be found in the general BIDS specification (https://bids-specification.readthedocs.io/en/stable/).

PET-BIDS Summary

The extension of BIDS to include PET data largely aligns with the BIDS specification for other imaging modalities, describing a way for organizing PET data and specifying metadata for PET experiments. A detailed overview of the directory structure is given in Fig. 2. Each subject’s data corresponds to a directory which contains subdirectories for each acquisition session (e.g., baseline, rescan or intervention) and imaging modality (e.g., pet). The subdirectories are also accompanied by a dataset_description.json file containing generic information about the dataset, providing full credit to the authors sharing the data. Furthermore, a free form text file (README) describing the dataset in more detail such as issues with the data (e.g. excessive head movement or a lesion in the brain for a given subject or scan) should be provided (Fig. 2). Within each subject directory, the /pet subdirectory contains the PET imaging data and the corresponding metadata. PET imaging data are stored in 4D (or 3D if only one volume was acquired) NIfTI files with a _pet suffix. When acquiring several volumes (frames in PET terminology) these should be stored in 4D in chronological order (the order they were acquired in). The imaging data are accompanied by a _pet.json sidecar file which details the metadata of the PET acquisition. All the metadata are in accordance with the guidelines for the content and format of brain PET data[8].

Fig. 2

Exemplary PET-BIDS dataset with a dataset description. This includes the adequate acknowledgements (1), previews of PET files (2,3), including blood (4,5) and MRI data (6). The left side shows a directory tree of a common PET-BIDS dataset, with files in the root directory describing the dataset (README and data description.json), a file with participant-specific information (participants.tsv), and a JSON sidecar file describing the metadata needed to understand the corresponding TSV file. Next to the files in the root directory, there are subject directories named sub- for each study participant. In the subject directory lies all acquired data divided into modalities (anat and pet, for the structural MRI and PET, respectively). The content of the pet directory are displayed in the right side of the figure, including the metadata of the raw PET data (2), and the associated imaging data (3). The metadata for the blood data acquired using manual sampling is stored in a JSON sidecar file (4), with the corresponding blood specified in the TSV file (5). The columns in the TSV file contain (1) time, (2) plasma radioactivity, (3) whole blood radioactivity, (4) metabolite parent fraction, and (5) metabolite polar fraction. Blood data acquired using an autosampler is also available following a similar structure as (4) and (5). The PET data may be accompanied with MRI data for co-registration and region definition (6). In this case, it is required to specify if the MRI data has been corrected for gradient non-linearities (NonLinearGradientCorrection) to allow for correct co-registration with the PET data.

The file naming structure for PET data closely follows the general BIDS guidelines[9], using specified key-value pairs joined by hyphens and separated by underscores. Multiple sessions (typically visits) are encoded by adding an extra layer of directories in the form of ses-. Hence, a single session study, sub- would have a subdirectory /pet which contains PET files using the naming pattern sub- ses-pet.nii[.gz] corresponding to several acquisitions of PET data. The session label should be used in cases such as test-retest or baseline-intervention setups. Additionally, a task- can be inserted in a similar way as task-based and resting state Blood Oxygen Level Dependent (BOLD) fMRI data in the existing BIDS standard. For example, in the case of studies using combined PET/fMRI, subject-specific tasks may be carried out during the acquisition within the same session. Therefore, it is possible to specify task- in accordance with the fMRI data. Specifically for PET, multiple acquisitions per subject using different tracers during the same session are possible and the trc- must be used to distinguish between different tracers. Please keep in mind that the label used is arbitrary and each file requires a separate JSON sidecar file with details of the tracer used. Also, specifically for PET a reconstruction key rec- can be used to distinguish different types of reconstructions of the PET data. The rec- has four reserved values: acdyn, for reconstructions with attenuation correction of dynamic data; acstat, for reconstructions with attenuation correction of static data; nacdyn, for reconstructions without attenuation correction of dynamic data; nacstat, for reconstructions without attenuation correction of static data. Further information about the reconstruction should be added to the accompanying _pet.json metadata file. Finally, the run- can be used if one scan type/contrast is repeated multiple times within the same scan session/visit. For example, for dynamic PET acquisitions, subjects may have to leave the scanner to use the bathroom. While leaving the scanner would interrupt an MR acquisition, in PET this disruption is more appropriately considered missing data during a run, and the acquisition would still be considered the same session/run. However, there are also cases of acquisitions where this definition might not be entirely clear, and it will be up to the researcher to decide what makes most sense. For example, dual-time-window acquisitions[15] could be considered two runs within the same session, but it could also be considered a single run with missing data between the two time windows.

If blood data are available, such as arterial or venous samples acquired during the PET experiment, they are stored in the /pet folder alongside the corresponding PET data (Fig. 2). Blood can be sampled by an autosampler, for continuous monitoring of whole blood radioactivity, and/or manually drawn for discrete blood samples. Therefore, the recording key recording- for blood data has two reserved values: (1) autosampler, and (2) manual. The blood metadata should be stored in a JSON sidecar file with a recording- and a _blood suffix, containing information about what blood data are available (e.g. radioactivity in plasma and/or whole blood and parent compound). The blood JSON sidecar file should be accompanied by a tabular TSV file with similar naming convention, containing all the values of the available blood data. All blood data should be reported according to a unique reference time-scale in relation to a predefined time zero defined by the PET data (Fig. 1). The definition of time zero will be further explained below.

Specific PET-BIDS Considerations

In order to construct the _pet.json sidecar file which details the PET experiment metadata, a description of a common PET experiment is necessary. In Fig. 1 we present an overview of a common PET experiment (which includes the sampling of blood data, including plasma, whole blood and metabolite data). The experiment is defined on a single time scale relative to a predefined “time zero”. Notably, “time zero” will often be defined as time of injection or scan start, and the injected dose, the PET data, and blood data should optimally all be decay-corrected to time zero. However, because the time of injection does not always coincide with scan start, PET data may not be decay-corrected to the time of injection. Whether the image has been decay-corrected may be indicated in its metadata (using the fields: ImageDecayCorrected and ImageDecayCorrectionTime). The flexibility in choice of time zero to either scan start or injection time was chosen to maximize ease of use and adoption of PET-BIDS by the broadest possible spectrum of the PET community, due to potentially large differences in experimental design between PET studies. For example, scan start and injection time may not always coincide, and due to radioactive decay of the radiotracer, it is important to be aware of post-hoc decay correction. Importantly, the injected dose should always be decay corrected to the time of injection.

Across the diverse set of radiotracers and experimental designs in PET, it will not always be possible to enter the required metadata following the guidelines for sharing of PET data[8]. For example, while the injected mass and specific radioactivity are required metadata according to the guidelines[8], this is not possible to measure for certain radiotracers such as [18F]FDG due to its mass being too low to measure. In these cases, the values for injected mass and specific activity may be set to “n/a” to indicate missing values. We note that for required metadata, this is currently only valid for injected mass and specific activity, although future releases of the PET-BIDS specification may tackle further challenges in use cases that deviate from the current guidelines.

In the case of including MRI data with PET data, it is necessary to pay specific attention to the format the MR images are in, such as whether they have been unwarped to correct for gradient non-linearities. There is a specific metadata field in the BIDS specification for MRI[9] named NonlinearGradientCorrection which indicates this (please see Fig. 2 for an example). The main reason for the importance of this is that the MRI needs to be corrected for nonlinear gradients causing spatial distortions in order to have the same shape as the accompanying PET scans for co-registration[8]. Therefore, it is required to specify whether the corresponding MR images have been corrected for gradient non-linearities, using the NonLinearGradientCorrection metadata field, if PET data are present.

In general, SI units must be used to describe the data such as “Bq/mL” for radioligand concentration, and seconds for time measurements relative to either scan start or injection time (“time zero”).

Public PET-BIDS Datasets

Several example datasets (with zero-byte, i.e., empty, NIfTI files) are publicly available in the BIDS-examples GitHub repository (https://github.com/bids-standard/bids-examples). The first two of these datasets (full version) are also openly available on OpenNeuro formatted using the PET-BIDS standard:

The CIMBI Database [11C]DASB PET Example Dataset consists of test and retest measurements from two individuals to measure serotonin transporter availability[16]. No blood data are available for this dataset. It was collected as a part of the CIMBI database (10.18112/openneuro.ds001420.v1.0.1)[17,18].

The NRM2018 Grand Challenge dataset consists of baseline and intervention data from five individuals. No blood data are available for this dataset (10.18112/openneuro.ds001705.v1.0.1)[19,20].

The CIMBI Database [11C]CIMBI-36 PET Example Dataset consists of a single dynamic PET measurement of a pig to measure serotonin 2A receptor availability. Blood and metabolite data are available for this dataset. It was collected as a part of the CIMBI database (https://github.com/bids-standard/bids-examples/tree/master/pet001).

The CIMBI Database [11C]CIMBI-36 Intervention Example Dataset consists of a single dynamic PET measurement of a pig using bolus-infusion, and with a pharmacological ketanserin intervention during the scan. Blood and metabolite data are available for this dataset (https://github.com/bids-standard/bids-examples/tree/master/pet004).

The CIMBI Database [11C]AZ10419369 Visual Stimuli Example Dataset consists of two dynamic PET measurements of a single subject using combined PET/MRI. The first scan is a baseline scan, whereas the second scan includes a visual stimuli task during the scan. No blood data are available for this dataset. It was included in Hansen et al.[21] (https://github.com/bids-standard/bids-examples/tree/master/pet005).

Community Tools For Data Sharing and Analysis

Data curated into the PET-BIDS standard can be validated for BIDS compliance by using the “bids-validator”[22], a JavaScript application checking for the completeness and consistency of the data. The BIDS validator runs locally as a command line version (via Node.js), as a Docker container, or as a browser-based application (https://bids-standard.github.io/bids-validator/). Using this important validation software, PET researchers are provided with feedback about incompatibility errors as well as warnings if important pieces of metadata are missing. In addition to providing a version of the validator as open source software, we are collaborating with the software developers of major PET analysis tools (PMOD, SPM, MIAKAT and PETsurfer[23]) to facilitate rapid adoption and support of this standard, and working with major PET centres to help PET researchers convert their data into the PET-BIDS standard. Several software tools already exist to convert dicom files and ECAT data into BIDS, such as dcm2niix[24] (https://github.com/rordenlab/dcm2niix), however, the output may need post-hoc editing if required metadata are not available in the imaging header files.

The BIDS starter kit is a tool to help researchers get started with the BIDS data structure (https://github.com/bids-standard/ bids-starter-kit). It consists of a collection of community-driven guides, tutorials, helper scripts, and wiki resources. A tutorial that describes how to create a BIDS-compatible PET data set has been provided on the starter-kit wiki (https://github.com/bids-standard/bids-starter-kit/wiki), and MATLAB (bids-matlab; https://github.com/bids-standard/bids-matlab) and Python (pybids; https://github.com/bids-standard/pybids) packages are also available to produce and/or work with PET sidecar JSON and TSV files. These packages are freely available on GitHub.

According to the guidelines for the content and format of brain PET data[8], acquired PET data are defined as PET data after reconstruction into 3D or 4D frames. These data may be shared in repositories such as OpenNeuro (https://openneuro.org), which is an open archive for analysis and sharing of public neuroimaging data spearheading the movement of best practices within MRI, MEG, EEG, iEEG, ECoG, ASL and now PET[25]. The OpenNeuro platform is the successor of OpenfMRI (established in 2011, https://openfmri.org/), and the project enjoys relatively wide acceptance by the field and capitalizes on the BIDS standard. It has been running for almost a decade and is one of the fastest growing image databases, with about 12 new datasets being added per month. All datasets on OpenNeuro are validated for BIDS compliance prior to upload.

Data Analysis Pipelines and Sharing of Derived PET Data

BIDS also offers the possibility to build fully reproducible analysis workflows using the concept of BIDS applications[26]. A BIDS application is a software container capturing all the dependencies of a neuroimaging analysis pipeline (e.g., fMRIprep) that takes a BIDS formatted data set as input. Each BIDS application has the same core set of command line arguments, making them easy to run and integrate into automated platforms, allowing for full computational reproducibility. Several open-source initiatives for PET are currently under development to BIDS applications, including PETSurfer (https://surfer.nmr.mgh.harvard.edu/fswiki/PetSurfer), APPIAN (https://github.com/APPIAN-PET/APPIAN) and kinfitr (https://github.com/mathesong/kinfitr) providing tools to carry out preprocessing and/or pharmacokinetic modeling of PET data. We also highly recommend the great resources by the TURKU PET centre (http://www.turkupetcentre.net/), which have been providing thorough documentation and analysis tools for PET for several decades. The resulting outputs from BIDS applications can be shared using BIDS derivatives standards, describing the outputs of common preprocessing pipelines and pharmacokinetic models. The specification for PET-BIDS derivatives (BIDS Extension Proposal 23) is currently under development, and will capture data and metadata sufficient for a researcher to understand and reuse the output of a common PET analysis pipeline, including preprocessing and pharmacokinetic modeling.

This article describes the new PET extension to BIDS specifying a structured way of storing PET data and metadata. BIDS is a community-driven project, and the PET-BIDS specification was created in a joint effort made by the PET community (open and with community peer review) aligning with the consensus guidelines for the content and format of PET brain data in publications and archives[8]. PET-BIDS will make data sharing and software development easier, facilitate rapid development, adoption and application of new tools and methods, and ultimately foster collaboration between researchers to combine data sets from different centers to achieve larger sample sizes and improved statistical power to test hypotheses.