A phenome-wide examination of neural and cognitive function.

This data descriptor outlines a shared neuroimaging dataset from the UCLA Consortium for Neuropsychiatric Phenomics, which focused on understanding the dimensional structure of memory and cognitive control (response inhibition) functions in both healthy individuals (130 subjects) and individuals with neuropsychiatric disorders including schizophrenia (50 subjects), bipolar disorder (49 subjects), and attention deficit/hyperactivity disorder (43 subjects). The dataset includes an extensive set of task-based fMRI assessments, resting fMRI, structural MRI, and high angular resolution diffusion MRI. The dataset is shared through the OpenfMRI project, and is formatted according to the Brain Imaging Data Structure (BIDS) standard.

Background and Summary

The Consortium for Neuropsychiatric Phenomics (CNP) was established on the principle that discovery of genetic mechanisms for complex mental disorders will ultimately demand integration of research on large numbers of phenotypes spanning multiple biological and behavioral scales, from genome to syndrome, and will require both novel informatics and data analytic strategies. The CNP focused on two phenotype domains—memory/working memory and response inhibition—because these neurocognitive domains have been examined across multiple levels of analysis and are relevant to multiple neuropsychiatric syndromes. To examine these domains the CNP collected interviews and rating scales, self-report measures, neurocognitive exams using both paper-pencil and computerized tests, and neuroimaging data comprising structural MRI (sMRI), high angular resolution diffusion imaging (HARDI), and functional MRI (fMRI) both at rest (rsfMRI) and during five different cognitive activation paradigms. The CNP was committed to the perspective that most phenotypes should be considered dimensional until evidence of categorical structure was clear, and for this reason focused most on assessing phenotypes across a broad range of healthy people rather than specific groups representing individuals with neuropsychiatric syndromes. The CNP sample additionally included smaller groups of people who had diagnoses of schizophrenia, bipolar disorder, and attention-deficit/hyperactivity disorder (following DSM-IV), to determine where these individuals’ scores would fall with respect to the distributions observed in healthy individuals. The CNP was one of 9 Interdisciplinary Research Consortia supported by the NIH Roadmap Initiative from 2007–2012. The CNP comprised 8 linked awards, including a Coordinating Center (UL1DE019580, UL1RR024911), five linked R01 awards (RL1MH083268, RL1MH083269, RL1DA024853, RL1MH083270, and RL1LM009833) and two center grant awards (PL1MH083271 and PL1NS062410). To date the neuroimaging data from the CNP study have been used in 6 publications[1-6], and another 53 publications are linked to CNP awards.

This Data Descriptor provides a description of the neuroimaging data and related data from the CNP, which have been shared via the OpenFMRI project (http://openfmri.org)[7,8]. The shared data are described according to the Brain Imaging Data Structure (BIDS: http://bids.neuroimaging.io)[9] which is a data organization format designed to facilitate the transfer, description, and storage of neuroimaging experiment data. BIDS provides a disciplined way to arrange the different data types that comprise a neuroimaging experiment.

Methods

Participants

Healthy adults were recruited by community advertisements from the Los Angeles area. Participants with adult ADHD, bipolar disorder, and schizophrenia were recruited using a patient-oriented strategy involving outreach to local clinics and online portals (separate from the methods used to recruit healthy volunteers). After receiving a verbal explanation of the study, participants gave written informed consent following procedures approved by the Institutional Review Boards at UCLA and the Los Angeles County Department of Mental Health.

For both healthy and patient groups, participants were men or women ages 21–50 years; NIH racial/ethnic category either White, not Hispanic or Latino; or Hispanic or Latino, of any racial group; primary language either English or Spanish; completed at least 8 years of formal education; no significant medical illness; adequate cooperation to complete assessments; visual acuity 20/60 or better; and urinalysis negative for drugs of abuse (Cocaine; Methamphetamine; Morphine; THC; and Benzodiazepines). Participants in the healthy group were excluded if they had lifetime diagnoses of Schizophrenia or Other Psychotic Disorder, Bipolar I or II Disorder, or Substance Abuse or Dependence (not counting caffeine or nicotine); or current Major Depressive Disorder; suicidality; Anxiety Disorder (Obsessive Compulsive Disorder, Panic Disorder, Generalized Anxiety Disorder, Post-Traumatic Stress Disorder), Attention Deficit Hyperactivity Disorder (ADHD). ADHD criteria were assessed using Adult ADHD Interview; healthy participants were screened for sub-threshold ADHD, defined as 4 or more ADHD inattentive or hyperactive/impulsive symptoms in either childhood and adulthood; in addition, they could not have had medication treatment for ADHD within the prior 12 months. Each of the patient groups (Schizophrenia, Bipolar Disorder, and ADHD) excluded anyone with one of these other diagnoses; stable medications were permitted for the patients. For MRI studies we excluded participants who were left handed, who believed they might be pregnant, or had other contraindications to scanning (e.g., claustrophobia, metal in body, body too large to fit in scanner).

Diagnoses followed the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition—Text Revision[10], and were based on the Structured Clinical Interview for DSM-IV (SCID-I[11]) supplemented by the Adult ADHD Interview (a structured interview form derived from the Kiddie Schedule for Affective Disorders and Schizophrenia, Present and Lifetime Version (KSADS-PL)[12]), in order to enable a more detailed characterization of lifetime history of ADHD in adults. Interviewers/raters were trained to criteria as described elsewhere[13]; in brief minimum standards of acceptable symptom agreement were overall kappa of .75, a kappa specificity of .75 and sensitivity of .75, and .85 kappa for diagnostic accuracy. Diagnostic and Symptom elicitation skill was also assessed with the SCID Checklist of Interviewer Behaviors[13] and the Symptom Checklist of Interviewer Behaviors[14]. Ongoing quality assurance checks documented kappa above .75 for each rater annually during the course of the study.

A subset of the larger healthy sample and patient sample took part in two separate fMRI sessions, which each included one hour of behavioral testing and a one-hour scan on the same day. Eligible English-speaking participants between the ages of 21 and 50 were recruited from the parent study if they successfully completed all previous testing sessions, and they did not meet the following additional exclusion criteria: left-handedness, pregnancy, history of head injury with loss of consciousness or cognitive sequelae, or other contraindications to scanning (e.g., claustrophobia, metal in body).

The final sample included imaging data for healthy individuals from the community (130 subjects), as well as samples of individuals diagnosed with Schizophrenia (50), Bipolar Disorder (49), and ADHD (43). Diagnostic information is included in Table 1. 81% of patients had at least one comorbid diagnosis; information on the prevalence comorbid conditions is included in Table 2, and additional details are available from the scid.tsv file in the phenotype directory of the datasets.

Table 1

Primary diagnostic information for patient groups (SZ: schizophrenia, BP: bipolar, ADHD: attention deficit/hyperactivity disorder).

Diagnosis (DSM code)	Number
SZ, Disorganized Type (295.10)	1
SZ, Paranoid Type (295.30)	21
SZ, Residual Type (295.60)	6
SZ, Undifferentiated Type (295.90)	11
Schizoaffective Disorder (295.70)	11
BP I, Most Recent Episode Hypomanic (296.40)	4
BP I, Most Recent Episode Manic, Mild (296.41)	2
BP I, Most Recent Episode Manic, Moderate (296.42)	1
BP I, Most Recent Episode Manic, In Partial Remission (296.45)	3
BP I, Most Recent Episode Manic, In Full Remission (296.46)	6
BP I, Most Recent Episode Depressed, Mild (296.51)	2
BP I, Most Recent Episode Depressed, Moderate (296.52)	4
BP I, Most Recent Episode Depressed, Severe Without Psychotic Features (296.53)	5
BP I, Most Recent Episode Depressed, In Partial Remission (296.55)	8
BP I, Most Recent Episode Depressed, In Full Remission (296.56)	5
BP I, Most Recent Episode Mixed, Moderate (296.62)	1
BP I, Most Recent Episode Mixed, Severe With Psychotic Features (296.64)	3
BP I, Most Recent Episode Unspecified (296.70)	5
ADHD, Predominantly Inattentive Type (314.00)	17
ADHD, Combined Type (314.01)	19
ADHD, Predominantly Hyperactive-Impulsive Type (314.01)	7

Table 2

Comorbid conditions present for patients.

Diagnosis (DSM code)	Number
292.12 Hallucinogen-Induced Psychotic Disorder, With Hallucinations	1
292.84 Hallucinogen-Induced Mood Disorder	1
296.26 Major Depressive Disorder, Single Episode, In Full Remission	10
296.32 Major Depressive Disorder, Recurrent, Moderate	2
296.35 Major Depressive Disorder, Recurrent, In Partial Remission	5
296.36 Major Depressive Disorder, Recurrent, In Full Remission	4
296.80 Bipolar Disorder NOS	3
300.02 Generalized Anxiety Disorder	1
300.3 Obsessive-Compulsive Disorder	6
300.4 Dysthymic Disorder	2
303.90 Alcohol Dependence	50
304.00 Opioid Dependence	10
304.10 Sedative, Hypnotic, or Anxiolytic Dependence	7
304.20 Cocaine Dependence	16
304.30 Cannabis Dependence	40
304.40 Amphetamine Dependence	21
304.50 Hallucinogen Dependence	4
304.60 Inhalant Dependence	1
304.80 Polysubstance Dependence	8
304.90 Other (or Unknown) Substance Dependence	3
305.00 Alcohol Abuse	22
305.20 Cannabis Abuse	17
305.30 Hallucinogen Abuse	5
305.40 Sedative, Hypnotic, or Anxiolytic Abuse	3
305.50 Opioid Abuse	2
305.60 Cocaine Abuse	9
305.70 Amphetamine Abuse	8
305.90 Other (or Unknown) Substance Abuse	4
309.28 Adjustment Disorder With Mixed Anxiety and Depressed Mood	1
311 Depressive Disorder NOS	16
314.00 ADHD, Predominantly Inattentive Type	10
314.01 ADHD, Combined Type	12

Behavioral assessment

Enrolled participants completed extensive neuropsychological testing; these data are included in the phenotype subdirectory, with each data file accompanied by an associated JSON metadata file. The list of tests performed during the behavioral session is presented in Table 3. Participants were debriefed and compensated for their time at the end of the final testing session.

Table 3

Behavioral session assessments.

Domain	Measure
Demographics & General Health	Study-specific demographic form
	General Health Questionnaire
	Smoking Status
Symptoms	Young Mania Rating Scale-C (YMRS)*
	Hamilton Psychiatric Rating Scale for Depression (HAMD-28)*
	Scale for the Assessment of Negative Symptoms (SANS), Scale for the Assessment of Positive Symptoms (SAPS)*
	Brief Psychiatric Rating Scale (BPRS)*
	Hopkins Symptom Checklist (HSCL)
	Adult Self-Report Scale v1.1 Screener (ASRS)
Traits	Barratt Impulsiveness Scale (BIS-11)
	Dickman Functional and Dysfunctional Impulsivity Scale
	Multidimensional Personality Questionnaire (MPQ)—Control subscale
	Impulsiveness, Venturesomeness and Empathy Scale (IVE)
	Scale for Traits that Increase Risk for Bipolar II Disorder
	Golden & Meehl’s Seven MMPI Items Selected by Taxonomic Method
	Hypomanic Personality Scale (HPS)
	Chapman Scales (Perceptual Aberrations, Social Anhedonia, Physical Anhedonia)
	Temperament and Character Inventory (TCI)
	Munich Chronotype Questionnaire (MCTQ)
Neurocognitive Tasks†	Task-switching Task
	Spatial Capacity Task (SCAP)
	Verbal Capacity Task (VCAP)
	Delay Discounting Task (DDT)
	Balloon Analog Risk Task (BART)
	Attention Network Task (ANT)
	Continuous Performance Go/NoGo Task (CPT)
	Stroop Color Word Task (SCWT)
	Stop Signal Task (SST)
	Scene Recognition Task
	Remember-Know Task
	Spatial Maintenance and Manipulation Task (SMNM)
	Verbal Maintenance and Manipulation Task (VMNM)
	Probabilistic Reversal Learning Task (PRLT)
	Deterministic Reversal Learning Task (DRLT)
Neuropsychological Assessment	California Verbal Learning Test (CVLT-II)
	WMS-IV Symbol Span
	WMS-IV Visual Reproduction
	WAIS-IV Letter Number Sequencing
	WMS-IV Digit Span
	WAIS-IV Vocabulary
	WAIS-IV Matrix Reasoning
	Color Trails Test

*In patients only.

†See the following wiki for further details on neurocognitive tasks: http://lcni-3.uoregon.edu/phenowiki/index.php/HTAC.

Neuroimaging assessment

Imaging acquisition

MRI data were acquired on one of two 3T Siemens Trio scanners, located at the Ahmanson-Lovelace Brain Mapping Center (Siemens version syngo MR B15) and the Staglin Center for Cognitive Neuroscience (Siemens version syngo MR B17) at UCLA. Functional MRI data were collected using a T2*-weighted echoplanar imaging (EPI) sequence with the following parameters: slice thickness=4 mm, 34 slices, TR=2 s, TE=30 ms, flip angle=90°, matrix 64×64, FOV=192 mm, oblique slice orientation. Additionally, a T2-weighted matched-bandwidth high-resolution anatomical scan (with the same slice prescription as the fMRI scan) and MPRAGE were collected; the former is not included in the current data distribution. The parameters for the high-resolution scan were: 4 mm slices, TR/TE=5000/34 ms, 4 averages, matrix=128×128, 90 degree flip angle. The parameters for MPRAGE were the following: TR=1.9 s, TE=2.26 ms, FOV=250 mm, matrix=256×256, sagittal plane, slice thickness=1 mm, 176 slices. Diffusion weighted imaging (DWI) data were collected using an echo-planar sequence with parameters: 64 directions, 2 mm slices, TR/TE=9000/93 ms, 1 average, 96×96 matrix, 90 degree flip angle, axial slices, b=1000 s/mm2.

Subjects participated in two scanning sessions (‘A’ and ‘B’) in a counterbalanced fashion. Session A included: localizer, MPRAGE, DWI, Matched Bandwidth Hires, Balloon Analog Risk Task, and Paired-associate Memory Task. Session B included: localizer, Matched Bandwidth Hires, Resting State, Breath Hold Task, Stop Signal Task, Spatial Capacity working memory task (SCAP), and Task Switching Task.

Resting fMRI

The resting fMRI scan lasted 304 s. Participants were asked to remain relaxed and keep their eyes open; they were not presented any stimuli or asked to respond during the scan.

Task fMRI

For each task described below, all participants received brief training on each task immediately before scanning. Each participant viewed the task stimuli through MRI-compatible goggles, responded with his or her right hand on an MR-compatible button box in the scanner, and performed one run of each task while in the scanner. The presentation and timing of all stimuli and response events were achieved using Matlab (Mathworks) and the Psychtoolbox (www.psychtoolbox.org) on an Apple Powerbook. Code for all tasks (except for Paired Associate Memory) is available at https://poldracklab.stanford.edu/software.

Balloon analog risk task

In this task[15], participants were allowed to pump a series of virtual green (experimental) and white (control) balloons. On each trial, participants chose to pump the balloon or cash out and collect their accumulated earnings for that round. For experimental balloons, after a trial in which the participant successfully pumped the balloon (meaning it did not result in an explosion), an image of a larger balloon was presented, the participant earned 5 points, and was able to pump again or cash out. After a trial in which the participant chose to cash out, the participant’s accumulated earnings for that round were displayed and the task moved onto the next round. On an explosion trial (necessarily following a pump trial), an exploded balloon was presented, the participant received no points for that round, and the task moved onto the next round. In this version of the BART, balloons exploded randomly on a number drawn from a uniform distribution over numbers of pumps, with 12 maximum pumps possible before an explosion or end of a round. Thus, participants experienced the probability as non-stationary, as the likelihood of a loss event increased with each trial in a round and as no information was provided to subjects about the probability of explosion. Subjects were not informed regarding the maximum pumps possible before an explosion or end of the round. Participants also responded to control balloons, which increased in size on successive trials, but which neither resulted in points nor exploded. For both balloons (green and white), the balloon would disappear from the screen once the participant responded, and each balloon trial was separated by a jittered delay. An outcome trial (following a cash-out trial or an explosion) was displayed for a fixed duration of 2 s. Each trial was separated by a blank screen that was presented for a variable duration (1–2 s, average 1.5 s); each round was separated by a blank screen that was presented for variable duration (1–12 s, average 4 s).

The task performed in the scanner was self-paced, but the task was programmed such that participants saw approximately 30 virtual balloons, with an approximate run time of 9 min. Each successful pump was worth 5 points, but participants did not collect their earnings at the end of the scan.

Paired-associate memory task

In this task, two scans were performed to assess declarative memory encoding and retrieval, respectively. The first scan consisted of a block of 64 encoding trials. Forty of the encoding trials were ‘memory’ trials and 24 were ‘control’ trials. During memory trials, first, two words appeared for 1 s, one on each side of the screen. Then, line drawings of two objects that matched the words appeared above the words and they were presented together for 3 more seconds, for a total of 4 s per trial. One of the objects was drawn in black and white, and one object drawn using a single color (e.g., orange). For control trials, pairs of scrambled stimuli, one black and white and one colored, appeared for 2 s. The subject indicated by button press which side the colored object was on. Subjects were instructed to remember the objects and the relationship between the objects. The trial types were intermixed. In total, the encoding run included 242 TRs of which 104 (42.9%) were active task, and the ITI between trials was jittered. The encoding run, in total, was 8.07 min long.

During the retrieval scan there were 104 total trials: 24 control trials, 40 correct trials, and 40 incorrect trials, each trial was 3 s long. The retrieval task required the subjects to look at a pair of objects and rate their confidence in their memory of the pairing. There were 4 possible response options ranging from ‘Sure correct’ to ‘Sure incorrect’, allowing the responses to be analyzed as a spectrum or binarized into yes/no type responses. During control trials, on one side of the screen was one of the four retrieval confidence response options ‘sure correct’, ‘maybe correct’, ‘maybe incorrect’, or ‘sure incorrect’. On the other side of the screen was ‘xxxx’. Subjects were asked to press the button (1–4) that corresponded to the response option displayed (for example, a response of 1 if ‘sure correct’ appeared). In the 40 correct trials, items were shown paired as they had been at encoding. During the 40 incorrect trials items were shown paired differently than they were at encoding; some individual objects had been displayed as part of pairs during encoding, but at retrieval were paired with different objects. The retrieval scan consisted of 268 TRs, of which 156 (58.2%) were active task, and was 8.93 min long.

Spatial working memory task

During the spatial delayed response task (or, spatial capacity task—SCAP), subjects were shown a target array of 1, 3, 5 or 7 yellow circles positioned pseudorandomly around a central fixation cross[16,17]. After a delay, subjects were shown a single green circle and were required to indicate whether that circle was in the same position as one of the target circles had been. A relatively long stimulus presentation time of two seconds was used to allow subjects to fully encode the target array, minimizing a potential encoding bias on the basis of set size interaction. Likewise, decision or selection requirements were kept constant across set sizes to reduce possible effects of set size on response processes. In addition to load, delay period was manipulated, with delays of 1.5, 3 or 4.5 s. Trial events included a 2-sec target-array presentation, a 1.5, 3 or 4.5 s delay period, and a 3-sec fixed response interval. A central fixation was visible throughout each of the 48 trials (12 per memory set size, with 4 at each delay length for each memory set). Half the trials were true-positive, and half were true-negative. Before starting the in-scanner task, subjects underwent a supervised instruction and training period outside of the scanner, and once in the scanner were again reminded of the instructions.

Stop signal task

In this task[18], participants were instructed to respond quickly when a ‘go’ stimulus was presented on the computer screen, except on the subset of trials where the ‘go’ stimulus was paired with a ‘stop’ signal. Specifically, participants were shown a series of go stimuli (left- and right-wards pointing arrows), to which participants were told to respond with left and right button presses, respectively (Go trials). On a subset of trials (25%), a stop-signal (a 500 Hz tone presented through headphones) was presented a short delay after the go stimulus appeared and lasted for 250 ms (Stop trials). Participants were instructed to respond as quickly and accurately as possible on all trials, but to withhold their response on Stop trials (on trials with the tone). They also were instructed that stopping and going were equally important.

On Stop trials, the delay of the onset of the stop-signal, or stop-signal delay (SSD), was varied, such that it was increased after the participant successfully inhibited in response to a stop-signal (making the next stop trial more difficult), and decreased after the participant failed to inhibit in response to a stop-signal (making the next stop trial less difficult). Each SSD increase or decrease was in 50-ms intervals. The SSD values were drawn from two interleaved staircases per block, resulting in 16 trials from each staircase for a total of 32 Stop trials per block. On the testing day, participants completed two experimental runs (one run outside of the scanner and one while inside of the scanner). In the first task run completed outside of the scanner, SSD values started at 250 and 350 ms for staircase 1 and 2, respectively. At the end of the behavioral run, the last SSD time from each staircase was then carried over to be the initial SSD for the scan run. This one-up/one-down tracking procedure ensured that subjects successfully inhibited on approximately 50% of inhibition trials. Also as a result, difficulty level is individualized across subjects and both behavioral performance and numbers of successful stop trials are equated across subjects.

Each experiment run contained 128 trials, 96 of which were Go trials and 32 of which were Stop trials, each presented randomly. All trials were preceded by a 500 ms fixation cross in the center of the screen, then each trial began with the appearance of an arrow and ended after 1000 ms, followed by the null period. Jittered null events separated every trial (with a blank screen), with the duration of null events sampled from an exponential distribution (null events ranged from 0.5 to 4 s, with a mean of 1 s).

Task-switching task

In this task[19], participants were presented with a series of one of four possible stimuli and asked to respond to the stimulus based on the task cue presented prior to, and above, the image. The four stimuli included a red triangle, red circle, green triangle, and green circle. Participants switched between responding to the image’s color (i.e., red or green) or the shape (i.e., triangle or circle). Cues presented included either ‘SHAPE’ or ‘S’ on trials where participants were expected to respond to the shape of the stimulus; cues presented included either ‘COLOR’ or ‘C’ on trials where participants were expected to respond to the color of the stimulus. On 33% of trials the instructions switched, such that participants were instructed to switch from responding from shape to color, or vice versa. On 67% of trials, the instructions remained the same but the cues changed. This task is designed to measure the changes in reaction time between trials requiring versus not requiring a switch in responding. Participants completed a total of 96 trials, for a total run time of 6 min 52 s.

Breath holding task

In this task, participants were asked to alternate between holding their breath and breathing regularly while resting. Participants were presented with bars on the screen in order to pace breath holding. When the bar was green (lasting 13.5 s), they were asked to breathe regularly; when the bar was yellow (lasting 3 s), they were asked to prepare to hold their breath; and when the bar was red, they were asked to hold their breath for 13.5 s. This cycle was repeated five times, for a total run time of 2 min and 46 s. Participants performed the BHT in the scanner while wearing a respiratory belt in order to measure breathing. The purpose of this task was to assess the contribution of respiratory rhythms to changes in the BOLD signal.

Physiological recording

During the breath holding and resting fMRI scans, physiological data were collected using a BIOPAC MP150 with Pulse Oximeter and Respiration modules. Data were sampled at 1000 Hz and recorded using BIOPAC Acknowledge software.

Data Records

The data set is hosted on OpenfMRI under the accession number ds000030 (Data Citation 1). The files are organized in the Brain Imaging Data Structure (BIDS) version 1.0.0rc3 format, and provided in multiple gzip-compressed tar archives[9]. BIDS provides a consistent naming convention and folder structure that is amenable to automated processing.

At the top-level, machine-readable text files in javascript object notation (JSON) are provided describing each fMRI task. A demographics file (participants.tsv; tab-delimited format) contains the unique participant number, age group, diagnosis, and an indication of which scan data are present for each subject.

The neuroimaging data collected for each participant are organized in a sub- data folder. The data for each imaging series are further organized into subfolders reflecting their intended use or purpose:

func: BOLD contrast fMRI image data (both task-based and resting state), fMRI task event files, physiological monitoring logs;

anat: high resolution T1-weighted image data used for structural analysis, cortical parcellation, and co-registration to standard templates;

dwi: Diffusion-weighted image data used for experimental white matter characterization methods such as Diffusion Tensor Imaging;

beh: Stop Signal task training run collected outside of the MRI scanner.

Each of these folders contain the following data files that collectively represent a participant’s study visit:

func/sub-_task-_bold.nii.gz: fMRI time series image data collected from the MRI in NIfTI format and compressed with gzip;

func/sub_task-_events.tsv: Tab-separated text file containing the timing of the events presented during the fMRI task and the responses given by the participant;

func/sub-_task-_bold.json: a JSON-structured file containing the MRI scanner parameters used to collect the image data;

func/sub-_task_physio.tsv.gz: A gzip-compressed tab-separated text file containing the output from physiological monitoring that occurred during the fMRI resting state and breath-holding fMRI experiments;

dwi/sub-_dwi.bval: Diffusion weighting (s/mm2) applied to each volume in the diffusion weighted image series;

dwi/sub-_dwi.bvec: Diffusion weighting gradient orientations applied to each volume in the diffusion weighted image series;

dwi/sub-_dwi.nii.gz: Diffusion image data in NIfTI format compressed with gzip;

dwi/sub-_dwi.json: Diffusion weighted imaging scan parameters and metadata in a JSON text file;

anat/sub-_T1w.nii.gz: High-resolution (1×1×1 mm resolution) T1-weighted image data for structural/anatomical use;

anat/sub-_T1w.json: MR scanning parameters and metadata in a JSON text file.

The label refers to data from one of the tasks performed during the visit:

bart: Balloon analog risk task

stopsignal: Stop signal task

pamenc: Paired memory encoding task

pamret: Paired memory retrieval task

scap: Spatial working memory task

taskswitch: Task switching task

bht: Breath-holding task

rest: Resting state (eyes open)

An additional directory named ‘phenotype’ contains additional phenotypic information (such as neuropsychological testing results); for each measurement, there is a tab-delimited text file along with a JSON file containing the data dictionary and measurement-level metadata.

Technical Validation

Imaging files were converted from primary DICOM data to Neuroimaging Informatics Technology Initiative (NIfTI) version 1.1 format using the dcm2niix (https://github.com/neurolabusc/dcm2niix) program. dcm2niix extracts the image pixel data and pertinent metadata parameters from the DICOM files to populate the NIfTI file header. For diffusion weighted images, dcm2niix also extracts the diffusion gradient strengths and orientations to separate files. The dcm2niix program output was saved for each conversion and inspected for errors. Additional metadata were extracted from the DICOM files using the gdcm library (http://gdcm.sourceforge.net/wiki/index.php/Main_Page) and converted to javascript object notation (JSON) format to accompany the respective image files.

All relevant metadata were combined and plotted over the timeframe of the study to reveal instances where individual scans were performed with MRI parameter settings that deviated from the standard protocol and to highlight scans with non-standard dimensions. There were very few deviations, and most were minor in nature that did not warrant data exclusion; however the table of parameters and plots are included with the data set as reference material.

After processing the source data, BOLD contrast (fMRI) and T1-weighted anatomical imaging data were processed by the MRI Quality Control protocol (MRIQC https://github.com/poldracklab/mriqc). MRIQC computes several quality control metrics available in the published literature.

Anatomical T1w scans

All T1-weighted images were skull-stripped[20] [AFNI 3dSkullStrip], corrected for intensity inhomogeneity due to B1 variations[21] [ANTS N4], and normalized to MNI-152 2 mm template space[22] [ANTS]. The background, gray matter, white matter, and cerebrospinal fluid were segmented[23] [FSL 5.0.8/FAST], and the resulting segmentations were used to compute the a set of quality control measures (See Table 4).

Table 4

Quality control measures for Anatomical MRI.

Measure	Description
Signal-to-noise ratio (snr)	Relative measure of image quality and MR system performance. This measure can be sensitive to acquisition parameters and hardware configurations[28].
Contrast-to-noise ratio (cnr)	A measure of how much contrast between the gray and white matter is captured by the scan[28].
Foreground-background energy ratio (fber)	The mean energy of image values within the head relative to outside the head. Higher values are better.
Entropy focus criterion (efc)	Sensitive to the ghosting and blurring artifacts caused by head motion during the acquisition. Lower values are better[29].
Artifact detection (qi1)	Detects structure in the background noise to determine if any artifacts from motion, flow, wrap-around, or ghosting are present[30].

In addition to the quality metrics listed above, summary information about the mean and standard deviation of background, foreground, gray matter, white matter, CSF, and average inhomogeneity bias are provided for comparison across the subject population (see Fig. 1).

Figure 1

Distribution of selected QC measures for T1w anatomical scans included in dataset.

 

After the MRIQC was completed, each high-resolution T1-weighted anatomical scan was scrubbed of facial features using the mri_deface program[24]. Visual inspection of all outputs was performed to ensure that the facial features were properly removed to preserve subject privacy.

BOLD/T2*-weighted functional scans

Each BOLD-contrast fMRI scan was corrected for head motion during the acquisition using AFNI[20], and a brain mask was computed to separate the brain from the skull and outside air[25]. Table 5 lists the QC metrics that were were computed from these images, and Fig. 2 shows the distribution of these measures.

Table 5

Functional MRI quality control measures.

Ghost-to-signal ratio (gsr)	An indication of how much nyquist (N/2) ghosting is present due to poor EPI phase correction[31].
Temporal derivative of RMS variance across all voxels within the brain (DVARS)	A measure of how much the signal intensity changes from volume to volume. Lower values are desirable[32].
Framewise displacement (fd)	Measures the amount of subject motion that occurred during the acquisition. The framewise displacement values are provided for each subject and task on a per-TR basis as well[26].
Mean temporal signal-to-noise ratio (m_tsnr)	Quantifies the influence of the physiological and thermal background noise in the overall voxel-wise time course.
Signal-to-noise ratio (snr)	Relative measure of image quality and MR system performance. Can be sensitive to acquisition parameters and hardware configurations[28].

Figure 2

Distribution of selected QC measures for BOLD contrast functional scans included in dataset.

 

fMRI task event and physiological monitoring data

fMRI task and event onsets were extracted from the raw matlab log files using R-matlab and are provided as tabular text files. These event files contain the onset, duration, type of trial, reaction time, and task-specific details. The onset time entries in each event file were corrected for the scanner trigger delay to account for extra T1 saturation equalization pre-scans performed by the scanner. In some cases the trigger delay was not explicitly available from the source data. Plots of all event files are provided with the data set for reference.

Similarly, physiological recordings collected during resting state and breath-holding fMRI scans were converted from raw data (Acknowledge format, BioPac) using the Bioread python package (https://github.com/njvack/bioread). The individual recording traces were saved as gzip-compressed tab-separated value files, and the raw data header information (Channel Name, Units, Sampling Rate) was extracted to accompanying JSON files. Plots of physiological traces are provided overlaid with breath-hold instruction timings (ready, hold, rest) for reference.

Diffusion-weighted imaging

The diffusion-weighted data were corrected for eddy currents and head motion using the B0 image as the reference, and the motion parameters were saved for later analysis[26]. The corrected volumes were then skull-stripped to remove the background and other non-brain scan regions[25]. Diffusion tensor estimation was performed for each subject, and the mean fractional anisotropy (FA) and average diffusivity (AD) were computed for all brain voxels[27]. Quality assurance was performed using a semi-automated method, including the following steps: confirmed that the b-values and applied directions were the same as expected, calculated mean in-mask FA and average or mean diffusivity (AD, MD), calculated and plotted motion, created a standard deviation image based on regular and motion corrected files, and used regional maps to calculate the percentage of cropped voxels in the occipital lobe, frontal lobe, temporal lobes, cerebellum, and the most superior portion of the brain. Trained individuals then used the results of this script to evaluate scan quality. This included a visual inspection of the FA map, visual inspection of the color map to ensure tracts were coded correctly, a check for cropping in which if more than 10 percent of voxels in a region were cropped, the map was visually inspected to ensure that the cropped region did not encroach on a major tract (i.e., regions that would be included in the FSL DTI skeleton), and the raw data were watched as a movie. Data were flagged for coverage (0=no cropping, 1=minor cropping, 2=severe unuseable cropping), motion (based on watching raw data as movie, and on motion plots), tensor direction flags (based on b-value and direction calculations, and on observation of the color map), artifact flags (including noise, striping, and the vibration artifact that affected many Siemens Trio systems during this time period). An overall quality score from 1–4 was generated from these measures, in which 1=excellent, 2=good, 3=fair (both 2 and 3 may be useable, but depending on analysis may want to consider the reason for the decreased score), and 4=unuseable (all individuals with vibration artifacts are in this category, along with other irreconcilable problems), 5=not evaluated; the distribution of these ratings is shown in Fig. 3. Ratings are included with the data download.

Figure 3

Results from the semi-automated DWI quality assessment, by subject group.

 

Usage Notes

All data are made available under the Public Domain Dedication and License v1.0 (http://www.opendatacommons.org/licenses/pddl/1.0/), which places no restrictions on the usage of the data. We expect that users of the data will follow the ODC Attribution-Sharealike Community Norms (http://opendatacommons.org/norms/odc-by-sa/), which state that any derivative works should be shared via an equally open license, and that the creators of this dataset should be credited by citation of the present data descriptor.

We ask that users with questions please use the NeuroStars Forum (https://neurostars.org) and attach the tags ‘openfmri’ and ‘ds000030’ in order to discuss and comment on this dataset.

Additional Information

How to cite this article: Poldrack, R. A. et al. A phenome-wide examination of neural and cognitive function. Sci. Data 3:160110 doi: 10.1038/sdata.2016.110 (2016).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.