Autoantibody discovery across monogenic, acquired, and COVID-19-associated autoimmunity with scalable PhIP-seq
Abstract
Phage immunoprecipitation sequencing (PhIP-seq) allows for unbiased, proteome-wide autoantibody discovery across a variety of disease settings, with identification of disease-specific autoantigens providing new insight into previously poorly understood forms of immune dysregulation. Despite several successful implementations of PhIP-seq for autoantigen discovery, including our previous work (Vazquez et al., 2020), current protocols are inherently difficult to scale to accommodate large cohorts of cases and importantly, healthy controls. Here, we develop and validate a high throughput extension of PhIP-seq in various etiologies of autoimmune and inflammatory diseases, including APS1, IPEX, RAG1/2 deficiency, Kawasaki disease (KD), multisystem inflammatory syndrome in children (MIS-C), and finally, mild and severe forms of COVID-19. We demonstrate that these scaled datasets enable machine-learning approaches that result in robust prediction of disease status, as well as the ability to detect both known and novel autoantigens, such as prodynorphin (PDYN) in APS1 patients, and intestinally expressed proteins BEST4 and BTNL8 in IPEX patients. Remarkably, BEST4 antibodies were also found in two patients with RAG1/2 deficiency, one of whom had very early onset IBD. Scaled PhIP-seq examination of both MIS-C and KD demonstrated rare, overlapping antigens, including CGNL1, as well as several strongly enriched putative pneumonia-associated antigens in severe COVID-19, including the endosomal protein EEA1. Together, scaled PhIP-seq provides a valuable tool for broadly assessing both rare and common autoantigen overlap between autoimmune diseases of varying origins and etiologies.
Editor's evaluation
This is an important paper that is methodologically compelling. The work presents a series of enhancements to the PhIP-seq method of autoantibody discovery, improving scaling to larger cohorts and control populations, and increasing the ability to discover disease-specific immune responses. The approach is used to discover a novel, frequent autoantibody response (BTNL8) in IPEX patients, and will be an accessible approach to investigate the presence and specificity of autoantibodies in diseases where these have been difficult to define.
https://doi.org/10.7554/eLife.78550.sa0Introduction
Autoantibodies provide critical insight into autoimmunity by informing specific protein targets of an aberrant immune response and serving as predictors of disease. Previously, disease-associated autoantigens have been discovered through candidate-based approaches or by screening tissue specific expression libraries. By contrast, the development of proteome-wide screening approaches has enabled the systematic and unbiased discovery of autoantigens.
Two complementary approaches for proteome-wide autoantibody discovery include printed protein arrays and phage-immunoprecipitation sequencing (PhIP-seq) (Jeong et al., 2012; Larman et al., 2011; Sharon and Snyder, 2014; Zhu et al., 2001). While powerful, printed protein arrays can be cost- and volume-prohibitive and are not flexible to adapting or generating new antigen libraries. On the other hand, PhIP-seq, originally developed by Larman et al., 2011, uses the economy of scale of arrayed oligonucleotide synthesis to enable very large libraries at comparatively low cost. However, phage-based techniques have remained hindered by labor-intensive protocols that prevent broader accessibility and scaling.
Recently, we and others have adapted and applied PhIP-seq to detect novel, disease-associated autoantibodies targeting autoantigens across a variety of autoimmune contexts (Larman et al., 2011; Mandel-Brehm et al., 2019; O’Donovan et al., 2020; Vazquez et al., 2020). However, both technical and biological limitations exist. From a technical standpoint, PhIP-seq libraries express programmed sets of linear peptides, and thus this technique is inherently less sensitive to detect reactivity to conformational antigens, such as Type I interferons (IFNs) and GAD65 (Vazquez et al., 2020; Wang et al., 2021). Nonetheless, the technique allows hundreds of thousands to millions of discrete peptide sequences to be represented in a deterministic manner, including the multiplicity of protein isoforms and variants that are known to exist in vivo. In this sense, PhIP-seq is complementary to mass spectrometry and other techniques that leverage fully native proteins. Given that many forms of autoimmunity exhibit significant phenotypic heterogeneity, the true number of patients with shared disease-associated autoantibodies may be low (Ferre et al., 2016). Therefore, the screening of large cohorts is an essential step for identifying shared antigens and would benefit from a scaled PhIP-seq approach for high throughput testing.
Beyond the benefits of detecting low-frequency or low-sensitivity antigens, a scaled approach to PhIP-seq would also facilitate increased size of healthy control cohorts. Recently, PhIP-seq has been deployed to explore emerging forms of autoimmune and inflammatory disease, including COVID-19-associated multisystem inflammatory syndrome in children (MIS-C) (Gruber et al., 2020). However, these studies suffer from a relative paucity of control samples, resulting in low confidence in the disease specificity of the suggested autoantigens. Questions of disease specificity, rare antigens, and antigen overlap can be addressed in larger, scaled experiments.
Here, we develop a high-throughput PhIP-seq protocol with markedly increased accessibility (not requiring robotics) and scale (enabling 600–800 samples to be run in parallel), with minimal plate-to-plate variability and low contamination potential, and without sacrificing data quality. We demonstrate the utility of this protocol in the context of an expanded, multi-cohort study, including APS1, patients with immune dysregulation, polyendocrinopathy, X-linked (IPEX), RAG1/2 deficiency with autoimmune features, a KD patient cohort, and emerging COVID-19 patient phenotypes with possible autoimmune underpinnings. In the future, scaled PhIP-seq cohort studies could be used across additional syndromic and sporadic autoimmune diseases to develop an atlas of linear B cell autoantigens.
Results and discussion
Design and implementation of a scaled PhIP-seq protocol
The ability to process large numbers of patient samples for PhIP-seq in a highly uniform manner has several important benefits, including reduction of batch effects between samples from the same cohorts as well as between disease and control cohorts, detection of lower-frequency autoantigens, and the ability to simultaneously include large numbers of control samples.
In creating a scaled protocol, we searched for a method that would allow 600–800 samples to be run fully in parallel to reduce any batch or plate-to-plate variability. Thus, each wash or transfer step needed to be performed in rapid succession for all plates. We also prioritized reduction of any well-to-well contamination, particularly given that small, early contamination can amplify across subsequent rounds of immunoprecipitation. Finally, we required our protocol to minimize consumable waste and maximize benchtop accessibility without robotics or other specialized equipment.
A benchtop vacuum-based protocol (rapid, consistent wash times) in deep-well 96-well filter plates with single-use foil seals (no well-to-well contamination) met our requirements (see schematic in Figure 1A). The data for APS1 samples on our moderate-throughput manual multichannel protocols were closely correlated with vacuum-based output, including identification of previously validated antigens within each sample (Figure 1B). Additional protocol improvements included: 3-D printing template for vacuum plate adaptors (for easy centrifugation and incubation steps); direct addition of protein A/G beads to Escherichia coli lysis without a preceding elution step; shortened lysis step by using square-well plates with semi-permeable membrane for aeration; and options for smaller volume and decreased reagent library preparation in both 96-well and 384-well formats. A detailed protocol, including custom part designs, for both high-throughput vacuum-based and moderate throughput multichannel-based protocols is available at protocols.io.
Large control cohorts are critical for identifying disease-specific autoantibodies
Some assays for autoantigen discovery, such as protein arrays, are often used as a quantitative measure of how autoreactive an individual serum sample may be. In contrast, PhIP-seq is an enrichment-based assay in which binders are serially enriched and amplified. A practical limitation of this technique is that non-specific phage may also be amplified, in addition to a wide array of autoreactive, but non-disease-related peptides. We tested whether we could detect global differences between case and control cohorts as a measure of autoreactivity. As each APS1 patient is known to have multiple, high-affinity antibodies to self-proteins (Fishman et al., 2017; Landegren et al., 2016; Meyer et al., 2016; Vazquez et al., 2020) we reasoned that this would be an ideal cohort to determine whether a global autoreactive state was discernible. As expected, each individual sample exhibits a spectrum of enriched genes, regardless of disease status (Figure 2A), indicating that measures of simple enrichment are inadequate for discrimination of cases from controls.
We and others have shown that PhIP-seq can robustly detect disease-associated antigens by comparing antigen-specific signal between disease and control cohorts (Larman et al., 2011; Mandel-Brehm et al., 2019; O’Donovan et al., 2020; Vazquez et al., 2020). In this dataset, encompassing 186 control samples and 128 APS1 samples, we further evaluated the importance of control cohort size. We iteratively downsampled the number of healthy control samples in our dataset to 5, 10, 25, 50, 100, or 150 (out of n=186 total control samples). The number of apparent hits was determined in each condition, where a gene-level hit was called when the following criteria were met: (1) at least 10% of APS1 samples and less than two control samples with a Z-score >10, (2) no control sample exhibiting higher signal than the highest patient signal, and (3) a minimum of one strong positive patient sample (50-fold enrichment over mock-IP). Genes that failed to meet these conditions were considered non-specific. Using these conservative criteria, a control set of 10 samples resulted in (on average) 404 apparent hits, while increasing the control set to 50 samples removed an additional 388 non-specific hits, leaving 16 apparent hits (Figure 2B). Further increasing the number of control samples to 150 samples had diminishing returns, although 4–5 more autoantigen candidates were removed as being non-specific, reducing the frequency of false positives, and ultimately leaving only ~1% of the original candidate list for further investigation. In sum, to improve downstream analysis aimed at detecting disease-associated hits, PhIP-seq experimental design should include a large and appropriate number of non-disease controls.
Scaled PhIP-seq replicates and expands autoantigen repertoires across multiple independent cohorts of APS1
We previously identified and validated PhIP-seq hits based on shared positivity of a given hit among 3 (out of 39) or more patients with APS1 (Vazquez et al., 2020). While this enabled us to robustly detect frequently shared antigens within a small disease cohort, antigens with lower frequencies – or with low detection sensitivity – would likely fall below this conservative detection threshold.
To improve both sensitivity and specificity, we utilized scaled PhIP-seq to explore expanded cohorts of disease, including a much larger (n=128) APS1 cohort spanning two North American cohorts and one Swedish cohort. All hits present in at least 10% of APS1 patients also spanned all three cohorts, thus further validating broad prevalence of antigens that were previously described by us (RFX6, ACPT, TRIM50, CROCC2, GIP, NKX6-3, KHDC3L) and others (SOX10, LCN1) (Figure 2C; Fishman et al., 2017; Hedstrand et al., 2001; Vazquez et al., 2020). At the gene level, we detected 39 candidate hits that were present in 6/128 (>4%) of APS1 and in 2 or fewer controls (2/186, 1%) (Figure 3A and B). As expected, the larger cohort yielded new candidate antigens that had not been detected previously. For example, PDYN is a secreted opioid peptide thought to be involved in regulation of addiction-related behaviors (reviewed in Fricker et al., 2020). PDYN was enriched by 6/128 (4.7%) of patient samples and subsequently was validated by Radioligand Binding Assay (RLBA) (Figure 3C). Indeed, this validated antigen was present in our previous investigations (Vazquez et al., 2020); however, it was enriched in too few samples to qualify for follow-up.
Other notable hits with relevant tissue-restricted expression were also observed. For example, SPAG17 is closely related to the known APS1 autoantigen SPAG16 and is expressed primarily in male germ cells and in the lung, with murine genetic mutations resulting in ciliary dyskinesis with pulmonary phenotypes (Fishman et al., 2017; Teves et al., 2013). Also potentially related to ciliary and/or pulmonary autoimmunity is CROCC2/Rootletin, a protein expressed in ciliated cells, which we previously observed, and now recognize at a high frequency across multiple cohorts (Uhlén et al., 2015; Yang et al., 2002). Similarly, GAS2L2 is a ciliary protein expressed in the airway, with genetic loss of function in mice resulting in impaired mucociliary clearance, and clustered closely with CROCC2 (Bustamante-Marin et al., 2019; Uhlén et al., 2015) in this dataset. These novel putative antigens together hint at potential mucociliary airway autoreactivity. CT45A10 and GPR64 are both proteins with expression restricted primarily to male gonadal tissues (Uhlén et al., 2015). GABRR1 is a GABA receptor expressed primarily in the central nervous system as well as on platelets (Ge et al., 2006; Zhu et al., 2019), and TRIM2 is implicated in genetic disorders of demyelination within the peripheral nervous system, and therefore may be of interest to the chronic inflammatory demyelinating polyneuropathy phenotype that can be seen in some patients with APS1 (Li et al., 2020; Valenzise et al., 2017). In addition to our previously described intestinally expressed autoantigens RFX6, NKX6-3, and GIP, we also identify CDHR5, a transmembrane cadherin-family protein expressed on the enterocyte cell surface, as a putative autoantigen in APS1 (Crawley et al., 2014; Uhlén et al., 2015).
APS1 disease prediction by machine learning
APS1 is a clinically heterogeneous disease, and it is also heterogeneous with respect to autoantibodies (Ferre et al., 2016; Fishman et al., 2017; Landegren et al., 2016; Meyer et al., 2016; Vazquez et al., 2020). Because PhIP-seq simultaneously interrogates autoreactivity to hundreds of thousands of peptides, we hypothesized that unsupervised machine learning techniques could be used to create a classifier that would distinguish APS1 cases from healthy controls. We applied a simple logistic regression classifier to our full gene-level APS1 (n=128) and control (n=186) datasets, resulting in excellent prediction of disease status (AUC = 0.95, Figure 4A) using fivefold cross-validation. Moreover, we found that the classification model was driven strongly by many of the previously identified autoantigens, including RFX6, KHDC3L, and others (Figure 4A), in addition to some targets that had not been previously examined (Vazquez et al., 2020). These results demonstrate that PhIP-seq autoreactive antigen enrichment profiles are amenable to machine learning techniques, and further suggest that such data could be used to derive diagnostic signatures with strong clinical predictive value.
Antoantibody discovery in IPEX
IPEX syndrome is characterized by defective peripheral immune tolerance secondary to impaired T regulatory cell (Treg) function. In IPEX, peripheral tolerance rather than central tolerance is impaired, resulting in a phenotypic constellation of autoimmunity that partially overlaps with APS1 (Bacchetta et al., 2006; Powell et al., 1982). Notably, the majority of IPEX patients exhibit severe enteropathy, with early-onset severe diarrhea and failure to thrive, with many of these children harboring anti-enterocyte antibodies detected by indirect immunofluorescence (Bacchetta et al., 2006; Gambineri et al., 2018; Powell et al., 1982). We hypothesized that the same PhIP-seq approach that was successful for APS1 would also yield informative hits for IPEX. A total of 27 patient samples were analyzed using scaled PhIP-seq, and the data processed in the same manner as for APS1.
A handful of IPEX autoantibodies targeting antigens expressed in the intestine have been described, including harmonin (USH1C) and ANKS4B (Eriksson et al., 2019; Kobayashi et al., 2011). In our data, enrichment of USH1C was observed in two IPEX patients, and this signal was fully correlated with anti-ANKS5B as previously described (Figure 5—figure supplement 1; Eriksson et al., 2019).
Several novel putative autoantigens were observed that were shared among three or more IPEX patients (Figure 5A). Among these were several with expression restricted to the intestine, including BEST4, a protein expressed by a specific subset of enterocytes, BTNL8, a butyrophilin-like molecule highly expressed in the gut epithelium, ST6GALNAC1, and ITGA4 (Figure 5A; Mayassi et al., 2019; Schaum et al., 2018; Uhlén et al., 2015). BEST4 and BTNL8 were selected for validation by whole protein immunoprecipitation. A total of 4/26 (15%) of IPEX patients were positive for anti-BEST4 autoantibodies (Figure 5B). In the case of BTLN8, orthogonal validation identified 11/26 (42%) of IPEX patients who were positive for anti-BTNL8 antibodies (Figure 5B). Of note, all patients with anti-BTNL8 and/or BEST4 antibodies also had clinical evidence of enteropathy (Supplementary file 1). In addition, to broaden the validity of our findings, we tested a second, independent cohort of IPEX patients for anti-BTNL8 antibodies and found 3/15 (20%) patients to be positive (Figure 5C). Taken together, these results suggest that anti-BEST4 and anti-BTNL8 are associated with IPEX enteropathy.
Overlap of intestinal autoantigen BEST4 in the setting of hypomorphic RAG1/2 mutations
Hypomorphic RAG1/2 mutations represent an additional and notoriously phenotypically heterogeneous form of monogenic immune dysregulation. Absent RAG complex activity leads to lack of peripheral T and B cells, therefore causing severe combined immunodeficiency (SCID). However, patients with hypomorphic RAG1/2 have residual capacity to generate T and B cells. Depending on the severity of the defect, these patients can present with Omenn syndrome, atypical SCID (AS), or combined immune deficiency with granulomas and autoimmunity (CID-G/AI) (Delmonte et al., 2018; Delmonte et al., 2020). Autoimmune manifestations are particularly common in patients with AS and CID-G/AI. Cytopenias are the most frequent form of autoimmunity in patients with RAG deficiency, but cutaneous, neuromuscular, and intestinal manifestations have been also reported. While anti-cytokine antibodies have been described, other disease-associated autoantibody targets remain to be identified (Delmonte et al., 2020).
62 patients with hypomorphic RAG1/2 mutations were screened by PhIP-seq to assess for overlap with APS1 and IPEX antigens, as well as for novel autoantigen specificities. Minimal overlap was observed between RAG1/2 deficiency and APS1 and IPEX. However, two samples from RAG1/2 patients indicated the presence of anti-BEST4 antibodies (Figure 5D), which were confirmed through orthogonal validation using whole protein (Figure 5E). Both positive patients had CID-G/AI, indicating the presence of autoimmune features. Remarkably – given that enteropathy in the setting of hypomorphic RAG1/2 deficiency is rare – one of the two individuals harboring anti-BEST4 antibodies also had very early-onset inflammatory bowel disease.
Several other putative antigens among the larger RAG1/2 deficiency cohort were revealed. Nearly half of the cohort sera were enriched for peptides derived from ZNF365 (Figure 5—figure supplement 1B). ZNF365 is a protein associated with multiple autoimmune diseases, as evidenced by genome-wide association studies (GWAS) show that variants in ZNF365 are associated with Crohn’s disease and autoimmune uveitis (Haritunians et al., 2011; Hou et al., 2020). Many patients also had evidence of autoantibodies targeting RELL2, a TNF receptor binding partner, and CEACAM3, a phagocyte receptor that recognizes human-specific pathogens and is important for opsonin-independent phagocytosis of bacteria (Bonsignore et al., 2019; Moua et al., 2017). Autoantibodies targeting these antigens could potentially play a role in the autoinflammation seen in certain cases of hypomorphic RAG1/2 and/or increased susceptibility to particular infections and will require additional future follow-up.
PhIP-seq identifies rare, shared candidate autoantigens in MIS-C
MIS-C leads to critical illness in ~70% of affected children and shares some common clinical features with KD, the most common cause of acquired pediatric heart disease in the United States. Despite hints for a role of abnormal adaptive immunity and autoantibodies in the pathogenesis of KD and MIS-C, the etiologies of both diseases remain enigmatic (Feldstein et al., 2020; Newburger et al., 2016). Recently, PhIP-seq has been deployed to explore COVID-19-associated MIS-C (Gruber et al., 2020). However, this study included only four healthy controls and nine MIS-C patients, and as our results have shown, removal of false-positive PhIP-seq hits requires the use of substantial numbers of unaffected controls (Figure 2C). Furthermore, these previously published hits lacked orthogonal validation. Therefore, we sought to examine an MIS-C cohort in light of these results, as well as to explore for possible autoantibody overlap between KD and MIS-C.
First, 20 MIS-C subjects, 20 pediatric febrile controls, and 20 COVID-19 controls were examined by PhIP-seq, each of which was compared to a cohort of adult healthy controls (n=87). No evidence for specific enrichment was observed for any of the previously reported candidate antigens that overlapped with our PhIP-seq library (Figure 6A). Methodologic and sample differences could, in part, account for differences in our results; however, these results suggest that PhIP-seq hits should be subjected to external replication and/or validation. Additionally, MIS-C patients are treated with IVIG, which has abundant low-affinity, low-avidity autoantibodies, so caution must be exercised in the timing of sample collection (Burbelo et al., 2021; Sacco et al., 2021). While it remains unclear whether high-affinity, novel autoantibodies are likely to be present in IVIG, the majority of our MIS-C samples are confirmed to be collected pre-IVIG. Though several of our samples are of unknown IVIG status, each of the candidate autoantibodies discussed below is present in at least one sample known to be IVIG-free.
Analysis of our MIS-C cohort for shared candidate hits yielded only three candidate hits, each in 2/20 patient sera. These were CD34, RPS6KB1, and CAPZB (Figure 6B). While these targets may be of interest, disease-association remains uncertain. These results suggest that a much larger MIS-C cohort, controlled by an equally large set of healthy controls, will be required to detect rare, shared antigens with confidence by PhIP-seq.
PhIP-seq screen of a cohort of KD patients
To screen for possible KD-specific autoantibodies, we analyzed a large cohort of 70 KD subjects by PhIP-seq. KD patients are also often treated with IVIG, so care was taken to ensure that each of these samples was collected prior to IVIG administration. Using the same hit selection criteria as previously, we detected 25 shared hits among 3 or more of the 70 KD samples, which were specific to KD relative to adult healthy controls (Figure 6C). Of these 25 shared KD hits, 17 were absent from additional control groups including the febrile pediatric patients. Each of these hits was present in only a small subset of KD samples, suggesting significant heterogeneity among samples.
Some of the candidate antigens have possible connections to the systemic inflammation seen in KD, including SPATA2 (6/70, 8.6%) and ALOX5B (4/70, 5.7%). SPATA2 is a protein known to regulate the TNF receptor response, with murine knockout of SPATA2 resulting in increased activation of NFkB and MAPK (Schlicher and Maurer, 2017). Similarly, ALOX15 family of lipoxygenases is known to be responsive to Th2-induced anti-inflammatory cytokine IL-4 and IL-13 in human macrophages and thereby likely plays a role in suppressing inflammatory responses, and polymorphisms have been linked to the development of coronary artery disease (Snodgrass and Brüne, 2019; Wuest et al., 2014).
Other KD candidate autoantigens exhibited tissue expression patterns that suggest possible relationships to sub-phenotypes of KD. For example, pancreatitis and psoriasis have been reported as rare complications of KD (Asano et al., 2005; Haddock et al., 2016). We found that 7/70 of the KD patients (10%) have increased signal for autoantibodies targeting FBXL19, a protein with associations with both acute pancreatitis (Ma et al., 2020) and psoriasis (Stuart et al., 2010) and 10/70 for pancreas-expressed protein RESP18 (Zhang et al., 2007).
Taken together, our dataset did not uncover the presence of common autoantibodies in KD or MIS-C. Nonetheless, our findings leave open the possibility of lower-frequency disease-associated autoantibodies, and future studies with increased cohort scaling, diverse representation of clinical subphenotypes, and high-sensitivity follow-up assays may shed further light on the role of autoantibodies in KD.
Autoantigen overlap between KD and MIS-C
Given the partial clinical overlap of KD with MIS-C, we also searched for the low-frequency, shared KD hits within MIS-C and found that 6 of the 17 KD hits were present in one or more MIS-C samples (Figure 6D). Of these, CGNL1 was of particular interest given the very high enrichment values in patients and absence of signal from all controls, as well as its anatomic expression pattern. CGNL1 is an endothelial junction protein and is highly expressed in the cardiac endothelium (Chrifi et al., 2017; Schaum et al., 2018). We confirmed anti-CGNL1 autoantibodies using a radioligand binding assay with whole CGNL1 protein (Figure 6D). We also identified an additional positive KD patient that had not been detected by PhIP-seq. These data suggest that anti-CGNL1 antibodies, while rare, may be associated with KD and/or MIS-C.
Application of scaled PhIP-seq to severe COVID-19
Recently, it was reported that over 10% of severe COVID-19 pneumonia is characterized by the presence of anti-Type I IFN autoantibodies, a specificity that overlaps with APS1 (Bastard et al., 2020; Bastard et al., 2021b; Bastard et al., 2021a; Meager et al., 2006; Meyer et al., 2016; van der Wijst et al., 2021). We therefore looked for possible overlap of additional antigens between APS1 and COVID-19. As expected, we had low sensitivity for the known anti-type I IFN autoantibodies by PhIP-seq, likely due to the conformational nature of these antigens; however, this cohort had been previously assessed for anti-Type I IFN antibodies by other techniques (Bastard et al., 2020; Bastard et al., 2021c). We did not detect substantial overlap of any of the other antigens that were found in 5% or more of APS1 samples, suggesting that autoantibody commonalities between APS1 and severe COVID-19 may be limited (Figure 7A).
PhIP-seq was then used to investigate patients with severe to critical COVID-19 pneumonia (n=117) relative to patients with mild to moderate COVID-19 (n=55) and non-COVID-19 healthy controls (n=125) (Figure 7A). A small number of putative autoantigens were identified in four or more of the severe COVID-19 group (>4%), but not within our control group, including MCAM, CHRM5, and EEA1. Of these, only EEA1, an early endosomal protein, had a frequency of >10% (Mu et al., 1995). Notably, the set of anti-EEA1 positive patients was almost fully distinct from the anti-Type I IFN antibody positive group, with only one patient shared between the groups (Figure 7B). Given the importance of the finding of autoantibodies to Type 1 IFNs in >15% of critical COVID-19 pneumonia, our data suggest that anti-EEA1 antibodies may also have predictive potential worth further investigation.
Discussion
PhIP-seq is a powerful tool for antigen discovery due to its throughput and scalability. Continually declining costs in sequencing paired with scaled protocols such as the ones presented here result in a low per-sample cost and experimental time, and declining costs of custom oligonucleotide based-libraries allow for extensive adaptation (Román-Meléndez et al., 2021; Vogl et al., 2021). Yet, as the technology is relatively new, increased discussion around best practices in experimental design, methodology, and data interpretation would be beneficial. In this work, we contribute to these efforts by presenting an accessible, carefully documented scaled lab protocol. We anticipate that scaling and low cost will facilitate cross-platform comparisons of autoimmune sera, enabling a more systematic understanding of each platform’s abilities and shortcomings, as well as the best approaches to combine platforms for comprehensive autoantibody profiling, including fixed protein arrays, yeast display (REAP), and PhIP-seq (Wang et al., 2021).
Our scaled data reinforces the importance of including large control cohorts to improve specificity and orthogonal validations to evaluate sensitivity. The absence of appropriately sized control cohorts can result in false-positive associations between enriched protein sequences and disease which can lead to misinterpretation. Furthermore, PhIP-seq is optimized for linear antigens and thus is inherently less robust for detection of confirmational or modified epitopes. For example, our previous work highlights that known anti-IFN and anti-GAD65 antibodies can be detected only in a handful of APS1 samples by PhIP-seq, while orthogonal assays using whole conformation protein demonstrate increased sensitivity (Vazquez et al., 2020). Similarly, the increased sensitivity of whole protein validation relative to PhIP-seq for BTNL8 autoantibodies highlights the need for rapid and sensitive secondary assays to confirm or even increase the importance of a given candidate antigen. These principles may be extended to other related PhIP-seq modalities (Mina et al., 2019; Román-Meléndez et al., 2021; Vogl et al., 2021).
Several novel APS1 antigens were identified in this work by virtue of the large cohort size and large control set, enabling the detection of low-frequency hits such as PDYN. PDYN is a secreted opioid precursor of the central nervous system, and while PDYN-knockout mice are largely phenotypically normal, they do experience hyperalgesia and altered anxiety-related behaviors relative to wildtype mice (McLaughlin et al., 2003; Sharifi et al., 2001; Wang et al., 2001). Future studies will be required to determine if anti-PDYN antibodies may themselves mediate these phenotypes. In addition to uncovering lower-frequency autoantigens, antigen overlap between APS1 and other syndromic autoimmune diseases, including IPEX and RAG1/2 deficiency, were evaluated. Detection of common APS1 antigens, such as RFX6 and KHDC3L, was rare among other disease contexts, suggesting APS1 displays limited sharing with respect to other multiorgan autoimmune syndromes, including IPEX and RAG deficiency.
In contrast, antigen overlap between IPEX and RAG-hypomorphic patients was detected in the form of anti-BEST4 antibodies. BEST4 is a well-characterized, cell-surface intestinal ion channel (Qu and Hartzell, 2008). Recently, BEST4 has become a standard and specific marker for a subset of enterocytes in the duodenum and colon, including CFTR+ enterocytes, which are involved in fluid homeostasis (Busslinger et al., 2021; Elmentaite et al., 2021; Smillie et al., 2019). The presence of an IPEX antigen BEST4 within CID-G/AI suggests a possible etiologic link between the two diseases, possibly relating to Treg dysfunction. Furthermore, given that BEST4 autoantibodies were found across two distinct etiologies of IBD, in addition to the intestinal localization of BEST4 expression, future experiments specifically searching for the presence of anti-BEST4 antibodies in IBD and other forms of autoimmune enteropathy are warranted.
Another novel IPEX antigen, BTNL8, is also of particular interest given its high frequency in IPEX (11/26) as well as its biological functions. BTN and BTNL-family members belong to the family of B7 co-stimulatory receptors known to modulate T-cell responses, with structural similarity to CD80 and PD-L1 (Chapoval et al., 2013; Rhodes et al., 2016). Broadly, BTNL-family members are thought to participate primarily in regulation of gamma-delta T cells, with BTNL8/BTNL3 implicated in gut epithelial immune homeostasis (Di Marco Barros et al., 2016; Chapoval et al., 2013). Given the cell-surface expression pattern of BTNL8, it is conceivable that antibodies to these proteins could play a functional role in immune checkpoints, rather than only the bystander role implicated for most autoantibodies to intracellular antigens. Interestingly, recent studies in patients with celiac disease have demonstrated a reduction in BTNL3/BTNL8 expression following acute episodes of inflammation, with an associated loss of the physiological normal gamma/delta T-cell subset of gut intraepithelial lymphocytes (Mayassi et al., 2019).
While some disease cohorts yield strongly enriched common antigens, such as BTNL8 and BEST4, other cohorts appear heterogeneous. Such was the case for MIS-C and KD, where our results suggested only rare, shared putative antigens. Nonetheless, a number of observations point to a possible role for abnormal adaptive immunity and autoantibodies (reviewed in Sakurai, 2019). Both diseases respond to varying degrees to IVIG therapy, and genetic studies suggest some children with polymorphisms in genes involving B-cell and antibody function can be predisposed to developing KD (reviewed in Onouchi, 2018). Autoantibodies targeting clinically relevant tissues such as the heart or endothelium have been detected in KD, though their role in disease pathogenesis is uncertain (Cunningham et al., 1999; Fujieda et al., 1997; Grunebaum et al., 2002). Although anti-CGNL1 antibodies were only found at low frequency, it is possible that these antibodies may represent a subset of specificities within anti-endothelial cell antibodies, given the endothelial expression of CGNL1 as well as its implications in cardiovascular disease (Chrifi et al., 2017). Particularly in the setting of MIS-C, expanded cohorts both to further validate CGNL1 as well as for as-of-yet undescribed antigen specificities may add additional information to these otherwise poorly understood diseases.
There is increasing evidence for a role for autoantibodies in other acquired disease states that have not been classified per se as autoimmune diseases, including severe to critical COVID-19 pneumonia (Bastard et al., 2020; Bastard et al., 2021c; van der Wijst et al., 2021). Peptides derived from EEA1 were enriched in 11% of patients with severe COVID-19 pneumonia but not in patients with mild COVID-19. Interestingly, similarly to anti-IFN antibodies, anti-EEA1 antibodies have also previously been reported in the setting of systemic lupus erythematosus (Selak et al., 2000; Stinton et al., 2004; Waite et al., 1998). While intriguing, future studies will be needed to determine whether the specificity of EEA1 autoantibodies to severe-critical COVID-19 can be replicated in independent cohorts, as well as in orthogonal assays.
This work suggests that PhIP-seq data may have predictive or diagnostic value when the cohorts are sufficiently large. Implementation of a logistic regression machine learning model with cross-validation successfully discriminates APS1 samples from healthy control samples with high sensitivity and specificity (Figure 4A). Thus, scaled PhIP-seq represents an avenue to not only discover and compare novel autoantigens within and across cohorts of autoimmune disease, but also suggests that machine learning techniques will be effective tools for determining, without supervision, the underlying proteins that drive successful classification using this data type. This, in turn, may have diagnostic or prognostic value in a wide spectrum of immune dysregulatory contexts.
Finally, we emphasize that to broadly understand autoantibody profiles, full PhIP-seq data availability is critical for enabling expanded analysis and comparison across datasets, cohorts, and research groups. Importantly, groups with expertise in specific disease areas may add value by re-analyzing and evaluating candidate hits using additional metrics for hit prioritization, including increased weighting of orthogonal expression and/or genetic data. Full PhIP-seq data for all cohorts presented here will be linked to this publication and will be available for download at Dryad.
Materials and methods
PhIP-seq protocols
Request a detailed protocolAll PhIP-seq protocols described in detail are available at protocols.io, located at the links below.
Vacuum-based scaled: https://www.protocols.io/view/scaled-high-throughput-vacuum-phip-protocol-ewov1459kvr2/v1
Multichannel-based scaled: https://www.protocols.io/view/scaled-moderate-throughput-multichannel-phip-proto-8epv5zp6dv1b/v1
Library Preparation: https://www.protocols.io/view/phage-display-library-prep-method-rm7vz3945gx1/v1
PhIP-seq data alignment and normalization
Request a detailed protocolFastq files were aligned at the level of amino acids using RAPSearch2 (Zhao et al., 2012). For gene-level analysis, all peptide counts mapping to the same gene were summed. 0.5 reads were added to all genes, and raw reads were normalized by converting to percentage of total reads per sample. Fold change (FC) over mock-IP was calculated on a gene-by-gene basis by dividing sample read percentage by mean read percentage in corresponding AG bead-only samples. Z-scores were calculated using FC values; for each disease sample by using all corresponding healthy controls, and for each healthy control samples by using all other healthy controls. The positive threshold used for detection of shared candidate antigens was a Z-score ≥10. Shared hits were then determined by positive rate in the specified percentage of patient samples and under a specified percentage (<2% unless otherwise specified) in controls. In addition, at least one positive sample was required to have a minimum FC of 50 or above. Finally, no candidate antigens were allowed where any positive control samples signal fell above the highest patient sample.
Radioligand binding assay
Request a detailed protocolRadioligand binding assay was performed as described in Vazquez et al., 2020 for each gene of interest (GOI). Briefly, GOI-myc-flag construct was used to in vitro transcribe and translate [35 S]-methionine-labeled protein, which was subsequently column purified, immunoprecipitated with patient serum on Sephadex protein A/G beads, and counts per million (cpm) were read out on a 96-well liquid scintillation reader (MicroBeta Trilux, Perkin Elmer). Constructs: PDYN (Origene, #RC205331), BEST4 (Origene, #RC211033), BTNL8(Origene, #RC215370), and CGNL1 (Origene, #RC223121). Anti-MYC (Cell Signaling Technologies, #2272 S) and/or anti-FLAG (Cell Signaling Technologies, #1493 S) antibodies were used as positive controls for immunoprecipitations.
Statistics
Statistics for radioligand binding assay were performed as described in Vazquez et al., 2020. Briefly, the antibody index for each sample was calculated as: (sample value – mean blank value)/(positive control antibody value – mean blank value).
We applied a logistic regression classifier on log-transformed PhIP-seq RPK values from APS1 patients (n=128) versus healthy controls (n=186) using the scikit-learn package (Pedregosa et al., 2011) with a liblinear solver and L1 regularization. The model was evaluated with fivefold cross-validation.
Ethics/human subjects research
Available clinical metadata varies across cohorts; all available clinical data is included in Supplementary file 1.
APS1
Request a detailed protocolNorth America - 1: All patient cohort data was collected and evaluated at the NIH, and all APECED/ APS1 patients were enrolled in a research study protocol approved by the NIH Institutional Review Board Committee and provided with written informed consent for study participation (protocol #11-I-0187, NCT01386437). All NIH patients gave consent for passive use of their medical record for research purposes. The majority of this human cohort data was previously published Ferre et al., 2016; Ferré et al., 2019.
Sweden: Serum samples were obtained from Finnish and Swedish patients with APS1. All individuals met the clinical diagnostic criteria for APS1, requiring two of the hallmark components: chronic mucocutaneous candidiasis, hypoparathyroidism, and adrenal insufficiency, or at least one of the hallmark components in siblings or children of APS1 patients. The project was approved by local ethical boards and was performed in accordance with the declarations of Helsinki. All patients and healthy subjects had given their informed consent for participation.
North America - 2: All patients underwent informed consent with research study protocols approved by the UCSF Human Research Protection Program (IRB# 10-02467).
RAG1/2 deficiency
Request a detailed protocolPatients with RAG1/2 deficiency were enrolled in research study protocols approved by the NIAID, NIH Institutional Review Board (protocols 05-I0213, 06-I-0015, NCT03394053, and NCT03610802). Peripheral blood samples were obtained upon written informed consent.
IPEX
Request a detailed protocolIPEX patient cohort data was collected and evaluated at Seattle Children’s Hospital under an IRB-approved protocol. Clinical data on this human cohort was previously published by Gambineri et al., 2018. For the independent validation cohort (Figure 5C), all patient cohort data was collected and evaluated at Stanford under an IRB-approved protocol and was previously published by Narula et al., 2022.
KD/MIS-C/febrile controls
Request a detailed protocolUCSD cohort (MIS-C/KD/febrile controls): The study was reviewed and approved by the institutional review board at the University of California, San Diego. Written informed consent from the parents or legal guardians and assent from patients were obtained as appropriate.
Rockefeller cohort (MIS-C): All individuals were recruited according to protocols approved by local Institutional Review Boards.
Severe COVID-19
Request a detailed protocolAll COVID-19 (non-MIS-C) patients were collected between March 01, 2020 and July 21, 2020 and had positive results by SARS-CoV-2 RT-PCR in nasopharyngeal swabs. ZSFG remnant specimens (institutional review board [IRB] number 20-30387) were approved by the IRB of the University of California, San Francisco. The committee judged that written consent was not required for use of remnant specimens.
Adult healthy controls
Request a detailed protocolNew York Blood Center & Vitalant Research Institute: Healthy, pre-COVID control plasma were obtained as deidentified samples. These samples were part of retention tubes collected at the time of blood donations from volunteer donors who provided informed consent for their samples to be used for research.
Rockefeller cohort: All individuals were recruited according to protocols approved by local Institutional Review Boards.
Data availability
Full PhIP-Seq data for all cohorts presented is available for download at Dryad at https://doi.org/10.5061/dryad.qfttdz0k4. All available deidentified clinical data for this study is available in Supplementary file 1.
-
Dryad Digital RepositoryAutoantibody discovery across monogenic, acquired, and COVID19-associated autoimmunity with scalable PhIP-Seq.https://doi.org/10.5061/dryad.qfttdz0k4
References
-
Acute pancreatitis with kawasaki disease: analysis of cases with elevated serum amylase levelsEuropean Journal of Pediatrics 164:180–181.https://doi.org/10.1007/s00431-004-1589-4
-
Defective regulatory and effector T cell functions in patients with FOXP3 mutationsThe Journal of Clinical Investigation 116:1713–1722.https://doi.org/10.1172/JCI25112
-
Auto-antibodies to type I ifns can underlie adverse reactions to yellow fever live attenuated vaccineThe Journal of Experimental Medicine 218:e20202486.https://doi.org/10.1084/jem.20202486
-
Preexisting autoantibodies to type I ifns underlie critical COVID-19 pneumonia in patients with APS-1The Journal of Experimental Medicine 218:e20210554.https://doi.org/10.1084/jem.20210554
-
CEACAM3-A prim(at)E invention for opsonin-independent phagocytosis of bacteriaFrontiers in Immunology 10:3160.https://doi.org/10.3389/fimmu.2019.03160
-
Lack of GAS2L2 causes pcd by impairing cilia orientation and mucociliary clearanceAmerican Journal of Human Genetics 104:229–245.https://doi.org/10.1016/j.ajhg.2018.12.009
-
Cgnl1, an endothelial junction complex protein, regulates GTPase mediated angiogenesisCardiovascular Research 113:1776–1788.https://doi.org/10.1093/cvr/cvx175
-
Anti-human cardiac myosin autoantibodies in kawasaki syndromeJournal of Immunology 163:1060–1065.
-
Rag deficiency: two genes, many diseasesJournal of Clinical Immunology 38:646–655.https://doi.org/10.1007/s10875-018-0537-4
-
The autoimmune targets in IPEX are dominated by gut epithelial proteinsThe Journal of Allergy and Clinical Immunology 144:327–330.https://doi.org/10.1016/j.jaci.2019.02.031
-
Multisystem inflammatory syndrome in U.SChildren and Adolescents. New Engl J Medicine 383:EJMoa2021680.https://doi.org/10.1056/NEJMoa2021680
-
Lymphocyte-driven regional immunopathology in pneumonitis caused by impaired central immune toleranceScience Translational Medicine 11:eaav5597.https://doi.org/10.1126/scitranslmed.aav5597
-
Five decades of research on opioid peptides: current knowledge and unanswered questionsMolecular Pharmacology 98:96–108.https://doi.org/10.1124/mol.120.119388
-
Antibodies to endothelial cells in Kawasaki disease lysE endothelial cells without cytokine pretreatmentClinical and Experimental Immunology 107:120–126.https://doi.org/10.1046/j.1365-2249.1997.d01-894.x
-
The role of anti-endothelial cell antibodies in Kawasaki disease-in vitro and in vivo studiesClinical and Experimental Immunology 130:233–240.https://doi.org/10.1046/j.1365-2249.2002.02000.x
-
Psoriasiform eruptions during Kawasaki disease (kD): a distinct phenotypeJournal of the American Academy of Dermatology 75:69–76.https://doi.org/10.1016/j.jaad.2016.02.1146
-
The transcription factors SOX9 and SOX10 are vitiligo autoantigens in autoimmune polyendocrine syndrome type IThe Journal of Biological Chemistry 276:35390–35395.https://doi.org/10.1074/jbc.M102391200
-
Rapid identification of monospecific monoclonal antibodies using a human proteome microarrayMolecular & Cellular Proteomics 11:111.https://doi.org/10.1074/mcp.O111.016253
-
Autoantigen discovery with a synthetic human peptidomeNature Biotechnology 29:535–541.https://doi.org/10.1038/nbt.1856
-
A recessive TRIM2 mutation causes an axonal neuropathy in miceNeurobiology of Disease 140:104845.https://doi.org/10.1016/j.nbd.2020.104845
-
Analysis of associations of FBXL19-AS1 with occurrence, development and prognosis of acute pancreatitisEuropean Review for Medical and Pharmacological Sciences 24:12763–12769.https://doi.org/10.26355/eurrev_202012_24176
-
Kelch-Like protein 11 antibodies in seminoma-associated paraneoplastic encephalitisThe New England Journal of Medicine 381:47–54.https://doi.org/10.1056/NEJMoa1816721
-
Kappa opioid receptor antagonism and prodynorphin gene disruption block stress-induced behavioral responsesThe Journal of Neuroscience 23:5674–5683.
-
RELT family members activate p38 and induce apoptosis by a mechanism distinct from TNFR1Biochemical and Biophysical Research Communications 491:25–32.https://doi.org/10.1016/j.bbrc.2017.07.022
-
EEA1, an early endosome-associated proteinJournal of Biological Chemistry 270:13503–13511.https://doi.org/10.1074/jbc.270.22.13503
-
Epigenetic and immunological indicators of IPEX disease in subjects with FOXP3 gene mutationThe Journal of Allergy and Clinical Immunology S0091-6749(22)01221-0.https://doi.org/10.1016/j.jaci.2022.09.013
-
Kawasaki diseaseJournal of the American College of Cardiology 67:1738–1749.https://doi.org/10.1016/j.jacc.2015.12.073
-
The genetics of Kawasaki diseaseInternational Journal of Rheumatic Diseases 21:26–30.https://doi.org/10.1111/1756-185X.13218
-
PedregosaScikit-learn: machine learning in pythonJournal of Machine Learning Research 12:2825–2830.
-
An X-linked syndrome of diarrhea, polyendocrinopathy, and fatal infection in infancyThe Journal of Pediatrics 100:731–737.https://doi.org/10.1016/s0022-3476(82)80573-8
-
Bestrophin cl- channels are highly permeable to HCO3-American Journal of Physiology. Cell Physiology 294:C1371–C1377.https://doi.org/10.1152/ajpcell.00398.2007
-
Regulation of immunity by butyrophilinsAnnual Review of Immunology 34:151–172.https://doi.org/10.1146/annurev-immunol-041015-055435
-
Autoimmune aspects of Kawasaki diseaseJournal of Investigational Allergology & Clinical Immunology 29:251–261.https://doi.org/10.18176/jiaci.0300
-
Generation of dynorphin knockout miceBrain Research. Molecular Brain Research 86:70–75.https://doi.org/10.1016/s0169-328x(00)00264-3
-
Serum profiling using protein microarrays to identify disease related antigensMethods in Molecular Biology 1176:169–178.https://doi.org/10.1007/978-1-4939-0992-6_14
-
Regulation and functions of 15-lipoxygenases in human macrophagesFrontiers in Pharmacology 10:719.https://doi.org/10.3389/fphar.2019.00719
-
Genome-Wide association analysis identifies three psoriasis susceptibility lociNature Genetics 42:1000–1004.https://doi.org/10.1038/ng.693
-
Sperm-Associated antigen–17 gene is essential for motile cilia function and neonatal survivalAmerican Journal of Respiratory Cell and Molecular Biology 48:765–772.https://doi.org/10.1165/rcmb.2012-0362OC
-
Type I interferon autoantibodies are associated with systemic immune alterations in patients with COVID-19Science Translational Medicine 13:eabh2624.https://doi.org/10.1126/scitranslmed.abh2624
-
Autoantibodies to a novel early endosome antigen 1Clinical Immunology and Immunopathology 86:81–87.https://doi.org/10.1006/clin.1997.4455
-
Pronociceptive actions of dynorphin maintain chronic neuropathic painThe Journal of Neuroscience 21:1779–1786.
-
Association of polymorphisms in the ALOX15B gene with coronary artery diseaseClinical Biochemistry 47:349–355.https://doi.org/10.1016/j.clinbiochem.2013.12.013
-
Rootletin, a novel coiled-coil protein, is a structural component of the ciliary rootletThe Journal of Cell Biology 159:431–440.https://doi.org/10.1083/jcb.200207153
Article and author information
Author details
Funding
National Institute of Allergy and Infectious Diseases (5P01AI118688)
- Mark S Anderson
National Institute of Allergy and Infectious Diseases (1ZIAAI001175)
- Michail S Lionakis
National Institute of Diabetes and Digestive and Kidney Diseases (1F30DK123915)
- Sara E Vazquez
Chan Zuckerberg Biohub
- Joseph L DeRisi
Parker Institute for Cancer Immunotherapy
- Mark S Anderson
Juvenile Diabetes Research Foundation United States of America
- Mark S Anderson
Helmsley Charitable Trust
- Mark S Anderson
National Institute of General Medical Sciences (5T32GM007618)
- Mark S Anderson
American Diabetes Association (1-19-PDF-131)
- Zoe Quandt
UCSF-CTSI TL1 (TR001871)
- Zoe Quandt
Division of Intramural Research, National Institute of Allergy and Infectious Diseases (1 ZIA AI001222)
- Luigi D Notarangelo
Eunice Kennedy Shriver National Institute of Child Health and Human Development (1R61HD105590)
- Adriana Tremoulet
Laboratory of Human Genetics of Infectious Diseases
- Jean-Laurent Casanova
FRM (EA20170638020)
- Paul Bastard
Imagine Institute (MD-PhD program)
- Paul Bastard
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank members of the DeRisi lab, members of the Anderson lab, Chukwuka Didigu, and Joseph M Replogle for helpful discussions. We thank the New York Blood Center and Vitalant Research Institute for providing deidentified healthy control plasma. The Laboratory of Human Genetics of Infectious Diseases is supported by the Howard Hughes Medical Institute, the Rockefeller University, the St. Giles Foundation, the National Institutes of Health (NIH) (R01AI088364 and R01AI163029), the National Center for Advancing Translational Sciences (NCATS), NIH Clinical and Translational Science Award (CTSA) program (UL1TR001866), the Fisher Center for Alzheimer's Research Foundation, the Meyer Foundation, the JPB Foundation, the French National Research Agency (ANR) under the "Investments for the Future" program (ANR-10-IAHU-01), the Integrative Biology of Emerging Infectious Diseases Laboratory of Excellence (ANR-10-LABX-62-IBEID), the French Foundation for Medical Research (FRM) (EQU201903007798), ANRS Nord-Sud (ANRS-COV05), ANR grants GENVIR (ANR-20-CE93-003), AABIFNCOV (ANR-20-CO11-0001), and GenMIS-C (ANR-21-COVR-0039), the European Union's Horizon 2020 research and innovation programme under grant agreement no. 824110 (EASI-Genomics), the Square Foundation, Grandir - Fonds de solidarité pour l'enfance, the SCOR Corporate Foundation for Science, Fondation du Souffle, Institut National de la Santé et de la Recherche Médicale (INSERM), REACTing-INSERM, and the University of Paris.
Ethics
Human subjects: Detailed information on consent, where applicable, is available in the methods section of the manuscript.
Copyright
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Metrics
-
- 2,145
- views
-
- 327
- downloads
-
- 29
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
- Immunology and Inflammation
The immunosuppressive microenvironment in pancreatic ductal adenocarcinoma (PDAC) prevents tumor control and strategies to restore anti-cancer immunity (i.e. by increasing CD8 T-cell activity) have had limited success. Here, we demonstrate how inducing localized physical damage using ionizing radiation (IR) unmasks the benefit of immunotherapy by increasing tissue-resident natural killer (trNK) cells that support CD8 T activity. Our data confirms that targeting mouse orthotopic PDAC tumors with IR together with CCR5 inhibition and PD1 blockade reduces E-cadherin positive tumor cells by recruiting a hypoactive NKG2D-ve NK population, phenotypically reminiscent of trNK cells, that supports CD8 T-cell involvement. We show an equivalent population in human single-cell RNA sequencing (scRNA-seq) PDAC cohorts that represents immunomodulatory trNK cells that could similarly support CD8 T-cell levels in a cDC1-dependent manner. Importantly, a trNK signature associates with survival in PDAC and other solid malignancies revealing a potential beneficial role for trNK in improving adaptive anti-tumor responses and supporting CCR5 inhibitor (CCR5i)/αPD1 and IR-induced damage as a novel therapeutic approach.
-
- Immunology and Inflammation
During thymic development, most γδ T cells acquire innate-like characteristics that are critical for their function in tumor surveillance, infectious disease, and tissue repair. The mechanisms, however, that regulate γδ T cell developmental programming remain unclear. Recently, we demonstrated that the SLAM/SAP signaling pathway regulates the development and function of multiple innate-like γδ T cell subsets. Here, we used a single-cell proteogenomics approach to identify SAP-dependent developmental checkpoints and to define the SAP-dependent γδ TCR repertoire in mice. SAP deficiency resulted in both a significant loss of an immature Gzma+Blk+Etv5+Tox2+ γδT17 precursor population and a significant increase in Cd4+Cd8+Rorc+Ptcra+Rag1+ thymic γδ T cells. SAP-dependent diversion of embryonic day 17 thymic γδ T cell clonotypes into the αβ T cell developmental pathway was associated with a decreased frequency of mature clonotypes in neonatal thymus, and an altered γδ TCR repertoire in the periphery. Finally, we identify TRGV4/TRAV13-4(DV7)-expressing T cells as a novel, SAP-dependent Vγ4 γδT1 subset. Together, the data support a model in which SAP-dependent γδ/αβ T cell lineage commitment regulates γδ T cell developmental programming and shapes the γδ TCR repertoire.