ABSTRACT
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.
ABSTRACT
BackgroundIncreased inflammation is a hallmark of COVID-19, with pulmonary and systemic inflammation identified in multiple cohorts of patients. Definitive cellular and molecular pathways driving severe forms of this disease remain uncertain. Neutrophils, the most numerous leukocytes in blood circulation, can contribute to immunopathology in infections, inflammatory diseases and acute respiratory distress syndrome (ARDS), a primary cause of morbidity and mortality in COVID-19. Changes in multiple neutrophil functions and circulating cytokine levels over time during COVID-19 may help define disease severity and guide care and decision making. MethodsBlood was obtained serially from critically ill COVID-19 patients for 11 days. Neutrophil oxidative burst, neutrophil extracellular trap formation (NETosis), phagocytosis and cytokine levels were assessed ex vivo. Lung tissue was obtained immediately post-mortem for immunostaining. ResultsElevations in neutrophil-associated cytokines IL-8 and IL-6, and general inflammatory cytokines IP-10, GM-CSF, IL-1b, IL-10 and TNF, were identified in COVID-19 plasma both at the first measurement and at multiple timepoints across hospitalization (p < 0.0001). Neutrophils had exaggerated oxidative burst (p < 0.0001), NETosis (p < 0.0001) and phagocytosis (p < 0.0001) relative to controls. Increased NETosis correlated with both leukocytosis and neutrophilia. Neutrophils and NETs were identified within airways and alveoli in the lung parenchyma of 40% of SARS-CoV-2 infected lungs. While elevations in IL-8 and ANC correlated to COVID-19 disease severity, plasma IL-8 levels alone correlated with death. ConclusionsCirculating neutrophils in COVID-19 exhibit an activated phenotype with increased oxidative burst, NETosis and phagocytosis. Readily accessible and dynamic, plasma IL-8 and circulating neutrophil function may be potential COVID-19 disease biomarkers.