Dysregulated transcriptional responses to SARS-CoV-2 in the periphery

Deep-HPI-pred: An R-Shiny applet for network-based classification and prediction of Host-Pathogen protein-protein interactions

Host-pathogen interactions (HPIs) are vital in numerous biological activities and are intrinsically linked to the onset and progression of infectious diseases. HPIs are pivotal in the entire lifecycle of diseases: from the onset of pathogen introduction, navigating through the mechanisms that bypass host cellular defenses, to its subsequent proliferation inside the host. At the heart of these stages lies the synergy of proteins from both the host and the pathogen. By understanding these interlinking protein dynamics, we can gain crucial insights into how diseases progress and pave the way for stronger plant defenses and the swift formulation of countermeasures. In the framework of current study, we developed a web-based R/Shiny app, Deep-HPI-pred, that uses network-driven feature learning method to predict the yet unmapped interactions between pathogen and host proteins. Leveraging citrus and CLas bacteria training datasets as case study, we spotlight the effectiveness of Deep-HPI-pred in discerning Protein-protein interaction (PPIs) between them. Deep-HPI-pred use Multilayer Perceptron (MLP) models for HPI prediction, which is based on a comprehensive evaluation of topological features and neural network architectures. When subjected to independent validation datasets, the predicted models consistently surpassed a Matthews correlation coefficient (MCC) of 0.80 in host-pathogen interactions. Remarkably, the use of Eigenvector Centrality as the leading topological feature further enhanced this performance. Further, Deep-HPI-pred also offers relevant gene ontology (GO) term information for each pathogen and host protein within the system. This protein annotation data contributes an additional layer to our understanding of the intricate dynamics within host-pathogen interactions. In the additional benchmarking studies, the Deep-HPI-pred model has proven its robustness by consistently delivering reliable results across different host-pathogen systems, including plant-pathogens (accuracy of 98.4% and 97.9%), human-virus (accuracy of 94.3%), and animal-bacteria (accuracy of 96.6%) interactomes. These results not only demonstrate the model's versatility but also pave the way for gaining comprehensive insights into the molecular underpinnings of complex host-pathogen interactions. Taken together, the Deep-HPI-pred applet offers a unified web service for both identifying and illustrating interaction networks. Deep-HPI-pred applet is freely accessible at its homepage: https://cbi.gxu.edu.cn/shiny-apps/Deep-HPI-pred/ and at github: https://github.com/tahirulqamar/Deep-HPI-pred.

Exploring COVID-19 Pandemic Disparities with Transcriptomic Meta-analysis from the Perspective of Personalized Medicine

Infection with SARS-CoV2, which is responsible for COVID-19, can lead to differences in disease development, severity and mortality rates depending on gender, age or the presence of certain diseases. Considering that existing studies ignore these differences, this study aims to uncover potential differences attributable to gender, age and source of sampling as well as viral load using bioinformatics and multi-omics approaches. Differential gene expression analyses were used to analyse the phenotypic differences between SARS-CoV-2 patients and controls at the mRNA level. Pathway enrichment analyses were performed at the gene set level to identify the activated pathways corresponding to the differences in the samples. Drug repurposing analysis was performed at the protein level, focusing on host-mediated drug candidates to uncover potential therapeutic differences. Significant differences (i.e. the number of differentially expressed genes and their characteristics) were observed for COVID-19 at the mRNA level depending on the sample source, gender and age of the samples. The results of the pathway enrichment show that SARS-CoV-2 can be combated more effectively in the respiratory tract than in the blood samples. Taking into account the different sample sources and their characteristics, different drug candidates were identified. Evaluating disease prediction, prevention and/or treatment strategies from a personalised perspective is crucial. In this study, we not only evaluated the differences in COVID-19 from a personalised perspective, but also provided valuable data for further experimental and clinical efforts. Our findings could shed light on potential pandemics.

Whole blood transcriptome signature predicts severe forms of COVID-19: Results from the COVIDeF cohort study

COVID-19 is associated with heterogeneous outcome. Early identification of a severe progression of the disease is essential to properly manage the patients and improve their outcome. Biomarkers reflecting an increased inflammatory response, as well as individual features including advanced age, male gender, and pre-existing comorbidities, are risk factors of severe COVID-19. Yet, these features show limited accuracy for outcome prediction. The aim was to evaluate the prognostic value of whole blood transcriptome at an early stage of the disease. Blood transcriptome of patients with mild pneumonia was profiled. Patients with subsequent severe COVID-19 were compared to those with favourable outcome, and a molecular predictor based on gene expression was built. Unsupervised classification discriminated patients who would later develop a COVID-19-related severe pneumonia. The corresponding gene expression signature reflected the immune response to the viral infection dominated by a prominent type I interferon, with IFI27 among the most over-expressed genes. A 48-genes transcriptome signature predicting the risk of severe COVID-19 was built on a training cohort, then validated on an external independent cohort, showing an accuracy of 81% for predicting severe outcome. These results identify an early transcriptome signature of severe COVID-19 pneumonia, with a possible relevance to improve COVID-19 patient management.

Quantitative and Longitudinal Assessment of Systemic Innate Immunity in Health and Disease Using a 2D Gene Model

Dysregulation of innate immunity is deeply involved in infectious and autoimmune diseases. For a better understanding of pathogenesis and improved management of these diseases, it is of vital importance to implement convenient monitoring of systemic innate immunity. Built upon our previous works on the host transcriptional response to infection in peripheral blood, we proposed a 2D gene model for the simultaneous assessment of two major components of systemic innate immunity, including VirSig as the signature of the host response to viral infection and BacSig as the signature of the host response to bacterial infection. The revelation of dysregulation in innate immunity by this 2D gene model was demonstrated with a wide variety of transcriptome datasets. In acute infection, distinctive patterns of VirSig and BacSig activation were observed in viral and bacterial infection. In comparison, both signatures were restricted to a defined range in the vast majority of healthy adults, regardless of age. In addition, BacSig showed significant elevation during pregnancy and an upward trend during development. In tuberculosis (TB), elevation of BacSig and VirSig was observed in a significant portion of active TB patients, and abnormal BacSig was also associated with a longer treatment course. In cystic fibrosis (CF), abnormal BacSig was observed in a subset of patients, and no overall change in BacSig abnormality was observed after the drug treatment. In systemic sclerosis-associated interstitial lung disease (SSc-ILD), significant elevation of VirSig and BacSig was observed in some patients, and treatment with a drug led to the further deviation of BacSig from the control level. In systemic lupus erythematosus (SLE), positivity for the anti-Ro autoantibody was associated with significant elevation of VirSig in SLE patients, and the additive effect of VirSig/BacSig activation was also observed in SLE patients during pregnancy. Overall, these data demonstrated that the 2D gene model can be used to assess systemic innate immunity in health and disease, with the potential clinical applications including patient stratification, prescription of antibiotics, understanding of pathogenesis, and longitudinal monitoring of treatment response.

Integrated analysis of RNA-seq datasets reveals novel targets and regulators of COVID-19 severity

During the COVID-19 pandemic, RNA-seq datasets were produced to investigate the virus-host relationship. However, much of these data remains underexplored. To improve the search for molecular targets and biomarkers, we performed an integrated analysis of multiple RNA-seq datasets, expanding the cohort and including patients from different countries, encompassing severe and mild COVID-19 patients. Our analysis revealed that severe COVID-19 patients exhibit overexpression of genes coding for proteins of extracellular exosomes, endomembrane system, and neutrophil granules (e.g., S100A9, LY96, and RAB1B), which may play an essential role in the cellular response to infection. Concurrently, these patients exhibit down-regulation of genes encoding components of the T cell receptor complex and nucleolus, including TP53, IL2RB, and NCL Finally, SPI1 may emerge as a central transcriptional factor associated with the up-regulated genes, whereas TP53, MYC, and MAX were associated with the down-regulated genes during COVID-19. This study identified targets and transcriptional factors, lighting on the molecular pathophysiology of syndrome coronavirus 2 infection.

Consequences of COVID-19 on Adipose Tissue Signatures

Since the emergence of coronavirus disease-19 (COVID-19) in 2019, it has been crucial to investigate the causes of severe cases, particularly the higher rates of hospitalization and mortality in individuals with obesity. Previous findings suggest that adipocytes may play a role in adverse COVID-19 outcomes in people with obesity. The impact of COVID-19 vaccination and infection on adipose tissue (AT) is currently unclear. We therefore analyzed 27 paired biopsies of visceral and subcutaneous AT from donors of the Leipzig Obesity BioBank that have been categorized into three groups (1: no infection/no vaccination; 2: no infection but vaccinated; 3: infected and vaccinated) based on COVID-19 antibodies to spike (indicating vaccination) and/or nucleocapsid proteins. We provide additional insights into the impact of COVID-19 on AT biology through a comprehensive histological transcriptome and serum proteome analysis. This study demonstrates that COVID-19 infection is associated with smaller average adipocyte size. The impact of infection on gene expression was significantly more pronounced in subcutaneous than in visceral AT and mainly due to immune system-related processes. Serum proteome analysis revealed the effects of the infection on circulating adiponectin, interleukin 6 (IL-6), and carbonic anhydrase 5A (CA5A), which are all related to obesity and blood glucose abnormalities.

Hypoxia and Activation of Neutrophil Degranulation-Related Genes in the Peripheral Blood of COVID-19 Patients

Severe COVID-19 is characterized by systematic hyper-inflammation and subsequent damage to various organs. Therefore, it is critical to trace this cascade of hyper-inflammation. Blood transcriptome has been routinely utilized in the interrogation of host immune response in COVID-19 and other infectious conditions. In this study, consensus gene dysregulation in the blood was obtained from 13 independent transcriptome studies on COVID-19. Among the up-regulated genes, the most prominent functional categories were neutrophil degranulation and cell cycle, which is clearly different from the classical activation of interferon signaling pathway in seasonal flu. As for the potential upstream causal factors of the atypical gene dysregulation, systemic hypoxia was further examined because it is much more widely reported in COVID-19 than that in seasonal flu. It was found that both physiological and pathological hypoxia can induce activation of neutrophil degranulation-related genes in the blood. Furthermore, COVID-19 patients with different requirement for oxygen intervention showed distinctive levels of gene expression related to neutrophil degranulation in the whole blood, which was validated in isolated neutrophils. Thus, activation of neutrophil degranulation-related genes in the blood of COVID-19 could be partially attributed to hypoxia. Interestingly, similar pattern was also observed in H1N1 infection (the cause of Spanish flu) and several other severe respiratory viral infections. As for the molecular mechanism, both HIF-dependent and HIF-independent pathways have been examined. Since the activation of neutrophil degranulation-related genes is highly correlated with disease severity in COVID-19, early detection of hypoxia and active intervention may prevent further activation of neutrophil degranulation-related genes and other harmful downstream hyper-inflammation. This common mechanism is applicable to current and future pandemic as well as the severe form of common respiratory infection.

The 'analysis of gene expression and biomarkers for point-of-care decision support in Sepsis' study; temporal clinical parameter analysis and validation of early diagnostic biomarker signatures for severe inflammation andsepsis-SIRS discrimination

Introduction

Early diagnosis of sepsis and discrimination from SIRS is crucial for clinicians to provide appropriate care, management and treatment to critically ill patients. We describe identification of mRNA biomarkers from peripheral blood leukocytes, able to identify severe, systemic inflammation (irrespective of origin) and differentiate Sepsis from SIRS, in adult patients within a multi-center clinical study.

Methods

Participants were recruited in Intensive Care Units (ICUs) from multiple UK hospitals, including fifty-nine patients with abdominal sepsis, eighty-four patients with pulmonary sepsis, forty-two SIRS patients with Out-of-Hospital Cardiac Arrest (OOHCA), sampled at four time points, in addition to thirty healthy control donors. Multiple clinical parameters were measured, including SOFA score, with many differences observed between SIRS and sepsis groups. Differential gene expression analyses were performed using microarray hybridization and data analyzed using a combination of parametric and non-parametric statistical tools.

Results

Nineteen high-performance, differentially expressed mRNA biomarkers were identified between control and combined SIRS/Sepsis groups (FC>20.0, p<0.05), termed 'indicators of inflammation' (I°I), including CD177, FAM20A and OLAH. Best-performing minimal signatures e.g. FAM20A/OLAH showed good accuracy for determination of severe, systemic inflammation (AUC>0.99). Twenty entities, termed 'SIRS or Sepsis' (S°S) biomarkers, were differentially expressed between sepsis and SIRS (FC>2·0, p-value<0.05).

Discussion

The best performing signature for discriminating sepsis from SIRS was CMTM5/CETP/PLA2G7/MIA/MPP3 (AUC=0.9758). The I°I and S°S signatures performed variably in other independent gene expression datasets, this may be due to technical variation in the study/assay platform.

Epigenetic and transcriptional responses in circulating leukocytes are associated with future decompensation during SARS-CoV-2 infection

To elucidate host response elements that define impending decompensation during SARS-CoV-2 infection, we enrolled subjects hospitalized with COVID-19 who were matched for disease severity and comorbidities at the time of admission. We performed combined single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) on peripheral blood mononuclear cells (PBMCs) at admission and compared subjects who improved from their moderate disease with those who later clinically decompensated and required invasive mechanical ventilation or died. Chromatin accessibility and transcriptomic immune profiles were markedly altered between the two groups, with strong signals in CD4+ T cells, inflammatory T cells, dendritic cells, and NK cells. Multiomic signature scores at admission were tightly associated with future clinical deterioration (auROC 1.0). Epigenetic and transcriptional changes in PBMCs reveal early, broad immune dysregulation before typical clinical signs of decompensation are apparent and thus may act as biomarkers to predict future severity in COVID-19.

Interferon response and profiling of interferon response genes in peripheral blood of vaccine-naive COVID-19 patients

Introduction

There is insufficient understanding on systemic interferon (IFN) responses during COVID-19 infection. Early reports indicated that interferon responses were suppressed by the coronavirus (SARS-CoV-2) and clinical trials of administration of various kinds of interferons had been disappointing. Expression of interferon-stimulated genes (ISGs) in peripheral blood (better known as interferon score) has been a well-established bioassay marker of systemic IFN responses in autoimmune diseases. Therefore, with archival samples of a cohort of COVID-19 patients collected before the availability of vaccination, we aimed to better understand this innate immune response by studying the IFN score and related ISGs expression in bulk and single cell RNAs sequencing expression datasets.

Methods

In this study, we recruited 105 patients with COVID-19 and 30 healthy controls in Hong Kong. Clinical risk factors, disease course, and blood sampling times were recovered. Based on a set of five commonly used ISGs (IFIT1, IFIT2, IFI27, SIGLEC1, IFI44L), the IFN score was determined in blood leukocytes collected within 10 days after onset. The analysis was confined to those blood samples collected within 10 days after disease onset. Additional public datasets of bulk gene and single cell RNA sequencing of blood samples were used for the validation of IFN score results.

Results

Compared to the healthy controls, we showed that ISGs expression and IFN score were significantly increased during the first 10 days after COVID infection in majority of patients (71%). Among those low IFN responders, they were more commonly asymptomatic patients (71% vs 25%). 22 patients did not mount an overall significant IFN response and were classified as low IFN responders (IFN score < 1). However, early IFN score or ISGs level was not a prognostic biomarker and could not predict subsequent disease severity. Both IFI27 and SIGLEC1 were monocyte-predominant expressing ISGs and IFI27 were activated even among those low IFN responders as defined by IFN score. In conclusion, a substantial IFN response was documented in this cohort of COVID-19 patients who experience a natural infection before the vaccination era. Like innate immunity towards other virus, the ISGs activation was observed largely during the early course of infection (before day 10). Single-cell RNA sequencing data suggested monocytes were the cell-type that primarily accounted for the activation of two highly responsive ISGs (IFI44L and IFI27).

Discussion

As sampling time and age were two major confounders of ISG expression, they may account for contradicting observations among previous studies. On the other hand, the IFN score was not associated with the severity of the disease.

Analysis of blood and nasal epithelial transcriptomes to identify mechanisms associated with control of SARS-CoV-2 viral load in the upper respiratory tract

OBJECTIVES

The amount of SARS-CoV-2 detected in the upper respiratory tract (URT viral load) is a key driver of transmission of infection. Current evidence suggests that mechanisms constraining URT viral load are different from those controlling lower respiratory tract viral load and disease severity. Understanding such mechanisms may help to develop treatments and vaccine strategies to reduce transmission. Combining mathematical modelling of URT viral load dynamics with transcriptome analyses we aimed to identify mechanisms controlling URT viral load.

METHODS

COVID-19 patients were recruited in Spain during the first wave of the pandemic. RNA sequencing of peripheral blood and targeted NanoString nCounter transcriptome analysis of nasal epithelium were performed and gene expression analysed in relation to paired URT viral load samples collected within 15 days of symptom onset. Proportions of major immune cells in blood were estimated from transcriptional data using computational differential estimation. Weighted correlation network analysis (adjusted for cell proportions) and fixed transcriptional repertoire analysis were used to identify associations with URT viral load, quantified as standard deviations (z-scores) from an expected trajectory over time.

RESULTS

Eighty-two subjects (50% female, median age 54 years (range 3-73)) with COVID-19 were recruited. Paired URT viral load samples were available for 16 blood transcriptome samples, and 17 respiratory epithelial transcriptome samples. Natural Killer (NK) cells were the only blood cell type significantly correlated with URT viral load z-scores (r = -0.62, P = 0.010). Twenty-four blood gene expression modules were significantly correlated with URT viral load z-score, the most significant being a module of genes connected around IFNA14 (Interferon Alpha-14) expression (r = -0.60, P = 1e-10). In fixed repertoire analysis, prostanoid-related gene expression was significantly associated with higher viral load. In nasal epithelium, only GNLY (granulysin) gene expression showed significant negative correlation with viral load.

CONCLUSIONS

Correlations between the transcriptional host response and inter-individual variations in SARS-CoV-2 URT viral load, revealed many molecular mechanisms plausibly favouring or constraining viral replication. Existing evidence corroborates many of these mechanisms, including likely roles for NK cells, granulysin, prostanoids and interferon alpha-14. Inhibition of prostanoid production and administration of interferon alpha-14 may be attractive transmission-blocking interventions.

The two-stage molecular scenery of SARS-CoV-2 infection with implications to disease severity: An in-silico quest

Introduction

The two-stage molecular profile of the progression of SARS-CoV-2 (SCOV2) infection is explored in terms of five key biological/clinical questions: (a) does SCOV2 exhibits a two-stage infection profile? (b) SARS-CoV-1 (SCOV1) vs. SCOV2: do they differ? (c) does and how SCOV2 differs from Influenza/INFL infection? (d) does low viral-load and (e) does COVID-19 early host response relate to the two-stage SCOV2 infection profile? We provide positive answers to the above questions by analyzing the time-series gene-expression profiles of preserved cell-lines infected with SCOV1/2 or, the gene-expression profiles of infected individuals with different viral-loads levels and different host-response phenotypes.

Methods

Our analytical methodology follows an in-silico quest organized around an elaborate multi-step analysis pipeline including: (a) utilization of fifteen gene-expression datasets from NCBI's gene expression omnibus/GEO repository; (b) thorough designation of SCOV1/2 and INFL progression stages and COVID-19 phenotypes; (c) identification of differentially expressed genes (DEGs) and enriched biological processes and pathways that contrast and differentiate between different infection stages and phenotypes; (d) employment of a graph-based clustering process for the induction of coherent groups of networked genes as the representative core molecular fingerprints that characterize the different SCOV2 progression stages and the different COVID-19 phenotypes. In addition, relying on a sensibly selected set of induced fingerprint genes and following a Machine Learning approach, we devised and assessed the performance of different classifier models for the differentiation of acute respiratory illness/ARI caused by SCOV2 or other infections (diagnostic classifiers), as well as for the prediction of COVID-19 disease severity (prognostic classifiers), with quite encouraging results.

Results

The central finding of our experiments demonstrates the down-regulation of type-I interferon genes (IFN-1), interferon induced genes (ISGs) and fundamental innate immune and defense biological processes and molecular pathways during the early SCOV2 infection stages, with the inverse to hold during the later ones. It is highlighted that upregulation of these genes and pathways early after infection may prove beneficial in preventing subsequent uncontrolled hyperinflammatory and potentially lethal events.

Discussion

The basic aim of our study was to utilize in an intuitive, efficient and productive way the most relevant and state-of-the-art bioinformatics methods to reveal the core molecular mechanisms which govern the progression of SCOV2 infection and the different COVID-19 phenotypes.

Transcriptome-based analysis of human peripheral blood reveals regulators of immune response in different viral infections

Introduction

There are difficulties in creating direct antiviral drugs for all viruses, including new, suddenly arising infections, such as COVID-19. Therefore, pathogenesis-directed therapy is often necessary to treat severe viral infections and comorbidities associated with them. Despite significant differences in the etiopathogenesis of viral diseases, in general, they are associated with significant dysfunction of the immune system. Study of common mechanisms of immune dysfunction caused by different viral infections can help develop novel therapeutic strategies to combat infections and associated comorbidities.

Methods

To identify common mechanisms of immune functions disruption during infection by nine different viruses (cytomegalovirus, Ebstein-Barr virus, human T-cell leukemia virus type 1, Hepatitis B and C viruses, human immunodeficiency virus, Dengue virus, SARS-CoV, and SARS-CoV-2), we analyzed the corresponding transcription profiles from peripheral blood mononuclear cells (PBMC) using the originally developed pipeline that include transcriptome data collection, processing, normalization, analysis and search for master regulators of several viral infections. The ten datasets containing transcription data from patients infected by nine viruses and healthy people were obtained from Gene Expression Omnibus. The analysis of the data was performed by Genome Enhancer pipeline.

Results

We revealed common pathways, cellular processes, and master regulators for studied viral infections. We found that all nine viral infections cause immune activation, exhaustion, cell proliferation disruption, and increased susceptibility to apoptosis. Using network analysis, we identified PBMC receptors, representing proteins at the top of signaling pathways that may be responsible for the observed transcriptional changes and maintain the current functional state of cells.

Discussion

The identified relationships between some of them and virus-induced alteration of immune functions are new and have not been found earlier, e.g., receptors for autocrine motility factor, insulin, prolactin, angiotensin II, and immunoglobulin epsilon. Modulation of the identified receptors can be investigated as one of therapeutic strategies for the treatment of severe viral infections.

Differential host responses within the upper respiratory tract and peripheral blood of children and adults with SARS-CoV-2 infection

Age is among the strongest risk factors for severe outcomes from SARS-CoV-2 infection. We sought to evaluate associations between age and both mucosal and systemic host responses to SARS-CoV-2 infection. We profiled the upper respiratory tract (URT) and peripheral blood transcriptomes of 201 participants (age range of 1 week to 83 years), including 137 non-hospitalized individuals with mild SARS-CoV-2 infection and 64 uninfected individuals. Among uninfected children and adolescents, young age was associated with upregulation of innate and adaptive immune pathways within the URT, suggesting that young children are primed to mount robust mucosal immune responses to exogeneous respiratory pathogens. SARS-CoV-2 infection was associated with broad induction of innate and adaptive immune responses within the URT of children and adolescents. Peripheral blood responses among SARS-CoV-2-infected children and adolescents were dominated by interferon pathways, while upregulation of myeloid activation, inflammatory, and coagulation pathways was observed only in adults. Systemic symptoms among SARS-CoV-2-infected subjects were associated with blunted innate and adaptive immune responses in the URT and upregulation of many of these same pathways within peripheral blood. Finally, within individuals, robust URT immune responses were correlated with decreased peripheral immune activation, suggesting that effective immune responses in the URT may promote local viral control and limit systemic immune activation and symptoms. These findings demonstrate that there are differences in immune responses to SARS-CoV-2 across the lifespan, including between young children and adolescents, and suggest that these varied host responses contribute to observed differences in the clinical presentation of SARS-CoV-2 infection by age.

One Sentence Summary

Age is associated with distinct upper respiratory and peripheral blood transcriptional responses among children and adults with SARS-CoV-2 infection.

Systems-level profiling of early peripheral host-response landscape variations across COVID-19 severity states in an Indian cohort

Host immune response to COVID-19 plays a significant role in regulating disease severity. Although big data analysis has provided significant insights into the host biology of COVID-19 across the world, very few such studies have been performed in the Indian population. This study utilizes a transcriptome-integrated network analysis approach to compare the immune responses between asymptomatic or mild and moderate-severe COVID-19 patients in an Indian cohort. An immune suppression phenotype is observed in the early stages of moderate-severe COVID-19 manifestation. A number of pathways are identified that play crucial roles in the host control of the disease such as the type I interferon response and classical complement pathway which show different activity levels across the severity spectrum. This study also identifies two transcription factors, IRF7 and ESR1, to be important in regulating the severity of COVID-19. Overall this study provides a deep understanding of the peripheral immune landscape in the COVID-19 severity spectrum in the Indian genetic background and opens up future research avenues to compare immune responses across global populations.

Impact of the SARS-CoV-2 nucleocapsid 203K/204R mutations on the inflammatory immune response in COVID-19 severity

BACKGROUND

The excessive inflammatory responses provoked by SARS-CoV-2 infection are critical factors affecting the severity and mortality of COVID-19. Previous work found that two adjacent co-occurring mutations R203K and G204R (KR) on the nucleocapsid (N) protein correlate with increased disease severity in COVID-19 patients. However, links with the host immune response remain unclear.

METHODS

Here, we grouped nasopharyngeal swab samples of COVID-19 patients into two cohorts based on the presence and absence of SARS-CoV-2 nucleocapsid KR mutations. We performed nasopharyngeal transcriptome analysis of age, gender, and ethnicity-matched COVID-19 patients infected with either SARS-CoV-2 with KR mutations in the N protein (KR patients n = 39) or with the wild-type N protein (RG patients n = 39) and compared to healthy controls (n = 34). The impact of KR mutation on immune response was further characterized experimentally by transcriptomic and proteomic profiling of virus-like-particle (VLP) incubated cells.

RESULTS

We observed markedly elevated expression of proinflammatory cytokines, chemokines, and interferon-stimulated (ISGs) genes in the KR patients compared to RG patients. Using nasopharyngeal transcriptome data, we found significantly higher levels of neutrophils and neutrophil-to-lymphocyte (NLR) ratio in KR patients than in the RG patients. Furthermore, transcriptomic and proteomic profiling of VLP incubated cells confirmed a similar hyper-inflammatory response mediated by the KR variant.

CONCLUSIONS

Our data demonstrate an unforeseen connection between nucleocapsid KR mutations and augmented inflammatory immune response in severe COVID-19 patients. These findings provide insights into how mutations in SARS-CoV-2 modulate host immune output and pathogenesis and may contribute to more efficient therapeutics and vaccine development.

Repurposing Niclosamide as a Novel Anti-SARS-CoV-2 Drug by Restricting Entry Protein CD147

The outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to the global coronavirus disease 2019 (COVID-19) pandemic, and the search for effective treatments has been limited. Furthermore, the rapid mutations of SARS-CoV-2 have posed challenges to existing vaccines and neutralizing antibodies, as they struggle to keep up with the increased viral transmissibility and immune evasion. However, there is hope in targeting the CD147-spike protein, which serves as an alternative point for the entry of SARS-CoV-2 into host cells. This protein has emerged as a promising therapeutic target for the development of drugs against COVID-19. Here, we demonstrate that the RNA-binding protein Human-antigen R (HuR) plays a crucial role in the post-transcriptional regulation of CD147 by directly binding to its 3'-untranslated region (UTR). We observed a decrease in CD147 levels across multiple cell lines upon HuR depletion. Furthermore, we identified that niclosamide can reduce CD147 by lowering the cytoplasmic translocation of HuR and reducing CD147 glycosylation. Moreover, our investigation revealed that SARS-CoV-2 infection induces an upregulation of CD147 in ACE2-expressing A549 cells, which can be effectively neutralized by niclosamide in a dose-dependent manner. Overall, our study unveils a novel regulatory mechanism of regulating CD147 through HuR and suggests niclosamide as a promising therapeutic option against COVID-19.

A CRISPR-enhanced metagenomic NGS test to improve pandemic preparedness

The lack of preparedness for detecting and responding to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pathogen (i.e., COVID-19) has caused enormous harm to public health and the economy. Testing strategies deployed on a population scale at day zero, i.e., the time of the first reported case, would be of significant value. Next-generation sequencing (NGS) has such capabilities; however, it has limited detection sensitivity for low-copy-number pathogens. Here, we leverage the CRISPR-Cas9 system to effectively remove abundant sequences not contributing to pathogen detection and show that NGS detection sensitivity of SARS-CoV-2 approaches that of RT-qPCR. The resulting sequence data can also be used for variant strain typing, co-infection detection, and individual human host response assessment, all in a single molecular and analysis workflow. This NGS work flow is pathogen agnostic and, therefore, has the potential to transform how large-scale pandemic response and focused clinical infectious disease testing are pursued in the future.

Dynamics of Gene Expression Profiling and Identification of High-Risk Patients for Severe COVID-19

The clinical manifestations of SARS-CoV-2 infection vary widely, from asymptomatic infection to the development of acute respiratory distress syndrome (ARDS) and death. The host response elicited by SARS-CoV-2 plays a key role in determining the clinical outcome. We hypothesized that determining the dynamic whole blood transcriptomic profile of hospitalized adult COVID-19 patients and characterizing the subgroup that develops severe disease and ARDS would broaden our understanding of the heterogeneity in clinical outcomes. We recruited 60 hospitalized patients with RT-PCR-confirmed SARS-CoV-2 infection, among whom 19 developed ARDS. Peripheral blood was collected using PAXGene RNA tubes within 24 h of admission and on day 7. There were 2572 differently expressed genes in patients with ARDS at baseline and 1149 at day 7. We found a dysregulated inflammatory response in COVID-19 ARDS patients, with an increased expression of genes related to pro-inflammatory molecules and neutrophil and macrophage activation at admission, in addition to an immune regulation loss. This led, in turn, to a higher expression of genes related to reactive oxygen species, protein polyubiquitination, and metalloproteinases in the latter stages. Some of the most significant differences in gene expression found between patients with and without ARDS corresponded to long non-coding RNA involved in epigenetic control.

Repurposing niclosamide as a novel anti-SARS-Cov-2 drug by restricting entry protein CD147

Background The burst of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is causing the global COVID-19 pandemic. But until today only limited numbers of drugs are discovered to treat COVID-19 patients. Even worse, the rapid mutations of SARS-CoV-2 compromise the effectiveness of existing vaccines and neutralizing antibodies due to the increased viral transmissibility and immune escape. CD147-spike protein, one of the entries of SRAR-CoV-2 into host cells, has been reported as a promising therapeutic target for developing drugs against COVID-19. Methods CRISPR-Cas9 induced gene knockout, western blotting, tet-off protein overexpression, ribonucleoprotein IP and RNA-IP were used to confirm the regulation of HuR on mRNA of CD147. Regulation of niclosamide on HuR nucleo-translocation was assessed by immunofluorescence staining of cell lines, IHC staining of tissue of mouse model and western blotting. Finally, the suppression of niclosamide on SARS-CoV-2 infection induced CD147 was evaluated by ACE2-expressing A549 cells and western blotting. Results We first discovered a novel regulation mechanism of CD147 via the RNA-binding protein HuR. We found that HuR regulates CD147 post-transcription by directly bound to its 3'-UTR. The loss of HuR reduced CD147 in multiple cell lines. Niclosamide inhibited CD147 function by blocking HuR cytoplasmic translocation and diminishing CD147 glycosylation. SARS-CoV-2 infection induced CD147 in ACE2-expressing A549 cells, which could be neutralized by niclosamide in a dose-dependent manner. Conclusion Together, our study reveals a novel regulation mechanism of CD147 and niclosamide can be repurposed as an effective COVID-19 drug by targeting the virus entry, CD147-spike protein.

Classification of COVID-19 Patients into Clinically Relevant Subsets by a Novel Machine Learning Pipeline Using Transcriptomic Features

The persistent impact of the COVID-19 pandemic and heterogeneity in disease manifestations point to a need for innovative approaches to identify drivers of immune pathology and predict whether infected patients will present with mild/moderate or severe disease. We have developed a novel iterative machine learning pipeline that utilizes gene enrichment profiles from blood transcriptome data to stratify COVID-19 patients based on disease severity and differentiate severe COVID cases from other patients with acute hypoxic respiratory failure. The pattern of gene module enrichment in COVID-19 patients overall reflected broad cellular expansion and metabolic dysfunction, whereas increased neutrophils, activated B cells, T-cell lymphopenia, and proinflammatory cytokine production were specific to severe COVID patients. Using this pipeline, we also identified small blood gene signatures indicative of COVID-19 diagnosis and severity that could be used as biomarker panels in the clinical setting.

Early peripheral blood MCEMP1 and HLA-DRA expression predicts COVID-19 prognosis

BACKGROUND

Mass vaccination has dramatically reduced the incidence of severe COVID-19, with most cases now presenting as self-limiting upper respiratory tract infections. However, those with co-morbidities, the elderly and immunocompromised, as well as the unvaccinated, remain disproportionately vulnerable to severe COVID-19 and its sequelae. Furthermore, as the effectiveness of vaccination wanes with time, immune escape SARS-CoV-2 variants could emerge to cause severe COVID-19. Reliable prognostic biomarkers for severe disease could be used as early indicator of re-emergence of severe COVID-19 as well as for triaging of patients for antiviral therapy.

METHODS

We performed a systematic review and re-analysis of 7 publicly available datasets, analysing a total of 140 severe and 181 mild COVID-19 patients, to determine the most consistent differentially regulated genes in peripheral blood of severe COVID-19 patients. In addition, we included an independent cohort where blood transcriptomics of COVID-19 patients were prospectively and longitudinally monitored previously, to track the time in which these gene expression changes occur before nadir of respiratory function. Single cell RNA-sequencing of peripheral blood mononuclear cells from publicly available datasets was then used to determine the immune cell subsets involved.

FINDINGS

The most consistent differentially regulated genes in peripheral blood of severe COVID-19 patients were MCEMP1, HLA-DRA and ETS1 across the 7 transcriptomics datasets. Moreover, we found significantly heightened MCEMP1 and reduced HLA-DRA expression as early as four days before the nadir of respiratory function, and the differential expression of MCEMP1 and HLA-DRA occurred predominantly in CD14+ cells. The online platform which we developed is publicly available at https://kuanrongchan-covid19-severity-app-t7l38g.streamlitapp.com/, for users to query gene expression differences between severe and mild COVID-19 patients in these datasets.

INTERPRETATION

Elevated MCEMP1 and reduced HLA-DRA gene expression in CD14+ cells during the early phase of disease are prognostic of severe COVID-19.

FUNDING

K.R.C is funded by the National Medical Research Council (NMRC) of Singapore under the Open Fund Individual Research Grant (MOH-000610). E.E.O. is funded by the NMRC Senior Clinician-Scientist Award (MOH-000135-00). J.G.H.L. is funded by the NMRC under the Clinician-Scientist Award (NMRC/CSAINV/013/2016-01). S.K. is funded by the NMRC under the Transition Award. This study was sponsored in part by a generous gift from The Hour Glass.

Blood RNA alternative splicing events as diagnostic biomarkers for infectious disease

Assays detecting blood transcriptome changes are studied for infectious disease diagnosis. Blood-based RNA alternative splicing (AS) events, which have not been well characterized in pathogen infection, have potential normalization and assay platform stability advantages over gene expression for diagnosis. Here, we present a computational framework for developing AS diagnostic biomarkers. Leveraging a large prospective cohort of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and whole-blood RNA sequencing (RNA-seq) data, we identify a major functional AS program switch upon viral infection. Using an independent cohort, we demonstrate the improved accuracy of AS biomarkers for SARS-CoV-2 diagnosis compared with six reported transcriptome signatures. We then optimize a subset of AS-based biomarkers to develop microfluidic PCR diagnostic assays. This assay achieves nearly perfect test accuracy (61/62 = 98.4%) using a naive principal component classifier, significantly more accurate than a gene expression PCR assay in the same cohort. Therefore, our RNA splicing computational framework enables a promising avenue for host-response diagnosis of infection.

A 2-Gene Host Signature for Improved Accuracy of COVID-19 Diagnosis Agnostic to Viral Variants

The continued emergence of SARS-CoV-2 variants is one of several factors that may cause false-negative viral PCR test results. Such tests are also susceptible to false-positive results due to trace contamination from high viral titer samples. Host immune response markers provide an orthogonal indication of infection that can mitigate these concerns when combined with direct viral detection. Here, we leverage nasopharyngeal swab RNA-seq data from patients with COVID-19, other viral acute respiratory illnesses, and nonviral conditions (n = 318) to develop support vector machine classifiers that rely on a parsimonious 2-gene host signature to diagnose COVID-19. We find that optimal classifiers include an interferon-stimulated gene that is strongly induced in COVID-19 compared with nonviral conditions, such as IFI6, and a second immune-response gene that is more strongly induced in other viral infections, such as GBP5. The IFI6+GBP5 classifier achieves an area under the receiver operating characteristic curve (AUC) greater than 0.9 when evaluated on an independent RNA-seq cohort (n = 553). We further provide proof-of-concept demonstration that the classifier can be implemented in a clinically relevant RT-qPCR assay. Finally, we show that its performance is robust across common SARS-CoV-2 variants and is unaffected by cross-contamination, demonstrating its utility for improved accuracy of COVID-19 diagnostics. IMPORTANCE In this work, we study upper respiratory tract gene expression to develop and validate a 2-gene host-based COVID-19 diagnostic classifier and then demonstrate its implementation in a clinically practical qPCR assay. We find that the host classifier has utility for mitigating false-negative results, for example due to SARS-CoV-2 variants harboring mutations at primer target sites, and for mitigating false-positive viral PCR results due to laboratory cross-contamination. Both types of error carry serious consequences of either unrecognized viral transmission or unnecessary isolation and contact tracing. This work is directly relevant to the ongoing COVID-19 pandemic given the continued emergence of viral variants and the continued challenges of false-positive PCR assays. It also suggests the feasibility of pan-respiratory virus host-based diagnostics that would have value in congregate settings, such as hospitals and nursing homes, where unrecognized respiratory viral transmission is of particular concern.

The role of peroxisome proliferator-activated receptors in the modulation of hyperinflammation induced by SARS-CoV-2 infection: A perspective for COVID-19 therapy

Coronavirus disease 2019 (COVID-19) is a severe respiratory disease caused by infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that affects the lower and upper respiratory tract in humans. SARS-CoV-2 infection is associated with the induction of a cascade of uncontrolled inflammatory responses in the host, ultimately leading to hyperinflammation or cytokine storm. Indeed, cytokine storm is a hallmark of SARS-CoV-2 immunopathogenesis, directly related to the severity of the disease and mortality in COVID-19 patients. Considering the lack of any definitive treatment for COVID-19, targeting key inflammatory factors to regulate the inflammatory response in COVID-19 patients could be a fundamental step to developing effective therapeutic strategies against SARS-CoV-2 infection. Currently, in addition to well-defined metabolic actions, especially lipid metabolism and glucose utilization, there is growing evidence of a central role of the ligand-dependent nuclear receptors and peroxisome proliferator-activated receptors (PPARs) including PPARα, PPARβ/δ, and PPARγ in the control of inflammatory signals in various human inflammatory diseases. This makes them attractive targets for developing therapeutic approaches to control/suppress the hyperinflammatory response in patients with severe COVID-19. In this review, we (1) investigate the anti-inflammatory mechanisms mediated by PPARs and their ligands during SARS-CoV-2 infection, and (2) on the basis of the recent literature, highlight the importance of PPAR subtypes for the development of promising therapeutic approaches against the cytokine storm in severe COVID-19 patients.

Plasma Proteomic Variables Related to COVID-19 Severity: An Untargeted nLC-MS/MS Investigation

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) infection leads to a wide range of clinical manifestations and determines the need for personalized and precision medicine. To better understand the biological determinants of this heterogeneity, we explored the plasma proteome of 43 COVID-19 patients with different outcomes by an untargeted liquid chromatography-mass spectrometry approach. The comparison between asymptomatic or pauci-symptomatic subjects (MILDs), and hospitalised patients in need of oxygen support therapy (SEVEREs) highlighted 29 proteins emerged as differentially expressed: 12 overexpressed in MILDs and 17 in SEVEREs. Moreover, a supervised analysis based on a decision-tree recognised three proteins (Fetuin-A, Ig lambda-2chain-C-region, Vitronectin) that are able to robustly discriminate between the two classes independently from the infection stage. In silico functional annotation of the 29 deregulated proteins pinpointed several functions possibly related to the severity; no pathway was associated exclusively to MILDs, while several only to SEVEREs, and some associated to both MILDs and SEVEREs; SARS-CoV-2 signalling pathway was significantly enriched by proteins up-expressed in SEVEREs (SAA1/2, CRP, HP, LRG1) and in MILDs (GSN, HRG). In conclusion, our analysis could provide key information for 'proteomically' defining possible upstream mechanisms and mediators triggering or limiting the domino effect of the immune-related response and characterizing severe exacerbations.

Machine learning-driven blood transcriptome-based discovery of SARS-CoV-2 specific severity biomarkers

The Coronavirus disease 2019 (COVID-19) pandemic, caused by rapidly evolving variants of severe acute respiratory syndrome coronavirus (SARS-CoV-2), continues to be a global health threat. SARS-CoV-2 infection symptoms often intersect with other nonsevere respiratory infections, making early diagnosis challenging. There is an urgent need for early diagnostic and prognostic biomarkers to predict severity and reduce mortality when a sudden outbreak occurs. This study implemented a novel approach of integrating bioinformatics and machine learning algorithms over publicly available clinical COVID-19 transcriptome data sets. The robust 7-gene biomarker identified through this analysis can not only discriminate SARS-CoV-2 associated acute respiratory illness (ARI) from other types of ARIs but also can discriminate severe COVID-19 patients from nonsevere COVID-19 patients. Validation of the 7-gene biomarker in an independent blood transcriptome data set of longitudinal analysis of COVID-19 patients across various stages of the disease showed that the dysregulation of the identified biomarkers during severe disease is restored during recovery, showing their prognostic potential. The blood biomarkers identified in this study can serve as potential diagnostic candidates and help reduce COVID-19-associated mortality.

Increase of mast cells in COVID-19 pneumonia may contribute to pulmonary fibrosis and thrombosis

AIMS

Lung tissue from COVID-19 patients shares similar histomorphological features with chronic lung allograft disease, also suggesting activation of autoimmune-related pathways in COVID-19. To more clearly understand the underlying spectrum of pathophysiology in COVID-19 pneumonia, we analysed mRNA expression of autoimmune-related genes in post-mortem lung tissue from COVID-19 patients.

METHODS AND RESULTS

Formalin-fixed, paraffin-embedded lung tissue samples of 18 COVID-19 patients and eight influenza patients were used for targeted gene expression profiling using NanoString technology. Multiplex immunofluorescence for tryptase and chymase was applied for validation. Genes related to mast cells were significantly increased in COVID-19. This finding was strengthened by multiplex immunofluorescence also showing a significant increase of tryptase- and chymase-positive cells in COVID-19. Furthermore, receptors for advanced glycation end-products (RAGE) and pro-platelet basic protein (PPBP) were up-regulated in COVID-19 compared to influenza. Genes associated with Type I interferon signalling showed a significant correlation to detected SARS-CoV2 pathway-related genes. The comparison of lung tissue samples from both groups based on the presence of histomorphological features indicative of acute respiratory distress syndrome did not result in finding any specific gene or pathways.

CONCLUSION

Two separate means of measuring show a significant increase of mast cells in SARS-CoV-2-infected lung tissue compared to influenza. Additionally, several genes involved in fibrosis and thrombosis, among which are RAGE and PPBP, are up-regulated in COVID-19. As mast cells are able to induce thrombosis and fibrosis, they may play an important role in the pathogenesis of COVID-19.

Longitudinal home self-collection of capillary blood using homeRNA correlates interferon and innate viral defense pathways with SARS-CoV-2 viral clearance

Blood transcriptional profiling is a powerful tool to evaluate immune responses to infection; however, blood collection via traditional phlebotomy remains a barrier to precise characterization of the immune response in dynamic infections (e.g., respiratory viruses). Here we present an at-home self-collection methodology, homeRNA, to study the host transcriptional response during acute SARS-CoV-2 infections. This method uniquely enables high frequency measurement of the host immune kinetics in non-hospitalized adults during the acute and most dynamic stage of their infection. COVID-19+ and healthy participants self-collected blood every other day for two weeks with daily nasal swabs and symptom surveys to track viral load kinetics and symptom burden, respectively. While healthy uninfected participants showed remarkably stable immune kinetics with no significant dynamic genes, COVID-19+ participants, on the contrary, depicted a robust response with over 418 dynamic genes associated with interferon and innate viral defense pathways. When stratified by vaccination status, we detected distinct response signatures between unvaccinated and breakthrough (vaccinated) infection subgroups; unvaccinated individuals portrayed a response repertoire characterized by higher innate antiviral responses, interferon signaling, and cytotoxic lymphocyte responses while breakthrough infections portrayed lower levels of interferon signaling and enhanced early cell-mediated response. Leveraging cross-platform longitudinal sampling (nasal swabs and blood), we observed that IFI27, a key viral response gene, tracked closely with SARS-CoV-2 viral clearance in individual participants. Taken together, these results demonstrate that at-home sampling can capture key host antiviral responses and facilitate frequent longitudinal sampling to detect transient host immune kinetics during dynamic immune states.

Integrated Data Analysis Uncovers New COVID-19 Related Genes and Potential Drug Re-Purposing Candidates

The COVID-19 pandemic is an acute and rapidly evolving global health crisis. To better understand this disease's molecular basis and design therapeutic strategies, we built upon the recently proposed concept of an integrated cell, iCell, fusing three omics, tissue-specific human molecular interaction networks. We applied this methodology to construct infected and control iCells using gene expression data from patient samples and three cell lines. We found large differences between patient-based and cell line-based iCells (both infected and control), suggesting that cell lines are ill-suited to studying this disease. We compared patient-based infected and control iCells and uncovered genes whose functioning (wiring patterns in iCells) is altered by the disease. We validated in the literature that 18 out of the top 20 of the most rewired genes are indeed COVID-19-related. Since only three of these genes are targets of approved drugs, we applied another data fusion step to predict drugs for re-purposing. We confirmed with molecular docking that the predicted drugs can bind to their predicted targets. Our most interesting prediction is artenimol, an antimalarial agent targeting ZFP62, one of our newly identified COVID-19-related genes. This drug is a derivative of artemisinin drugs that are already under clinical investigation for their potential role in the treatment of COVID-19. Our results demonstrate further applicability of the iCell framework for integrative comparative studies of human diseases.

A two-gene marker for the two-tiered innate immune response in COVID-19 patients

For coronavirus disease 2019 (COVID-19), a pandemic disease characterized by strong immune dysregulation in severe patients, convenient and efficient monitoring of the host immune response is critical. Human hosts respond to viral and bacterial infections in different ways, the former is characterized by the activation of interferon stimulated genes (ISGs) such as IFI27, while the latter is characterized by the activation of anti-bacterial associated genes (ABGs) such as S100A12. This two-tiered innate immune response has not been examined in COVID-19. In this study, the activation patterns of this two-tiered innate immune response represented by IFI27 and S100A12 were explored based on 1421 samples from 17 transcriptome datasets derived from the blood of COVID-19 patients and relevant controls. It was found that IFI27 activation occurred in most of the symptomatic patients and displayed no correlation with disease severity, while S100A12 activation was more restricted to patients under severe and critical conditions with a stepwise activation pattern. In addition, most of the S100A12 activation was accompanied by IFI27 activation. Furthermore, the activation of IFI27 was most pronounced within the first week of symptom onset, but generally waned after 2-3 weeks. On the other hand, the activation of S100A12 displayed no apparent correlation with disease duration and could last for several months in certain patients. These features of the two-tiered innate immune response can further our understanding on the disease mechanism of COVID-19 and may have implications to the clinical triage. Development of a convenient two-gene protocol for the routine serial monitoring of this two-tiered immune response will be a valuable addition to the existing laboratory tests.

Temporal Transcriptome Analysis of SARS-CoV-2-Infected Lung and Spleen in Human ACE2-Transgenic Mice

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmissible and potentially fatal virus. So far, most comprehensive analyses encompassing clinical and transcriptional manifestation have concentrated on the lungs. Here, we confirmed evident signs of viral infection in the lungs and spleen of SARS-CoV-2-infected K18-hACE2 mice, which replicate the phenotype and infection symptoms in hospitalized humans. Seven days post viral detection in organs, infected mice showed decreased vital signs, leading to death. Bronchopneumonia due to infiltration of leukocytes in the lungs and reduction in the spleen lymphocyte region were observed. Transcriptome profiling implicated the meticulous regulation of distress and recovery from cytokine-mediated immunity by distinct immune cell types in a time-dependent manner. In lungs, the chemokine-driven response to viral invasion was highly elevated at 2 days post infection (dpi). In late infection, diseased lungs, post the innate immune process, showed recovery signs. The spleen established an even more immediate line of defense than the lungs, and the cytokine expression profile dropped at 7 dpi. At 5 dpi, spleen samples diverged into two distinct groups with different transcriptome profile and pathophysiology. Inhibition of consecutive host cell viral entry and massive immunoglobulin production and proteolysis inhibition seemed that one group endeavored to survive, while the other group struggled with developmental regeneration against consistent viral intrusion through the replication cycle. Our results may contribute to improved understanding of the longitudinal response to viral infection and development of potential therapeutics for hospitalized patients affected by SARS-CoV-2.

Genome-Wide Pleiotropy Study Identifies Association of PDGFB with Age-Related Macular Degeneration and COVID-19 Infection Outcomes

Age-related macular degeneration (AMD) has been implicated as a risk factor for severe consequences from COVID-19. We evaluated the genetic architecture shared between AMD and COVID-19 (critical illness, hospitalization, and infections) using analyses of genetic correlations and pleiotropy (i.e., cross-phenotype meta-analysis) of AMD (n = 33,976) and COVID-19 (n ≥ 1,388,342) and subsequent analyses including expression quantitative trait locus (eQTL), differential gene expression, and Mendelian randomization (MR). We observed a significant genetic correlation between AMD and COVID-19 infection (r_G = 0.10, p = 0.02) and identified novel genome-wide significant associations near PDGFB (best SNP: rs130651; p = 2.4 × 10^-8) in the pleiotropy analysis of the two diseases. The disease-risk allele of rs130651 was significantly associated with increased gene expression levels of PDGFB in multiple tissues (best eQTL p = 1.8 × 10^-11 in whole blood) and immune cells (best eQTL p = 7.1 × 10^-20 in T-cells). PDGFB expression was observed to be higher in AMD cases than AMD controls {fold change (FC) = 1.02; p = 0.067}, as well as in the peak COVID-19 symptom stage (11-20 days after the symptom onset) compared to early/progressive stage (0-10 days) among COVID-19 patients over age 40 (FC = 2.17; p = 0.03) and age 50 (FC = 2.15; p = 0.04). Our MR analysis found that the liability of AMD risk derived from complement system dysfunction {OR (95% CI); hospitalization = 1.02 (1.01-1.03), infection = 1.02 (1.01-1.03) and increased levels of serum cytokine PDGF-BB {β (95% CI); critical illness = 0.07 (0.02-0.11)} are significantly associated with COVID-19 outcomes. Our study demonstrated that the liability of AMD is associated with an increased risk of COVID-19, and PDGFB may be responsible for the severe COVID-19 outcomes among AMD patients.

Multi-objective optimization identifies a specific and interpretable COVID-19 host response signature

The identification of a COVID-19 host response signature in blood can increase the understanding of SARS-CoV-2 pathogenesis and improve diagnostic tools. Applying a multi-objective optimization framework to both massive public and new multi-omics data, we identified a COVID-19 signature regulated at both transcriptional and epigenetic levels. We validated the signature's robustness in multiple independent COVID-19 cohorts. Using public data from 8,630 subjects and 53 conditions, we demonstrated no cross-reactivity with other viral and bacterial infections, COVID-19 comorbidities, or confounders. In contrast, previously reported COVID-19 signatures were associated with significant cross-reactivity. The signature's interpretation, based on cell-type deconvolution and single-cell data analysis, revealed prominent yet complementary roles for plasmablasts and memory T cells. Although the signal from plasmablasts mediated COVID-19 detection, the signal from memory T cells controlled against cross-reactivity with other viral infections. This framework identified a robust, interpretable COVID-19 signature and is broadly applicable in other disease contexts. A record of this paper's transparent peer review process is included in the supplemental information.

DNA Repair Mechanisms are Activated in Circulating Lymphocytes of Hospitalized Covid-19 Patients

Purpose

Reactive oxygen species (ROS) are an important part of the inflammatory response during infection but can also promote DNA damage. Due to the sustained inflammation in severe Covid-19, we hypothesized that hospitalized Covid-19 patients would be characterized by increased levels of oxidative DNA damage and dysregulation of the DNA repair machinery.

Patients and Methods

Levels of the oxidative DNA lesion 8-oxoG and levels of base excision repair (BER) proteins were measured in peripheral blood mononuclear cells (PBMC) from patients (8-oxoG, n = 22; BER, n = 17) and healthy controls (n = 10) (Cohort 1). Gene expression related to DNA repair was investigated in two independent cohorts of hospitalized Covid-19 patients (Cohort 1; 15 patents and 5 controls, Cohort 2; 15 patients and 6 controls), and by publicly available datasets.

Results

Patients and healthy controls showed comparable amounts of oxidative DNA damage as assessed by 8-oxoG while levels of several BER proteins were increased in Covid-19 patients, indicating enhanced DNA repair in acute Covid-19 disease. Furthermore, gene expression analysis demonstrated regulation of genes involved in BER and double strand break repair (DSBR) in PBMC of Covid-19 patients and expression level of several DSBR genes correlated with the degree of respiratory failure. Finally, by re-analyzing publicly available data, we found that the pathway Hallmark DNA repair was significantly more regulated in circulating immune cells during Covid-19 compared to influenza virus infection, bacterial pneumonia or acute respiratory infection due to seasonal coronavirus.

Conclusion

Although beneficial by protecting against DNA damage, long-term activation of the DNA repair machinery could also contribute to persistent inflammation, potentially through mechanisms such as the induction of cellular senescence. However, further studies that also include measurements of additional markers of DNA damage are required to determine the role and precise molecular mechanisms for DNA repair in SARS-CoV-2 infection.

Endometrial gene expression differences in women with coronavirus disease 2019

OBJECTIVE

To study the potential effect of coronavirus disease (COVID-19) on the endometrial transcriptome of affected, symptomatic women for the detection of altered gene expression.

DESIGN

Pilot study of the endometrial transcriptomes of women manifesting COVID-19 compared with those of women without COVID-19 undergoing hysteroscopic procedures for benign gynecologic disorders using RNA sequencing.

SETTING

Hospital and university laboratories.

PATIENT(S)

Women with (n = 14) and without a COVID-19 (n = 10) diagnosis based on a nasopharyngeal swab analysis using quantitative reverse-transcription polymerase chain reaction. The endometrium of the patients with COVID-19 had previously been tested for severe acute respiratory syndrome coronavirus 2 infection, revealing the absence of the virus in this tissue.

INTERVENTION(S)

Endometrial biopsy sample collection.

MAIN OUTCOMES MEASURE(S)

Endometrial gene expression and functional analysis of symptomatic patients with COVID-19 vs. individuals without the infection.

RESULT(S)

The systemic disease COVID-19 altered endometrial gene expression in 75% of the women, with the patients exhibiting a preponderance of 163 up-regulated (e.g., UTS2, IFI6, IFIH1, and BNIP3) and 72 down-regulated genes (e.g., CPZ, CDH3, and IRF4) (false discovery rate<0.05). A total of 161 dysregulated functions (36 up-regulated and 125 down-regulated) were typically enriched in the endometria of the patients with COVID-19, including up-regulation in pathways involved in the development of immune responses to viruses and cytokine inflammation, reflecting elicitation of a COVID-19 response pathway.

CONCLUSION(S)

Coronavirus disease 2019 affects endometrial gene expression despite the absence of severe acute respiratory syndrome coronavirus 2 RNA in endometrial tissues.

α2,6-Sialylation Is Upregulated in Severe COVID-19, Implicating the Complement Cascade

Better understanding of the molecular mechanisms underlying COVID-19 severity is desperately needed in current times. Although hyper-inflammation drives severe COVID-19, precise mechanisms triggering this cascade and what role glycosylation might play therein are unknown. Here we report the first high-throughput glycomic analysis of COVID-19 plasma samples and autopsy tissues. We find that α2,6-sialylation is upregulated in the plasma of patients with severe COVID-19 and in autopsied lung tissue. This glycan motif is enriched on members of the complement cascade (e.g., C5, C9), which show higher levels of sialylation in severe COVID-19. In the lung tissue, we observe increased complement deposition, associated with elevated α2,6-sialylation levels, corresponding to elevated markers of poor prognosis (IL-6) and fibrotic response. We also observe upregulation of the α2,6-sialylation enzyme ST6GAL1 in patients who succumbed to COVID-19. Our work identifies a heretofore undescribed relationship between sialylation and complement in severe COVID-19, potentially informing future therapeutic development.

Exploration of blood−derived coding and non-coding RNA diagnostic immunological panels for COVID-19 through a co-expressed-based machine learning procedure

Severe acute respiratory syndrome coronavirus 2 (SARS- CoV-2) is the causative virus of the pandemic coronavirus disease 2019 (COVID-19). Evaluating the immunological factors and other implicated processes underlying the progression of COVID-19 is essential for the recognition and then the design of efficacious therapies. Therefore, we analyzed RNAseq data obtained from PBMCs of the COVID-19 patients to explore coding and non-coding RNA diagnostic immunological panels. For this purpose, we integrated multiple RNAseq data and analyzed them overall as well as by considering the state of disease including severe and non-severe conditions. Afterward, we utilized a co-expressed-based machine learning procedure comprising weighted-gene co-expression analysis and differential expression gene as filter phase and recursive feature elimination-support vector machine as wrapper phase. This procedure led to the identification of two modules containing 5 and 84 genes which are mostly involved in cell dysregulation and innate immune suppression, respectively. Moreover, the role of vitamin D in regulating some classifiers was highlighted. Further analysis disclosed the role of discriminant miRNAs including miR-197-3p, miR-150-5p, miR-340-5p, miR-122-5p, miR-1307-3p, miR-34a-5p, miR-98-5p and their target genes comprising GAN, VWC2, TNFRSF6B, and CHST3 in the metabolic pathways. These classifiers differentiate the final fate of infection toward severe or non-severe COVID-19. The identified classifier genes and miRNAs may help in the proper design of therapeutic procedures considering their involvement in the immune and metabolic pathways.

Unpacking COVID-19 Systems Biology in Lung and Whole Blood with Transcriptomics and miRNA Regulators

COVID-19 is a systemic disease affecting tissues and organs, including and beyond the lung. Apart from the current pandemic context, we also have vastly inadequate knowledge of consequences of repeated exposures to SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), the virus causing COVID-19, in multiple organ systems and the whole organism scales when the disease evolves from a pandemic to an endemic state. This calls for a systems biology and systems medicine approach and unpacking the effects of COVID-19 in lung as well as other tissues. We report here original findings from transcriptomics analyses and differentially expressed genes (DEGs) in lung samples from 60 patients and 27 healthy controls, and in whole blood samples from 255 patients and 103 healthy individuals. A total of 11 datasets with RNA-seq transcriptomic data were obtained from the Gene Expression Omnibus and the European Nucleotide Archive. The identified DEGs were used to construct protein interaction and functional networks and to identify related pathways and miRNAs. We found 35 DEGs common between lung and the whole blood, and importantly, 2 novel genes, namely CYP1B1 and TNFAIP6, which have not been previously implicated with COVID-19. We also identified four novel miRNA potential regulators, hsa-mir-192-5p, hsa-mir-221-3p, hsa-mir-4756-3p, and hsa-mir-10a-5p, implicated in lung or other diseases induced by coronaviruses. In summary, these findings offer new molecular leads and insights to unpack COVID-19 systems biology in a whole organism context and might inform future antiviral drug, diagnostics, and vaccine discovery efforts.

Deciphering COVID-19 host transcriptomic complexity and variations for therapeutic discovery against new variants

The molecular manifestations of host cells responding to SARS-CoV-2 and its evolving variants of infection are vastly different across the studied models and conditions, imposing challenges for host-based antiviral drug discovery. Based on the postulation that antiviral drugs tend to reverse the global host gene expression induced by viral infection, we retrospectively evaluated hundreds of signatures derived from 1,700 published host transcriptomic profiles of SARS/MERS/SARS-CoV-2 infection using an iterative data-driven approach. A few of these signatures could be reversed by known anti-SARS-CoV-2 inhibitors, suggesting the potential of extrapolating the biology for new variant research. We discovered IMD-0354 as a promising candidate to reverse the signatures globally with nanomolar IC50 against SARS-CoV-2 and its five variants. IMD-0354 stimulated type I interferon antiviral response, inhibited viral entry, and down-regulated hijacked proteins. This study demonstrates that the conserved coronavirus signatures and the transcriptomic reversal approach that leverages polypharmacological effects could guide new variant therapeutic discovery.

Biomarker screening in preeclampsia: an RNA-sequencing approach based on data from multiple studies

OBJECTIVE

Biomarkers have become important in the prognosis and diagnosis of various diseases. High-throughput methods, such as RNA sequencing facilitate the detection of differentially expressed genes (DEGs), hence potential biomarker candidates. Individual studies suggest long lists of DEGs, hampering the identification of clinically relevant ones. Concerning preeclampsia - a major obstetric burden with high risk for adverse maternal and/or neonatal outcomes - limitations in diagnosis and prediction are still important issues. We, therefore, developed a workflow to facilitate the screening for biomarkers.

METHODS

On the basis of the tool DESeq2, a comprehensive workflow for identifying DEGs was established, analyzing data from several publicly available RNA-sequencing studies. We applied it to four RNA-sequencing datasets (one blood, three placenta) analyzing patients with preeclampsia and normotensive controls. We compared our results with other published approaches and evaluated their performance.

RESULTS

We identified 110 genes that are dysregulated in preeclampsia, observed in at least three of the studies analyzed, six even in all four studies. These included FLT-1, TREM-1, and FN1, which either represent established biomarkers at protein level, or promising candidates based on recent studies. For comparison, using a published meta-analysis approach, 5240 DEGs were obtained.

CONCLUSION

This study presents a data analysis workflow for preeclampsia biomarker screening, capable of identifying promising biomarker candidates, while drastically reducing the numbers of candidates. Moreover, we were also able to confirm its performance for heart failure. This approach can be applied to additional diseases for biomarker identification, and the set of DEGs identified in preeclampsia represents a resource for further studies.

Comparative Upper Respiratory Tract Transcriptomic Profiling Reveals a Potential Role of Early Activation of Interferon Pathway in Severe COVID-19

Infection with SARS-CoV-2 results in Coronavirus disease 2019 (COVID-19) is known to cause mild to acute respiratory infection and sometimes progress towards respiratory failure and death. The mechanisms driving the progression of the disease and accumulation of high viral load in the lungs without initial symptoms remain elusive. In this study, we evaluated the upper respiratory tract host transcriptional response in COVID-19 patients with mild to severe symptoms and compared it with the control COVID-19 negative group using RNA-sequencing (RNA-Seq). Our results reveal an upregulated early type I interferon response in severe COVID-19 patients as compared to mild or negative COVID-19 patients. Moreover, severely symptomatic patients have pronounced induction of interferon stimulated genes (ISGs), particularly the oligoadenylate synthetase (OAS) family of genes. Our results are in concurrence with other studies depicting the early induction of IFN-I response in severe COVID-19 patients, providing novel insights about the ISGs involved.

The Variation of Transcriptomic Perturbations is Associated with the Development and Progression of Various Diseases

Background

Although transcriptomic data have been widely applied to explore various diseases, few studies have investigated the association between transcriptomic perturbations and disease development in a wide variety of diseases.

Methods

Based on a previously developed algorithm for quantifying intratumor heterogeneity at the transcriptomic level, we defined the variation of transcriptomic perturbations (VTP) of a disease relative to the health status. Based on publicly available transcriptome datasets, we compared VTP values between the disease and health status and analyzed correlations between VTP values and disease progression or severity in various diseases, including neurological disorders, infectious diseases, cardiovascular diseases, respiratory diseases, liver diseases, kidney diseases, digestive diseases, and endocrine diseases. We also identified the genes and pathways whose expression perturbations correlated positively with VTP across diverse diseases.

Results

VTP values were upregulated in various diseases relative to their normal controls. VTP values were significantly greater in define than in possible or probable Alzheimer's disease. VTP values were significantly larger in intensive care unit (ICU) COVID-19 patients than in non-ICU patients, and in COVID-19 patients requiring mechanical ventilatory support (MVS) than in those not requiring MVS. VTP correlated positively with viral loads in acquired immune deficiency syndrome (AIDS) patients. Moreover, the AIDS patients treated with abacavir or zidovudine had lower VTP values than those without such therapies. In pulmonary tuberculosis (TB) patients, VTP values followed the pattern: active TB > latent TB > normal controls. VTP values were greater in clinically apparent than in presymptomatic malaria. VTP correlated negatively with the cardiac index of left ventricular ejection fraction (LVEF). In chronic obstructive pulmonary disease (COPD), VTP showed a negative correlation with forced expiratory volume in the first second (FEV1). VTP values increased with H. pylori infection and were upregulated in atrophic gastritis caused by H. pylori infection. The genes and pathways whose expression perturbations correlated positively with VTP scores across diseases were mainly involved in the regulation of immune, metabolic, and cellular activities.

Conclusions

VTP is upregulated in the disease versus health status, and its upregulation is associated with disease progression and severity in various diseases. Thus, VTP has potential clinical implications for disease diagnosis and prognosis.

PD-1/PD-L1 blockade abrogates a dysfunctional innate-adaptive immune axis in critical β-coronavirus disease

Severe COVID-19 is associated with hyperinflammation and weak T cell responses against SARS-CoV-2. However, the links between those processes remain partially characterized. Moreover, whether and how therapeutically manipulating T cells may benefit patients are unknown. Our genetic and pharmacological evidence demonstrates that the ion channel TMEM176B inhibited inflammasome activation triggered by SARS-CoV-2 and SARS-CoV-2-related murine β-coronavirus. Tmem176b^-/- mice infected with murine β-coronavirus developed inflammasome-dependent T cell dysfunction and critical disease, which was controlled by modulating dysfunctional T cells with PD-1 blockers. In critical COVID-19, inflammasome activation correlated with dysfunctional T cells and low monocytic TMEM176B expression, whereas PD-L1 blockade rescued T cell functionality. Here, we mechanistically link T cell dysfunction and inflammation, supporting a cancer immunotherapy to reinforce T cell immunity in critical β-coronavirus disease.

COVID-19 patients exhibit unique transcriptional signatures indicative of disease severity

COVID-19 manifests a spectrum of respiratory symptoms, with the more severe often requiring hospitalization. To identify markers for disease progression, we analyzed longitudinal gene expression data from patients with confirmed SARS-CoV-2 infection admitted to the intensive care unit (ICU) for acute hypoxic respiratory failure (AHRF) as well as other ICU patients with or without AHRF and correlated results of gene set enrichment analysis with clinical features. The results were then compared with a second dataset of COVID-19 patients separated by disease stage and severity. Transcriptomic analysis revealed that enrichment of plasma cells (PCs) was characteristic of all COVID-19 patients whereas enrichment of interferon (IFN) and neutrophil gene signatures was specific to patients requiring hospitalization. Furthermore, gene expression results were used to divide AHRF COVID-19 patients into 2 groups with differences in immune profiles and clinical features indicative of severe disease. Thus, transcriptomic analysis reveals gene signatures unique to COVID-19 patients and provides opportunities for identification of the most at-risk individuals.

COVID-19-specific transcriptomic signature detectable in blood across multiple cohorts

The coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is spreading across the world despite vast global vaccination efforts. Consequently, many studies have looked for potential human host factors and immune mechanisms associated with the disease. However, most studies have focused on comparing COVID-19 patients to healthy controls, while fewer have elucidated the specific host factors distinguishing COVID-19 from other infections. To discover genes specifically related to COVID-19, we reanalyzed transcriptome data from nine independent cohort studies, covering multiple infections, including COVID-19, influenza, seasonal coronaviruses, and bacterial pneumonia. The identified COVID-19-specific signature consisted of 149 genes, involving many signals previously associated with the disease, such as induction of a strong immunoglobulin response and hemostasis, as well as dysregulation of cell cycle-related processes. Additionally, potential new gene candidates related to COVID-19 were discovered. To facilitate exploration of the signature with respect to disease severity, disease progression, and different cell types, we also offer an online tool for easy visualization of the selected genes across multiple datasets at both bulk and single-cell levels.

High Circulating Levels of the Homeostatic Chemokines CCL19 and CCL21 Predict Mortality and Disease Severity in COVID-19

BACKGROUND

Immune dysregulation is a major factor in the development of severe coronavirus disease 2019 (COVID-19). The homeostatic chemokines CCL19 and CCL21 have been implicated as mediators of tissue inflammation, but data on their regulation in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is limited. We thus investigated the levels of these chemokines in COVID-19 patients.

METHODS

Serial blood samples were obtained from patients hospitalized with COVID-19 (n = 414). Circulating CCL19 and CCL21 levels during hospitalization and 3-month follow-up were analyzed. In vitro assays and analysis of RNAseq data from public repositories were performed to further explore possible regulatory mechanisms.

RESULTS

A consistent increase in circulating levels of CCL19 and CCL21 was observed, with high levels correlating with disease severity measures, including respiratory failure, need for intensive care, and 60-day all-cause mortality. High levels of CCL21 at admission were associated with persisting impairment of pulmonary function at the 3-month follow-up.

CONCLUSIONS

Our findings highlight CCL19 and CCL21 as markers of immune dysregulation in COVID-19. This may reflect aberrant regulation triggered by tissue inflammation, as observed in other chronic inflammatory and autoimmune conditions. Determination of the source and regulation of these chemokines and their effects on lung tissue is warranted to further clarify their role in COVID-19.

CLINICAL TRIALS REGISTRATION

NCT04321616 and NCT04381819.

α2,6-Sialylation is Upregulated in Severe COVID-19 Implicating the Complement Cascade

Better understanding of the mechanisms of COVID-19 severity is desperately needed in current times. Although hyper-inflammation drives severe COVID-19, precise mechanisms triggering this cascade and what role glycosylation might play therein is unknown. Here we report the first high-throughput glycomic analysis of COVID-19 plasma samples and autopsy tissues. We find α2,6-sialylation is upregulated in plasma of patients with severe COVID-19 and in the lung. This glycan motif is enriched on members of the complement cascade, which show higher levels of sialylation in severe COVID-19. In the lung tissue, we observe increased complement deposition, associated with elevated α2,6-sialylation levels, corresponding to elevated markers of poor prognosis (IL-6) and fibrotic response. We also observe upregulation of the α2,6-sialylation enzyme ST6GAL1 in patients who succumbed to COVID-19. Our work identifies a heretofore undescribed relationship between sialylation and complement in severe COVID-19, potentially informing future therapeutic development.

Unraveling T Cell Responses for Long Term Protection of SARS-CoV-2 Infection

Due to the COVID-19 pandemic, the global need for vaccines to prevent the disease is imperative. To date, several manufacturers have made efforts to develop vaccines against SARS-CoV-2. In spite of the success of developing many useful vaccines so far, it will be helpful for future vaccine designs, targetting long-term disease protection. For this, we need to know more details of the mechanism of T cell responses to SARS-CoV-2. In this study, we first detected pairwise differentially expressed genes among the healthy, mild, and severe COVID-19 groups of patients based on the expression of CD4+ T cells and CD8+ T cells, respectively. The CD4+ T cells dataset contains 6 mild COVID-19 patients, 8 severe COVID-19 patients, and 6 healthy donors, while the CD8+ T cells dataset has 15 mild COVID-19 patients, 22 severe COVID-19 patients, and 4 healthy donors. Furthermore, we utilized the deep learning algorithm to investigate the potential of differentially expressed genes in distinguishing different disease states. Finally, we built co-expression networks among those genes separately. For CD4+ T cells, we identified 6 modules for the healthy network, 4 modules for the mild network, and 1 module for the severe network; for CD8+ T cells, we detected 6 modules for the healthy network, 4 modules for the mild network, and 3 modules for the severe network. We also obtained hub genes for each module and evaluated the differential connectivity of each gene between pairs of networks constructed on different disease states. Summarizing the results, we find that the following genes TNF, CCL4, XCL1, and IFITM1 can be highly identified with SARS-CoV-2. It is interesting to see that IFITM1 has already been known to inhibit multiple infections with other enveloped viruses, including coronavirus. In addition, our networks show some specific patterns of connectivity among genes and some meaningful clusters related to COVID-19. The results might improve the insight of gene expression mechanisms associated with both CD4+ and CD8+ T cells, expand our understanding of COVID-19 and help develop vaccines with long-term protection.