Epigenetic memory of coronavirus infection in innate immune cells and their progenitors.

Inflammation can trigger lasting phenotypes in immune and non-immune cells. Whether and how human infections and associated inflammation can form innate immune memory in hematopoietic stem and progenitor cells (HSPC) has remained unclear. We found that circulating HSPC, enriched from peripheral blood, captured the diversity of bone marrow HSPC, enabling investigation of their epigenomic reprogramming following coronavirus disease 2019 (COVID-19). Alterations in innate immune phenotypes and epigenetic programs of HSPC persisted for months to 1 year following severe COVID-19 and were associated with distinct transcription factor (TF) activities, altered regulation of inflammatory programs, and durable increases in myelopoiesis. HSPC epigenomic alterations were conveyed, through differentiation, to progeny innate immune cells. Early activity of IL-6 contributed to these persistent phenotypes in human COVID-19 and a mouse coronavirus infection model. Epigenetic reprogramming of HSPC may underlie altered immune function following infection and be broadly relevant, especially for millions of COVID-19 survivors.

Antiviral innate immune memory in alveolar macrophages following SARS-CoV-2 infection ameliorates secondary influenza A virus disease.

Pathogen encounter can result in epigenetic remodeling that shapes disease caused by heterologous pathogens. Here, we examined innate immune memory in the context of commonly circulating respiratory viruses. Single-cell analyses of airway-resident immune cells in a disease-relevant murine model of SARS-CoV-2 recovery revealed epigenetic reprogramming in alveolar macrophages following infection. Post-COVID-19 human monocytes exhibited similar epigenetic signatures. In airway-resident macrophages, past SARS-CoV-2 infection increased activity of type I interferon (IFN-I)-related transcription factors and epigenetic poising of antiviral genes. Viral pattern recognition and canonical IFN-I signaling were required for the establishment of this innate immune memory and augmented secondary antiviral responses. Antiviral innate immune memory mounted by airway-resident macrophages post-SARS-CoV-2 was necessary and sufficient to ameliorate secondary disease caused by influenza A virus and curtailed hyperinflammatory dysregulation and mortality. Our findings provide insights into antiviral innate immune memory in the airway that may facilitate the development of broadly effective therapeutic strategies.

HIV-1 control in vivo is related to the number but not the fraction of infected cells with viral unspliced RNA.

In the absence of antiretroviral therapy (ART), a subset of individuals, termed HIV controllers, have levels of plasma viremia that are orders of magnitude lower than non-controllers (NC) who are at higher risk for HIV disease progression. In addition to having fewer infected cells resulting in fewer cells with HIV RNA, it is possible that lower levels of plasma viremia in controllers are due to a lower fraction of the infected cells having HIV-1 unspliced RNA (HIV usRNA) compared with NC. To directly test this possibility, we used sensitive and quantitative single-cell sequencing methods to compare the fraction of infected cells that contain one or more copies of HIV usRNA in peripheral blood mononuclear cells (PBMC) obtained from controllers and NC. The fraction of infected cells containing HIV usRNA did not differ between the two groups. Rather, the levels of viremia were strongly associated with the total number of infected cells that had HIV usRNA, as reported by others, with controllers having 34-fold fewer infected cells per million PBMC. These results reveal that viremic control is not associated with a lower fraction of proviruses expressing HIV usRNA, unlike what is reported for elite controllers, but is only related to having fewer infected cells overall, maybe reflecting greater immune clearance of infected cells. Our findings show that proviral silencing is not a key mechanism for viremic control and will help to refine strategies toward achieving HIV remission without ART.

Integration in oncogenes plays only a minor role in determining the in vivo distribution of HIV integration sites before or during suppressive antiretroviral therapy.

HIV persists during antiretroviral therapy (ART) as integrated proviruses in cells descended from a small fraction of the CD4+ T cells infected prior to the initiation of ART. To better understand what controls HIV persistence and the distribution of integration sites (IS), we compared about 15,000 and 54,000 IS from individuals pre-ART and on ART, respectively, with approximately 395,000 IS from PBMC infected in vitro. The distribution of IS in vivo is quite similar to the distribution in PBMC, but modified by selection against proviruses in expressed genes, by selection for proviruses integrated into one of 7 specific genes, and by clonal expansion. Clones in which a provirus integrated in an oncogene contributed to cell survival comprised only a small fraction of the clones persisting in on ART. Mechanisms that do not involve the provirus, or its location in the host genome, are more important in determining which clones expand and persist.

Combined HIV-1 sequence and integration site analysis informs viral dynamics and allows reconstruction of replicating viral ancestors.

Understanding HIV-1 persistence despite antiretroviral therapy (ART) is of paramount importance. Both single-genome sequencing (SGS) and integration site analysis (ISA) provide useful information regarding the structure of persistent HIV DNA populations; however, until recently, there was no way to link integration sites to their cognate proviral sequences. Here, we used multiple-displacement amplification (MDA) of cellular DNA diluted to a proviral endpoint to obtain full-length proviral sequences and their corresponding sites of integration. We applied this method to lymph node and peripheral blood mononuclear cells from 5 ART-treated donors to determine whether groups of identical subgenomic sequences in the 2 compartments are the result of clonal expansion of infected cells or a viral genetic bottleneck. We found that identical proviral sequences can result from both cellular expansion and viral genetic bottlenecks occurring prior to ART initiation and following ART failure. We identified an expanded T cell clone carrying an intact provirus that matched a variant previously detected by viral outgrowth assays and expanded clones with wild-type and drug-resistant defective proviruses. We also found 2 clones from 1 donor that carried identical proviruses except for nonoverlapping deletions, from which we could infer the sequence of the intact parental virus. Thus, MDA-SGS can be used for "viral reconstruction" to better understand intrapatient HIV-1 evolution and to determine the clonality and structure of proviruses within expanded clones, including those with drug-resistant mutations. Importantly, we demonstrate that identical sequences observed by standard SGS are not always sufficient to establish proviral clonality.

HIVIntact: a python-based tool for HIV-1 genome intactness inference.

The characterisation of the HIV-1 reservoir, which consists of replication-competent integrated proviruses that persist on antiretroviral therapy (ART), is made difficult by the rarity of intact proviruses relative to those that are defective. While the only conclusive test for the replication-competence of HIV-1 proviruses is carried out in cell culture, genetic characterization of genomes by near full-length (NFL) PCR and sequencing can be used to determine whether particular proviruses have insertions, deletions, or substitutions that render them defective. Proviruses that are not excluded by having such defects can be classified as genetically intact and, possibly, replication competent. Identifying and quantifying proviruses that are potentially replication-competent is important for the development of strategies towards a functional cure. However, to date, there are no programs that can be incorporated into deep-sequencing pipelines for the automated characterization and annotation of HIV genomes. Existing programs that perform this work require manual intervention, cannot be widely installed, and do not have easily adjustable settings. Here, we present HIVIntact, a python-based software tool that characterises genomic defects in NFL HIV-1 sequences, allowing putative intact genomes to be identified in-silico. Unlike other applications that assess the genetic intactness of HIV genomes, this tool can be incorporated into existing sequence-analysis pipelines and applied to large next-generation sequencing datasets.

Intact HIV Proviruses Persist in Children Seven to Nine Years after Initiation of Antiretroviral Therapy in the First Year of Life.

In adults starting antiretroviral therapy (ART) during acute infection, 2% of proviruses that persist on ART are genetically intact by sequence analysis. In contrast, a recent report in children treated early failed to detect sequence-intact proviruses. In another cohort of children treated early, we sought to detect and characterize proviral sequences after 6 to 9 years on suppressive ART. Peripheral blood mononuclear cells (PBMC) from perinatally infected children from the Children with HIV Early antiRetroviral (CHER) study were analyzed. Nearly full-length proviral amplification and sequencing (NFL-PAS) were performed at one time point after 6 to 9 years on ART. Amplicons with large internal deletions were excluded (<9 kb). All amplicons of ≥9 kb were sequenced and analyzed through a bioinformatic pipeline to detect indels, frameshifts, or hypermutations that would render them defective. In eight children who started ART at a median age of 5.4 months (range, 2.0 to 11.1 months), 733 single NFL-PAS amplicons were generated. Of these, 534 (72.9%) had large internal deletions, 174 (23.7%) had hypermutations, 15 (1.4%) had small internal deletions, 3 (1.0%) had deletions in the packaging signal/major splice donor site, and 7 (1.0%) were sequence intact. These 7 intact sequences were from three children who initiated ART after 2.3 months of age, one of whom had two identical intact sequences, suggestive of a cell clone harboring a replication-competent provirus. No intact proviruses were detected in four children who initiated ART before 2.3 months of age. Rare, intact proviruses can be detected in children who initiate ART after 2.3 months of age and are probably, as in adults, maintained by clonal expansion of cells infected before ART initiation.IMPORTANCE There are limited data about the proviral landscape in children exhibiting long-term suppression after early treatment, particularly in Sub-Saharan Africa where HIV-1 subtype C predominates. Investigating the sequence-intact reservoir could provide insight on the mechanisms by which intact proviruses persist and inform ongoing cure efforts. Through nearly full-length proviral amplification and sequencing (NFL-PAS), we generated 733 NFL-PAS amplicons from eight children. We showed that rare, genetically intact proviruses could be detected in children who initiated ART after 2.3 months of age. The frequency of intact proviruses was lower (P < 0.05) than that reported for HIV subtype B-infected adults treated during early HIV infection. We show that cells harboring genetically intact HIV proviruses are rare in children exhibiting long-term suppression after early treatment and may require the processing of a large number of cells to assess reservoir size. This points to the need for efficient methods to accurately quantify latent reservoirs, particularly in pediatric studies where sample availability is limited.

Effector memory differentiation increases detection of replication-competent HIV-l in resting CD4+ T cells from virally suppressed individuals.

Studies have demonstrated that intensive ART alone is not capable of eradicating HIV-1, as the virus rebounds within a few weeks upon treatment interruption. Viral rebound may be induced from several cellular subsets; however, the majority of proviral DNA has been found in antigen experienced resting CD4+ T cells. To achieve a cure for HIV-1, eradication strategies depend upon both understanding mechanisms that drive HIV-1 persistence as well as sensitive assays to measure the frequency of infected cells after therapeutic interventions. Assays such as the quantitative viral outgrowth assay (QVOA) measure HIV-1 persistence during ART by ex vivo activation of resting CD4+ T cells to induce latency reversal; however, recent studies have shown that only a fraction of replication-competent viruses are inducible by primary mitogen stimulation. Previous studies have shown a correlation between the acquisition of effector memory phenotype and HIV-1 latency reversal in quiescent CD4+ T cell subsets that harbor the reservoir. Here, we apply our mechanistic understanding that differentiation into effector memory CD4+ T cells more effectively promotes HIV-1 latency reversal to significantly improve proviral measurements in the QVOA, termed differentiation QVOA (dQVOA), which reveals a significantly higher frequency of the inducible HIV-1 replication-competent reservoir in resting CD4+ T cells.

Correction to: An analytical pipeline for identifying and mapping the integration sites of HIV and other retroviruses.

An amendment to this paper has been published and can be accessed via the original article.

An analytical pipeline for identifying and mapping the integration sites of HIV and other retroviruses.

BACKGROUND: All retroviruses, including human immunodeficiency virus (HIV), must integrate a DNA copy of their genomes into the genome of the infected host cell to replicate. Although integrated retroviral DNA, known as a provirus, can be found at many sites in the host genome, integration is not random. The adaption of linker-mediated PCR (LM-PCR) protocols for high-throughput integration site mapping, using randomly-sheared genomic DNA and Illumina paired-end sequencing, has dramatically increased the number of mapped integration sites. Analysis of samples from human donors has shown that there is clonal expansion of HIV infected cells and that clonal expansion makes an important contribution to HIV persistence. However, analysis of HIV integration sites in samples taken from patients requires extensive PCR amplification and high-throughput sequencing, which makes the methodology prone to certain specific artifacts. RESULTS: To address the problems with artifacts, we use a comprehensive approach involving experimental procedures linked to a bioinformatics analysis pipeline. Using this combined approach, we are able to reduce the number of PCR/sequencing artifacts that arise and identify the ones that remain. Our streamlined workflow combines random cleavage of the DNA in the samples, end repair, and linker ligation in a single step. We provide guidance on primer and linker design that reduces some of the common artifacts. We also discuss how to identify and remove some of the common artifacts, including the products of PCR mispriming and PCR recombination, that have appeared in some published studies. Our improved bioinformatics pipeline rapidly parses the sequencing data and identifies bona fide integration sites in clonally expanded cells, producing an Excel-formatted report that can be used for additional data processing. CONCLUSIONS: We provide a detailed protocol that reduces the prevalence of artifacts that arise in the analysis of retroviral integration site data generated from in vivo samples and a bioinformatics pipeline that is able to remove the artifacts that remain.

HIV evolution and diversity in ART-treated patients.

Characterizing HIV genetic diversity and evolution during antiretroviral therapy (ART) provides insights into the mechanisms that maintain the viral reservoir during ART. This review describes common methods used to obtain and analyze intra-patient HIV sequence data, the accumulation of diversity prior to ART and how it is affected by suppressive ART, the debate on viral replication and evolution in the presence of ART, HIV compartmentalization across various tissues, and mechanisms for the emergence of drug resistance. It also describes how CD4+ T cells that were likely infected with latent proviruses prior to initiating treatment can proliferate before and during ART, providing a renewable source of infected cells despite therapy. Some expanded cell clones carry intact and replication-competent proviruses with a small fraction of the clonal siblings being transcriptionally active and a source for residual viremia on ART. Such cells may also be the source for viral rebound after interrupting ART. The identical viral sequences observed for many years in both the plasma and infected cells of patients on long-term ART are likely due to the proliferation of infected cells both prior to and during treatment. Studies on HIV diversity may reveal targets that can be exploited in efforts to eradicate or control the infection without ART.

Mechanisms underlying a thalamocortical transformation during active tactile sensation.

During active somatosensation, neural signals expected from movement of the sensors are suppressed in the cortex, whereas information related to touch is enhanced. This tactile suppression underlies low-noise encoding of relevant tactile features and the brain's ability to make fine tactile discriminations. Layer (L) 4 excitatory neurons in the barrel cortex, the major target of the somatosensory thalamus (VPM), respond to touch, but have low spike rates and low sensitivity to the movement of whiskers. Most neurons in VPM respond to touch and also show an increase in spike rate with whisker movement. Therefore, signals related to self-movement are suppressed in L4. Fast-spiking (FS) interneurons in L4 show similar dynamics to VPM neurons. Stimulation of halorhodopsin in FS interneurons causes a reduction in FS neuron activity and an increase in L4 excitatory neuron activity. This decrease of activity of L4 FS neurons contradicts the "paradoxical effect" predicted in networks stabilized by inhibition and in strongly-coupled networks. To explain these observations, we constructed a model of the L4 circuit, with connectivity constrained by in vitro measurements. The model explores the various synaptic conductance strengths for which L4 FS neurons actively suppress baseline and movement-related activity in layer 4 excitatory neurons. Feedforward inhibition, in concert with recurrent intracortical circuitry, produces tactile suppression. Synaptic delays in feedforward inhibition allow transmission of temporally brief volleys of activity associated with touch. Our model provides a mechanistic explanation of a behavior-related computation implemented by the thalamocortical circuit.

Microsecond-scale timing precision in rodent trigeminal primary afferents.

Communication in the nervous system occurs by spikes: the timing precision with which spikes are fired is a fundamental limit on neural information processing. In sensory systems, spike-timing precision is constrained by first-order neurons. We found that spike-timing precision of trigeminal primary afferents in rats and mice is limited both by stimulus speed and by electrophysiological sampling rate. High-speed video of behaving mice revealed whisker velocities of at least 17,000°/s, so we delivered an ultrafast "ping" (>50,000°/s) to single whiskers and sampled primary afferent activity at 500 kHz. Median spike jitter was 17.4 µs; 29% of neurons had spike jitter < 10 µs. These results indicate that the input stage of the trigeminal pathway has extraordinary spike-timing precision and very high potential information capacity. This timing precision ranks among the highest in biology.

Low-dimensional sensory feature representation by trigeminal primary afferents.

In any sensory system, the primary afferents constitute the first level of sensory representation and fundamentally constrain all subsequent information processing. Here, we show that the spike timing, reliability, and stimulus selectivity of primary afferents in the whisker system can be accurately described by a simple model consisting of linear stimulus filtering combined with spike feedback. We fitted the parameters of the model by recording the responses of primary afferents to filtered, white noise whisker motion in anesthetized rats. The model accurately predicted not only the response of primary afferents to white noise whisker motion (median correlation coefficient 0.92) but also to naturalistic, texture-induced whisker motion. The model accounted both for submillisecond spike-timing precision and for non-Poisson spike train structure. We found substantial diversity in the responses of the afferent population, but this diversity was accurately captured by the model: a 2D filter subspace, corresponding to different mixtures of position and velocity sensitivity, captured 94% of the variance in the stimulus selectivity. Our results suggest that the first stage of the whisker system can be well approximated as a bank of linear filters, forming an overcomplete representation of a low-dimensional feature space.

Mechanisms of Epigenomic and Functional Convergence Between Glucocorticoid- and IL4-Driven Macrophage Programming.

Macrophages adopt distinct phenotypes in response to environmental cues, with type-2 cytokine interleukin-4 promoting a tissue-repair homeostatic state (M2IL4). Glucocorticoids, widely used anti-inflammatory therapeutics, reportedly impart a similar phenotype (M2GC), but how such disparate pathways may functionally converge is unknown. We show using integrative functional genomics that M2IL4 and M2GC transcriptomes share a striking overlap mirrored by a shift in chromatin landscape in both common and signal-specific gene subsets. This core homeostatic program is enacted by transcriptional effectors KLF4 and the GC receptor, whose genome-wide occupancy and actions are integrated in a stimulus-specific manner by the nuclear receptor cofactor GRIP1. Indeed, many of the M2IL4:M2GC-shared transcriptomic changes were GRIP1-dependent. Consistently, GRIP1 loss attenuated phagocytic activity of both populations in vitro and macrophage tissue-repair properties in the murine colitis model in vivo. These findings provide a mechanistic framework for homeostatic macrophage programming by distinct signals, to better inform anti-inflammatory drug design.

Microbial cancer immunotherapy reprograms hematopoietic stem cells to enhance anti-tumor immunity.

Mycobacterium bovis BCG is the vaccine against tuberculosis and an immunotherapy for bladder cancer. When administered intravenously, BCG reprograms bone marrow hematopoietic stem and progenitor cells (HSPCs), leading to heterologous protection against infections. Whether HSPC-reprogramming contributes to the anti-tumor effects of BCG administered into the bladder is unknown. We demonstrate that BCG administered in the bladder in both mice and humans reprograms HSPCs to amplify myelopoiesis and functionally enhance myeloid cell antigen presentation pathways. Reconstitution of naive mice with HSPCs from bladder BCG-treated mice enhances anti-tumor immunity and tumor control, increases intratumor dendritic cell infiltration, reprograms pro-tumorigenic neutrophils, and synergizes with checkpoint blockade. We conclude that bladder BCG acts systemically, reprogramming HSPC-encoded innate immunity, highlighting the broad potential of modulating HSPC phenotypes to improve tumor immunity.

Myeloid progenitor dysregulation fuels immunosuppressive macrophages in tumors.

Monocyte-derived macrophages (mo-macs) drive immunosuppression in the tumor microenvironment (TME) and tumor-enhanced myelopoiesis in the bone marrow (BM) fuels these populations. Here, we performed paired transcriptome and chromatin analysis over the continuum of BM myeloid progenitors, circulating monocytes, and tumor-infiltrating mo-macs in mice and in patients with lung cancer to identify myeloid progenitor programs that fuel pro-tumorigenic mo-macs. Analyzing chromatin accessibility and histone mark changes, we show that lung tumors prime accessibility for Nfe2l2 (NRF2) in BM myeloid progenitors as a cytoprotective response to oxidative stress. NRF2 activity is sustained and increased during monocyte differentiation into mo-macs in the lung TME to regulate oxidative stress, in turn promoting metabolic adaptation, resistance to cell death, and contributing to immunosuppressive phenotype. NRF2 genetic deletion and pharmacological inhibition significantly reduced mo-macs' survival and immunosuppression in the TME, enabling NK and T cell therapeutic antitumor immunity and synergizing with checkpoint blockade strategies. Altogether, our study identifies a targetable epigenetic node of myeloid progenitor dysregulation that sustains immunoregulatory mo-macs in the TME.

The largest HIV-1-infected T cell clones in children on long-term combination antiretroviral therapy contain solo LTRs.

Combination antiretroviral therapy (cART) suppresses viral replication but does not cure HIV infection because a reservoir of infectious (intact) HIV proviruses persists in long-lived CD4+T cells. However, a large majority (>95%) of HIV-infected cells that persist on effective cART carry defective (non-infectious) proviruses. Defective proviruses consisting of only a single LTR (solo long terminal repeat) are commonly found as endogenous retroviruses in many animal species, but the frequency of solo-LTR HIV proviruses has not been well defined. Here we show that, in five pediatric donors whose viremia was suppressed on cART for at least 5 years, the proviruses in the nine largest clones of HIV-infected cells were solo LTRs. The sizes of five of these clones were assayed longitudinally by integration site-specific quantitative PCR. Minor waxing and waning of the clones was observed, suggesting that these clones are generally stable over time. Our findings show that solo LTRs comprise a large fraction of the proviruses in infected cell clones that persist in children on long-term cART. IMPORTANCE This work highlights that severely deleted HIV-1 proviruses comprise a significant proportion of the proviral landscape and are often overlooked.

Application of ultrasensitive digital ELISA for p24 enables improved evaluation of HIV-1 reservoir diversity and growth kinetics in viral outgrowth assays.

The advent of combined antiretroviral therapy (cART) has been instrumental in controlling HIV-1 replication and transmission and decreasing associated morbidity and mortality. However, cART alone is not able to cure HIV-1 due to the presence of long-lived, latently infected immune cells, which re-seed plasma viremia when cART is interrupted. Assessment of HIV-cure strategies using ex vivo culture methods for further understanding of the diversity of reactivated HIV, viral outgrowth, and replication dynamics are enhanced using ultrasensitive digital ELISA based on single-molecule array (Simoa) technology to increase the sensitivity of endpoint detection. In viral outgrowth assays (VOA), exponential HIV-1 outgrowth has been shown to be dependent upon initial virus burst size surpassing a critical growth threshold of 5100 HIV-1 RNA copies. Here, we show an association between ultrasensitive HIV-1 Gag p24 concentrations and HIV-1 RNA copy number that characterize viral dynamics below the exponential replication threshold. Single-genome sequencing (SGS) revealed the presence of multiple identical HIV-1 sequences, indicative of low-level replication occurring below the threshold of exponential outgrowth early during a VOA. However, SGS further revealed diverse related HIV variants detectable by ultrasensitive methods that failed to establish exponential outgrowth. Overall, our data suggest that viral outgrowth occurring below the threshold necessary for establishing exponential growth in culture does not preclude replication competence of reactivated HIV, and ultrasensitive detection of HIV-1 p24 may provide a method to detect previously unquantifiable variants. These data strongly support the use of the Simoa platform in a multi-prong approach to measuring latent viral burden and efficacy of therapeutic interventions aimed at an HIV-1 cure.

Comparison of latency and rate coding for the direction of whisker deflection in the subcortical somatosensory pathway.

The response of many neurons in the whisker somatosensory system depends on the direction in which a whisker is deflected. Although it is known that the spike count conveys information about this parameter, it is not known how important spike timing might be. The aim of this study was to compare neural codes based on spike count and first-spike latency, respectively. We extracellularly recorded single units from either the rat trigeminal ganglion (primary sensory afferents) or ventroposteromedial (VPM) thalamic nucleus in response to deflection in different directions and quantified alternative neural codes using mutual information. We found that neurons were diverse: some (58% in ganglion, 32% in VPM) conveyed information only by spike count; others conveyed additional information by latency. An issue with latency coding is that latency is measured with respect to the time of stimulus onset, a quantity known to the experimenter but not directly to the subject's brain. We found a potential solution using the integrated population activity as an internal timing signal: in this way, 91% of the first-spike latency information could be recovered. Finally, we asked how well direction could be decoded. For large populations, spike count and latency codes performed similarly; for small ones, decoding was more accurate using the latency code. Our findings indicate that whisker deflection direction is more efficiently encoded by spike timing than by spike count. Spike timing decreases the population size necessary for reliable information transmission and may thereby bring significant advantages in both wiring and metabolic efficiency.