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

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.

Transcriptomics-inferred dynamics of SARS-CoV-2 interactions with host epithelial cells

Virus-host interactions can reveal potentially effective and selective therapeutic targets for treating infection. Here, we performed an integrated analysis of the dynamics of virus replication and the host cell transcriptional response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection using human Caco-2 colon cancer cells as a model. Time-resolved RNA sequencing revealed that, upon infection, cells immediately transcriptionally activated genes associated with inflammatory pathways that mediate the antiviral response, which was followed by an increase in the expression of genes involved in ribosome and mitochondria function, thus suggesting rapid alterations in protein production and cellular energy supply. At later stages, between 24 and 48 hours after infection, the expression of genes involved in metabolic processes-in particular, those related to xenobiotic metabolism-was decreased. Mathematical modeling incorporating SARS-CoV-2 replication suggested that SARS-CoV-2 proteins inhibited the host antiviral response and that virus transcripts exceeded the translation capacity of the host cells. Targeting kinase-dependent pathways that exhibited increases in transcription in host cells was as effective as a virus-targeted inhibitor at repressing viral replication. Our findings in this model system delineate a sequence of SARS-CoV-2 virus-host interactions that may facilitate the identification of druggable host pathways to suppress infection.

MAP3K19 regulatory variation in populations with African ancestry may increase COVID-19 severity

To identify ancestry-linked genetic risk variants associated with COVID-19 hospitalization, we performed an integrative analysis of two genome-wide association studies and resolved four single nucleotide polymorphisms more frequent in COVID-19-hospitalized patients with non-European ancestry. Among them, the COVID-19 risk SNP rs16831827 shows the largest difference in minor allele frequency (MAF) between populations with African and European ancestry and also shows higher MAF in hospitalized COVID-19 patients among cohorts of mixed ancestry (odds ratio [OR] = 1.20, 95% CI: 1.10-1.30) and entirely African ancestry (OR = 1.30, 95% CI: 1.02-1.67). rs16831827 is an expression quantitative trait locus of MAP3K19. MAP3K19 expression is induced during ciliogenesis and most abundant in ciliated tissues including lungs. Single-cell RNA sequencing analyses revealed that MAP3K19 is highly expressed in multiple ciliated cell types. As rs16831827∗T is associated with reduced MAP3K19 expression, it may increase the risk of severe COVID-19 by reducing MAP3K19 expression.

SARS-CoV-2 infection induces a long-lived pro-inflammatory transcriptional profile

BACKGROUND

The immune response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in COVID-19 patients has been extensively investigated. However, much less is known about the long-term effects of infection in patients and how it could affect the immune system and its capacity to respond to future perturbations.

METHODS

Using a targeted single-cell multiomics approach, we have recently identified a prolonged anti-inflammatory gene expression signature in T and NK cells in type 1 diabetes patients treated with low-dose IL-2. Here, we investigated the dynamics of this signature in three independent cohorts of COVID-19 patients: (i) the Oxford COVID-19 Multi-omics Blood Atlas (COMBAT) dataset, a cross-sectional cohort including 77 COVID-19 patients and ten healthy donors; (ii) the INCOV dataset, consisting of 525 samples taken from 209 COVID-19 patients during and after infection; and (iii) a longitudinal dataset consisting of 269 whole-blood samples taken from 139 COVID-19 patients followed for a period of up to 7 months after the onset of symptoms using a bulk transcriptomic approach.

RESULTS

We discovered that SARS-CoV-2 infection leads to a prolonged alteration of the gene expression profile of circulating T, B and NK cells and monocytes. Some of the genes affected were the same as those present in the IL-2-induced anti-inflammatory gene expression signature but were regulated in the opposite direction, implying a pro-inflammatory status. The altered transcriptional profile was detected in COVID-19 patients for at least 2 months after the onset of the disease symptoms but was not observed in response to influenza infection or sepsis. Gene network analysis suggested a central role for the transcriptional factor NF-κB in the regulation of the observed transcriptional alterations.

CONCLUSIONS

SARS-CoV-2 infection causes a prolonged increase in the pro-inflammatory transcriptional status that could predispose post-acute patients to the development of long-term health consequences, including autoimmune disease, reactivation of other viruses and disruption of the host immune system-microbiome ecosystem.

Exploring the Immune Response against RSV and SARS-CoV-2 Infection in Children

Viral respiratory tract infections are a significant public health concern, particularly in children. RSV is a prominent cause of lower respiratory tract infections among infants, whereas SARS-CoV-2 has caused a global pandemic with lower overall severity in children than in adults. In this review, we aimed to compare the innate and adaptive immune responses induced by RSV and SARS-CoV-2 to better understand differences in the pathogenesis of infection. Some studies have demonstrated that children present a more robust immune response against SARS-CoV-2 than adults; however, this response is dissimilar to that of RSV. Each virus has a distinctive mechanism to escape the immune response. Understanding the mechanisms underlying these differences is crucial for developing effective treatments and improving the management of pediatric respiratory infections.

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

SARS-CoV-2 initiates infection in the conducting airways, which rely on mucocilliary clearance (MCC) to minimize pathogen penetration. However, it is unclear how MCC impacts SARS-CoV-2 spread after infection is established. To understand viral spread at this site, we performed live imaging of SARS-CoV-2 infected differentiated primary human bronchial epithelium cultures for up to 9 days. Fluorescent markers for cilia and mucus allowed longitudinal monitoring of MCC, ciliary motion, and infection. The number of infected cells peaked at 4 days post-infection in characteristic foci that followed mucus movement. Inhibition of MCC using physical and genetic perturbations limited foci. Later in infection, MCC was diminished despite relatively subtle ciliary function defects. Resumption of MCC and infection spread after mucus removal suggests that mucus secretion mediates this effect. We show that MCC facilitates SARS-CoV-2 spread early in infection while later decreases in MCC inhibit spread, suggesting a complex interplay between SARS-CoV-2 and MCC.

Unambiguous detection of SARS-CoV-2 subgenomic mRNAs with single-cell RNA sequencing

Single-cell RNA sequencing (scRNA-Seq) studies have provided critical insight into the pathogenesis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19). scRNA-Seq library preparation methods and data processing workflows are generally designed for the detection and quantification of eukaryotic host mRNAs and not viral RNAs. Here, we compare different scRNA-Seq library preparation methods for their ability to quantify and detect SARS-CoV-2 RNAs with a focus on subgenomic mRNAs (sgmRNAs). We show that compared to 10X Genomics Chromium Next GEM Single Cell 3' (10X 3') libraries or 10X Genomics Chromium Next GEM Single Cell V(D)J (10X 5') libraries sequenced with standard read configurations, 10X 5' libraries sequenced with an extended length read 1 (R1) that covers both cell barcode and transcript sequence (termed "10X 5' with extended R1") increase the number of unambiguous reads spanning leader-sgmRNA junction sites. We further present a data processing workflow, single-cell coronavirus sequencing (scCoVseq), which quantifies reads unambiguously assigned to viral sgmRNAs or viral genomic RNA (gRNA). We find that combining 10X 5' with extended R1 library preparation/sequencing and scCoVseq data processing maximizes the number of viral UMIs per cell quantified by scRNA-Seq. Corresponding sgmRNA expression levels are highly correlated with expression in matched bulk RNA-Seq data sets quantified with established tools for SARS-CoV-2 analysis. Using this scRNA-Seq approach, we find that SARS-CoV-2 gene expression is highly correlated across individual infected cells, which suggests that the proportion of viral sgmRNAs remains generally consistent throughout infection. Taken together, these results and corresponding data processing workflow enable robust quantification of coronavirus sgmRNA expression at single-cell resolution, thereby supporting high-resolution studies of viral RNA processes in individual cells. IMPORTANCE Single-cell RNA sequencing (scRNA-Seq) has emerged as a valuable tool to study host-virus interactions, especially for coronavirus disease 2019 (COVID-19). Here we compare the performance of different scRNA-Seq library preparation methods and sequencing strategies to detect SARS-CoV-2 RNAs and develop a data processing workflow to quantify unambiguous sequence reads derived from SARS-CoV-2 genomic RNA and subgenomic mRNAs. After establishing a workflow that maximizes the detection of SARS-CoV-2 subgenomic mRNAs, we explore patterns of SARS-CoV-2 gene expression across cells with variable levels of total viral RNA, assess host gene expression differences between infected and bystander cells, and identify non-canonical and lowly abundant SARS-CoV-2 RNAs. The sequencing and data processing strategies developed here can enhance studies of coronavirus RNA biology at single-cell resolution and thereby contribute to our understanding of viral pathogenesis.

The Dual Role of the Airway Epithelium in Asthma: Active Barrier and Regulator of Inflammation

Chronic airway inflammation is the cornerstone on which bronchial asthma arises, and in turn, chronic inflammation arises from a complex interplay between environmental factors such as allergens and pathogens and immune cells as well as structural cells constituting the airway mucosa. Airway epithelial cells (AECs) are at the center of these processes. On the one hand, they represent the borderline separating the body from its environment in order to keep inner homeostasis. The airway epithelium forms a multi-tiered, self-cleaning barrier that involves an unstirred, discontinuous mucous layer, the dense and rigid mesh of the glycocalyx, and the cellular layer itself, consisting of multiple, densely interconnected cell types. On the other hand, the airway epithelium represents an immunologically highly active tissue once its barrier has been penetrated: AECs play a pivotal role in releasing protective immunoglobulin A. They express a broad spectrum of pattern recognition receptors, enabling them to react to environmental stressors that overcome the mucosal barrier. By releasing alarmins-proinflammatory and regulatory cytokines-AECs play an active role in the formation, strategic orientation, and control of the subsequent defense reaction. Consequently, the airway epithelium is of vital importance to chronic inflammatory diseases, such as asthma.

A comparative study of in vitro air-liquid interface culture models of the human airway epithelium evaluating cellular heterogeneity and gene expression at single cell resolution

BACKGROUND

The airway epithelium is composed of diverse cell types with specialized functions that mediate homeostasis and protect against respiratory pathogens. Human airway epithelial (HAE) cultures at air-liquid interface are a physiologically relevant in vitro model of this heterogeneous tissue and have enabled numerous studies of airway disease. HAE cultures are classically derived from primary epithelial cells, the relatively limited passage capacity of which can limit experimental methods and study designs. BCi-NS1.1, a previously described and widely used basal cell line engineered to express hTERT, exhibits extended passage lifespan while retaining the capacity for differentiation to HAE. However, gene expression and innate immune function in BCi-NS1.1-derived versus primary-derived HAE cultures have not been fully characterized.

METHODS

BCi-NS1.1-derived HAE cultures (n = 3 independent differentiations) and primary-derived HAE cultures (n = 3 distinct donors) were characterized by immunofluorescence and single cell RNA-Seq (scRNA-Seq). Innate immune functions were evaluated in response to interferon stimulation and to infection with viral and bacterial respiratory pathogens.

RESULTS

We confirm at high resolution that BCi-NS1.1- and primary-derived HAE cultures are largely similar in morphology, cell type composition, and overall gene expression patterns. While we observed cell-type specific expression differences of several interferon stimulated genes in BCi-NS1.1-derived HAE cultures, we did not observe significant differences in susceptibility to infection with influenza A virus and Staphylococcus aureus.

CONCLUSIONS

Taken together, our results further support BCi-NS1.1-derived HAE cultures as a valuable tool for the study of airway infectious disease.

LY6E is a pan-coronavirus restriction factor in the respiratory tract

LY6E is an antiviral restriction factor that inhibits coronavirus spike-mediated fusion, but the cell types in vivo that require LY6E for protection from respiratory coronavirus infection are unknown. Here we used a panel of seven conditional Ly6e knockout mice to define which Ly6e-expressing cells confer control of airway infection by murine coronavirus and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Loss of Ly6e in Lyz2-expressing cells, radioresistant Vav1-expressing cells and non-haematopoietic cells increased susceptibility to murine coronavirus. Global conditional loss of Ly6e expression resulted in clinical disease and higher viral burden after SARS-CoV-2 infection, but little evidence of immunopathology. We show that Ly6e expression protected secretory club and ciliated cells from SARS-CoV-2 infection and prevented virus-induced loss of an epithelial cell transcriptomic signature in the lung. Our study demonstrates that lineage confined rather than broad expression of Ly6e sufficiently confers resistance to disease caused by murine and human coronaviruses.