Integrated analysis of rectal mucosal microbiome and transcriptome reveals a distinct microenvironment among young MSM.

Cross-talk between the microbiome and gut mucosa-resident immune cells plays a pivotal role in modulating immune responses to pathogens, including responses to HIV infection. However, how these interactions may differ between young MSM (YMSM) disproportionately impacted by HIV, as compared to older adult MSM (AMSM), is not well understood. A broad analysis of associations between the microbiome and rectal transcriptome revealed 10 microbial families/genera correlated with immunologic gene pathways. Specifically, the rectal transcriptome of YMSM is characterized by upregulation of T cell activation/differentiation pathways and signaling from multiple cytokine families, compared to AMSM. The microbiome of YMSM is enriched with pathogenic genera including Peptostreptococcus, shown to be positively correlated with type I interferon pathways important for antiviral immunity. These findings demonstrate that YMSM have a unique immune phenotype and rectal microenvironment and support further evaluation of biological factors that influence rectal HIV transmission.

Distinct SIV-specific CD8+ T cells in the lymph node exhibit simultaneous effector and stem-like profiles and are associated with limited SIV persistence.

Human immunodeficiency virus (HIV) cure efforts are increasingly focused on harnessing CD8+ T cell functions, which requires a deeper understanding of CD8+ T cells promoting HIV control. Here we identifiy an antigen-responsive TOXhiTCF1+CD39+CD8+ T cell population with high expression of inhibitory receptors and low expression of canonical cytolytic molecules. Transcriptional analysis of simian immunodeficiency virus (SIV)-specific CD8+ T cells and proteomic analysis of purified CD8+ T cell subsets identified TOXhiTCF1+CD39+CD8+ T cells as intermediate effectors that retained stem-like features with a lineage relationship with terminal effector T cells. TOXhiTCF1+CD39+CD8+ T cells were found at higher frequency than TCF1-CD39+CD8+ T cells in follicular microenvironments and were preferentially located in proximity of SIV-RNA+ cells. Their frequency was associated with reduced plasma viremia and lower SIV reservoir size. Highly similar TOXhiTCF1+CD39+CD8+ T cells were detected in lymph nodes from antiretroviral therapy-naive and antiretroviral therapy-suppressed people living with HIV, suggesting this population of CD8+ T cells contributes to limiting SIV and HIV persistence.

Single-cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells.

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10× Genomics Chromium single-cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human monocyte-derived dendritic cells infected with ZIKV at the single-cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN-dependent and -independent genes (the antiviral module). We modeled the ZIKV-specific antiviral state at the protein level, leveraging experimentally derived protein interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per-cell basis with experimental protein interaction data. IMPORTANCE: Zika virus (ZIKV) remains a public health threat given its potential for re-emergence and the detrimental fetal outcomes associated with infection during pregnancy. Understanding the dynamics between ZIKV and its host is critical to understanding ZIKV pathogenesis. Through ZIKV-inclusive single-cell RNA sequencing (scRNA-seq), we demonstrate on the single-cell level the dynamic interplay between ZIKV and the host: the transcriptional program that restricts viral infection and ZIKV-mediated inhibition of that response. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool for gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

The CARD8 inflammasome dictates HIV/SIV pathogenesis and disease progression.

While CD4+ T cell depletion is key to disease progression in people living with HIV and SIV-infected macaques, the mechanisms underlying this depletion remain incompletely understood, with most cell death involving uninfected cells. In contrast, SIV infection of "natural" hosts such as sooty mangabeys does not cause CD4+ depletion and AIDS despite high-level viremia. Here, we report that the CARD8 inflammasome is activated immediately after HIV entry by the viral protease encapsulated in incoming virions. Sensing of HIV protease activity by CARD8 leads to rapid pyroptosis of quiescent cells without productive infection, while T cell activation abolishes CARD8 function and increases permissiveness to infection. In humanized mice reconstituted with CARD8-deficient cells, CD4+ depletion is delayed despite high viremia. Finally, we discovered loss-of-function mutations in CARD8 from "natural hosts," which may explain the peculiarly non-pathogenic nature of these infections. Our study suggests that CARD8 drives CD4+ T cell depletion during pathogenic HIV/SIV infections.

Single cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells.

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10x Genomics Chromium single cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human moDCs infected with ZIKV at the single cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN dependent and independent genes (antiviral module). We modeled the ZIKV specific antiviral state at the protein level leveraging experimentally derived protein-interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per cell basis with experimental protein interaction data. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool to gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

Vpr attenuates antiviral immune responses and is critical for full pathogenicity of SIVmac239 in rhesus macaques.

The accessory viral protein R (Vpr) is encoded by all primate lentiviruses. Vpr counteracts DNA repair pathways, modulates viral immune sensing, and induces cell-cycle arrest in cell culture. However, its impact in vivo is controversial. Here, we show that deletion of vpr is associated with delayed viral replication kinetics, rapid innate immune activation, development and maintenance of strong B and T cell responses, and increased neutralizing activity against SIVmac239 in rhesus macaques. All wild-type SIVmac239-infected animals maintained high viral loads, and five of six developed fatal immunodeficiency during ∼80 weeks of follow-up. Lack of Vpr was associated with better preservation of CD4+ T cells, lower viral loads, and an attenuated clinical course of infection in most animals. Our results show that Vpr contributes to efficient viral immune evasion and the full pathogenic potential of SIVmacin vivo. Inhibition of Vpr may improve humoral immune control of viral replication.

Systems analysis reveals differential expression of endocervical genes in African women randomized to DMPA-IM, LNG implant or cu-IUD.

Although effective contraceptives are crucial for preventing unintended pregnancies, evidence suggests that their use may perturb the female genital tract (FGT). A comparative analysis of the effects of the most common contraceptives on the FGT have not been evaluated in a randomized clinical trial setting. Here, we evaluated the effect of three long-acting contraceptive methods: depot medroxyprogesterone acetate(DMPA-IM), levonorgestrel(LNG) implant, and a copper intrauterine device (Cu-IUD), on the endocervical host transcriptome in 188 women from the Evidence for Contraceptive Options and HIV Outcomes Trial (ECHO) trial. Cu-IUD usage showed the most extensive transcriptomic changes, and was associated with inflammatory and anti-viral host responses. DMPA-IM usage was enriched for pathways associated with T cell responses. LNG implant had the mildest effect on endocervical gene expression, and was associated with growth factor signaling. These data provide a mechanistic basis for the diverse influence that varying contraceptives have on the FGT.

Intradermal but not intramuscular modified vaccinia Ankara immunizations protect against intravaginal tier2 simian-human immunodeficiency virus challenges in female macaques.

Route of immunization can markedly influence the quality of immune response. Here, we show that intradermal (ID) but not intramuscular (IM) modified vaccinia Ankara (MVA) vaccinations provide protection from acquisition of intravaginal tier2 simian-human immunodeficiency virus (SHIV) challenges in female macaques. Both routes of vaccination induce comparable levels of serum IgG with neutralizing and non-neutralizing activities. The protection in MVA-ID group correlates positively with serum neutralizing and antibody-dependent phagocytic activities, and envelope-specific vaginal IgA; while the limited protection in MVA-IM group correlates only with serum neutralizing activity. MVA-ID immunizations induce greater germinal center Tfh and B cell responses, reduced the ratio of Th1 to Tfh cells in blood and showed lower activation of intermediate monocytes and inflammasome compared to MVA-IM immunizations. This lower innate activation correlates negatively with induction of Tfh responses. These data demonstrate that the MVA-ID vaccinations protect against intravaginal SHIV challenges by modulating the innate and T helper responses.

Modulation of type I interferon responses potently inhibits SARS-CoV-2 replication and inflammation in rhesus macaques.

Type I interferons (IFN-I) are critical mediators of innate control of viral infections but also drive the recruitment of inflammatory cells to sites of infection, a key feature of severe coronavirus disease 2019. Here, IFN-I signaling was modulated in rhesus macaques (RMs) before and during acute SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) infection using a mutated IFN-α2 (IFN-modulator; IFNmod), which has previously been shown to reduce the binding and signaling of endogenous IFN-I. IFNmod treatment in uninfected RMs was observed to induce a modest up-regulation of only antiviral IFN-stimulated genes (ISGs); however, in SARS-CoV-2-infected RMs, IFNmod reduced both antiviral and inflammatory ISGs. IFNmod treatment resulted in a potent reduction in SARS-CoV-2 viral loads both in vitro in Calu-3 cells and in vivo in bronchoalveolar lavage (BAL), upper airways, lung, and hilar lymph nodes of RMs. Furthermore, in SARS-CoV-2-infected RMs, IFNmod treatment potently reduced inflammatory cytokines, chemokines, and CD163+ MRC1- inflammatory macrophages in BAL and expression of Siglec-1 on circulating monocytes. In the lung, IFNmod also reduced pathogenesis and attenuated pathways of inflammasome activation and stress response during acute SARS-CoV-2 infection. Using an intervention targeting both IFN-α and IFN-ß pathways, this study shows that, whereas early IFN-I restrains SARS-CoV-2 replication, uncontrolled IFN-I signaling critically contributes to SARS-CoV-2 inflammation and pathogenesis in the moderate disease model of RMs.

CD8+ T cells promote HIV latency by remodeling CD4+ T cell metabolism to enhance their survival, quiescence, and stemness.

HIV infection persists during antiretroviral therapy (ART) due to a reservoir of latently infected cells that harbor replication-competent virus and evade immunity. Previous ex vivo studies suggested that CD8+ T cells from people with HIV may suppress HIV expression via non-cytolytic mechanisms, but the mechanisms responsible for this effect remain unclear. Here, we used a primary cell-based in vitro latency model and demonstrated that co-culture of autologous activated CD8+ T cells with HIV-infected memory CD4+ T cells promoted specific changes in metabolic and/or signaling pathways resulting in increased CD4+ T cell survival, quiescence, and stemness. Collectively, these pathways negatively regulated HIV expression and ultimately promoted the establishment of latency. As shown previously, we observed that macrophages, but not B cells, promoted latency in CD4+ T cells. The identification of CD8-specific mechanisms of pro-latency activity may favor the development of approaches to eliminate the viral reservoir in people with HIV.

TREM2+ and interstitial-like macrophages orchestrate airway inflammation in SARS-CoV-2 infection in rhesus macaques.

The immunopathological mechanisms driving the development of severe COVID-19 remain poorly defined. Here, we utilize a rhesus macaque model of acute SARS-CoV-2 infection to delineate perturbations in the innate immune system. SARS-CoV-2 initiates a rapid infiltration of plasmacytoid dendritic cells into the lower airway, commensurate with IFNA production, natural killer cell activation, and a significant increase of blood CD14-CD16+ monocytes. To dissect the contribution of lung myeloid subsets to airway inflammation, we generate a longitudinal scRNA-Seq dataset of airway cells, and map these subsets to corresponding populations in the human lung. SARS-CoV-2 infection elicits a rapid recruitment of two macrophage subsets: CD163+MRC1-, and TREM2+ populations that are the predominant source of inflammatory cytokines. Treatment with baricitinib (Olumiant®), a JAK1/2 inhibitor is effective in eliminating the influx of non-alveolar macrophages, with a reduction of inflammatory cytokines. This study delineates the major lung macrophage subsets driving airway inflammation during SARS-CoV-2 infection.

Corrigendum: The Vacc-SeqQC project: Benchmarking RNA-Seq for clinical vaccine studies.

[This corrects the article DOI: 10.3389/fimmu.2022.1093242.].

Modulation of type I interferon responses potently inhibits SARS-CoV-2 replication and inflammation in rhesus macaques.

Type-I interferons (IFN-I) are critical mediators of innate control of viral infections, but also drive recruitment of inflammatory cells to sites of infection, a key feature of severe COVID-19. Here, and for the first time, IFN-I signaling was modulated in rhesus macaques (RMs) prior to and during acute SARS-CoV-2 infection using a mutated IFNα2 (IFN-modulator; IFNmod), which has previously been shown to reduce the binding and signaling of endogenous IFN-I. In SARS-CoV-2-infected RMs, IFNmod reduced both antiviral and inflammatory ISGs. Notably, IFNmod treatment resulted in a potent reduction in (i) SARS-CoV-2 viral load in Bronchoalveolar lavage (BAL), upper airways, lung, and hilar lymph nodes; (ii) inflammatory cytokines, chemokines, and CD163+MRC1-inflammatory macrophages in BAL; and (iii) expression of Siglec-1, which enhances SARS-CoV-2 infection and predicts disease severity, on circulating monocytes. In the lung, IFNmod also reduced pathogenesis and attenuated pathways of inflammasome activation and stress response during acute SARS-CoV-2 infection. This study, using an intervention targeting both IFN-α and IFN-ß pathways, shows that excessive inflammation driven by type 1 IFN critically contributes to SARS-CoV-2 pathogenesis in RMs, and demonstrates the potential of IFNmod to limit viral replication, SARS-CoV-2 induced inflammation, and COVID-19 severity.

Systems Analysis Reveals Contraceptive-Induced Alteration of Cervicovaginal Gene Expression in a Randomized Trial.

Hormonal contraceptives (HCs) are vital in managing the reproductive health of women. However, HC usage has been linked to perturbations in cervicovaginal immunity and increased risk of sexually transmitted infections. Here, we evaluated the impact of three HCs on the cervicovaginal environment using high-throughput transcriptomics. From 2015 to 2017, 130 adolescent females aged 15-19 years were enrolled into a substudy of UChoose, a single-site, open-label randomized, crossover trial (NCT02404038) and randomized to injectable norethisterone-enanthate (Net-En), combined oral contraceptives (COC), or etonorgesterol/ethinyl-estradiol-combined contraceptive vaginal ring (CCVR). Cervicovaginal samples were collected after 16 weeks of randomized HC use and analyzed by RNA-Seq, 16S rRNA gene sequencing, and Luminex analysis. Participants in the CCVR arm had a significant elevation of transcriptional networks driven by IL-6, IL-1, and NFKB, and lower expression of genes supporting epithelial barrier integrity. An integrated multivariate analysis demonstrated that networks of microbial dysbiosis and inflammation best discriminated the CCVR arm from the other contraceptive groups, while genes involved in epithelial cell differentiation were predictive of the Net-En and COC arms. Collectively, these data from a randomized trial represent the most comprehensive "omics" analyses of the cervicovaginal response to HCs and provide important mechanistic guidelines for the provision of HCs in sub-Saharan Africa.

T-cell activation and B-cell interaction signatures in rectal tissues are associated with HIV replication in ex-vivo model of infection.

OBJECTIVE: The rectal mucosa is a critical site of HIV vulnerability. We sought to identify transcriptomic features of rectal mucosal tissue prior to exposure associated with support or restriction of HIV replication. DESIGN: Rectal tissue from HIV-negative cis gender men ( n  = 57) underwent concurrent RNAseq transcriptomic analyses (two biopsies/participant) and challenge with HIV in the ex-vivo explant model of infection (three biopsies challenged/participant) as part of a larger cohort study to understand the rectal mucosal immune environment among MSM. METHODS: P24 was quantified in the explant supernatants over a culture period of 18 days via ELISA. Participant median p24 log area under the curve was correlated with bulk transcriptomic data (Illumina HiSeq3000) to identify associations between gene expression and p24 production. Significant differentially expressed genes (DEGs) were identified via DESeq2 analysis and analyzed with Reactome to identify pathways of interest. RESULTS: In total, 183 DEG (181 upregulated, two downregulated) were associated with higher p24 accumulation in the ex-vivo challenge model, including T-cell activation, B-cell function, and chemokine DEG. Reactome analysis of the upregulated genes identified 'Adaptive Immune System', 'Cytokine Signaling in Immune System', and 'Innate Immune System' as significantly upregulated pathways. CONCLUSION: For the first time, we identified rectal tissue transcriptomic signatures associated with increased p24 production utilizing an ex-vivo model. Our findings are highly relevant to HIV transmission and the early establishment of HIV reservoirs in humans, and future studies should examine the identified pathways as targets for new or improved biomedical prevention or treatment interventions.

The role of CD101-expressing CD4 T cells in HIV/SIV pathogenesis and persistence.

Despite the advent of effective antiretroviral therapy (ART), human immunodeficiency virus (HIV) continues to pose major challenges, with extensive pathogenesis during acute and chronic infection prior to ART initiation and continued persistence in a reservoir of infected CD4 T cells during long-term ART. CD101 has recently been characterized to play an important role in CD4 Treg potency. Using the simian immunodeficiency virus (SIV) model of HIV infection in rhesus macaques, we characterized the role and kinetics of CD101+ CD4 T cells in longitudinal SIV infection. Phenotypic analyses and single-cell RNAseq profiling revealed that CD101 marked CD4 Tregs with high immunosuppressive potential, distinct from CD101- Tregs, and these cells also were ideal target cells for HIV/SIV infection, with higher expression of CCR5 and α4ß7 in the gut mucosa. Notably, during acute SIV infection, CD101+ CD4 T cells were preferentially depleted across all CD4 subsets when compared with their CD101- counterpart, with a pronounced reduction within the Treg compartment, as well as significant depletion in mucosal tissue. Depletion of CD101+ CD4 was associated with increased viral burden in plasma and gut and elevated levels of inflammatory cytokines. While restored during long-term ART, the reconstituted CD101+ CD4 T cells display a phenotypic profile with high expression of inhibitory receptors (including PD-1 and CTLA-4), immunsuppressive cytokine production, and high levels of Ki-67, consistent with potential for homeostatic proliferation. Both the depletion of CD101+ cells and phenotypic profile of these cells found in the SIV model were confirmed in people with HIV on ART. Overall, these data suggest an important role for CD101-expressing CD4 T cells at all stages of HIV/SIV infection and a potential rationale for targeting CD101 to limit HIV pathogenesis and persistence, particularly at mucosal sites.

Altered Response Pattern following AZD5582 Treatment of SIV-Infected, ART-Suppressed Rhesus Macaque Infants.

The "shock and kill" strategy for HIV-1 cure incorporates latency-reversing agents (LRA) in combination with interventions that aid the host immune system in clearing virally reactivated cells. LRAs have not yet been investigated in pediatric clinical or preclinical studies. Here, we evaluated an inhibitor of apoptosis protein (IAP) inhibitor (IAPi), AZD5582, that activates the noncanonical NF-κB (ncNF-κB) signaling pathway to reverse latency. Ten weekly doses of AZD5582 were intravenously administered at 0.1 mg/kg to rhesus macaque (RM) infants orally infected with SIVmac251 at 4 weeks of age and treated with a triple ART regimen for over 1 year. During AZD5582 treatment, on-ART viremia above the limit of detection (LOD, 60 copies/mL) was observed in 5/8 infant RMs starting at 3 days post-dose 4 and peaking at 771 copies/mL. Of the 135 measurements during AZD5582 treatment in these 5 RM infants, only 8 were above the LOD (6%), lower than the 46% we have previously reported in adult RMs. Pharmacokinetic analysis of plasma AZD5582 levels revealed a lower Cmax in treated infants compared to adults (294 ng/mL versus 802 ng/mL). RNA-Sequencing of CD4+ T cells comparing pre- and post-AZD5582 dosing showed many genes that were similarly upregulated in infants and adults, but the expression of key ncNF-κB genes, including NFKB2 and RELB, was significantly higher in adult RMs. Our results suggest that dosing modifications for this latency reversal approach may be necessary to maximize virus reactivation in the pediatric setting for successful "shock and kill" strategies. IMPORTANCE While antiretroviral therapy (ART) has improved HIV-1 disease outcome and reduced transmission, interruption of ART results in rapid viral rebound due to the persistent latent reservoir. Interventions to reduce the viral reservoir are of critical importance, especially for children who must adhere to lifelong ART to prevent disease progression. Here, we used our previously established pediatric nonhuman primate model of oral SIV infection to evaluate AZD5582, identified as a potent latency-reversing agent in adult macaques, in the controlled setting of daily ART. We demonstrated the safety of the IAPi AZD5582 and evaluate the pharmacokinetics and pharmacodynamics of repeated dosing. The response to AZD5582 in macaque infants differed from what we previously showed in adult macaques with weaker latency reversal in infants, likely due to altered pharmacokinetics and less inducibility of infant CD4+ T cells. These data supported the contention that HIV-1 cure strategies for children are best evaluated using pediatric model systems.

Lymph node CXCR5+ NK cells associate with control of chronic SHIV infection.

The persistence of virally infected cells as reservoirs despite effective antiretroviral therapy is a major barrier to an HIV/SIV cure. These reservoirs are predominately contained within cells present in the B cell follicles (BCFs) of secondary lymphoid tissues, a site that is characteristically difficult for most cytolytic antiviral effector cells to penetrate. Here, we identified a population of NK cells in macaque lymph nodes that expressed BCF-homing receptor CXCR5 and accumulated within BCFs during chronic SHIV infection. These CXCR5+ follicular NK cells exhibited an activated phenotype coupled with heightened effector functions and a unique transcriptome characterized by elevated expression of cytolytic mediators (e.g., perforin and granzymes, LAMP-1). CXCR5+ NK cells exhibited high expression of FcγRIIa and FcγRIIIa, suggesting a potential for elevated antibody-dependent effector functionality. Consistently, accumulation of CXCR5+ NK cells showed a strong inverse association with plasma viral load and the frequency of germinal center follicular Th cells that comprise a significant fraction of the viral reservoir. Moreover, CXCR5+ NK cells showed increased expression of transcripts associated with IL-12 and IL-15 signaling compared with the CXCR5- subset. Indeed, in vitro treatment with IL-12 and IL-15 enhanced the proliferation of CXCR5+ granzyme B+ NK cells. Our findings suggest that follicular homing NK cells might be important in immune control of chronic SHIV infection, and this may have important implications for HIV cure strategies.

The Vacc-SeqQC project: Benchmarking RNA-Seq for clinical vaccine studies.

Introduction: Over the last decade, the field of systems vaccinology has emerged, in which high throughput transcriptomics and other omics assays are used to probe changes of the innate and adaptive immune system in response to vaccination. The goal of this study was to benchmark key technical and analytical parameters of RNA sequencing (RNA-seq) in the context of a multi-site, double-blind randomized vaccine clinical trial. Methods: We collected longitudinal peripheral blood mononuclear cell (PBMC) samples from 10 subjects before and after vaccination with a live attenuated Francisella tularensis vaccine and performed RNA-Seq at two different sites using aliquots from the same sample to generate two replicate datasets (5 time points for 50 samples each). We evaluated the impact of (i) filtering lowly-expressed genes, (ii) using external RNA controls, (iii) fold change and false discovery rate (FDR) filtering, (iv) read length, and (v) sequencing depth on differential expressed genes (DEGs) concordance between replicate datasets. Using synthetic mRNA spike-ins, we developed a method for empirically establishing minimal read-count thresholds for maintaining fold change accuracy on a per-experiment basis. We defined a reference PBMC transcriptome by pooling sequence data and established the impact of sequencing depth and gene filtering on transcriptome representation. Lastly, we modeled statistical power to detect DEGs for a range of sample sizes, effect sizes, and sequencing depths. Results and Discussion: Our results showed that (i) filtering lowly-expressed genes is recommended to improve fold-change accuracy and inter-site agreement, if possible guided by mRNA spike-ins (ii) read length did not have a major impact on DEG detection, (iii) applying fold-change cutoffs for DEG detection reduced inter-set agreement and should be used with caution, if at all, (iv) reduction in sequencing depth had a minimal impact on statistical power but reduced the identifiable fraction of the PBMC transcriptome, (v) after sample size, effect size (i.e. the magnitude of fold change) was the most important driver of statistical power to detect DEG. The results from this study provide RNA sequencing benchmarks and guidelines for planning future similar vaccine studies.

TREM2+ and interstitial macrophages orchestrate airway inflammation in SARS-CoV-2 infection in rhesus macaques.

The COVID-19 pandemic remains a global health crisis, yet, the immunopathological mechanisms driving the development of severe disease remain poorly defined. Here, we utilize a rhesus macaque (RM) model of SARS-CoV-2 infection to delineate perturbations in the innate immune system during acute infection using an integrated systems analysis. We found that SARS-CoV-2 initiated a rapid infiltration (two days post infection) of plasmacytoid dendritic cells into the lower airway, commensurate with IFNA production, natural killer cell activation, and induction of interferon-stimulated genes. At this early interval, we also observed a significant increase of blood CD14-CD16+ monocytes. To dissect the contribution of lung myeloid subsets to airway inflammation, we generated a novel compendium of RM-specific lung macrophage gene expression using a combination of sc-RNA-Seq data and bulk RNA-Seq of purified populations under steady state conditions. Using these tools, we generated a longitudinal sc-RNA-seq dataset of airway cells in SARS-CoV-2-infected RMs. We identified that SARS-CoV-2 infection elicited a rapid recruitment of two subsets of macrophages into the airway: a C206+MRC1-population resembling murine interstitial macrophages, and a TREM2+ population consistent with CCR2+ infiltrating monocytes, into the alveolar space. These subsets were the predominant source of inflammatory cytokines, accounting for ~75% of IL6 and TNF production, and >90% of IL10 production, whereas the contribution of CD206+MRC+ alveolar macrophages was significantly lower. Treatment of SARS-CoV-2 infected RMs with baricitinib (Olumiant ® ), a novel JAK1/2 inhibitor that recently received Emergency Use Authorization for the treatment of hospitalized COVID-19 patients, was remarkably effective in eliminating the influx of infiltrating, non-alveolar macrophages in the alveolar space, with a concomitant reduction of inflammatory cytokines. This study has delineated the major subsets of lung macrophages driving inflammatory and anti-inflammatory cytokine production within the alveolar space during SARS-CoV-2 infection. ONE SENTENCE SUMMARY: Multi-omic analyses of hyperacute SARS-CoV-2 infection in rhesus macaques identified two population of infiltrating macrophages, as the primary orchestrators of inflammation in the lower airway that can be successfully treated with baricitinib.