Validation of a combined image derived input function and venous sampling approach for the quantification of [18F]GE-179 PET binding in the brain

Individualized gray matter morphological abnormalities uncover two robust transdiagnostic biotypes

Psychiatric disorders exhibit a shared neuropathology, yet the diverse presentations among patients necessitate the identification of transdiagnostic subtypes to enhance diagnostic and treatment strategies. This study aims to unveil potential transdiagnostic subtypes based on personalized gray matter morphological abnormalities. A total of 496 patients with psychiatric disorders and 255 healthy controls (HCs) from three distinct datasets (one for discovery and two for validation) were enrolled. Individualized gray matter morphological abnormalities were determined using normative modeling to identify transdiagnostic subtypes. In the discovery dataset, two transdiagnostic subtypes with contrasting patterns of structural abnormalities compared to HCs were identified. Reproducibility and generalizability analyses demonstrated that these subtypes could be generalized to new patients and even to new disorders in the validation datasets. These subtypes were characterized by distinct disease epicenters. The gray matter abnormal pattern in subtype 1 was mainly linked to excitatory receptors, whereas subtype 2 showed a predominant association with inhibitory receptors. Furthermore, we observed that the gray matter abnormal pattern in subtype 2 was correlated with transcriptional profiles of inflammation-related genes, while subtype 1 did not show this association. Our findings reveal two robust transdiagnostic biotypes, offering novel insights into psychiatric nosology.

Assessment of [18F]PI-2620 Tau-PET Quantification via Non-Invasive Automatized Image Derived Input Function

PURPOSE

[18F]PI-2620 positron emission tomography (PET) detects misfolded tau in progressive supranuclear palsy (PSP) and Alzheimer's disease (AD). We questioned the feasibility and value of absolute [18F]PI-2620 PET quantification for assessing tau by regional distribution volumes (VT). Here, arterial input functions (AIF) represent the gold standard, but cannot be applied in routine clinical practice, whereas image-derived input functions (IDIF) represent a non-invasive alternative. We aimed to validate IDIF against AIF and we evaluated the potential to discriminate patients with PSP and AD from healthy controls by non-invasive quantification of [18F] PET.

METHODS

In the first part of the study, we validated AIF derived from radial artery whole blood against IDIF by investigating 20 subjects (ten controls and ten patients). IDIF were generated by manual extraction of the carotid artery using the average and the five highest (max5) voxel intensity values and by automated extraction of the carotid artery using the average and the maximum voxel intensity value. In the second part of the study, IDIF quantification using the IDIF with the closest match to the AIF was transferred to group comparison of a large independent cohort of 40 subjects (15 healthy controls, 15 PSP patients and 10 AD patients). We compared VT and VT ratios, both calculated by Logan plots, with distribution volume (DV) ratios using simplified reference tissue modelling and standardized uptake value (SUV) ratios.

RESULTS

AIF and IDIF showed highly correlated input curves for all applied IDIF extraction methods (0.78 < r < 0.83, all p < 0.0001; area under the curves (AUC): 0.73 < r ≤ 0.82, all p ≤ 0.0003). Regarding the VT values, correlations were mainly found between those generated by the AIF and by the IDIF methods using the maximum voxel intensity values. Lowest relative differences (RD) were observed by applying the manual method using the five highest voxel intensity values (max5) (AIF vs. IDIF manual, avg: RD = -82%; AIF vs. IDIF automated, avg: RD = -86%; AIF vs. IDIF manual, max5: RD = -6%; AIF vs. IDIF automated, max: RD = -26%). Regional VT values revealed considerable variance at group level, which was strongly reduced upon scaling by the inferior cerebellum. The resulting VT ratio values were adequate to detect group differences between patients with PSP or AD and healthy controls (HC) (PSP target region (globus pallidus): HC vs. PSP vs. AD: 1.18 vs. 1.32 vs. 1.16; AD target region (Braak region I): HC vs. PSP vs. AD: 1.00 vs. 1.00 vs. 1.22). VT ratios and DV ratios outperformed SUV ratios and VT in detecting differences between PSP and healthy controls, whereas all quantification approaches performed similarly in comparing AD and healthy controls.

CONCLUSION

Blood-free IDIF is a promising approach for quantification of [18F]PI-2620 PET, serving as correlating surrogate for invasive continuous arterial blood sampling. Regional [18F]PI-2620 VT show large variance, in contrast to regional [18F]PI-2620 VT ratios scaled with the inferior cerebellum, which are appropriate for discriminating PSP, AD and healthy controls. DV ratios obtained by simplified reference tissue modeling are similarly suitable for this purpose.

Contributions of network structure, chemoarchitecture and diagnostic categories to transitions between cognitive topographies

The mechanisms linking the brain's network structure to cognitively relevant activation patterns remain largely unknown. Here, by leveraging principles of network control, we show how the architecture of the human connectome shapes transitions between 123 experimentally defined cognitive activation maps (cognitive topographies) from the NeuroSynth meta-analytic database. Specifically, we systematically integrated large-scale multimodal neuroimaging data from functional magnetic resonance imaging, diffusion tractography, cortical morphometry and positron emission tomography to simulate how anatomically guided transitions between cognitive states can be reshaped by neurotransmitter engagement or by changes in cortical thickness. Our model incorporates neurotransmitter-receptor density maps (18 receptors and transporters) and maps of cortical thickness pertaining to a wide range of mental health, neurodegenerative, psychiatric and neurodevelopmental diagnostic categories (17,000 patients and 22,000 controls). The results provide a comprehensive look-up table charting how brain network organization and chemoarchitecture interact to manifest different cognitive topographies, and establish a principled foundation for the systematic identification of ways to promote selective transitions between cognitive topographies.

Convergent functional effects of antidepressants in major depressive disorder: a neuroimaging meta-analysis

Background

Neuroimaging studies have provided valuable insights into the macroscale impacts of antidepressants on brain functions in patients with major depressive disorder. However, the findings of individual studies are inconsistent. Here, we aimed to provide a quantitative synthesis of the literature to identify convergence of the reported findings at both regional and network levels and to examine their associations with neurotransmitter systems.

Methods

Through a comprehensive search in PubMed and Scopus databases, we reviewed 5,258 abstracts and identified 36 eligible functional neuroimaging studies on antidepressant effects in major depressive disorder. Activation likelihood estimation was used to investigate regional convergence of the reported foci of consistent antidepressant effects, followed by functional decoding and connectivity mapping of the convergent clusters. Additionally, utilizing group-averaged data from the Human Connectome Project, we assessed convergent resting-state functional connectivity patterns of the reported foci. Next, we compared the convergent circuit with the circuits targeted by transcranial magnetic stimulation (TMS) therapy. Last, we studied the association of regional and network-level convergence maps with selected neurotransmitter receptors/transporters maps.

Results

No regional convergence was found across foci of treatment-associated alterations in functional imaging. Subgroup analysis across the Treated > Untreated contrast revealed a convergent cluster in the left dorsolateral prefrontal cortex, which was associated with working memory and attention behavioral domains. Moreover, we found network-level convergence of the treatment-associated alterations in a circuit more prominent in the frontoparietal areas. This circuit was co-aligned with circuits targeted by "anti-subgenual" and "Beam F3" TMS therapy. We observed no significant correlations between our meta-analytic findings with the maps of neurotransmitter receptors/transporters.

Conclusion

Our findings highlight the importance of the frontoparietal network and the left dorsolateral prefrontal cortex in the therapeutic effects of antidepressants, which may relate to their role in improving executive functions and emotional processing.

A 2.5D Self-Training Strategy for Carotid Artery Segmentation in T1-Weighted Brain Magnetic Resonance Images

Precise annotations for large medical image datasets can be time-consuming. Additionally, when dealing with volumetric regions of interest, it is typical to apply segmentation techniques on 2D slices, compromising important information for accurately segmenting 3D structures. This study presents a deep learning pipeline that simultaneously tackles both challenges. Firstly, to streamline the annotation process, we employ a semi-automatic segmentation approach using bounding boxes as masks, which is less time-consuming than pixel-level delineation. Subsequently, recursive self-training is utilized to enhance annotation quality. Finally, a 2.5D segmentation technique is adopted, wherein a slice of a volumetric image is segmented using a pseudo-RGB image. The pipeline was applied to segment the carotid artery tree in T1-weighted brain magnetic resonance images. Utilizing 42 volumetric non-contrast T1-weighted brain scans from four datasets, we delineated bounding boxes around the carotid arteries in the axial slices. Pseudo-RGB images were generated from these slices, and recursive segmentation was conducted using a Res-Unet-based neural network architecture. The model's performance was tested on a separate dataset, with ground truth annotations provided by a radiologist. After recursive training, we achieved an Intersection over Union (IoU) score of (0.68 ± 0.08) on the unseen dataset, demonstrating commendable qualitative results.

Adolescent maturation of cortical excitation-inhibition balance based on individualized biophysical network modeling

The balance of excitation and inhibition is a key functional property of cortical microcircuits which changes through the lifespan. Adolescence is considered a crucial period for the maturation of excitation-inhibition balance. This has been primarily observed in animal studies, yet human in vivo evidence on adolescent maturation of the excitation-inhibition balance at the individual level is limited. Here, we developed an individualized in vivo marker of regional excitation-inhibition balance in human adolescents, estimated using large-scale simulations of biophysical network models fitted to resting-state functional magnetic resonance imaging data from two independent cross-sectional (N = 752) and longitudinal (N = 149) cohorts. We found a widespread relative increase of inhibition in association cortices paralleled by a relative age-related increase of excitation, or lack of change, in sensorimotor areas across both datasets. This developmental pattern co-aligned with multiscale markers of sensorimotor-association differentiation. The spatial pattern of excitation-inhibition development in adolescence was robust to inter-individual variability of structural connectomes and modeling configurations. Notably, we found that alternative simulation-based markers of excitation-inhibition balance show a variable sensitivity to maturational change. Taken together, our study highlights an increase of inhibition during adolescence in association areas using cross sectional and longitudinal data, and provides a robust computational framework to estimate microcircuit maturation in vivo at the individual level.

Disrupted topological properties of structural brain networks present a glutamatergic neuropathophysiology in people with narcolepsy

STUDY OBJECTIVES

Growing evidences have documented various abnormalities of the white matter bundles in people with narcolepsy. We sought to evaluate topological properties of brain structural networks, and their association with symptoms and neuropathophysiological features in people with narcolepsy.

METHODS

Diffusion tensor imaging was conducted for people with narcolepsy (n = 30) and matched healthy controls as well as symptoms assessment. Structural connectivity for each participant was generated to analyze global and regional topological properties and their correlations with narcoleptic features. Further human brain transcriptome was extracted and spatially registered for connectivity vulnerability. Genetic functional enrichment analysis was performed and further clarified using in vivo emission computed tomography data.

RESULTS

A wide and dramatic decrease in structural connectivities was observed in people with narcolepsy, with descending network degree and global efficiency. These metrics were not only correlated with sleep latency and awakening features, but also reflected alterations of sleep macrostructure in people with narcolepsy. Network-based statistics identified a small hyperenhanced subnetwork of cingulate gyrus that was closely related to rapid eye movement sleep behavior disorder (RBD) in narcolepsy. Further imaging genetics analysis suggested glutamatergic signatures were responsible for the preferential vulnerability of connectivity alterations in people with narcolepsy, while additional PET/SPECT data verified that structural alteration was significantly correlated with metabotropic glutamate receptor 5 (mGlutR5) and N-methyl-D-aspartate receptor (NMDA).

CONCLUSIONS

People with narcolepsy endured a remarkable decrease in the structural architecture, which was not only closely related to narcolepsy symptoms but also glutamatergic signatures.

Advances and challenges in immunoPET methodology

Immuno-positron emission tomography (immunoPET) enables imaging of specific targets that play a role in targeted therapy and immunotherapy, such as antigens on cell membranes, targets in the disease microenvironment, or immune cells. The most common immunoPET applications use a monoclonal antibody labeled with a relatively long-lived positron emitter such as 89Zr (T 1/2 = 78.4 h), but smaller antibody-based constructs labeled with various other positron emitting radionuclides are also being investigated. This molecular imaging technique can thus guide the development of new drugs and may have a pivotal role in selecting patients for a particular therapy. In early phase immunoPET trials, multiple imaging time points are used to examine the time-dependent biodistribution and to determine the optimal imaging time point, which may be several days after tracer injection due to the slow kinetics of larger molecules. Once this has been established, usually only one static scan is performed and semi-quantitative values are reported. However, total PET uptake of a tracer is the sum of specific and nonspecific uptake. In addition, uptake may be affected by other factors such as perfusion, pre-/co-administration of the unlabeled molecule, and the treatment schedule. This article reviews imaging methodologies used in immunoPET studies and is divided into two parts. The first part summarizes the vast majority of clinical immunoPET studies applying semi-quantitative methodologies. The second part focuses on a handful of studies applying pharmacokinetic models and includes preclinical and simulation studies. Finally, the potential and challenges of immunoPET quantification methodologies are discussed within the context of the recent technological advancements provided by long axial field of view PET/CT scanners.

Kinetic Modeling of Brain [18-F]FDG Positron Emission Tomography Time Activity Curves with Input Function Recovery (IR) Method

Accurate positron emission tomography (PET) data quantification relies on high-quality input plasma curves, but venous blood sampling may yield poor-quality data, jeopardizing modeling outcomes. In this study, we aimed to recover sub-optimal input functions by using information from the tail (5th-100th min) of curves obtained through the frequent sampling protocol and an input recovery (IR) model trained with reference curves of optimal shape. Initially, we included 170 plasma input curves from eight published studies with clamp [18F]-fluorodeoxyglucose PET exams. Model validation involved 78 brain PET studies for which compartmental model (CM) analysis was feasible (reference (ref) + training sets). Recovered curves were compared with original curves using area under curve (AUC), max peak standardized uptake value (maxSUV). CM parameters (ref + training sets) and fractional uptake rate (FUR) (all sets) were computed. Original and recovered curves from the ref set had comparable AUC (d = 0.02, not significant (NS)), maxSUV (d = 0.05, NS) and comparable brain CM results (NS). Recovered curves from the training set were different from the original according to maxSUV (d = 3) and biologically plausible according to the max theoretical K1 (53//56). Brain CM results were different in the training set (p < 0.05 for all CM parameters and brain regions) but not in the ref set. FUR showed reductions similarly in the recovered curves of the training and test sets compared to the original curves (p < 0.05 for all regions for both sets). The IR method successfully recovered the plasma inputs of poor quality, rescuing cases otherwise excluded from the kinetic modeling results. The validation approach proved useful and can be applied to different tracers and metabolic conditions.

An update on the use of image-derived input functions for human PET studies: new hopes or old illusions?

BACKGROUND

The need for arterial blood data in quantitative PET research limits the wider usability of this imaging method in clinical research settings. Image-derived input function (IDIF) approaches have been proposed as a cost-effective and non-invasive alternative to gold-standard arterial sampling. However, this approach comes with its own limitations-partial volume effects and radiometabolite correction among the most important-and varying rates of success, and the use of IDIF for brain PET has been particularly troublesome.

MAIN BODY

This paper summarizes the limitations of IDIF methods for quantitative PET imaging and discusses some of the advances that may make IDIF extraction more reliable. The introduction of automated pipelines (both commercial and open-source) for clinical PET scanners is discussed as a way to improve the reliability of IDIF approaches and their utility for quantitative purposes. Survey data gathered from the PET community are then presented to understand whether the field's opinion of the usefulness and validity of IDIF is improving. Finally, as the introduction of next-generation PET scanners with long axial fields of view, ultra-high sensitivity, and improved spatial and temporal resolution, has also brought IDIF methods back into the spotlight, a discussion of the possibilities offered by these state-of-the-art scanners-inclusion of large vessels, less partial volume in small vessels, better description of the full IDIF kinetics, whole-body modeling of radiometabolite production-is included, providing a pathway for future use of IDIF.

CONCLUSION

Improvements in PET scanner technology and software for automated IDIF extraction may allow to solve some of the major limitations associated with IDIF, such as partial volume effects and poor temporal sampling, with the exciting potential for accurate estimation of single kinetic rates. Nevertheless, until individualized radiometabolite correction can be performed effectively, IDIF approaches remain confined at best to a few tracers.

The biological role of local and global fMRI BOLD signal variability in human brain organization

Variability drives the organization and behavior of complex systems, including the human brain. Understanding the variability of brain signals is thus necessary to broaden our window into brain function and behavior. Few empirical investigations of macroscale brain signal variability have yet been undertaken, given the difficulty in separating biological sources of variance from artefactual noise. Here, we characterize the temporal variability of the most predominant macroscale brain signal, the fMRI BOLD signal, and systematically investigate its statistical, topographical and neurobiological properties. We contrast fMRI acquisition protocols, and integrate across histology, microstructure, transcriptomics, neurotransmitter receptor and metabolic data, fMRI static connectivity, and empirical and simulated magnetoencephalography data. We show that BOLD signal variability represents a spatially heterogeneous, central property of multi-scale multi-modal brain organization, distinct from noise. Our work establishes the biological relevance of BOLD signal variability and provides a lens on brain stochasticity across spatial and temporal scales.

Assortative mixing in micro-architecturally annotated brain connectomes

The wiring of the brain connects micro-architecturally diverse neuronal populations, but the conventional graph model, which encodes macroscale brain connectivity as a network of nodes and edges, abstracts away the rich biological detail of each regional node. Here, we annotate connectomes with multiple biological attributes and formally study assortative mixing in annotated connectomes. Namely, we quantify the tendency for regions to be connected based on the similarity of their micro-architectural attributes. We perform all experiments using four cortico-cortical connectome datasets from three different species, and consider a range of molecular, cellular, and laminar annotations. We show that mixing between micro-architecturally diverse neuronal populations is supported by long-distance connections and find that the arrangement of connections with respect to biological annotations is associated to patterns of regional functional specialization. By bridging scales of cortical organization, from microscale attributes to macroscale connectivity, this work lays the foundation for next-generation annotated connectomics.

Transitions between cognitive topographies: contributions of network structure, neuromodulation, and disease

Patterns of neural activity underlie human cognition. Transitions between these patterns are orchestrated by the brain's network architecture. What are the mechanisms linking network structure to cognitively relevant activation patterns? Here we implement principles of network control to investigate how the architecture of the human connectome shapes transitions between 123 experimentally defined cognitive activation maps (cognitive topographies) from the NeuroSynth meta-analytic engine. We also systematically incorporate neurotransmitter receptor density maps (18 receptors and transporters) and disease-related cortical abnormality maps (11 neurodegenerative, psychiatric and neurodevelopmental diseases; N = 17 000 patients, N = 22 000 controls). Integrating large-scale multimodal neuroimaging data from functional MRI, diffusion tractography, cortical morphometry, and positron emission tomography, we simulate how anatomically-guided transitions between cognitive states can be reshaped by pharmacological or pathological perturbation. Our results provide a comprehensive look-up table charting how brain network organisation and chemoarchitecture interact to manifest different cognitive topographies. This computational framework establishes a principled foundation for systematically identifying novel ways to promote selective transitions between desired cognitive topographies.

In Vivo N-Methyl-d-Aspartate Receptor (NMDAR) Density as Assessed Using Positron Emission Tomography During Recovery From NMDAR-Antibody Encephalitis

This case-control study uses a radiotracer and positron emission tomography to assess N-methyl-d-aspartate receptor (NMDAR) density changes during recovery from NMDAR-antibody encephalitis.

Non-invasive kinetic modelling approaches for quantitative analysis of brain PET studies

Pharmacokinetic modelling with arterial sampling is the gold standard for analysing dynamic PET data of the brain. However, the invasive character of arterial sampling prevents its widespread clinical application. Several methods have been developed to avoid arterial sampling, in particular reference region methods. Unfortunately, for some tracers or diseases, no suitable reference region can be defined. For these cases, other potentially non-invasive approaches have been proposed: (1) a population based input function (PBIF), (2) an image derived input function (IDIF), or (3) simultaneous estimation of the input function (SIME). This systematic review aims to assess the correspondence of these non-invasive methods with the gold standard. Studies comparing non-invasive pharmacokinetic modelling methods with the current gold standard methods using an input function derived from arterial blood samples were retrieved from PubMed/MEDLINE (until December 2021). Correlation measurements were extracted from the studies. The search yielded 30 studies that correlated outcome parameters (V_T, DVR, or BP_ND for reversible tracers; K_i or CMR_glu for irreversible tracers) from a potentially non-invasive method with those obtained from modelling using an arterial input function. Some studies provided similar results for PBIF, IDIF, and SIME-based methods as for modelling with an arterial input function (R² = 0.59-1.00, R² = 0.71-1.00, R² = 0.56-0.96, respectively), if the non-invasive input curve was calibrated with arterial blood samples. Even when the non-invasive input curve was calibrated with venous blood samples or when no calibration was applied, moderate to good correlations were reported, especially for the IDIF and SIME (R² = 0.71-1.00 and R² = 0.36-0.96, respectively). Overall, this systematic review illustrates that non-invasive methods to generate an input function are still in their infancy. Yet, IDIF and SIME performed well, not only with arterial blood calibration, but also with venous or no blood calibration, especially for some tracers without plasma metabolites, which would potentially make these methods better suited for clinical application. However, these methods should still be properly validated for each individual tracer and application before implementation.

Mapping neurotransmitter systems to the structural and functional organization of the human neocortex

Neurotransmitter receptors support the propagation of signals in the human brain. How receptor systems are situated within macro-scale neuroanatomy and how they shape emergent function remain poorly understood, and there exists no comprehensive atlas of receptors. Here we collate positron emission tomography data from more than 1,200 healthy individuals to construct a whole-brain three-dimensional normative atlas of 19 receptors and transporters across nine different neurotransmitter systems. We found that receptor profiles align with structural connectivity and mediate function, including neurophysiological oscillatory dynamics and resting-state hemodynamic functional connectivity. Using the Neurosynth cognitive atlas, we uncovered a topographic gradient of overlapping receptor distributions that separates extrinsic and intrinsic psychological processes. Finally, we found both expected and novel associations between receptor distributions and cortical abnormality patterns across 13 disorders. We replicated all findings in an independently collected autoradiography dataset. This work demonstrates how chemoarchitecture shapes brain structure and function, providing a new direction for studying multi-scale brain organization.

Cerebral Oxidative Stress in Early Alzheimer's Disease Evaluated by 64Cu-ATSM PET/MRI: A Preliminary Study

Oxidative stress imaging using diacetyl-bis (N4-methylthiosemicarbazone) (Cu-ATSM) was applied to the evaluation of patients with early Alzheimer's disease (eAD). Ten eAD patients (72 ± 9 years) and 10 age-matched healthy controls (HCs) (73 ± 9 years) participated in this study. They underwent dynamic PET/MRI using 11C-PiB and 64Cu-ATSM with multiple MRI sequences. To evaluate cerebral oxidative stress, three parameters of 64Cu-ATSM PET were compared: standardized uptake value (SUV), tracer influx rate (Kin), and a rate constant k3. The input functions were estimated by the image-derived input function method. The relative differences were analyzed by statistical parametric mapping (SPM) using SUV and Kin images. All eAD patients had positive and HC subjects had negative PiB accumulation, and MMSE scores were significantly different between them. The 64Cu-ATSM accumulation tended to be higher in eAD than in HCs for both SUV and Kin. When comparing absolute values, eAD patients had a greater Kin in the posterior cingulate cortex and a greater k3 in the hippocampus compared with lobar cortical values of HCs. In SPM analysis, eAD had an increased left operculum and decreased bilateral hippocampus and anterior cingulate cortex compared to HCs. 64Cu-ATSM PET/MRI and tracer kinetic analysis elucidated cerebral oxidative stress in the eAD patients, particularly in the cingulate cortex and hippocampus.

Profiling Transcription Initiation in Peripheral Leukocytes Reveals Severity-Associated Cis-Regulatory Elements in Critical COVID-19

The contribution of transcription factors (TFs) and gene regulatory programs in the immune response to COVID-19 and their relationship to disease outcome is not fully understood. Analysis of genome-wide changes in transcription at both promoter-proximal and distal cis-regulatory DNA elements, collectively termed the 'active cistrome,' offers an unbiased assessment of TF activity identifying key pathways regulated in homeostasis or disease. Here, we profiled the active cistrome from peripheral leukocytes of critically ill COVID-19 patients to identify major regulatory programs and their dynamics during SARS-CoV-2 associated acute respiratory distress syndrome (ARDS). We identified TF motifs that track the severity of COVID- 19 lung injury, disease resolution, and outcome. We used unbiased clustering to reveal distinct cistrome subsets delineating the regulation of pathways, cell types, and the combinatorial activity of TFs. We found critical roles for regulatory networks driven by stimulus and lineage determining TFs, showing that STAT and E2F/MYB regulatory programs targeting myeloid cells are activated in patients with poor disease outcomes and associated with single nucleotide genetic variants implicated in COVID-19 susceptibility. Integration with single-cell RNA-seq found that STAT and E2F/MYB activation converged in specific neutrophils subset found in patients with severe disease. Collectively we demonstrate that cistrome analysis facilitates insight into disease mechanisms and provides an unbiased approach to evaluate global changes in transcription factor activity and stratify patient disease severity.