Single Cell RNA-Seq Characterization of an Adaptive Population of NK Cells after Primary CMV Infection in Rhesus Macaques
Authors: Cordes Stefan , Mortlock Ryland D , Truitt Lauren , Yang Di , Espinoza Diego A , Fan Xing , Ram Daniel , Mostrom Matilda , Tran Dollnovan , Sprehe Lesli , Reeves R. Keith , Donahue Robert E. , Kelly Michael , Chen Jinguo , Hong So Gun , Kaur Amitinder , Wu Chuanfeng , Dunbar Cynthia E. 
Published: November 2021

Abstract
By virtue of their direct cytotoxicity to transformed and virus infected cells, Natural Killer (NK) cells play crucial roles in immunity. NK cells modulate and coordinate innate and adaptive responses through the release of chemokines and cytokines. Although NK cells are endowed only with germ-line encoded receptors, evidence has been accumulating, that subsets of NK cells can bestow adoptively transferable, long-lasting and antigen-specific immune responses to certain haptens and viruses. Growing evidence suggests that adaptive immune responses lie on a spectrum. Rechallenge of cells, canonically belonging to the innate immune system, can result in enhanced responsiveness - a process termed 'trained immunity' and thought to be maintained by epigenetic and metabolic reprogramming. In previous work, our lab studied the role of NK cell responses to rhesus cytomegalovirus (rhCMV) in a genetic barcoding model. We found that new clones arose in the CD16 + NK compartment after primary rhCMV infection. There was rapid clearance without the emergence of new clones in subsequent rechallenge with rhCMV.
In this study we used 3'-end single cell RNA-seq (3'-scRNA-seq) with CITE-seq to profile NK and T cells from an initially CMV-naïve rhesus macaque (RM) at four time points before and after primary and secondary infections with rhCMV. We immunophenotypically sorted NK and T cells from peripheral Blood (PB) samples at 'baseline', 30 days after initial rhCMV infection ('primary infection'), ca. 500 days after initial rhCMV infection ('steady state') and 10 days after rmCMV reinfection ('secondary infection'). Alongside the PB samples at 'steady state' and after 'secondary infection', we also sorted NK and T cells from lymph nodes (LN). We applied CD16 and CD56 CITE-seq antibodies to NK cells from all samples; NK cells from the 'steady state' and 'secondary infection' samples were also labeled with CX3CL1 CITE-seq antibodies. We multiplexed NK and T cells from each time point in 4:1 ratios before preparing 3'-scRNA-seq libraries. We used scanpy and scvi-tools as well as custom python code to demultiplex NK from T cells, harmonize 3'-scRNA-seq with CITE-seq data and integrate the 6 different samples. We used scvelo and cellrank to compute RNA velocities and infer trajectories, respectively.
We obtained a total of 35,523 high-quality cells. We identified 20 clusters of NK and T cells, on the basis of community detection via the Leiden algorithm. All clusters contained cells from both tissue sources. The 4 clusters characterized by expression of CD56 exhibit higher expression of KLRC1 (protein: NKG2A), IL7R and the transcription factors LEF1 and MYC. The 8 clusters of CD16 + cells are distinguished by high expression of the transcription factors ZEB2 and TBX21/T-BET, cytotoxicity markers, GZMB and PRF1, and activating receptors, KLRC2 (protein: NKG2C), KLRC3 (protein: NKG2E) and NCR3 (protein: NKp30). An adaptive population of NK cells is identified on the basis of high KLRC2 and low FCER1G expression. We analyzed changes in the proportions of cells in each cluster of the time course of CMV infection using a binomial generalized linear model. Clusters associated with proliferation and acute inflammation were increased in proportion after primary rhCMV infection; the proportion of the adaptive population did not significantly change during the acute phase of primary infection but increased markedly by the later 'steady state' samples. RNA velocity and inferred developmental trajectories suggest transitions between the adaptive, proliferating CD16 + and mature effector subsets; the predominant path into the adaptive population occurring from the proliferating CD16 + subset after primary infection. There is a notable paucity of inferred transitions between the CD56 + and CD16 +subpopulations under all the experimental conditions we observed.
We have characterized the single cell transcriptional states and dynamics of RM NK cells in response to rhCMV infection. We focus on a subset transcriptionally resembling a previously identified subset with adaptive function and find it arises from a proliferating population of effector cells after primary infection. This may be analogous to the dedifferentiation of effector CD8 T cells into memory T cells proposed by Youngblood et al. Confirmatory experiments to analyze the reconstitution of the CD16 + compartment after treatment with a depleting antibody are on-going.
Figure 1 Figure 1.

Disclosures
No relevant conflicts of interest to declare.