JAK-STAT signaling maintains homeostasis in T cells and macrophages.

Immune cells need to sustain a state of constant alertness over a lifetime. Yet, little is known about the regulatory processes that control the fluent and fragile balance that is called homeostasis. Here we demonstrate that JAK-STAT signaling, beyond its role in immune responses, is a major regulator of immune cell homeostasis. We investigated JAK-STAT-mediated transcription and chromatin accessibility across 12 mouse models, including knockouts of all STAT transcription factors and of the TYK2 kinase. Baseline JAK-STAT signaling was detected in CD8+ T cells and macrophages of unperturbed mice-but abrogated in the knockouts and in unstimulated immune cells deprived of their normal tissue context. We observed diverse gene-regulatory programs, including effects of STAT2 and IRF9 that were independent of STAT1. In summary, our large-scale dataset and integrative analysis of JAK-STAT mutant and wild-type mice uncovered a crucial role of JAK-STAT signaling in unstimulated immune cells, where it contributes to a poised epigenetic and transcriptional state and helps prepare these cells for rapid response to immune stimuli.

Artemisinins Target GABAA Receptor Signaling and Impair α Cell Identity.

Type 1 diabetes is characterized by the destruction of pancreatic ß cells, and generating new insulin-producing cells from other cell types is a major aim of regenerative medicine. One promising approach is transdifferentiation of developmentally related pancreatic cell types, including glucagon-producing α cells. In a genetic model, loss of the master regulatory transcription factor Arx is sufficient to induce the conversion of α cells to functional ß-like cells. Here, we identify artemisinins as small molecules that functionally repress Arx by causing its translocation to the cytoplasm. We show that the protein gephyrin is the mammalian target of these antimalarial drugs and that the mechanism of action of these molecules depends on the enhancement of GABAA receptor signaling. Our results in zebrafish, rodents, and primary human pancreatic islets identify gephyrin as a druggable target for the regeneration of pancreatic ß cell mass from α cells.

Temporal dissection of an enhancer cluster reveals distinct temporal and functional contributions of individual elements.

Many genes are regulated by multiple enhancers that often simultaneously activate their target gene. However, how individual enhancers collaborate to activate transcription is not well understood. Here, we dissect the functions and interdependencies of five enhancer elements that together activate Fgf5 expression during exit from naive murine pluripotency. Four intergenic elements form a super-enhancer, and most of the elements contribute to Fgf5 induction at distinct time points. A fifth, poised enhancer located in the first intron contributes to Fgf5 expression at every time point by amplifying overall Fgf5 expression levels. Despite low individual enhancer activity, together these elements strongly induce Fgf5 expression in a super-additive fashion that involves strong accumulation of RNA polymerase II at the intronic enhancer. Finally, we observe a strong anti-correlation between RNA polymerase II levels at enhancers and their distance to the closest promoter, and we identify candidate elements with properties similar to the intronic enhancer.

RANK links thymic regulatory T cells to fetal loss and gestational diabetes in pregnancy.

Successful pregnancies rely on adaptations within the mother1, including marked changes within the immune system2. It has long been known that the thymus, the central lymphoid organ, changes markedly during pregnancy3. However, the molecular basis and importance of this process remain largely obscure. Here we show that the osteoclast differentiation receptor RANK4,5 couples female sex hormones to the rewiring of the thymus during pregnancy. Genetic deletion of Rank (also known as Tnfrsf11a) in thymic epithelial cells results in impaired thymic involution and blunted expansion of natural regulatory T (Treg) cells in pregnant female mice. Sex hormones, in particular progesterone, drive the development of thymic Treg cells through RANK in a manner that depends on AIRE+ medullary thymic epithelial cells. The depletion of Rank in the mouse thymic epithelium results in reduced accumulation of natural Treg cells in the placenta, and an increase in the number of miscarriages. Thymic deletion of Rank also results in impaired accumulation of Treg cells in visceral adipose tissue, and is associated with enlarged adipocyte size, tissue inflammation, enhanced maternal glucose intolerance, fetal macrosomia, and a long-lasting transgenerational alteration in glucose homeostasis, which are all key hallmarks of gestational diabetes. Transplantation of Treg cells rescued fetal loss, maternal glucose intolerance and fetal macrosomia. In human pregnancies, we found that gestational diabetes also correlates with a reduced number of Treg cells in the placenta. Our findings show that RANK promotes the hormone-mediated development of thymic Treg cells during pregnancy, and expand the functional role of maternal Treg cells to the development of gestational diabetes and the transgenerational metabolic rewiring of glucose homeostasis.

Mitochondria Are a Subset of Extracellular Vesicles Released by Activated Monocytes and Induce Type I IFN and TNF Responses in Endothelial Cells.

RATIONALE: Extracellular vesicles, including microvesicles, are increasingly recognized as important mediators in cardiovascular disease. The cargo and surface proteins they carry are considered to define their biological activity, including their inflammatory properties. Monocyte to endothelial cell signaling is a prerequisite for the propagation of inflammatory responses. However, the contribution of microvesicles in this process is poorly understood. OBJECTIVE: To elucidate the mechanisms by which microvesicles derived from activated monocytic cells exert inflammatory effects on endothelial cells. METHODS AND RESULTS: LPS (lipopolysaccharide)-stimulated monocytic cells release free mitochondria and microvesicles with mitochondrial content as demonstrated by flow cytometry, quantitative polymerase chain reaction, Western Blot, and transmission electron microscopy. Using RNAseq analysis and quantitative reverse transcription-polymerase chain reaction, we demonstrated that both mitochondria directly isolated from and microvesicles released by LPS-activated monocytic cells, as well as circulating microvesicles isolated from volunteers receiving low-dose LPS-injections, induce type I IFN (interferon), and TNF (tumor necrosis factor) responses in endothelial cells. Depletion of free mitochondria significantly reduced the ability of these microvesicles to induce type I IFN and TNF-dependent genes. We identified mitochondria-associated TNFα and RNA from stressed mitochondria as major inducers of these responses. Finally, we demonstrated that the proinflammatory potential of microvesicles and directly isolated mitochondria were drastically reduced when they were derived from monocytic cells with nonrespiring mitochondria or monocytic cells cultured in the presence of pyruvate or the mitochondrial reactive oxygen species scavenger MitoTEMPO. CONCLUSIONS: Mitochondria and mitochondria embedded in microvesicles constitute a major subset of extracellular vesicles released by activated monocytes, and their proinflammatory activity on endothelial cells is determined by the activation status of their parental cells. Thus, mitochondria may represent critical intercellular mediators in cardiovascular disease and other inflammatory settings associated with type I IFN and TNF signaling.

A combinatorial screen of the CLOUD uncovers a synergy targeting the androgen receptor.

Approved drugs are invaluable tools to study biochemical pathways, and further characterization of these compounds may lead to repurposing of single drugs or combinations. Here we describe a collection of 308 small molecules representing the diversity of structures and molecular targets of all FDA-approved chemical entities. The CeMM Library of Unique Drugs (CLOUD) covers prodrugs and active forms at pharmacologically relevant concentrations and is ideally suited for combinatorial studies. We screened pairwise combinations of CLOUD drugs for impairment of cancer cell viability and discovered a synergistic interaction between flutamide and phenprocoumon (PPC). The combination of these drugs modulates the stability of the androgen receptor (AR) and resensitizes AR-mutant prostate cancer cells to flutamide. Mechanistically, we show that the AR is a substrate for γ-carboxylation, a post-translational modification inhibited by PPC. Collectively, our data suggest that PPC could be repurposed to tackle resistance to antiandrogens in prostate cancer patients.

A T cell-specific deletion of HDAC1 protects against experimental autoimmune encephalomyelitis.

Multiple sclerosis (MS) is a human neurodegenerative disease characterized by the invasion of autoreactive T cells from the periphery into the CNS. Application of pan-histone deacetylase inhibitors (HDACi) ameliorates experimental autoimmune encephalomyelitis (EAE), an animal model for MS, suggesting that HDACi might be a potential therapeutic strategy for MS. However, the function of individual HDAC members in the pathogenesis of EAE is not known. In this study we report that mice with a T cell-specific deletion of HDAC1 (using the Cd4-Cre deleter strain; HDAC1-cKO) were completely resistant to EAE despite the ability of HDAC1cKO CD4+ T cells to differentiate into Th17 cells. RNA sequencing revealed STAT1 as a prominent upstream regulator of differentially expressed genes in activated HDAC1-cKO CD4+ T cells and this was accompanied by a strong increase in phosphorylated STAT1 (pSTAT1). This suggests that HDAC1 controls STAT1 activity in activated CD4+ T cells. Increased pSTAT1 levels correlated with a reduced expression of the chemokine receptors Ccr4 and Ccr6, which are important for the migration of T cells into the CNS. Finally, EAE susceptibility was restored in WT:HDAC1-cKO mixed BM chimeric mice, indicating a cell-autonomous defect. Our data demonstrate a novel pathophysiological role for HDAC1 in EAE and provide evidence that selective inhibition of HDAC1 might be a promising strategy for the treatment of MS.

Mapping the chemical chromatin reactivation landscape identifies BRD4-TAF1 cross-talk.

Bromodomain-containing proteins of the BET family recognize histone lysine acetylation and mediate transcriptional activation of target genes such as the MYC oncogene. Pharmacological inhibitors of BET domains promise therapeutic benefits in a variety of cancers. We performed a high-diversity chemical compound screen for agents capable of modulating BRD4-dependent heterochromatization of a generic reporter in human cells. In addition to known and new compounds targeting BRD4, we identified small molecules that mimic BRD4 inhibition without direct engagement. One such compound was a potent inhibitor of the second bromodomain of TAF1. Using this inhibitor, we discovered that TAF1 synergizes with BRD4 to control proliferation of cancer cells, making TAF1 an attractive epigenetic target in cancers driven by MYC.

Single-cell transcriptomes reveal characteristic features of human pancreatic islet cell types.

Pancreatic islets of Langerhans contain several specialized endocrine cell types, which are commonly identified by the expression of single marker genes. However, the established marker genes cannot capture the complete spectrum of cellular heterogeneity in human pancreatic islets, and existing bulk transcriptome datasets provide averages across several cell populations. To dissect the cellular composition of the human pancreatic islet and to establish transcriptomes for all major cell types, we performed single-cell RNA sequencing on 70 cells sorted from human primary tissue. We used this dataset to validate previously described marker genes at the single-cell level and to identify specifically expressed transcription factors for all islet cell subtypes. All data are available for browsing and download, thus establishing a useful resource of single-cell expression profiles for endocrine cells in human pancreatic islets.

Parallel reverse genetic screening in mutant human cells using transcriptomics.

Reverse genetic screens have driven gene annotation and target discovery in model organisms. However, many disease-relevant genotypes and phenotypes cannot be studied in lower organisms. It is therefore essential to overcome technical hurdles associated with large-scale reverse genetics in human cells. Here, we establish a reverse genetic approach based on highly robust and sensitive multiplexed RNA sequencing of mutant human cells. We conduct 10 parallel screens using a collection of engineered haploid isogenic cell lines with knockouts covering tyrosine kinases and identify known and unexpected effects on signaling pathways. Our study provides proof of concept for a scalable approach to link genotype to phenotype in human cells, which has broad applications. In particular, it clears the way for systematic phenotyping of still poorly characterized human genes and for systematic study of uncharacterized genomic features associated with human disease.

Comparative genomics suggests an independent origin of cytoplasmic incompatibility in Cardinium hertigii.

Terrestrial arthropods are commonly infected with maternally inherited bacterial symbionts that cause cytoplasmic incompatibility (CI). In CI, the outcome of crosses between symbiont-infected males and uninfected females is reproductive failure, increasing the relative fitness of infected females and leading to spread of the symbiont in the host population. CI symbionts have profound impacts on host genetic structure and ecology and may lead to speciation and the rapid evolution of sex determination systems. Cardinium hertigii, a member of the Bacteroidetes and symbiont of the parasitic wasp Encarsia pergandiella, is the only known bacterium other than the Alphaproteobacteria Wolbachia to cause CI. Here we report the genome sequence of Cardinium hertigii cEper1. Comparison with the genomes of CI-inducing Wolbachia pipientis strains wMel, wRi, and wPip provides a unique opportunity to pinpoint shared proteins mediating host cell interaction, including some candidate proteins for CI that have not previously been investigated. The genome of Cardinium lacks all major biosynthetic pathways but harbors a complete biotin biosynthesis pathway, suggesting a potential role for Cardinium in host nutrition. Cardinium lacks known protein secretion systems but encodes a putative phage-derived secretion system distantly related to the antifeeding prophage of the entomopathogen Serratia entomophila. Lastly, while Cardinium and Wolbachia genomes show only a functional overlap of proteins, they show no evidence of laterally transferred elements that would suggest common ancestry of CI in both lineages. Instead, comparative genomics suggests an independent evolution of CI in Cardinium and Wolbachia and provides a novel context for understanding the mechanistic basis of CI.

The endosymbiont Amoebophilus asiaticus encodes an S-adenosylmethionine carrier that compensates for its missing methylation cycle.

All organisms require S-adenosylmethionine (SAM) as a methyl group donor and cofactor for various biologically important processes. However, certain obligate intracellular parasitic bacteria and also the amoeba symbiont Amoebophilus asiaticus have lost the capacity to synthesize this cofactor and hence rely on its uptake from host cells. Genome analyses revealed that A. asiaticus encodes a putative SAM transporter. The corresponding protein was functionally characterized in Escherichia coli: import studies demonstrated that it is specific for SAM and S-adenosylhomocysteine (SAH), the end product of methylation. SAM transport activity was shown to be highly dependent on the presence of a membrane potential, and by targeted analyses, we obtained direct evidence for a proton-driven SAM/SAH antiport mechanism. Sequence analyses suggest that SAM carriers from Rickettsiales might operate in a similar way, in contrast to chlamydial SAM transporters. SAM/SAH antiport is of high physiological importance, as it allows for compensation for the missing methylation cycle. The identification of a SAM transporter in A. asiaticus belonging to the Bacteroidetes phylum demonstrates that SAM transport is more widely spread than previously assumed and occurs in bacteria belonging to three different phyla (Proteobacteria, Chlamydiae, and Bacteroidetes).

Unity in variety--the pan-genome of the Chlamydiae.

Chlamydiae are evolutionarily well-separated bacteria that live exclusively within eukaryotic host cells. They include important human pathogens such as Chlamydia trachomatis as well as symbionts of protozoa. As these bacteria are experimentally challenging and genetically intractable, our knowledge about them is still limited. In this study, we obtained the genome sequences of Simkania negevensis Z, Waddlia chondrophila 2032/99, and Parachlamydia acanthamoebae UV-7. This enabled us to perform the first comprehensive comparative and phylogenomic analysis of representative members of four major families of the Chlamydiae, including the Chlamydiaceae. We identified a surprisingly large core gene set present in all genomes and a high number of diverse accessory genes in those Chlamydiae that do not primarily infect humans or animals, including a chemosensory system in P. acanthamoebae and a type IV secretion system. In S. negevensis, the type IV secretion system is encoded on a large conjugative plasmid (pSn, 132 kb). Phylogenetic analyses suggested that a plasmid similar to the S. negevensis plasmid was originally acquired by the last common ancestor of all four families and that it was subsequently reduced, integrated into the chromosome, or lost during diversification, ultimately giving rise to the extant virulence-associated plasmid of pathogenic chlamydiae. Other virulence factors, including a type III secretion system, are conserved among the Chlamydiae to variable degrees and together with differences in the composition of the cell wall reflect adaptation to different host cells including convergent evolution among the four chlamydial families. Phylogenomic analysis focusing on chlamydial proteins with homology to plant proteins provided evidence for the acquisition of 53 chlamydial genes by a plant progenitor, lending further support for the hypothesis of an early interaction between a chlamydial ancestor and the primary photosynthetic eukaryote.

Emergence of SARS-CoV-2 Alpha lineage and its correlation with quantitative wastewater-based epidemiology data.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) gave rise to an international public health emergency in 3 months after its emergence in Wuhan, China. Typically for an RNA virus, random mutations occur constantly leading to new lineages, incidental with a higher transmissibility. The highly infective alpha lineage, firstly discovered in the UK, led to elevated mortality and morbidity rates as a consequence of Covid-19, worldwide. Wastewater surveillance proved to be a powerful tool for early detection and subsequent monitoring of the dynamics of SARS-CoV-2 and its variants in a defined catchment. Using a combination of sequencing and RT-qPCR approaches, we investigated the total SARS-CoV-2 concentration and the emergence of the alpha lineage in wastewater samples in Vienna, Austria linking it to clinical data. Based on a non-linear regression model and occurrence of signature mutations, we conclude that the alpha variant was present in Vienna sewage samples already in December 2020, even one month before the first clinical case was officially confirmed and reported by the health authorities. This provides evidence that a well-designed wastewater monitoring approach can provide a fast snapshot and may detect the circulating lineages in wastewater weeks before they are detectable in the clinical samples. Furthermore, declining 14 days prevalence data with simultaneously increasing SARS-CoV-2 total concentration in wastewater indicate a different shedding behavior for the alpha variant. Overall, our results support wastewater surveillance to be a suitable approach to spot early circulating SARS-CoV-2 lineages based on whole genome sequencing and signature mutations analysis.

Viral variant-resolved wastewater surveillance of SARS-CoV-2 at national scale.

SARS-CoV-2 surveillance by wastewater-based epidemiology is poised to provide a complementary approach to sequencing individual cases. However, robust quantification of variants and de novo detection of emerging variants remains challenging for existing strategies. We deep sequenced 3,413 wastewater samples representing 94 municipal catchments, covering >59% of the population of Austria, from December 2020 to February 2022. Our system of variant quantification in sewage pipeline designed for robustness (termed VaQuERo) enabled us to deduce the spatiotemporal abundance of predefined variants from complex wastewater samples. These results were validated against epidemiological records of >311,000 individual cases. Furthermore, we describe elevated viral genetic diversity during the Delta variant period, provide a framework to predict emerging variants and measure the reproductive advantage of variants of concern by calculating variant-specific reproduction numbers from wastewater. Together, this study demonstrates the power of national-scale WBE to support public health and promises particular value for countries without extensive individual monitoring.

A bacterial genome in transition--an exceptional enrichment of IS elements but lack of evidence for recent transposition in the symbiont Amoebophilus asiaticus.

BACKGROUND: Insertion sequence (IS) elements are important mediators of genome plasticity and are widespread among bacterial and archaeal genomes. The 1.88 Mbp genome of the obligate intracellular amoeba symbiont Amoebophilus asiaticus contains an unusually large number of transposase genes (n = 354; 23% of all genes). RESULTS: The transposase genes in the A. asiaticus genome can be assigned to 16 different IS elements termed ISCaa1 to ISCaa16, which are represented by 2 to 24 full-length copies, respectively. Despite this high IS element load, the A. asiaticus genome displays a GC skew pattern typical for most bacterial genomes, indicating that no major rearrangements have occurred recently. Additionally, the high sequence divergence of some IS elements, the high number of truncated IS element copies (n = 143), as well as the absence of direct repeats in most IS elements suggest that the IS elements of A. asiaticus are transpositionally inactive. Although we could show transcription of 13 IS elements, we did not find experimental evidence for transpositional activity, corroborating our results from sequence analyses. However, we detected contiguous transcripts between IS elements and their downstream genes at nine loci in the A. asiaticus genome, indicating that some IS elements influence the transcription of downstream genes, some of which might be important for host cell interaction. CONCLUSIONS: Taken together, the IS elements in the A. asiaticus genome are currently in the process of degradation and largely represent reflections of the evolutionary past of A. asiaticus in which its genome was shaped by their activity.

Complex Interplay Between MAZR and Runx3 Regulates the Generation of Cytotoxic T Lymphocyte and Memory T Cells.

The BTB zinc finger transcription factor MAZR (also known as PATZ1) controls, partially in synergy with the transcription factor Runx3, the development of CD8 lineage T cells. Here we explored the role of MAZR as well as combined activities of MAZR/Runx3 during cytotoxic T lymphocyte (CTL) and memory CD8+ T cell differentiation. In contrast to the essential role of Runx3 for CTL effector function, the deletion of MAZR had a mild effect on the generation of CTLs in vitro. However, a transcriptome analysis demonstrated that the combined deletion of MAZR and Runx3 resulted in much more widespread downregulation of CTL signature genes compared to single Runx3 deletion, indicating that MAZR partially compensates for loss of Runx3 in CTLs. Moreover, in line with the findings made in vitro, the analysis of CTL responses to LCMV infection revealed that MAZR and Runx3 cooperatively regulate the expression of CD8α, Granzyme B and perforin in vivo. Interestingly, while memory T cell differentiation is severely impaired in Runx3-deficient mice, the deletion of MAZR leads to an enlargement of the long-lived memory subset and also partially restored the differentiation defect caused by loss of Runx3. This indicates distinct functions of MAZR and Runx3 in the generation of memory T cell subsets, which is in contrast to their cooperative roles in CTLs. Together, our study demonstrates complex interplay between MAZR and Runx3 during CTL and memory T cell differentiation, and provides further insight into the molecular mechanisms underlying the establishment of CTL and memory T cell pools.

Rapid, early and accurate SARS-CoV-2 detection using RT-qPCR in primary care: a prospective cohort study (REAP-1).

OBJECTIVES: We explore the importance of SARS-CoV-2 sentinel surveillance testing in primary care during a regional COVID-19 outbreak in Austria. DESIGN: Prospective cohort study. SETTING: A single sentinel practice serving 22 829 people in the ski-resort of Schladming-Dachstein. PARTICIPANTS: All 73 patients presenting with mild-to-moderate flu-like symptoms between 24 February and 03 April, 2020. INTERVENTION: Nasopharyngeal sampling to detect SARS-CoV-2 using real-time reverse transcriptase-quantitative PCR (RT-qPCR). OUTCOME MEASURES: We compared RT-qPCR at presentation with confirmed antibody status. We split the outbreak in two parts, by halving the period from the first to the last case, to characterise three cohorts of patients with confirmed infection: early acute (RT-qPCR reactive) in the first half; and late acute (reactive) and late convalescent (non-reactive) in the second half. For each cohort, we report the number of cases detected, the accuracy of RT-qPCR, the duration and variety of symptoms, and the number of viral clades present. RESULTS: Twenty-two patients were diagnosed with COVID-19 (eight early acute, seven late acute and seven late convalescent), 44 patients tested SARS-CoV-2 negative and 7 were excluded. The sensitivity of RT-qPCR was 100% among all acute cases, dropping to 68.1% when including convalescent. Test specificity was 100%. Mean duration of symptoms for each group were 2 days (range 1-4) among early acute, 4.4 days (1-7) among late acute and 8 days (2-12) among late convalescent. Confirmed infection was associated with loss of taste. Acute infection was associated with loss of taste, nausea/vomiting, breathlessness, sore throat and myalgia; but not anosmia, fever or cough. Transmission clusters of three viral clades (G, GR and L) were identified. CONCLUSIONS: RT-qPCR testing in primary care can rapidly and accurately detect SARS-CoV-2 among people with flu-like illness in a heterogeneous viral outbreak. Targeted testing in primary care can support national sentinel surveillance of COVID-19.

An inhibitor-mediated beta-cell dedifferentiation model reveals distinct roles for FoxO1 in glucagon repression and insulin maturation.

OBJECTIVE: The loss of forkhead box protein O1 (FoxO1) signaling in response to metabolic stress contributes to the etiology of type II diabetes, causing the dedifferentiation of pancreatic beta cells to a cell type reminiscent of endocrine progenitors. Lack of methods to easily model this process in vitro, however, have hindered progress into the identification of key downstream targets and potential inhibitors. We therefore aimed to establish such an in vitro cellular dedifferentiation model and apply it to identify novel agents involved in the maintenance of beta-cell identity. METHODS: The murine beta-cell line, Min6, was used for primary experiments and high-content screening. Screens encompassed a library of small-molecule drugs representing the chemical and target space of all FDA-approved small molecules with an automated immunofluorescence readout. Validation experiments were performed in a murine alpha-cell line as well as in primary murine and human diabetic islets. Developmental effects were studied in zebrafish and C. elegans models, while diabetic db/db mouse models were used to elucidate global glucose metabolism outcomes. RESULTS: We show that short-term pharmacological FoxO1 inhibition can model beta-cell dedifferentiation by downregulating beta-cell-specific transcription factors, resulting in the aberrant expression of progenitor genes and the alpha-cell marker glucagon. From a high-content screen, we identified loperamide as a small molecule that can prevent FoxO inhibitor-induced glucagon expression and further stimulate insulin protein processing and secretion by altering calcium levels, intracellular pH, and FoxO1 localization. CONCLUSIONS: Our study provides novel models, molecular targets, and drug candidates for studying and preventing beta-cell dedifferentiation.

SARS-CoV-2 mutations in MHC-I-restricted epitopes evade CD8+ T cell responses.

CD8+ T cell immunity to SARS-CoV-2 has been implicated in COVID-19 severity and virus control. Here, we identified nonsynonymous mutations in MHC-I-restricted CD8+ T cell epitopes after deep sequencing of 747 SARS-CoV-2 virus isolates. Mutant peptides exhibited diminished or abrogated MHC-I binding in a cell-free in vitro assay. Reduced MHC-I binding of mutant peptides was associated with decreased proliferation, IFN-γ production and cytotoxic activity of CD8+ T cells isolated from HLA-matched COVID-19 patients. Single cell RNA sequencing of ex vivo expanded, tetramer-sorted CD8+ T cells from COVID-19 patients further revealed qualitative differences in the transcriptional response to mutant peptides. Our findings highlight the capacity of SARS-CoV-2 to subvert CD8+ T cell surveillance through point mutations in MHC-I-restricted viral epitopes.