The immune response to COVID-19 infection is variable. How COVID-19 influences clinical outcomes in hospitalized patients needs to be understood through readily obtainable biological materials, such as blood. We hypothesized that a high-density analysis of host (and pathogen) blood RNA in hospitalized patients with SARS-CoV-2 would provide mechanistic insights into the heterogeneity of response amongst COVID-19 patients when combined with advanced multidimensional bioinformatics for RNA. We enrolled 36 hospitalized COVID-19 patients (11 died) and 15 controls, collecting 74 blood PAXgene RNA tubes at multiple timepoints, one early and in 23 patients after treatment with various therapies. Total RNAseq was performed at high-density, with >160 million paired-end, 150 base pair reads per sample, representing the most sequenced bases per sample for any publicly deposited blood PAXgene tube study. There are 770 genes significantly altered in the blood of COVID-19 patients associated with antiviral defense, mitotic cell cycle, type I interferon signaling, and severe viral infections. Immune genes activated include those associated with neutrophil mechanisms, secretory granules, and neutrophil extracellular traps (NETs), along with decreased gene expression in lymphocytes and clonal expansion of the acquired immune response. Therapies such as convalescent serum and dexamethasone reduced many of the blood expression signatures of COVID-19. Severely ill or deceased patients are marked by various secondary infections, unique gene patterns, dysregulated innate response, and peripheral organ damage not otherwise found in the cohort. High-density transcriptomic data offers shared gene expression signatures, providing unique insights into the immune system and individualized signatures of patients that could be used to understand the patient’s clinical condition. Whole blood transcriptomics provides patient-level insights for immune activation, immune repertoire, and secondary infections that can further guide precision treatment.

Ornithine decarboxylase 1 (ODC1 gene) has been linked through gain-of-function variants to a rare disease featuring developmental delay, alopecia, macrocephaly, and structural brain anomalies. ODC1 has been linked to additional diseases like cancer, with growing evidence for neurological contributions to schizophrenia, mood disorders, anxiety, epilepsy, learning, and suicidal behavior. The evidence of ODC1 connection to neural disorders highlights the need for a systematic analysis of ODC1 genotype-to-phenotype associations. An analysis of variants from ClinVar, Geno2MP, TOPMed, gnomAD, and COSMIC revealed an intellectual disability and seizure connected loss-of-function variant, ODC G84R (rs138359527, NC_000002.12:g.10444500C > T). The missense variant is found in ~1% of South Asian individuals and results in 2.5-fold decrease in enzyme function. Expression quantitative trait loci (eQTLs) reveal multiple functionally annotated, non-coding variants regulating ODC1 that associate with psychiatric/neurological phenotypes. Further dissection of RNA-Seq during fetal brain development and within cerebral organoids showed an association of ODC1 expression with cell proliferation of neural progenitor cells, suggesting gain-of-function variants with neural over-proliferation and loss-of-function variants with neural depletion. The linkage from the expression data of ODC1 in early neural progenitor proliferation to phenotypes of neurodevelopmental delay and to the connection of polyamine metabolites in brain function establish ODC1 as a bona fide neurodevelopmental disorder gene.

Precision medicine requires the translation of basic biological understanding to medical insights, mainly applied to characterization of each unique patient. In many clinical settings, this requires tools that can be broadly used to identify pathology and risks. Patients often present to the intensive care unit with broad phenotypes, including multiple organ dysfunction syndrome (MODS) resulting from infection, trauma, or other disease processes. Etiology and outcomes are unique to individuals, making it difficult to cohort patients with MODS, but presenting a prime target for testing/developing tools for precision medicine. Using multitime point whole blood (cellular/acellular) total transcriptomics in 27 patients, we highlight the promise of simultaneously mapping viral/bacterial load, cell composition, tissue damage biomarkers, balance between syndromic biology versus environmental response, and unique biological insights in each patient using a single platform measurement. Integration of a transcriptome workflow yielded unexpected insights into the complex interplay between host genetics and viral/bacterial specific mechanisms, highlighted by a unique case of virally induced genetics (VIG) within one of these 27 patients. The power of RNA-Seq to study unique patient biology while investigating environmental contributions can be a critical tool moving forward for translational sciences applied to precision medicine.

Multiple organ dysfunction syndrome (MODS) occurs in the setting of a variety of pathologies including infection and trauma. Some patients decompensate and require Veno-Arterial extra corporeal membrane oxygenation (ECMO) as a palliating manoeuvre for recovery of cardiopulmonary function. The molecular mechanisms driving progression from MODS to cardiopulmonary collapse remain incompletely understood, and no biomarkers have been defined to identify those MODS patients at highest risk for progression to requiring ECMO support. Whole blood RNA-seq profiling was performed for 23 MODS patients at three time points during their ICU stay (at diagnosis of MODS, 72 hours after, and 8 days later), as well as four healthy controls undergoing routine sedation. Of the 23 MODS patients, six required ECMO support (ECMO patients). The predictive power of conventional demographic and clinical features was quantified for differentiating the MODS and ECMO patients. We then compared the performance of markers derived from transcriptomic profiling including [1] transcriptomically imputed leukocyte subtype distribution, [2] relevant published gene signatures and [3] a novel differential gene expression signature computed from our data set. The predictive power of our novel gene expression signature was then validated using independently published datasets. None of the five demographic characteristics and 14 clinical features, including The Paediatric Logistic Organ Dysfunction (PELOD) score, could predict deterioration of MODS to ECMO at baseline. From previously published sepsis signatures, only the signatures positively associated with patient's mortality could differentiate ECMO patients from MODS patients, when applied to our transcriptomic dataset (P-value ranges from 0.01 to 0.04, Student's test). Deconvolution of bulk RNA-Seq samples suggested that lower neutrophil counts were associated with increased risk of progression from MODS to ECMO (P-value = 0.03, logistic regression, OR=2.82 [95% CI 0.63 - 12.45]). A total of 30 genes were differentially expressed between ECMO and MODS patients at baseline (log2 fold change ≥ 1 or ≤ -1 with false discovery rate ≤ 0.01). These genes are involved in protein maintenance and epigenetic-related processes. Further univariate analysis of these 30 genes suggested a signature of seven DE genes associated with ECMO (OR > 3.0, P-value ≤ 0.05, logistic regression). Notably, this contains a set of histone marker genes, including H1F0, HIST2H3C, HIST1H2AI, HIST1H4, HIST1H2BL and HIST1H1B, that were highly expressed in ECMO. A risk score derived from expression of these genes differentiated ECMO and MODS patients in our dataset (AUC = 0.91, 95% CI 0.79-1.00, P-value = 7e-04, logistic regression) as well as validation dataset (AUC= 0.73, 95% CI 0.53-0.93, P-value = 2e-02, logistic regression). This study demonstrates that transcriptomic features can serve as indicators of severity that could be superior to traditional methods of ascertaining acuity in MODS patients. Analysis of expression of signatures identified in this study could help clinicians in the diagnosis and prognostication of MODS patients after arrival to the Hospital.

Viral infections affecting the lower respiratory tract place enormous burdens on hospitals. As neither vaccines nor antiviral agents exist for many viruses, understanding risk factors and outcomes in each patient using minimally invasive analysis, such as blood, can lead to improved health care delivery. A cohort of PAXgene RNA sequencing of infants admitted with moderate or severe acute bronchiolitis and respiratory syncytial virus were compared with case-control statistical analysis and cohort-based outlier mapping for precision transcriptomics. Patients with severe bronchiolitis had signatures connected to the immune system, interferon signaling, and cytokine signaling, with marked sex differences in XIST, RPS4Y1, KDM5D, and LINC00278 for severity. Several patients had unique secondary infections, cytokine activation, immune responses, biological pathways, and immune cell activation, highlighting the need for defining patient-level transcriptomic signatures. Balancing relative contributions of cohort-based biomarker discoveries with patient’s biological responses is needed to understand the totality of mechanisms of adverse outcomes in viral bronchiolitis.