Single-cell longitudinal analysis of SARS-CoV-2 infection in human airway epithelium identifies target cells, alterations in gene expression, and cell state changes

Metabolic modeling elucidates phenformin and atpenin A5 as broad-spectrum antiviral drugs against RNA viruses

The SARS-CoV-2 pandemic has reemphasized the urgent need for broad-spectrum antiviral therapies. We developed a computational workflow using scRNA-Seq data to assess cellular metabolism during viral infection. With this workflow we predicted the capacity of cells to sustain SARS-CoV-2 virion production in patients and found a tissue-wide induction of metabolic pathways that support viral replication. Expanding our analysis to influenza A and dengue viruses, we identified metabolic targets and inhibitors for potential broad-spectrum antiviral treatment. These targets were highly enriched for known interaction partners of all analyzed viruses. Indeed, phenformin, an NADH:ubiquinone oxidoreductase inhibitor, suppressed SARS-CoV-2 and dengue virus replication. Atpenin A5, blocking succinate dehydrogenase, inhibited SARS-CoV-2, dengue virus, respiratory syncytial virus, and influenza A virus with high selectivity indices. In vivo, phenformin showed antiviral activity against SARS-CoV-2 in a Syrian hamster model. Our work establishes host metabolism as druggable for broad-spectrum antiviral strategies, providing invaluable tools for pandemic preparedness.

Identification and targeting of regulators of SARS-CoV-2-host interactions in the airway epithelium

The impact of SARS-CoV-2 in the lung has been extensively studied, yet the molecular regulators of host-cell programs hijacked by the virus in distinct human airway epithelial cell populations remain poorly understood. Some of the reasons include overreliance on transcriptomic profiling and use of nonprimary cell systems. Here we report a network-based analysis of single-cell transcriptomic profiles able to identify master regulator (MR) proteins controlling SARS-CoV-2-mediated reprogramming in pathophysiologically relevant human ciliated, secretory, and basal cells. This underscored chromatin remodeling, endosomal sorting, ubiquitin pathways, as well as proviral factors identified by CRISPR assays as components of the viral-host response in these cells. Large-scale drug perturbation screens revealed 11 candidate drugs able to invert the entire MR signature activated by SARS-CoV-2. Leveraging MR analysis and perturbational profiles of human primary cells represents an innovative approach to investigate pathogen-host interactions in multiple airway conditions for drug prioritization.

Assessment of repurposed compounds against coronaviruses highlights the antiviral broad-spectrum activity of host-targeting iminosugars and confirms the activity of potent directly acting antivirals

The COVID-19 pandemic highlights the need for novel antiviral drug discovery approaches that could dramatically shorten timelines from compound discovery to clinical development. At the beginning of the pandemic, repurposing approaches were at the forefront of early research efforts to screen for antiviral activity against SARS-CoV-2 in over 2500 compounds. Here, we report cellular screening results of 100 FDA-approved and experimental compounds against SARS-CoV-2 in the human Calu-3 cell line. We observed 13 compounds showing antiviral activity against SARS-CoV-2, including seven FDA-approved compounds (remdesivir, boceprevir, amiloride, nafamostat, cisplatin, silmitasertib, and miglustat), and six compounds in pre-clinical and clinical development (tarloxotinib, lucerastat (NB-DGJ), MON-DNJ, NAP-DNJ, NN-DGJ and NN-DNJ). Further, we observed that our screening hits include several host-targeting antivirals, namely iminosugars, that are largely non-toxic and offer a large therapeutic window. The most-developed iminosugar MON-DNJ (UV-4B), which has been evaluated in a Phase 1 clinical trial, shows antiviral activity against SARS-CoV-2 wild type as well as alpha, beta, gamma, delta, and Omicron variants. Its activity also extended to another betacoronavirus HCoV OC43, but not alphacoronavirus HCoV 229E. Our cellular screening results add to the body of knowledge on antivirals against coronaviruses and confirm the antiviral efficacy of iminosugars in cellular assays using the human lung-cell line Calu-3.

Unlocking the potential of remdesivir: innovative approaches to drug delivery

Given the recurrent waves of COVID-19 and the emergence of new viral infections, optimizing the potential of remdesivir as an antiviral agent is critical. While several reviews have explored the efficacy of remdesivir, few have comprehensively addressed its challenges, such as the necessity for intravenous infusion, suboptimal lung accumulation, and safety concerns related to its formulation. This review critically examines these challenges while proposing innovative solutions and effective combinations with other antiviral agents and repurposed drugs. By highlighting the role of complex generics, we aim to enhance therapeutic efficacy in ways not previously discussed in existing literature. Furthermore, we address the development of novel drug delivery systems which specifically aim to improve remdesivir's pharmacological profile. By analyzing recent findings, we assess both the successes and limitations of current approaches, providing insights into ongoing challenges and strategies for further optimization. This review uniquely focuses on targeted drug delivery systems and innovative formulations, thereby maximizing remdesivir's therapeutic benefits and broadening its application in combating emerging viral threats. In doing so, we fill a critical gap in literature, offering a comprehensive overview that informs future research and clinical strategies.

Alveolar epithelial type 2 cell specific loss of IGFBP2 activates inflammation in COVID-19

The coronavirus disease 2019 (COVID-19) global pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, our understanding of SARS-CoV-2-induced inflammation in alveolar epithelial cells remains very limited. The contributions of intracellular insulin-like growth factor binding protein-2 (IGFBP2) to SARS-CoV-2 pathogenesis are also unclear. In this study, we have uncovered a critical role for IGFBP2, specifically in alveolar epithelial type 2 cells (AEC2), in the immunopathogenesis of COVID-19. Using bulk RNA sequencing, we show that IGFBP2 mRNA expression is significantly downregulated in primary AEC2 cells isolated from fibrotic lung regions from patients with COVID-19-acute respiratory distress syndrome (ARDS) compared to those with idiopathic pulmonary fibrosis (IPF) alone or IPF with a history of COVID-19. Using multicolor immunohistochemistry, we demonstrated that IGFBP2 and its selective ligands IGF1 and IGF2 were significantly reduced in AEC2 cells from patients with COVID-ARDS, IPF alone, or IPF with COVID history than in those from age-matched donor controls. Further, we demonstrated that lentiviral expression of Igfbp2 significantly reduced mRNA expression of proinflammatory cytokines-Tnf-α, Il1β, Il6, Stat3, Stat6 and chemokine receptors-Ccr2 and Ccr5-in mouse lung epithelial cells challenged with SARS-CoV-2 spike protein injury (S2; 500 ng/mL). Finally, we demonstrated higher levels of cytokines-TNF-α; IL-6 and chemokine receptor-CCR5 in AEC2 cells from COVID-ARDS patients compared to the IPF alone and the IPF with COVID history patients. Altogether, these data suggest that anti-inflammatory properties of IGFBP2 in AEC2 cells and its localized delivery may serve as potential therapeutic strategy for patients with COVID-19.

Intrinsic OASL expression licenses interferon induction during influenza A virus infection

Effective control of viral infection requires rapid induction of the innate immune response, especially the type I and type III interferon (IFN) systems. Despite the critical role of IFN induction in host defense, numerous studies have established that most cells fail to produce IFNs in response to viral stimuli. The specific factors that govern cellular heterogeneity in IFN induction potential during infection are not understood. To identify specific host factors that license some cells but not others to mount an IFN response to viral infection, we developed an approach for analyzing temporal scRNA-seq data of influenza A virus (IAV)-infected cells. This approach identified the expression of several interferon stimulated genes (ISGs) within pre-infection cells as correlates of IFN induction potential of those cells, post-infection. Validation experiments confirmed that intrinsic expression of the ISG OASL is essential for robust IFNL induction during IAV infection. Altogether, our findings reveal an important role for IFN-independent, intrinsic expression of ISGs in promoting IFN induction and provide new insights into the mechanisms that regulate cell-to-cell heterogeneity in innate immune activation.

Dual-approach co-expression analysis framework (D-CAF) enables identification of novel circadian co-regulation from multi-omic timeseries data

BACKGROUND

The circadian clock is a central driver of many biological and behavioral processes, regulating the levels of many genes and proteins, termed clock controlled genes and proteins (CCGs/CCPs), to impart biological timing at the molecular level. While transcriptomic and proteomic data has been analyzed to find potential CCGs and CCPs, multi-omic modeling of circadian data, which has the potential to enhance the understanding of circadian control of biological timing, remains relatively rare due to several methodological hurdles. To address this gap, a dual-approach co-expression analysis framework (D-CAF) was created to perform co-expression analysis that is robust to Gaussian noise perturbations on time-series measurements of both transcripts and proteins.

RESULTS

Applying this D-CAF framework to previously gathered transcriptomic and proteomic data from mouse macrophages gathered over circadian time, we identified small, highly significant clusters of oscillating transcripts and proteins in the unweighted similarity matrices and larger, less significant clusters of of oscillating transcripts and proteins using the weighted similarity network. Functional enrichment analysis of these clusters identified novel immunological response pathways that appear to be under circadian control.

CONCLUSIONS

Overall, our findings suggest that D-CAF is a tool that can be used by the circadian community to integrate multi-omic circadian data to improve our understanding of the mechanisms of circadian regulation of molecular processes.

Single-cell transcriptomics of bronchoalveolar lavage during PRRSV infection with different virulence

Porcine reproductive and respiratory syndrome virus (PRRSV) causes significant economic losses in the global swine industry due to its high genetic diversity and different virulence levels, which complicate disease management and vaccine development. This study evaluated longitudinal changes in the immune cell composition of bronchoalveolar lavage fluid and the clinical outcomes across PRRSV strains with varying virulence, using techniques including single-cell transcriptomics. In highly virulent infection, faster viral replication results in an earlier peak lung-damage time point, marked by significant interstitial pneumonia, a significant decrease in macrophages, and an influx of lymphocytes. Viral tracking reveals less than 5% of macrophages are directly infected, and further analysis indicates bystander cell death, likely regulated by exosomal microRNAs as a significant factor. In contrast, the peak intermediate infection shows a delayed lung-damage time point with fewer cell population modifications. Furthermore, anti-inflammatory M2-like macrophages (SPP1-CXCL14high) are identified and their counts increase during the peak lung-damage time point, likely contributing to local defense and lung recovery, which is not observed in high virulent infection. These findings provide a comprehensive description of the immune cellular landscape and differential PRRSV virulence mechanisms, which will help build new hypotheses to understand PRRSV pathogenesis and other respiratory infections.

Airway Basal Stem Cells Inflammatory Alterations in COVID-19 and Mitigation by Mesenchymal Stem Cells

SARS-CoV-2 infection and the resultant COVID-19 pneumonia cause significant damage to the airway and lung epithelium. This damage manifests as mucus hypersecretion, pulmonary inflammation and fibrosis, which often lead to long-term complications collectively referred to as long COVID or post-acute sequelae of COVID-19 (PASC). The airway epithelium, as the first line of defence against respiratory pathogens, depends on airway basal stem cells (BSCs) for regeneration. Alterations in BSCs are associated with impaired epithelial repair and may contribute to the respiratory complications observed in PASC. Given the critical role of BSCs in maintaining epithelial integrity, understanding their alterations in COVID-19 is essential for developing effective therapeutic strategies. This study investigates the intrinsic properties of BSCs derived from COVID-19 patients and evaluates the modulatory effects of mesenchymal stem cells (MSCs). Through a combination of functional assessments and transcriptomic profiling, we identified key phenotypic and molecular deviations in COVID-19 patient-derived BSCs, including goblet cell hyperplasia, inflammation and fibrosis, which may underlie their contribution to PASC. Notably, MSC co-culture significantly mitigated these adverse effects, potentially through modulation of the interferon signalling pathway. This is the first study to isolate BSCs from COVID-19 patients in the Chinese population and establish a COVID-19 BSC-based xenograft model. Our findings reveal critical insights into the role of BSCs in epithelial repair and their inflammatory alterations in COVID-19 pathology, with potential relevance to PASC and virus-induced respiratory sequelae. Additionally, our study highlights MSC-based therapies as a promising strategy to address respiratory sequelae and persistent symptoms.

Identification of virus-rich intermediate cells as crucial players in SARS-CoV-2 infection and differentiation dynamics of human airway epithelium

Understanding the early interactions between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and human airway epithelial cells is essential for unraveling viral replication and spread mechanisms. In this study, we investigated the early dynamics of airway epithelial cells during SARS-CoV-2 infection using well-differentiated human nasal and tracheal epithelial cell cultures by incorporating three publicly available single-cell RNA sequencing datasets. We identified a previously uncharacterized cell population, termed virus-rich intermediate (VRI) cells, representing an intermediate differentiation stage between basal and ciliated cells. These VRI cells exhibited high viral loads at all infection time points, strong interferon and inflammatory responses, increased mRNA expression of microvilli-related genes (PAK1, PAK4, VIL1), and suppression of apoptosis markers (BAX, CASP3) alongside increased anti-apoptotic gene expression (BCL2). Cell-cell interaction analysis revealed that VRI cells send signals to basal cells via receptor-ligand pathways such as EPHA and VEGF, likely promoting basal cell differentiation and proliferation through MAPK signaling. These findings suggest that SARS-CoV-2 utilizes VRI cells as a primary site for replication and spread, leveraging these cells' unique differentiation state to evade host cell death and facilitate viral propagation. This study provides insights into the early cellular responses to SARS-CoV-2 infection and highlights potential therapeutic targets to limit viral spread.

Identification and Targeting of Regulators of SARS-CoV-2-Host Interactions in the Airway Epithelium

Although the impact of SARS-CoV-2 in the lung has been extensively studied, the molecular regulators and targets of the host-cell programs hijacked by the virus in distinct human airway epithelial cell populations remain poorly understood. This is in part ascribed to the use of nonprimary cell systems, overreliance on single-cell gene expression profiling that does not ultimately reflect protein activity, and bias toward the downstream effects rather than their mechanistic determinants. Here we address these issues by network-based analysis of single cell transcriptomic profiles of pathophysiologically relevant human adult basal, ciliated and secretory cells to identify master regulator (MR) protein modules controlling their SARS-CoV-2-mediated reprogramming. This uncovered chromatin remodeling, endosomal sorting, ubiquitin pathways, as well as proviral factors identified by CRISPR analyses as components of the host response collectively or selectively activated in these cells. Large-scale perturbation assays, using a clinically relevant drug library, identified 11 drugs able to invert the entire MR signature activated by SARS-CoV-2 in these cell types. Leveraging MR analysis and perturbational profiles of human primary cells represents a novel mechanism-based approach and resource that can be directly generalized to interrogate signatures of other airway conditions for drug prioritization.

Seasonal influenza a virus lineages exhibit divergent abilities to antagonize interferon induction and signaling

The circulation of seasonal influenza A viruses (IAVs) in humans relies on effective evasion and subversion of the host immune response. While the evolution of seasonal H1N1 and H3N2 viruses to avoid humoral immunity is well characterized, relatively little is known about the evolution of innate immune antagonism phenotypes in these viruses. Numerous studies have established that only a small subset of infected cells is responsible for initiating the type I and type III interferon (IFN) response during IAV infection, emphasizing the importance of single cell studies to accurately characterize the IFN response during infection. We developed a flow cytometry-based method to examine transcriptional changes in IFN and interferon stimulated gene (ISG) expression at the single cell level. We observed that NS segments derived from seasonal H3N2 viruses are more efficient at antagonizing IFN signaling but less effective at suppressing IFN induction, compared to the pdm2009 H1N1 lineage. We compared a collection of NS segments spanning the natural history of the current seasonal IAV lineages and demonstrate long periods of stability in IFN antagonism potential, punctuated by occasional phenotypic shifts. Altogether, our data reveal significant differences in how seasonal and pandemic H1N1 and H3N2 viruses antagonize the human IFN response at the single cell level.

Human placental cells are resistant to SARS-CoV-2 infection and replication

Background

Infection during pregnancy with SARS-CoV-2 can have a serious impact on both maternal and foetal health. Clinical studies have shown that SARS-CoV-2 transmission from the mother to the foetus typically does not occur. However, there is evidence that SARS-CoV-2 can infect the placenta in utero. Here we sought to quantify the permissiveness of placental cells to SARS-CoV-2 infection and to determine if they support viral release.

Methods

By using publicly available single-cell RNA sequencing (scRNAseq) data sets and confocal microscopy we compared ACE2 transcript and protein expression across human first trimester and term placental cells. We also used in vitro infection assays to quantify the infection rates of a range of placenta-derived cells. Finally, we quantified the viral egress from these cells.

Results

ACE2 transcripts are found in a range of placental cell types across gestation, including trophoblast. However, ACE2 protein expression does not significantly change across placental cell types from first trimester to term. We find that 0.5±0.15 % of term trophoblast cells can be infected with SARS-CoV-2 while primary placental fibroblasts and macrophages, and JEG-3, JAR and HUVEC cell lines are resistant to infection. Furthermore, primary trophoblast cells poorly support viral release while JEG-3 cells allow relatively high levels of viral release.

Conclusions

The low level of viral release by primary placental cells provides insight into how the virus is impaired from crossing the placenta to the foetus.

Live imaging of airway epithelium reveals that mucociliary clearance modulates SARS-CoV-2 spread

SARS-CoV-2 initiates infection in the conducting airways, where mucociliary clearance inhibits pathogen penetration. However, it is unclear how mucociliary clearance impacts SARS-CoV-2 spread after infection is established. To investigate viral spread at this site, we perform live imaging of SARS-CoV-2 infected differentiated primary human bronchial epithelium cultures for up to 12 days. Using a fluorescent reporter virus and markers for cilia and mucus, we longitudinally monitor mucus motion, ciliary motion, and infection. Infected cell numbers peak at 4 days post infection, forming characteristic foci that tracked mucus movement. Inhibition of MCC using physical and genetic perturbations limits foci. Later in infection, mucociliary clearance deteriorates. Increased mucus secretion accompanies ciliary motion defects, but mucociliary clearance and vectorial infection spread resume after mucus removal, suggesting that mucus secretion may mediate MCC deterioration. Our work shows that while MCC can facilitate SARS-CoV-2 spread after initial infection, subsequent MCC decreases inhibit spread, revealing a complex interplay between SARS-CoV-2 and MCC.

Method of moments framework for differential expression analysis of single-cell RNA sequencing data

Differential expression analysis of single-cell RNA sequencing (scRNA-seq) data is central for characterizing how experimental factors affect the distribution of gene expression. However, distinguishing between biological and technical sources of cell-cell variability and assessing the statistical significance of quantitative comparisons between cell groups remain challenging. We introduce Memento, a tool for robust and efficient differential analysis of mean expression, variability, and gene correlation from scRNA-seq data, scalable to millions of cells and thousands of samples. We applied Memento to 70,000 tracheal epithelial cells to identify interferon-responsive genes, 160,000 CRISPR-Cas9 perturbed T cells to reconstruct gene-regulatory networks, 1.2 million peripheral blood mononuclear cells (PBMCs) to map cell-type-specific quantitative trait loci (QTLs), and the 50-million-cell CELLxGENE Discover corpus to compare arbitrary cell groups. In all cases, Memento identified more significant and reproducible differences in mean expression compared with existing methods. It also identified differences in variability and gene correlation that suggest distinct transcriptional regulation mechanisms imparted by perturbations.

COVID-19: a multi-organ perspective

In this mini review, we explore the complex network of inflammatory reactions incited by SARS-CoV-2 infection, which extends its reach well beyond the respiratory domain to influence various organ systems. Synthesizing existing literature, it elucidates how the hyperinflammation observed in COVID-19 patients affects multiple organ systems leading to physiological impairments that can persist over long after the resolution of infection. By exploring the systemic manifestations of this inflammatory cascade, from acute respiratory distress syndrome (ARDS) to renal impairment and neurological sequelae, the review highlights the profound interplay between inflammation and organ dysfunction. By synthesizing recent research and clinical observations, this mini review aims to provide an overview of the systemic interactions and complications associated with COVID-19, underscoring the need for an integrated approach to treatment and management. Understanding these systemic effects is crucial for improving patient outcomes and preparing for future public health challenges.

Characterization and implementation of the MarathonRT template-switching reaction to expand the capabilities of RNA-seq

End-to-end RNA-sequencing methods that capture 5'-sequence content without cumbersome library manipulations are of great interest, particularly for analysis of long RNAs. While template-switching methods have been developed for RNA sequencing by distributive short-read RTs, such as the MMLV RTs used in SMART-Seq methods, they have not been adapted to leverage the power of ultraprocessive RTs, such as those derived from group II introns. To facilitate this transition, we dissected the individual processes that guide the enzymatic specificity and efficiency of the multistep template-switching reaction carried out by RTs, in this case, by MarathonRT. Remarkably, this is the first study of its kind, for any RT. First, we characterized the nucleotide specificity of nontemplated addition (NTA) reaction that occurs when the RT extends past the RNA 5'-terminus. We then evaluated the binding specificity of specialized template-switching oligonucleotides, optimizing their sequences and chemical properties to guide efficient template-switching reaction. Having dissected and optimized these individual steps, we then unified them into a procedure for performing RNA sequencing with MarathonRT enzymes, using a well-characterized RNA reference set. The resulting reads span a six-log range in transcript concentration and accurately represent the input RNA identities in both length and composition. We also performed RNA-seq from total human RNA and poly(A)-enriched RNA, with short- and long-read sequencing demonstrating that MarathonRT enhances the discovery of unseen RNA molecules by conventional RT. Altogether, we have generated a new pipeline for rapid, accurate sequencing of complex RNA libraries containing mixtures of long RNA transcripts.

Dual-Approach Co-expression Analysis Framework (D-CAF) Enables Identification of Novel Circadian Regulation From Multi-Omic Timeseries Data

The circadian clock is a central driver of many biological and behavioral processes, regulating the levels of many genes and proteins, termed clock controlled genes and proteins (CCGs/CCPs), to impart biological timing at the molecular level. While transcriptomic and proteomic data has been analyzed to find potential CCGs and CCPs, multi-omic modeling of circadian data, which has the potential to enhance the understanding of circadian control of biological timing, remains relatively rare due to several methodological hurdles. To address this gap, a Dual-approach Co-expression Analysis Framework (D-CAF) was created to perform perturbation-robust co-expression analysis on time-series measurements of both transcripts and proteins. Applying this D-CAF framework to previously gathered transcriptomic and proteomic data from mouse macrophages gathered over circadian time, we identified small, highly significant clusters of oscillating transcripts and proteins in the unweighted similarity matrices and larger, less significant clusters of of oscillating transcripts and proteins using the weighted similarity network. Functional enrichment analysis of these clusters identified novel immunological response pathways that appear to be under circadian control. Overall, our findings suggest that D-CAF is a tool that can be used by the circadian community to integrate multi-omic circadian data to improve our understanding of the mechanisms of circadian regulation of molecular processes.

Harnessing high-throughput OMICS in emerging zoonotic virus preparedness and response activities

Emerging and re-emerging viruses pose unpredictable and significant challenges to global health. Emerging zoonotic infectious diseases, which are transmitted between humans and non-human animals, have been estimated to be responsible for nearly two-thirds of emerging infectious disease events and emergence events attributed to these pathogens have been increasing in frequency with the potential for high global health and economic burdens. In this review we will focus on the application of highthroughput OMICS approaches to emerging zoonotic virus investigtations. We highlight the key contributions of transcriptome and proteome investigations to emerging zoonotic virus preparedness and response activities with a focus on SARS-CoV-2, avian influenza virus subtype H5N1, and Orthoebolavirus investigations.

A basally active cGAS-STING pathway limits SARS-CoV-2 replication in a subset of ACE2 positive airway cell models

Host factors that define the cellular tropism of SARS-CoV-2 beyond the cognate ACE2 receptor are poorly defined. Here we report that SARS-CoV-2 replication is restricted at a post-entry step in a number of ACE2-positive airway-derived cell lines due to tonic activation of the cGAS-STING pathway mediated by mitochondrial DNA leakage and naturally occurring cGAS and STING variants. Genetic and pharmacological inhibition of the cGAS-STING and type I/III IFN pathways as well as ACE2 overexpression overcome these blocks. SARS-CoV-2 replication in STING knockout cell lines and primary airway cultures induces ISG expression but only in uninfected bystander cells, demonstrating efficient antagonism of the type I/III IFN-pathway in productively infected cells. Pharmacological inhibition of STING in primary airway cells enhances SARS-CoV-2 replication and reduces virus-induced innate immune activation. Together, our study highlights that tonic activation of the cGAS-STING and IFN pathways can impact SARS-CoV-2 cellular tropism in a manner dependent on ACE2 expression levels.

Interactions of SARS-CoV-2 with Human Target Cells-A Metabolic View

Viruses are obligate intracellular parasites, and they exploit the cellular pathways and resources of their respective host cells to survive and successfully multiply. The strategies of viruses concerning how to take advantage of the metabolic capabilities of host cells for their own replication can vary considerably. The most common metabolic alterations triggered by viruses affect the central carbon metabolism of infected host cells, in particular glycolysis, the pentose phosphate pathway, and the tricarboxylic acid cycle. The upregulation of these processes is aimed to increase the supply of nucleotides, amino acids, and lipids since these metabolic products are crucial for efficient viral proliferation. In detail, however, this manipulation may affect multiple sites and regulatory mechanisms of host-cell metabolism, depending not only on the specific viruses but also on the type of infected host cells. In this review, we report metabolic situations and reprogramming in different human host cells, tissues, and organs that are favorable for acute and persistent SARS-CoV-2 infection. This knowledge may be fundamental for the development of host-directed therapies.

Rhinovirus infection of airway epithelial cells uncovers the non-ciliated subset as a likely driver of genetic risk to childhood-onset asthma

Asthma is a complex disease caused by genetic and environmental factors. Studies show that wheezing during rhinovirus infection correlates with childhood asthma development. Over 150 non-coding risk variants for asthma have been identified, many affecting gene regulation in T cells, but the effects of most risk variants remain unknown. We hypothesized that airway epithelial cells could also mediate genetic susceptibility to asthma given they are the first line of defense against respiratory viruses and allergens. We integrated genetic data with transcriptomics of airway epithelial cells subject to different stimuli. We demonstrate that rhinovirus infection significantly upregulates childhood-onset asthma-associated genes, particularly in non-ciliated cells. This enrichment is also observed with influenza infection but not with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) or cytokine activation. Overall, our results suggest that rhinovirus infection is an environmental factor that interacts with genetic risk factors through non-ciliated airway epithelial cells to drive childhood-onset asthma.

Cell Culture Differentiation and Proliferation Conditions Influence the In Vitro Regeneration of the Human Airway Epithelium

The human airway mucociliary epithelium can be recapitulated in vitro using primary cells cultured in an air-liquid interface (ALI), a reliable surrogate to perform pathophysiological studies. As tremendous variations exist among media used for ALI-cultured human airway epithelial cells, the aim of our study was to evaluate the impact of several media (BEGM, PneumaCult, Half & Half, and Clancy) on cell type distribution using single-cell RNA sequencing and imaging. Our work revealed the impact of these media on cell composition, gene expression profile, cell signaling, and epithelial morphology. We found higher proportions of multiciliated cells in PneumaCult-ALI and Half & Half, stronger EGF signaling from basal cells in BEGM-ALI, differential expression of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) entry factor ACE2, and distinct secretome transcripts depending on the media used. We also established that proliferation in PneumaCult-Ex Plus favored secretory cell fate, showing the key influence of proliferation media on late differentiation epithelial characteristics. Altogether, our data offer a comprehensive repertoire for evaluating the effects of culture conditions on airway epithelial differentiation and will aid in choosing the most relevant medium according to the processes to be investigated, such as cilia, mucus biology, or viral infection. We detail useful parameters that should be explored to document airway epithelial cell fate and morphology.

SARS-CoV-2-related bat viruses evade human intrinsic immunity but lack efficient transmission capacity

Circulating bat coronaviruses represent a pandemic threat. However, our understanding of bat coronavirus pathogenesis and transmission potential is limited by the lack of phenotypically characterized strains. We created molecular clones for the two closest known relatives of SARS-CoV-2, BANAL-52 and BANAL-236. We demonstrated that BANAL-CoVs and SARS-CoV-2 have similar replication kinetics in human bronchial epithelial cells. However, BANAL-CoVs have impaired replication in human nasal epithelial cells and in the upper airway of mice. We also observed reduced pathogenesis in mice and diminished transmission in hamsters. Further, we observed that diverse bat coronaviruses evade interferon and downregulate major histocompatibility complex class I. Collectively, our study demonstrates that despite high genetic similarity across bat coronaviruses, prediction of pandemic potential of a virus necessitates functional characterization. Finally, the restriction of bat coronavirus replication in the upper airway highlights that transmission potential and innate immune restriction can be uncoupled in this high-risk family of emerging viruses.

A therapy for suppressing canonical and noncanonical SARS-CoV-2 viral entry and an intrinsic intrapulmonary inflammatory response

The prevalence of "long COVID" is just one of the conundrums highlighting how little we know about the lung's response to viral infection, particularly to syndromecoronavirus-2 (SARS-CoV-2), for which the lung is the point of entry. We used an in vitro human lung system to enable a prospective, unbiased, sequential single-cell level analysis of pulmonary cell responses to infection by multiple SARS-CoV-2 strains. Starting with human induced pluripotent stem cells and emulating lung organogenesis, we generated and infected three-dimensional, multi-cell-type-containing lung organoids (LOs) and gained several unexpected insights. First, SARS-CoV-2 tropism is much broader than previously believed: Many lung cell types are infectable, if not through a canonical receptor-mediated route (e.g., via Angiotensin-converting encyme 2(ACE2)) then via a noncanonical "backdoor" route (via macropinocytosis, a form of endocytosis). Food and Drug Administration (FDA)-approved endocytosis blockers can abrogate such entry, suggesting adjunctive therapies. Regardless of the route of entry, the virus triggers a lung-autonomous, pulmonary epithelial cell-intrinsic, innate immune response involving interferons and cytokine/chemokine production in the absence of hematopoietic derivatives. The virus can spread rapidly throughout human LOs resulting in mitochondrial apoptosis mediated by the prosurvival protein Bcl-xL. This host cytopathic response to the virus may help explain persistent inflammatory signatures in a dysfunctional pulmonary environment of long COVID. The host response to the virus is, in significant part, dependent on pulmonary Surfactant Protein-B, which plays an unanticipated role in signal transduction, viral resistance, dampening of systemic inflammatory cytokine production, and minimizing apoptosis. Exogenous surfactant, in fact, can be broadly therapeutic.

Pseudomonas aeruginosa breaches respiratory epithelia through goblet cell invasion in a microtissue model

Pseudomonas aeruginosa, a leading cause of severe hospital-acquired pneumonia, causes infections with up to 50% mortality rates in mechanically ventilated patients. Despite some knowledge of virulence factors involved, it remains unclear how P. aeruginosa disseminates on mucosal surfaces and invades the tissue barrier. Using infection of human respiratory epithelium organoids, here we observed that P. aeruginosa colonization of apical surfaces is promoted by cyclic di-GMP-dependent asymmetric division. Infection with mutant strains revealed that Type 6 Secretion System activities promote preferential invasion of goblet cells. Type 3 Secretion System activity by intracellular bacteria induced goblet cell death and expulsion, leading to epithelial rupture which increased bacterial translocation and dissemination to the basolateral epithelium. These findings show that under physiological conditions, P. aeruginosa uses coordinated activity of a specific combination of virulence factors and behaviours to invade goblet cells and breach the epithelial barrier from within, revealing mechanistic insight into lung infection dynamics.

Self-inhibiting percolation and viral spreading in epithelial tissue

SARS-CoV-2 induces delayed type-I/III interferon production, allowing it to escape the early innate immune response. The delay has been attributed to a deficiency in the ability of cells to sense viral replication upon infection, which in turn hampers activation of the antiviral state in bystander cells. Here, we introduce a cellular automaton model to investigate the spatiotemporal spreading of viral infection as a function of virus and host-dependent parameters. The model suggests that the considerable person-to-person heterogeneity in SARS-CoV-2 infections is a consequence of high sensitivity to slight variations in biological parameters near a critical threshold. It further suggests that within-host viral proliferation can be curtailed by the presence of remarkably few cells that are primed for IFN production. Thus, the observed heterogeneity in defense readiness of cells reflects a remarkably cost-efficient strategy for protection.

Uncovering strain- and age-dependent innate immune responses to SARS-CoV-2 infection in air-liquid-interface cultured nasal epithelia

Continuous assessment of the impact of SARS-CoV-2 on the host at the cell-type level is crucial for understanding key mechanisms involved in host defense responses to viral infection. We investigated host response to ancestral-strain and Alpha-variant SARS-CoV-2 infections within air-liquid-interface human nasal epithelial cells from younger adults (26-32 Y) and older children (12-14 Y) using single-cell RNA-sequencing. Ciliated and secretory-ciliated cells formed the majority of highly infected cell-types, with the latter derived from ciliated lineages. Strong innate immune responses were observed across lowly infected and uninfected bystander cells and heightened in Alpha-infection. Alpha highly infected cells showed increased expression of protein-refolding genes compared with ancestral-strain-infected cells in children. Furthermore, oxidative phosphorylation-related genes were down-regulated in bystander cells versus infected and mock-control cells, underscoring the importance of these biological functions for viral replication. Overall, this study highlights the complexity of cell-type-, age- and viral strain-dependent host epithelial responses to SARS-CoV-2.

CASi: A framework for cross-timepoint analysis of single-cell RNA sequencing data

Single-cell RNA sequencing (scRNA-seq) technology has been widely used to study the differences in gene expression at the single cell level, providing insights into the research of cell development, differentiation, and functional heterogeneity. Various pipelines and workflows of scRNA-seq analysis have been developed but few considered multi-timepoint data specifically. In this study, we develop CASi, a comprehensive framework for analyzing multiple timepoints' scRNA-seq data, which provides users with: (1) cross-timepoint cell annotation, (2) detection of potentially novel cell types emerged over time, (3) visualization of cell population evolution, and (4) identification of temporal differentially expressed genes (tDEGs). Through comprehensive simulation studies and applications to a real multi-timepoint single cell dataset, we demonstrate the robust and favorable performance of the proposal versus existing methods serving similar purposes.

SARS-CoV-2 variant of concern fitness and adaptation in primary human airway epithelia

The severe acute respiratory syndrome coronavirus 2 pandemic is characterized by the emergence of novel variants of concern (VOCs) that replace ancestral strains. Here, we dissect the complex selective pressures by evaluating variant fitness and adaptation in human respiratory tissues. We evaluate viral properties and host responses to reconstruct forces behind D614G through Omicron (BA.1) emergence. We observe differential replication in airway epithelia, differences in cellular tropism, and virus-induced cytotoxicity. D614G accumulates the most mutations after infection, supporting zoonosis and adaptation to the human airway. We perform head-to-head competitions and observe the highest fitness for Gamma and Delta. Under these conditions, RNA recombination favors variants encoding the B.1.617.1 lineage 3' end. Based on viral growth kinetics, Alpha, Gamma, and Delta exhibit increased fitness compared to D614G. In contrast, the global success of Omicron likely derives from increased transmission and antigenic variation. Our data provide molecular evidence to support epidemiological observations of VOC emergence.

SARS-CoV-2-infected human airway epithelial cell cultures uniquely lack interferon and immediate early gene responses caused by other coronaviruses

Objectives

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a member of a class of highly pathogenic coronaviruses. The large family of coronaviruses, however, also includes members that cause only mild symptoms, like human coronavirus-229E (HCoV-229E) or OC43 (HCoV-OC43). Unravelling how molecular (and cellular) pathophysiology differs between highly and low pathogenic coronaviruses is important for the development of therapeutic strategies.

Methods

Here, we analysed the transcriptome of primary human bronchial epithelial cells (PBEC), differentiated at the air-liquid interface (ALI) after infection with SARS-CoV-2, SARS-CoV, Middle East Respiratory Syndrome (MERS)-CoV and HCoV-229E using bulk RNA sequencing.

Results

ALI-PBEC were efficiently infected by all viruses, and SARS-CoV, MERS-CoV and HCoV-229E infection resulted in a largely similar transcriptional response. The response to SARS-CoV-2 infection differed markedly as it uniquely lacked the increase in expression of immediate early genes, including FOS, FOSB and NR4A1 that was observed with all other coronaviruses. This finding was further confirmed in publicly available experimental and clinical datasets. Interfering with NR4A1 signalling in Calu-3 lung epithelial cells resulted in a 100-fold reduction in extracellular RNA copies of SARS-CoV-2 and MERS-CoV, suggesting an involvement in virus replication. Furthermore, a lack in induction of interferon-related gene expression characterised the main difference between the highly pathogenic coronaviruses and low pathogenic viruses HCoV-229E and HCoV-OC43.

Conclusion

Our results demonstrate a previously unknown suppression of a host response gene set by SARS-CoV-2 and confirm a difference in interferon-related gene expression between highly pathogenic and low pathogenic coronaviruses.

Overexpression of human ACE2 protein in mouse fibroblasts stably transfected with the intact ACE2 gene

Infection by SARS-CoV-2 is dependent on binding of the viral spike protein to angiotensin converting enzyme 2 (ACE2), a membrane glycoprotein expressed on epithelial cells in the human upper respiratory tract. Recombinant ACE2 protein has potential application for anti-viral therapy. Here we co-transfected mouse fibroblasts (A9 cells) with a cloned fragment of human genomic DNA containing the intact ACE2 gene and an unlinked neomycin phosphotransferase gene, and then selected stable neomycin-resistant transfectants. Transfectant clones expressed ACE2 protein at levels that were generally proportional to the number of ACE2 gene copies integrated in the cell genome, ranging up to approximately 50 times the level of ACE2 present of Vero-E6 cells. Cells overexpressing ACE2 were hypersensitive to infection by spike-pseudotyped vesicular stomatitis virus (VSV-S), and adsorption of VSV-S to these cells occurred at an accelerated rate compared to Vero-E6 cells. The transfectant cell clones described here therefore have favorable attributes as feedstocks for large-scale production of recombinant human ACE2 protein.

Differential beta-coronavirus infection dynamics in human bronchial epithelial organoids

The lower respiratory system serves as the target and barrier for beta-coronavirus (beta-CoV) infections. In this study, we explored beta-CoV infection dynamics in human bronchial epithelial (HBE) organoids, focusing on HCoV-OC43, SARS-CoV, MERS-CoV, and SARS-CoV-2. Utilizing advanced organoid culture techniques, we observed robust replication for all beta-CoVs, particularly noting that SARS-CoV-2 reached peak viral RNA levels at 72 h postinfection. Through comprehensive transcriptomic analysis, we identified significant shifts in cell population dynamics, marked by an increase in goblet cells and a concurrent decrease in ciliated cells. Furthermore, our cell tropism analysis unveiled distinct preferences in viral targeting: HCoV-OC43 predominantly infected club cells, while SARS-CoV had a dual tropism for goblet and ciliated cells. In contrast, SARS-CoV-2 primarily infected ciliated cells, and MERS-CoV showed a marked affinity for goblet cells. Host factor analysis revealed the upregulation of genes encoding viral receptors and proteases. Notably, HCoV-OC43 induced the unfolded protein response pathway, which may facilitate viral replication. Our study also reveals a complex interplay between inflammatory pathways and the suppression of interferon responses during beta-CoV infections. These findings provide insights into host-virus interactions and antiviral defense mechanisms, contributing to our understanding of beta-CoV infections in the respiratory tract.

Human Bronchial Epithelial Cell Transcriptome Changes in Response to Serum from Patients with Different Status of Inflammation

PURPOSE

To investigate the transcriptome of human bronchial epithelial cells (HBEC) in response to serum from patients with different degrees of inflammation.

METHODS

Serum from 19 COVID-19 patients obtained from the Hannover Unified Biobank was used. At the time of sampling, 5 patients had a WHO Clinical Progression Scale (WHO-CPS) score of 9 (severe illness). The remaining 14 patients had a WHO-CPS of below 9 (range 1-7), and lower illness. Multiplex immunoassay was used to assess serum inflammatory markers. The culture medium of HBEC was supplemented with 2% of the patient's serum, and the cells were cultured at 37 °C, 5% CO2 for 18 h. Subsequently, cellular RNA was used for RNA-Seq.

RESULTS

Patients with scores below 9 had significantly lower albumin and serum levels of E-selectin, IL-8, and MCP-1 than patients with scores of 9. Principal component analysis based on 500 "core genes" of RNA-seq segregated cells into two subsets: exposed to serum from 4 (I) and 15 (II) patients. Cells from a subset (I) treated with serum from 4 patients with a score of 9 showed 5566 differentially expressed genes of which 2793 were up- and 2773 downregulated in comparison with cells of subset II treated with serum from 14 patients with scores between 1 and 7 and one with score = 9. In subset I cells, a higher expression of TLR4 and CXCL8 but a lower CDH1, ACE2, and HMOX1, and greater effects on genes involved in metabolic regulation, cytoskeletal organization, and kinase activity pathways were observed.

CONCLUSION

This simple model could be useful to characterize patient serum and epithelial cell properties.

scPathoQuant: a tool for efficient alignment and quantification of pathogen sequence reads from 10× single cell sequencing datasets

MOTIVATION

Currently there is a lack of efficient computational pipelines/tools for conducting simultaneous genome mapping of pathogen-derived and host reads from single cell RNA sequencing (scRNAseq) output from pathogen-infected cells. Contemporary options include processes involving multiple steps and/or running multiple computational tools, increasing user operations time.

RESULTS

To address the need for new tools to directly map and quantify pathogen and host sequence reads from within an infected cell from scRNAseq datasets in a single operation, we have built a python package, called scPathoQuant. scPathoQuant extracts sequences that were not aligned to the primary host genome, maps them to a pathogen genome of interest (here as demonstrated for viral pathogens), quantifies total reads mapping to the entire pathogen, quantifies reads mapping to individual pathogen genes, and finally integrates pathogen sequence counts into matrix files that are used by standard single cell pipelines for downstream analyses with only one command. We demonstrate that scPathoQuant provides a scRNAseq viral and host genome-wide sequence read abundance analysis that can differentiate and define multiple viruses in a single sample scRNAseq output.

AVAILABILITY AND IMPLEMENTATION

The SPQ package is available software accessible at https://github.com/galelab/scPathoQuant (DOI 10.5281/zenodo.10463670) with test codes and datasets available https://github.com/galelab/Whitmore_scPathoQuant_testSets (DOI 10.5281/zenodo.10463677) to serve as a resource for the community.

Synthesis and biological evaluation of novel peptidomimetic inhibitors of the coronavirus 3C-like protease

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and related variants, are responsible for the devastating coronavirus disease 2019 (COVID-19) pandemic. The SARS-CoV-2 main protease (Mpro) plays a central role in the replication of the virus and represents an attractive drug target. Herein, we report the discovery of novel SARS-CoV-2 Mpro covalent inhibitors, including highly effective compound NIP-22c which displays high potency against several key variants and clinically relevant nirmatrelvir Mpro E166V mutants.

Comparative analysis of cell-cell communication at single-cell resolution

Inference of cell-cell communication from single-cell RNA sequencing data is a powerful technique to uncover intercellular communication pathways, yet existing methods perform this analysis at the level of the cell type or cluster, discarding single-cell-level information. Here we present Scriabin, a flexible and scalable framework for comparative analysis of cell-cell communication at single-cell resolution that is performed without cell aggregation or downsampling. We use multiple published atlas-scale datasets, genetic perturbation screens and direct experimental validation to show that Scriabin accurately recovers expected cell-cell communication edges and identifies communication networks that can be obscured by agglomerative methods. Additionally, we use spatial transcriptomic data to show that Scriabin can uncover spatial features of interaction from dissociated data alone. Finally, we demonstrate applications to longitudinal datasets to follow communication pathways operating between timepoints. Our approach represents a broadly applicable strategy to reveal the full structure of niche-phenotype relationships in health and disease.

Rhinovirus infection of airway epithelial cells uncovers the non-ciliated subset as a likely driver of genetic susceptibility to childhood-onset asthma

Asthma is a complex disease caused by genetic and environmental factors. Epidemiological studies have shown that in children, wheezing during rhinovirus infection (a cause of the common cold) is associated with asthma development during childhood. This has led scientists to hypothesize there could be a causal relationship between rhinovirus infection and asthma or that RV-induced wheezing identifies individuals at increased risk for asthma development. However, not all children who wheeze when they have a cold develop asthma. Genome-wide association studies (GWAS) have identified hundreds of genetic variants contributing to asthma susceptibility, with the vast majority of likely causal variants being non-coding. Integrative analyses with transcriptomic and epigenomic datasets have indicated that T cells drive asthma risk, which has been supported by mouse studies. However, the datasets ascertained in these integrative analyses lack airway epithelial cells. Furthermore, large-scale transcriptomic T cell studies have not identified the regulatory effects of most non-coding risk variants in asthma GWAS, indicating there could be additional cell types harboring these "missing regulatory effects". Given that airway epithelial cells are the first line of defense against rhinovirus, we hypothesized they could be mediators of genetic susceptibility to asthma. Here we integrate GWAS data with transcriptomic datasets of airway epithelial cells subject to stimuli that could induce activation states relevant to asthma. We demonstrate that epithelial cultures infected with rhinovirus significantly upregulate childhood-onset asthma-associated genes. We show that this upregulation occurs specifically in non-ciliated epithelial cells. This enrichment for genes in asthma risk loci, or 'asthma heritability enrichment' is also significant for epithelial genes upregulated with influenza infection, but not with SARS-CoV-2 infection or cytokine activation. Additionally, cells from patients with asthma showed a stronger heritability enrichment compared to cells from healthy individuals. Overall, our results suggest that rhinovirus infection is an environmental factor that interacts with genetic risk factors through non-ciliated airway epithelial cells to drive childhood-onset asthma.

Intestinal immunological events of acute and resolved SARS-CoV-2 infection in non-human primates

SARS-CoV-2 infection has been associated with intestinal mucosal barrier damage, leading to microbial and endotoxin translocation, heightened inflammatory responses, and aggravated disease outcomes. This study aimed to investigate the immunological mechanisms associated with impaired intestinal barrier function. We conducted a comprehensive analysis of gut damage and inflammation markers and phenotypic characterization of myeloid and lymphoid populations in the ileum and colon of SARS-CoV-2-exposed macaques during both the acute and resolved infection phases. Our findings revealed a significant accumulation of terminally differentiated and activated CD4+ and CD8+ T cells, along with memory B cells, within the gastrointestinal tract up to 43 days after exposure to SARS-CoV-2. This robust infection-induced immune response was accompanied by a notable depletion of plasmacytoid dendritic cells, myeloid dendritic cells, and macrophages, particularly affecting the colon during the resolved infection phase. Additionally, we identified a population of CX3CR1Low inflammatory macrophages associated with intestinal damage during active viral replication. Elevated levels of immune activation and gut damage markers, and perturbation of macrophage homeostasis, persisted even after the resolution of the infection, suggesting potential long-term clinical sequelae. These findings enhance our understanding of gastrointestinal immune pathology following SARS-CoV-2 infection and provide valuable information for developing and testing medical countermeasures.

Single-Cell RNA-Sequencing of RVFV Infection

On the RNA level, viral infections are characterized by perturbations in the host cell transcriptome as well as the development of viral genetic information. Investigating the abundance and dynamic of RNA molecules can provide ample information to understand many aspects of the infection, from viral replication to pathogenesis. A key aspect therein is the resolution of the data, as infections are generally highly heterogeneous. Even in simple model systems such as cell lines, viral infections happen in a very asynchronous way. Quantifying RNAs at single-cell resolution can therefore substantially increase our understanding of these processes.Whereas measuring the RNA in bulk, that is, in samples containing thousands to hundreds of thousands of cells, is established and widely used since many years, methods for studying not only just a few different RNAs in individual cells became widely available only recently. Here, I outline and compare current concepts and methodologies for using single-cell RNA-sequencing to study virus infections. This covers sample preparation, cell preservation, biosafety considerations, and various experimental methods, with a special focus on the aspects that are important for studying virus infections. Since there is not "the one" method for doing single-cell RNA-sequencing, I will not provide a detailed protocol. Rather, this chapter should serve as a primer for getting started with single-cell RNA-sequencing experiments of virus infections and discusses the criteria that allow readers to choose the best procedures for their specific research question.

Transcriptomic analysis of sorted lung cells revealed a proviral activity of the NF-κB pathway toward SARS-CoV-2

Investigations of cellular responses to viral infection are commonly performed on mixed populations of infected and uninfected cells or using single-cell RNA sequencing, leading to inaccurate and low-resolution gene expression interpretations. Here, we performed deep polyA+ transcriptome analyses and novel RNA profiling of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infected lung epithelial cells, sorted based on the expression of the viral spike (S) protein. Infection caused a massive reduction in mRNAs and long non-coding RNAs (lncRNAs), including transcripts coding for antiviral factors, such as interferons (IFNs). This absence of IFN signaling probably explained the poor transcriptomic response of bystander cells co-cultured with S+ ones. NF-κB pathway and the inflammatory response escaped the global shutoff in S+ cells. Functional investigations revealed the proviral function of the NF-κB pathway and the antiviral activity of CYLD, a negative regulator of the pathway. Thus, our transcriptomic analysis on sorted cells revealed additional genes that modulate SARS-CoV-2 replication in lung cells.

Ancestral, Delta, and Omicron (BA.1) SARS-CoV-2 strains are dependent on serine proteases for entry throughout the human respiratory tract

BACKGROUND

The SARS-CoV-2 Omicron BA.1 variant emerged in late 2021 and became the globally dominant variant by January 2022. Authentic virus and pseudovirus systems have shown Omicron spike has an increased dependence on the endosomal pathway for entry.

METHODS

We investigated the entry mechanisms of Omicron, Delta, and ancestral viruses in cell models that represent different parts of the human respiratory tract, including nasal epithelial cells (hNECs), large-airway epithelial cells (LAECs), small-airway epithelial cells, and embryonic stem cell-derived type II alveolar cells.

FINDINGS

Omicron had an early replication advantage in LAECs, while Delta grew to higher titers in all cells. Omicron maintained dependence on serine proteases for entry in all culture systems. While serine protease inhibition with camostat was less robust for Omicron in hNECs, endosomal entry was not enhanced.

CONCLUSIONS

Our findings demonstrate that entry of Omicron BA.1 SARS-CoV-2 is dependent on serine proteases for entry throughout the respiratory tract.

FUNDING

This work was supported by The Medical Research Future Fund (MRF9200007; K.S., J.M.P.) and the DHHS Victorian State Government grant (Victorian State Government; DJPR/COVID-19; K.S, J.M.P.). K.S. is supported by a National Health and Medical Research Council of Australia Investigator grant (APP1177174).

Targeting Cyclophilin A and CD147 to Inhibit Replication of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and SARS-CoV-2-Induced Inflammation

Identification and development of effective therapeutics for coronavirus disease 2019 (COVID-19) are still urgently needed. The CD147-spike interaction is involved in the severe acute respiratory syndrome coronavirus (SARS-CoV)-2 invasion process in addition to angiotensin-converting enzyme 2 (ACE2). Cyclophilin A (CyPA), the extracellular ligand of CD147, has been found to play a role in the infection and replication of coronaviruses. In this study, our results show that CyPA inhibitors such as cyclosporine A (CsA) and STG-175 can suppress the intracellular replication of SARS-CoV-2 by inhibiting the binding of CyPA to the SARS-CoV-2 nucleocapsid C-terminal domain (N-CTD), and the IC50 is 0.23 μM and 0.17 μM, respectively. Due to high homology, CsA also had inhibitory effects on SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV), and the IC50 is 3.2 μM and 2.8 μM, respectively. Finally, we generated a formulation of phosphatidylserine (PS)-liposome-CsA for pulmonary drug delivery. These findings provide a scientific basis for identifying CyPA as a potential drug target for the treatment of COVID-19 as well as for the development of broad-spectrum inhibitors for coronavirus via targeting CyPA. Highlights: 1) SARS-CoV-2 infects cells via the binding of its S protein and CD147; 2) binding of SARS-CoV-2 N protein and CyPA is essential for viral replication; 3) CD147 and CyPA are potential therapeutic targets for SARS-CoV-2; and 4) CsA is a potential therapeutic strategy by interrupting CD147/CyPA interactions. SIGNIFICANCE STATEMENT: New severe acute respiratory syndrome coronavirus (SARS-CoV)-2 variants and other pathogenic coronaviruses (CoVs) are continually emerging, and new broad-spectrum anti-CoV therapy is urgently needed. We found that binding sites of cyclophilin A/cyclosporin A (CyPA/CsA) overlap with CyPA/N-CTD (nucleocapsid C-terminal domain), which shows the potential to target CyPA during SARS-CoV-2 infection. Here, we provide new evidence for targeting CyPA in the treatment of coronavirus disease 2019 (COVID-19) as well as the potential of developing CyPA inhibitors for broad-spectrum inhibition of CoVs.

Nasal spray of an IgM-like ACE2 fusion protein HH-120 prevents SARS-CoV-2 infection: Two investigator-initiated postexposure prophylaxis trials

HH-120, an IgM-like angiotensin converting enzyme 2 (ACE2) fusion protein, has been developed as a nasal spray against Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and is currently undergoing human trials. HH-120 nasal spray was assessed for postexposure prophylaxis (PEP) in two investigator-initiated (NS01 and NS02) trials with different risk levels of SARS-CoV-2 exposure. NS01 enrolled family caregiver participants who had continuous contacts with laboratory-confirmed index cases; NS02 enrolled participants who had general contacts (Part 1) or close contacts (Part 2) with index cases. The primary endpoints were safety and laboratory-confirmed and/or symptomatic SARS-CoV-2 infection. In NS01 trial (14 participants), the SARS-CoV-2 infection rates were 25% in the HH-120 group and 83.3% in the external control group (relative risk reduction [RRR]: 70.0%). In NS02-Part 1 (193 participants), the infection rates were 4% (HH-120) versus 11.3% (placebo), symptomatic infection rates were 0.8% versus 3.5%, hence with a RRR of 64.6% and 77.1%, respectively. In Part 2 (76 participants), the infection rates were 17.1% (HH-120) versus 30.4% (placebo), symptomatic infection rates were 7.5% versus 27.3%, with a RRR of 43.8% and 72.5%, respectively. No HH-120-related serious adverse effects were observed. The HH-120 nasal spray used as PEP was safe and effective in preventing laboratory-confirmed and symptomatic SARS-CoV-2 infection.

A High-Throughput, High-Containment Human Primary Epithelial Airway Organ-on-Chip Platform for SARS-CoV-2 Therapeutic Screening

COVID-19 emerged as a worldwide pandemic in early 2020, and while the rapid development of safe and efficacious vaccines stands as an extraordinary achievement, the identification of effective therapeutics has been less successful. This process has been limited in part by a lack of human-relevant preclinical models compatible with therapeutic screening on the native virus, which requires a high-containment environment. Here, we report SARS-CoV-2 infection and robust viral replication in PREDICT96-ALI, a high-throughput, human primary cell-based organ-on-chip platform. We evaluate unique infection kinetic profiles across lung tissue from three human donors by immunofluorescence, RT-qPCR, and plaque assays over a 6-day infection period. Enabled by the 96 devices/plate throughput of PREDICT96-ALI, we also investigate the efficacy of Remdesivir and MPro61 in a proof-of-concept antiviral study. Both compounds exhibit an antiviral effect against SARS-CoV-2 in the platform. This demonstration of SARS-CoV-2 infection and antiviral dosing in a high-throughput organ-on-chip platform presents a critical capability for disease modeling and therapeutic screening applications in a human physiology-relevant in vitro system.

Causal identification of single-cell experimental perturbation effects with CINEMA-OT

Recent advancements in single-cell technologies allow characterization of experimental perturbations at single-cell resolution. While methods have been developed to analyze such experiments, the application of a strict causal framework has not yet been explored for the inference of treatment effects at the single-cell level. Here we present a causal-inference-based approach to single-cell perturbation analysis, termed CINEMA-OT (causal independent effect module attribution + optimal transport). CINEMA-OT separates confounding sources of variation from perturbation effects to obtain an optimal transport matching that reflects counterfactual cell pairs. These cell pairs represent causal perturbation responses permitting a number of novel analyses, such as individual treatment-effect analysis, response clustering, attribution analysis, and synergy analysis. We benchmark CINEMA-OT on an array of treatment-effect estimation tasks for several simulated and real datasets and show that it outperforms other single-cell perturbation analysis methods. Finally, we perform CINEMA-OT analysis of two newly generated datasets: (1) rhinovirus and cigarette-smoke-exposed airway organoids, and (2) combinatorial cytokine stimulation of immune cells. In these experiments, CINEMA-OT reveals potential mechanisms by which cigarette-smoke exposure dulls the airway antiviral response, as well as the logic that governs chemokine secretion and peripheral immune cell recruitment.

Chronic alcohol use primes bronchial cells for altered inflammatory response and barrier dysfunction during SARS-CoV-2 infection

Alcohol use disorder (AUD) is a significant public health concern and people with AUD are more likely to develop severe acute respiratory distress syndrome (ARDS) in response to respiratory infections. To examine whether AUD was a risk factor for more severe outcome in response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, we examined early responses to infection using cultured differentiated bronchial epithelial cells derived from brushings obtained from people with AUD or without AUD. RNA-seq analysis of uninfected cells determined that AUD cells were enriched for expression of epidermal genes as compared with non-AUD cells. Bronchial epithelial cells from patients with AUD showed a significant decrease in barrier function 72 h postinfection, as determined by transepithelial electrical resistance. In contrast, barrier function of non-AUD cells was enhanced 72 h after SARS-CoV-2 infection. AUD cells showed claudin-7 that did not colocalize with zonula occludens-1 (ZO-1), indicative of disorganized tight junctions. However, both AUD and non-AUD cells showed decreased β-catenin expression following SARS-CoV-2 infection. To determine the impact of AUD on the inflammatory response to SARS-CoV-2 infection, cytokine secretion was measured by multiplex analysis. SARS-CoV-2-infected AUD bronchial cells had enhanced secretion of multiple proinflammatory cytokines including TNFα, IL-1β, and IFNγ as opposed to non-AUD cells. In contrast, secretion of the barrier-protective cytokines epidermal growth factor (EGF) and granulocyte macrophage-colony stimulating factor (GM-CSF) was enhanced for non-AUD bronchial cells. Taken together, these data support the hypothesis that AUD is a risk factor for COVID-19, where alcohol primes airway epithelial cells for increased inflammation and increased barrier dysfunction and increased inflammation in response to infection by SARS-CoV-2.NEW & NOTEWORTHY Alcohol use disorder (AUD) is a significant risk factor for severe acute respiratory distress syndrome. We found that AUD causes a phenotypic shift in gene expression in human bronchial epithelial cells, enhancing expression of epidermal genes. AUD cells infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) had higher levels of proinflammatory cytokine secretion and barrier dysfunction not present in infected non-AUD cells, consistent with increased early COVID-19 severity due to AUD.

Single-cell RNA-seq public data reveal the gene regulatory network landscape of respiratory epithelial and peripheral immune cells in COVID-19 patients

Introduction

Infection with SARS-CoV-2 leads to coronavirus disease 2019 (COVID-19), which can result in acute respiratory distress syndrome and multiple organ failure. However, its comprehensive influence on pathological immune responses in the respiratory epithelium and peripheral immune cells is not yet fully understood.

Methods

In this study, we analyzed multiple public scRNA-seq datasets of nasopharyngeal swabs and peripheral blood to investigate the gene regulatory networks (GRNs) of healthy individuals and COVID-19 patients with mild/moderate and severe disease, respectively. Cell-cell communication networks among cell types were also inferred. Finally, validations were conducted using bulk RNA-seq and proteome data.

Results

Similar and dissimilar regulons were identified within or between epithelial and immune cells during COVID-19 severity progression. The relative transcription factors (TFs) and their targets were used to construct GRNs among different infection sites and conditions. Between respiratory epithelial and peripheral immune cells, different TFs tended to be used to regulate the activity of a cell between healthy individuals and COVID-19 patients, although they had some TFs in common. For example, XBP1, FOS, STAT1, and STAT2 were activated in both the epithelial and immune cells of virus-infected individuals. In contrast, severe COVID-19 cases exhibited activation of CEBPD in peripheral immune cells, while CEBPB was exclusively activated in respiratory epithelial cells. Moreover, in patients with severe COVID-19, although some inflammatory genes, such as S100A8/A9, were found to be upregulated in both respiratory epithelial and peripheral immune cells, their relative regulators can differ in terms of cell types. The cell-cell communication analysis suggested that epidermal growth factor receptor signaling among epithelia contributes to mild/moderate disease, and chemokine signaling among immune cells contributes to severe disease.

Conclusion

This study identified cell type- and condition-specific regulons in a wide range of cell types from the initial infection site to the peripheral blood, and clarified the diverse mechanisms of maladaptive responses to SARS-CoV-2 infection.

Systematic functional interrogation of SARS-CoV-2 host factors using Perturb-seq

Genomic and proteomic screens have identified numerous host factors of SARS-CoV-2, but efficient delineation of their molecular roles during infection remains a challenge. Here we use Perturb-seq, combining genetic perturbations with a single-cell readout, to investigate how inactivation of host factors changes the course of SARS-CoV-2 infection and the host response in human lung epithelial cells. Our high-dimensional data resolve complex phenotypes such as shifts in the stages of infection and modulations of the interferon response. However, only a small percentage of host factors showed such phenotypes upon perturbation. We further identified the NF-κB inhibitor IκBα (NFKBIA), as well as the translation factors EIF4E2 and EIF4H as strong host dependency factors acting early in infection. Overall, our study provides massively parallel functional characterization of host factors of SARS-CoV-2 and quantitatively defines their roles both in virus-infected and bystander cells.

The International Virus Bioinformatics Meeting 2023

The 2023 International Virus Bioinformatics Meeting was held in Valencia, Spain, from 24-26 May 2023, attracting approximately 180 participants worldwide. The primary objective of the conference was to establish a dynamic scientific environment conducive to discussion, collaboration, and the generation of novel research ideas. As the first in-person event following the SARS-CoV-2 pandemic, the meeting facilitated highly interactive exchanges among attendees. It served as a pivotal gathering for gaining insights into the current status of virus bioinformatics research and engaging with leading researchers and emerging scientists. The event comprised eight invited talks, 19 contributed talks, and 74 poster presentations across eleven sessions spanning three days. Topics covered included machine learning, bacteriophages, virus discovery, virus classification, virus visualization, viral infection, viromics, molecular epidemiology, phylodynamic analysis, RNA viruses, viral sequence analysis, viral surveillance, and metagenomics. This report provides rewritten abstracts of the presentations, a summary of the key research findings, and highlights shared during the meeting.