Transcriptome network analysis implicates CX3CR1-positive type 3 dendritic cells in non-infectious uveitis
Abstract
Background:
Type I interferons (IFNs) promote the expansion of subsets of CD1c+ conventional dendritic cells (CD1c+ DCs), but the molecular basis of CD1c+ DCs involvement in conditions not associated without elevated type I IFNs remains unclear.
Methods:
We analyzed CD1c+ DCs from two cohorts of non-infectious uveitis patients and healthy donors using RNA-sequencing followed by high-dimensional flow cytometry to characterize the CD1c+ DC populations.
Results:
We report that the CD1c+ DCs pool from patients with non-infectious uveitis is skewed toward a gene module with the chemokine receptor CX3CR1 as the key hub gene. We confirmed these results in an independent case–control cohort and show that the disease-associated gene module is not mediated by type I IFNs. An analysis of peripheral blood using flow cytometry revealed that CX3CR1+ DC3s were diminished, whereas CX3CR1− DC3s were not. Stimulated CX3CR1+ DC3s secrete high levels of inflammatory cytokines, including TNF-alpha, and CX3CR1+ DC3 like cells can be detected in inflamed eyes of patients.
Conclusions:
These results show that CX3CR1+ DC3s are implicated in non-infectious uveitis and can secrete proinflammatory mediators implicated in its pathophysiology.
Funding:
The presented work is supported by UitZicht (project number #2014-4, #2019-10, and #2021-4). The funders had no role in the design, execution, interpretation, or writing of the study.
Editor's evaluation
These findings are valuable to ocular immunologists who the study pathophysiologic mechanisms driving inflammation in human uveitis, and for future identification of novel therapeutic targets. The authors convincingly perform high dimensional multi-omic analysis of testing and replication cohorts, followed by characterization of a disease-specific cell type using comparative analysis with previously validated experimental datasets. The analysis will be of particular interest to basic and translational ocular immunologists, as well as dendritic cell biologists.
https://doi.org/10.7554/eLife.74913.sa0Introduction
Non-infectious uveitis refers to a group of chronic inflammatory eye diseases that are among the leading causes of preventable vision loss in western countries (Thorne et al., 2016; Suttorp-Schulten and Rothova, 1996). Currently, little is known about the disease mechanism of non-infectious uveitis. A large body of mechanistic studies using experimental autoimmune uveitis (EAU) in rodents suggest that T cells play a role in non-infectious uveitis (Lee et al., 2014; Caspi, 2010). Human non-infectious uveitis is characterized by inflammation and T cells infiltrating the eyes through unknown mechanisms. The genetic association between non-infectious uveitis and MHC, and ERAP1, ERAP2 genes indicates that antigen presentation is central to the etiology (Kuiper and Venema, 2020; Huang et al., 2020; Kuiper et al., 2018; Márquez et al., 2017). A key antigen presenting cell type are dendritic cells. Despite their important role in EAU, dendritic cells have yet to be fully investigated in human non-infectious uveitis (Chen et al., 2015b; Fu et al., 2019; Wang et al., 2021). CD1c-positive conventional dendritic cells (CD1c+ DCs) have been found to be associated with disease activity (Chen et al., 2016; Chen et al., 2015a; Chen et al., 2014), and are abundant in eye fluid of patients (O’Rourke et al., 2018). In order to understand the role of CD1c+ DCs in non-infectious uveitis, it is necessary to understand their functions.
Several single-cell studies have revealed that the CD1c+ DCs (and its murine equivalent, termed ‘cDC2s’) consists of multiple subsets derived from distinct progenitors (Villani et al., 2017; Dutertre et al., 2019; Cytlak et al., 2020). The type I IFN family of cytokines promotes the expansion of a subset of CD1c+ DCs called ‘DC3’ (Dutertre et al., 2019; Girard et al., 2020; Bourdely et al., 2020). DC3s are increased in blood of type I IFN-driven systemic lupus erythematosus (SLE) patients (Dutertre et al., 2019). A significant difference between non-infectious uveitis and SLE is that active uveitis is accompanied by lower levels of type I IFN (Wang et al., 2019; Kuiper et al., 2022). It should be noted that although type I IFNs can induce lupus-like disease, they can also suppress non-infectious uveitis (Wang et al., 2019; Rönnblom et al., 1991; Rönnblom et al., 1990), pointing to an alternative disease mechanism implicating CD1c+ DCs in non-infectious uveitis. Hence, we do not fully understand the characteristics of CD1c+ DC during autoimmunity, especially in conditions not driven by type I IFNs.
For the purpose of characterizing the core transcriptional features and subset composition of CD1c+ DCs in autoimmunity of the eye, we used whole transcriptome profiling by bulk RNA-sequencing of peripheral blood CD1c+ DCs and multiparameter flow cytometry of two cohorts of non-infectious uveitis patients and healthy donors. We constructed co-expression networks that identified a robust gene module associated with non-infectious uveitis in patients that helped identify a CX3CR1-positive CD1c+ DC subset.
Materials and methods
Patients and patient material
Request a detailed protocolThis study was conducted in compliance with the Helsinki principles. Ethical approval was requested and obtained from the Medical Ethical Research Committee in Utrecht (METC protocol number #14-065/M). All patients signed written informed consent before participation. We collected blood from a discovery cohort of 23 and a replication cohort of 28 adult patients (Table 1) with HLA-B27-associated acute anterior uveitis (AU), idiopathic intermediate uveitis (IU), or HLA-A29-associated birdshot uveitis (BU). Patients were recruited at the outbound patient clinic of the Department of Ophthalmology of the University Medical Center Utrecht between July 2014 and January 2017. We recruited 27 age- and sex-matched anonymous blood donors of European Ancestry with no history of ocular inflammatory disease at the same institute to serve as unaffected controls (Table 1). Uveitis was classified and graded in accordance with the SUN classification (Jabs et al., 2005). Each patient underwent a full ophthalmological examination by an ophthalmologist experienced in uveitis, routine laboratory screening, and an X-ray of the lungs. Laboratory screening included erythrocyte sedimentation rate, renal and liver function tests, angiotensin-converting enzyme, and screening for infectious agents in the serum and an Interferon-Gamma Release Assay (IGRA) was obtained for all patients. All patients with AU and BU were HLA-B27 or HLA-A29-positive, respectively (confirmed by HLA typing). All patients had active uveitis (new onset or relapse) and there was no clinical evidence for uveitis-associated systemic inflammatory disease (e.g., rheumatic condition) till the time of sampling. None of the patients received systemic immunomodulatory treatment in the last 3 months, other than low dose (≤10 mg) oral prednisolone in one BU patient of cohort II and one AU patient of cohort I.
CD1c+ DC purification
Request a detailed protocolPeripheral blood mononuclear cells (PBMCs) were isolated by standard ficoll density gradient centrifugation from 70 ml heparinized blood immediately after blood withdrawal (GE Healthcare, Uppsala, Sweden). For the first cohort, 10 batches (individual days) of 4–5 randomly selected patient and control samples of nitrogen stored PBMCs (mean storage time of 11 [range 0–31] months) were carefully thawed and subjected to sorting by the BD FACSAria III sorter after incubation with a panel of surface antibodies (Supplementary file 1A) and fluorescent-activated cell sorting (FACS) buffer (1% bovine serum albumin and 0.1% sodium azide in phosphate-buffered saline [PBS]). CD3−CD19−CD56−CD14−HLA-DR+CD11c+CD1c cells were sorted. The average number of collected cells by sorting was 56,881 (range 6669–243,385). For the second cohort, fresh PBMCs were immediately subjected to magnetic-activated cell sorting (MACS) for the removal (positive selection) of CD304+ cells (pDC), followed by the removal of CD19+ cells (B cell), and subsequently isolation of CD1c+ cells by using the CD1c+ (BDCA1) isolation kit (Miltenyi Biotec, Germany) according to the manufacturer’s instructions. The isolated CD1c+ fraction contained on average 147,114 cells (range 46,000–773,000) and purity was determined by flow cytometry (Supplementary file 1B) measured on the BD LSRFortessa Cell analyzer (Figure 1—figure supplement 1). Data were analyzed using FlowJo software (TreeStar Inc). MACS or FACS purified CD1c+ cells were immediately taken up in a lysis buffer (RLT plus, QIAGEN) containing 1% β-mercaptoethanol, snap frozen on dry ice, and stored at −80°C until RNA extraction was performed. Isolation of CD1c+ DC for functional experiments was done by MACS as described above. Purification of CD1c+DC subsets based on CD36 and CX3CR1 or CD14 expression from freshly isolated PBMCs was conducted by flow cytometry using the panel in Supplementary file 1C and shown in Figure 3—figure supplement 2B.
CD1c+ DC cultures and secretome analysis
Request a detailed protocolPurified CD1c+ DCs were cultured in RPMI Glutamax (Thermo Fisher Scientific) supplemented with 10% heat-inactivated fetal bovine serum (Biowest Riverside) and 1% penicillin/streptomycin (Thermo Fisher Scientific). CD1c+ DCs were cultured at a concentration of 0.5 × 106 cells/ml in a 96-well round-bottom plate (100 μl/well). Cells were stimulated overnight (18 hr) with multiple stimuli listed in Supplementary file 1D. After stimulation, cells were lysed in an RLT plus lysis buffer (QIAGEN) and stored at −80°C until RNA extraction was performed. Cell lysates were stored at −80°C until RNA extraction was performed for qPCR. In separate cultures, CD1c+ DC subsets (sorted based on CD36 and CX3CR1 expression) were cultured in the presence of 1 µg/ml lipoteichoic acid (LTA). After 18 hr of stimulation, supernatants were harvested and IL-23 cytokine production was analyzed by ELISA (R&D Systems). The levels of IL-2, IL-5, IL-6, IL-10, IL-12p70, IL-13, IL-17, IL-22, IL-27, TNF-alpha, IFN-alpha, IFN-beta, CCL1, CXCL10, CXCL13, VEGF, CD40L, FAS, TNFR1, TNFR2, Elastase, and Granzyme B were simultaneously measured in supernatant of CD1c+ DC cultures using the in-house multiplex immunoassay based on Luminex technology, as described previously (Bakker et al., 2022). Protein concentrations that were out of range were replaced with the LLOQ (lower limit of quantification) and ULOQ (upper limit of quantification) for each analyte and divided by 2 for the proteins detected below the range of detection or multiplied by 2 for values above the detection range (Supplementary file 1E).
Real-time quantitative PCR
Request a detailed protocolFirst-strand cDNA was synthesized from total RNA using Superscript IV kit (Thermo Fisher Scientific), and quantitative real-time PCR (RT-qPCR) was performed on the QuantStudio 12k flex System (LifeTechnologies), following the manufacturer’s instructions. Sequences of the primers used are listed in Supplementary file 1F and the Key Resource Table. RT-qPCR data were normalized to the expression of the selected housekeeping gene GUSB (ENSG00000169919). CT values were normalized to GUSB by subtracting the CT mean of GUSB (measured in duplo) from the CT mean of the target mRNA (e.g., RUNX3) = ΔCT. The fold change (FC) of each sample was calculated compared to ΔCt of the medium control using the formula FC = 2−ΔΔCt, where ΔΔCt = ΔCt sample − ΔCt reference.
Flow cytometry of CD1c+ DC populations
Request a detailed protocolPBMC samples from the two cohorts (HC = 11 samples; AU = 9 samples; IU = 6 samples; BU = 11 samples) were randomly selected and measured by flow cytometry in batches of 9–10 mixed samples per run, divided over 4 days. Per batch, 10 million PBMCs per sample were quickly thawed, washed with ice cold PBS and stained with the antibody panel depicted in Supplementary file 1C. PBMCs were incubated with Fixable Viability Dye eF780 (eBioscience) at room temperature for 10 min. Cells were then plated in V-bottomed plates (Greiner Bioone), washed with PBS and incubated for 30 min at 4°C in the dark with Brilliant Stain Buffer (BD) and the fluorescently conjugated antibodies. Next, the cells were washed and taken up in the FACS buffer. Flow cytometric analyses were performed on the BD FACSAria III sorter. Manual gating of data was done using FlowJo software (TreeStar inc San Carlos, CA, USA). FlowSOM v1.18.0 analysis was done as described previously (Laban et al., 2020). Lineage- (negative for CD3/CD56/CD19) HLA-DR+ data were transformed using the logicleTransform function of the flowCore v1.52.1 R package, using default parameters (Ellis et al., 2020). The SOM was trained for a 7 × 7 grid (49 clusters) with 2000 iterations. Consensus hierarchical clustering was used to annotate clusters, based on the ConsensusClusterPlus v1.50.0 R package (Wilkerson and Hayes, 2010). Principal component analysis (PCA) analysis was done on normalized expression data from flowSOM using the factoextra v 1.0.7.999 R package.
RNA isolation and RNA-sequencing
Request a detailed protocolTotal RNA from CD1c+ DC cell lysates from patients and controls was isolated using the AllPrep Universal Kit (QIAGEN) on the QIAcube (QIAGEN) according to the manufacturer’s instructions. For cohort I, RNA-seq libraries were generated by GenomeScan (Leiden, The Netherlands) with the TruSeq RNAseq RNA Library Prep Kit (Illumina Inc, Ipswich, MA, USA), and were sequenced using Illumina HiSeq 4000 generating ~20 million 150 bp paired end reads for each sample. Library preparation and Illumina sequencing was performed on samples of cohort II at BGI (Hong Kong). RNA-seq libraries were generated with the TruSeq RNAseq RNA Library Prep Kit (Illumina Inc, Ipswich, MA, USA) and were sequenced using Illumina NextSeq 500 generating approximately 20 million 100 bp paired end reads for each sample.
Power analysis
Request a detailed protocolWe conducted power analysis using the PROPER R package v 1.22.0 (Wu et al., 2015) with 100 simulations of the build-in RNA-seq count data from antigen presenting (B) cells from a cohort of 41 individuals (i.e., large biological variation as expected in our study) (Cheung et al., 2010). Simulation parameters used the default of 20,000 genes and an estimated 10% of genes being differentially expressed. We detected 0.8 power to detect differentially expressed genes (p < 0.05) at a log2(fold change) >1 for the smallest patient group (9 cases) and we considered the sample size reasonable for analysis.
Differential gene expression and statistical analysis
Request a detailed protocolQuality check of the raw sequences was performed using the FastQC tool. Reads were aligned to the human genome (GRCh38 build 79) using STAR aligner (Dobin et al., 2013) and the Python package HTSeq v0.6.1 was used to count the number of reads overlapping each annotated gene (Anders et al., 2015). We aligned the reads of the RNA-sequencing datasets to 65,217 annotated Ensemble Gene IDs. Raw count data were fed into DESeq2 v1.30.1 (Love et al., 2014) to identify differentially expressed genes between the four disease groups (AU, IU, BU, and HC). Using DESeq2, we modeled the biological variability and overdispersion in expression data following a negative binomial distribution. We used Wald’s test in each disease group versus control pairwise comparison and p values were corrected by the DESeq2 package using the Benjamini–Hochberg Procedure. We constructed co-expression gene networks with the WGCNA v 1.70-3 R package (Langfelder and Horvath, 2008) using the cumulative uveitis-associated genes from all pairwise comparisons and a soft power of 5. Module membership (MM) represents the intramodular connectivity of genes in a gene module. Gene Significance (GS >0.25) indicates a strong correlation between genes and non-infectious uveitis, whereas MM (MM >0.8) indicates a strong correlation with the EigenGene value of the modules. We calculated the intersection between the modules constructed from the two cohorts and used Fisher’s exact test to identify modules that exhibited significant overlap in genes. Gene expression data from runx3-knockout (KO) cDC2s, notch2-KO cDC2s, and ‘inflammatory’ cDC2s were obtained from the NCBI Gene Expression Omnibus (accession numbers GSE48590 [2 wild-type [WT] CD11b+ESAM+ splenic cDC2s vs. 2 CD11b+ESAM+ cDC2s from CD11c-DC-Runx3Δ mice], GSE119242 [2 untreated cDC2 vs. untreated cDC2 from CD11c-Cre notch2f/f mice], GSE149619 [5 CD172+MAR1− cDC2s in mock condition vs 3 CD172+MAR1+ cDC2 in virus condition]) using GEO2R in the GEO database, which builds on the GEOquery v2.58.0 and limma R v 3.46 packages (Davis and Meltzer, 2007; Ritchie et al., 2015). RNA-seq data from the mouse bone marrow stromal cell line OP9 expressing NOTCH ligand DLL1 (OP9-DLL1)-driven cDC2 cultures (GSE110577 [2 sorted CD11c+MHCII+B220−CD11b+ cDC2 from bone marrow cultures with FLT3L for 7 days vs. 2 sorted CD11c+MHCII+B220−CD11b +cDC2 from bone marrow cultures with FLT3L+OP9-DLL1 cells for 7 days]) were analyzed using DESeq2 and normalized count data plotted using the plotCounts function. RNA-seq count data from CD14+/− DC3 subsets from patients with SLE and systemic sclerosis were obtained via GEO (accession number: GSE136731) (Dutertre et al., 2019) and differential expression analysis was conducted using DESeq2 v1.30.1. Gene set enrichment analysis was done using the fgsea R package v1.16.0 and data plotted using the GSEA.barplot function from the PPInfer v 1.16.0 R package (Jung and Ge, 2020). Gene sets for runx3-KO, notch2-KO, inflammatory cDC2s, and cDC2s from OP9-DLL1 bone marrow cultures were generated by taking the top or bottom percentiles of ranked [−log10(p) × sign(log2(fold change))] genes from each dataset as indicated. Genes in the modules of interest that encode cell-surface proteins were identified according to surfaceome predictor SURFY (Bausch-Fluck et al., 2018).
Single-cell RNA-seq analysis of aqueous humor
Request a detailed protocolSingle-cell RNA-seq (scRNA-seq) data from a previous study of as reported by Kasper et al., 2021 of aqueous humor of four HLA-B27-positive anterior uveitis (identical to the AU group in this study) patients were obtained and downloaded via Gene Expression Omnibus (GEO) repository with the accession code GSE178833. Data were processed using the R package Seurat v4.1.0 (Stuart et al., 2019) using R v4.0.3. We removed low-quality cells (<200 or >2500 genes and mitochondrial percentages <5%) and normalized the data using the SCTransform() function accounting for mitochondrial percentage and cell cycle score (Hafemeister and Satija, 2019). Dimensionality reduction for all cells was achieved by adapting the original UMAP coordinates for each barcode as reported by Kasper et al., 2021 (see GSE178833). Data were subjected to scGate v1.0.0 (Andreatta et al., 2022) using CLEC10A+ and C5AR1 (CD88)− cells in our gating model to purify CD1c+ DCs in the scRNAseq dataset. Dimensionality reduction for CD1c+ DCs was conducted using the R package Seurat, and cells clustered using the FindNeighbors and FindClusters functions from Seurat. After clustering and visualization with UMAP, we used the DotPlot function from the Seurat package to visualize the average expression of genes in each cluster.
Results
A CX3CR1 gene module is associated with non-infectious uveitis
We characterized the transcriptome of primary CD1c+ DCs from patients with non-infectious uveitis (Figure 1A). RNA-seq analysis (RNA-seq) was performed on lineage (CD3−CD19−CD56−CD14−)-negative, and HLA-DR-positive, CD11c and CD1c-positive DCs purified from frozen PBMCs by flow cytometry from 36 patients with anterior (AU), intermediate (IU), or posterior non-infectious uveitis (BU) and healthy controls. A co-expression network was constructed using uveitis-associated genes identified by differential expression analyses (n = 2,016 genes at P<0.05, Figure 1B), which identified six modules, of which 3 were associated with non-infectious uveitis (GS >0.25, Figure 1C). The blue module was most associated with non-infectious uveitis (Figure 1C). Based on Module Membership, CX3CR1 was the top hub gene of the blue module (Figure 1D, Supplementary file 1G). Since CX3CR1 was previously associated with a distinct subset cDC2s that may also express CD14 (Brown et al., 2019; Fujita et al., 2019), we attempted to validate and expand the gene set associated with non-infectious uveitis by MACS-isolating CD1c+ DC cells from fresh blood of 28 patients and 14 healthy controls, followed by RNA-seq analysis of the highly purified CD1c+ DCs (median [interquartile range]% = 96 [3]% pure, Figure 1—figure supplement 1). We also constructed a co-expression network for uveitis-associated genes (n = 6794, p < 0.05) in the second cohort (Figure 1E), which revealed 24 gene modules (Supplementary file 1H). Note that patient samples did not cluster according to clinical parameters of disease activity (e.g., cell grade in eye fluid, macular thickness) (Figure 1—figure supplement 2). The three uveitis-associated modules in cohort I shared a significant number of co-expressed genes with one module in cohort II, the black module (Figure 1F). The black module was associated with non-infectious uveitis in cohort II (GS for uveitis >0.25) and CX3CR1 was also the hub gene for this module (Black Module Membership, p = 5.9 × 10−22; Supplementary file 1H, Figure 1G). According to these findings, the overlapping disease-associated gene modules appear to represent a single gene module. In cohort I, the separation of genes into three modules was possibly due to low sensitivity to detect disease-associated genes with low expression, as replicated genes of the black module were typically higher expressed (Figure 1—figure supplement 3). In total, we replicated 147 co-expressed genes between the two cohorts (which we will refer to as the ‘black module’), of which 94% also showed consistent direction of effect (e.g., upregulated in both cohorts) (Figure 1H, Supplementary file 1I). The black module was enriched for the GO term ‘positive regulation of cytokine production’ (GO:0001819, padj = 6.9 × 10−5). In addition to CX3CR1, the black module comprised CD36, CCR2, TLR-6,-7,-8, CD180, and transcription factors RUNX3, IRF8, and NFKB1 (Figure 1I), but not CD14. In summary, these results show that a gene module characterized by CX3CR1 in blood CD1c+ DCs is associated with non-infectious uveitis.
CX3CR1+ DC3 are diminished in peripheral blood of non-infectious uveitis patients
Type I IFN cytokines promote differentiation of CD1c+ DCs (Dutertre et al., 2019; Girard et al., 2020), but patients with active non-infectious uveitis have reduced blood levels of type I IFN cytokines (Wang et al., 2019; Kuiper et al., 2022). Assessment of the transcriptome of CD1c+ DCs from patients, found no enrichment for genes associated with murine type I IFN-dependent cDC2s (Bosteels et al., 2020; Figure 2A,B). Furthermore, while RUNX3 was downregulated in RNA-seq data from CD1c+ DCs from non-infectious uveitis patients (Figure 1I), stimulation of CD1c+ DCs with from healthy human donors with IFN-alpha resulted in upregulation of RUNX3 (Figure 2—figure supplement 1A). In contrast, the transcriptome of CD1c+ DCs from patients overlapped significantly with murine cDC2s knocked out for Runx3, or its upstream regulator Notch2 (Figure 2, Figure 2—figure supplement 1B ,C; Lewis et al., 2011; Fasnacht et al., 2014; Briseño et al., 2018). Given that cDC2 subsets differ by their dependence on NOTCH signaling (Lewis et al., 2011; Fasnacht et al., 2014; Briseño et al., 2018; Kirkling et al., 2018), we hypothesized that the transcriptomic signatures of the CD1c+ DC pool in patients might reflect changes in their proportions.
Therefore, we used flow-cytometry analysis to identify CD1c+ DC clusters in PBMCs samples from 26 cases and 11 controls. We designed a panel based on the black module (CX3CR1, CD36, CCR2, and CD180), other CD1c+ DC markers that were not in the black module (CD1c, CD11c, CD14, CD5, and CD163) (Dutertre et al., 2019; Korenfeld et al., 2017). FlowSOM (Van Gassen et al., 2015) was used on HLA-DR+ and lineage (CD3/CD19/CD56)− PBMCs to cluster cells into a predetermined number of 49 clusters (7 × 7 grid) to facilitate detection of CD1c+ DC phenotypes in blood. The analysis with flowSOM clearly distinguished four CD1c+ DC clusters (cluster number 22, 37, 44, and 45) (Figure 3A and Figure 3—figure supplement 1A and B). We extracted the data for these four CD1c+ DC clusters and conducted principal component analysis (PCA). The PCA biplot identified CD5 and CD163 as top loadings (Figure 3—figure supplement 1C), which defines the DC2s (cluster 45), CD5−CD163− DC3s (cluster 37), and CD5−CD163+ DC3s (clusters 22 and 44) (Dutertre et al., 2019; Figure 3B, C). Among the identified clusters, we detected a significant reduction in the frequency of cluster 44 in patients compared to controls (Welch t-test, p = 0.03, Figure 3D and Figure 3—figure supplement 1D). Clusters 44 as well as cluster 22 were CD36 and CD14 positive, which indicates these clusters may represent CD14+ DC3s in human blood (Dutertre et al., 2019). However, cluster 44 had relatively higher levels of CX3CR1 than cluster 22 (Figure 3E, F). This suggests that DC3s may be phenotypically bifurcated by CX3CR1 independently of CD14. This is supported by weak correlation between CD14 and CX3CR1 in our RNA-seq data from bulk CD1c+ DCs (Pearson correlation coefficient = 0.35, Figure 3—figure supplement 2A). In addition, we sorted CD14-positive and -negative fractions from CD1c+ DCs of six healthy donors (Figure 3—figure supplement 2B) which showed no significant difference in expression levels for CX3CR1 (Figure 3—figure supplement 2C). CX3CR1 levels were also not significantly different in sorted CD5−CD163+CD14-positive and CD14-negative DC3s from patients with autoimmune diseases, further indicating that CX3CR1 expression in CD1c+ DCs may be independent from CD14 expression in CD1c+ DCs (Figure 3—figure supplement 2D).
We validated by manual gating that cluster 44 represents the CD36+CX3CR1+ fraction of CD1c+ DCs in peripheral blood (~25% of total CD1c+ DCs) (Figure 3—figure supplement 1E). Comparison between patients and controls corroborated that the frequency of manual gated CD36+CX3CR1+ DC3s were decreased in the blood of non-infectious uveitis patients (Welch t-test, p = 0.029, Figure 3G). In detail, we show that CD14+CD1c+ DCs double positive for CD36+ and CX3CR1 were significantly decreased (p = 0.026), while CD14+CD1c+ DCs not positive for CX3CR1 were not (p = 0.43) (Figure 3G). This supports that CX3CR1 discerns a phenotypic subpopulation of CD14+ DC3s (Figure 3—figure supplement 3), that was diminished in the blood of patients with non-infectious uveitis.
CX3CR1+ DC3s can secrete pro-inflammatory cytokines upon stimulation
We compared the cytokine-producing abilities of CX3CR1+ DC3s to their negative counterparts since the gene module associated with CX3CR1 was enriched for genes involved in cytokine regulation. To this end, we freshly sorted primary human CD1c+ DC subsets based on the surface expression of CX3CR1 and CD36, of which double-positive and -negative subsets could be sorted from the selected healthy subjects in sufficient numbers for analysis (Figure 4—figure supplement 1). Since CD36 is involved in LTA-induced cytokine production (Jimenez-Dalmaroni et al., 2009), we overnight stimulated the CD1c+ subsets with LTA. Both subsets of CD1c+ DCs secreted IL-23 equally strongly (Figure 4A). To assess the secretome of the CD1c+ DC subsets in more detail, we profiled the supernatants of LTA-stimulated CD1c+ DC subsets for additional soluble immune mediators (Supplementary file 1E): The CD1c+ DC subsets could be distinguished based on the secreted protein profile (Figure 4B), of which the levels of TNF-alpha, IL-6, VEGF-A, and TNFR1 showed significant differences between the subsets (Figure 4C). These results show that CD1c+ DC subsets defined on the basis of surface co-expression of CD36 and CX3CR1 have the capacity to secrete pro-inflammatory mediators that participate in the pathophysiology of human non-infectious uveitis.
CX3CR1+ DC3s are detectable in the inflamed eye during non-infectious uveitis
We speculated that CX3CR1+ DC3s are important in the disease mechanisms of uveitis and may be found at increased abundance in the eye during active uveitis. We used single-cell RNA-sequencing data (scRNA-seq) of eye fluid biopsies of four noninfectious patients from Kasper et al., 2021. Cells positive for the CD1c+ DC-specific tissue-marker CLEC10A and negative for the monocyte marker C5AR1 (CD88) (Dutertre et al., 2019; Bourdely et al., 2020; Heger et al., 2018) were used to identify CD1c+ DCs among other immune cells in the scRNA-seq data (Figure 5A). Unsupervised clustering identified three clusters (1, 2, and 3) of different cells within the CD1c+ DC population (Figure 5B). We identified that cluster 1 expressed the gene profile associated with CX3CR1+ DC3s, including relatively higher levels of CX3CR1, CD36, CCR2, and lower levels of RUNX3 compared to the other two CD1c+ DC clusters (Figure 5C), which is in line with the gene profile identified by our bulk RNA-seq analysis. In summary, we conclude that CD1c+ DCs with a gene expression profile similar to CX3CR1+ DC3s can be detected in the eyes of patients during active non-infectious uveitis.
Discussion
In this study of non-infectious uveitis patients and controls, we identified and replicated a CX3CR1-associated gene module in CD1c+ DCs. We were able to track back the gene module to a CX3CR1+ DC3 subset that was diminished in peripheral blood of patients with non-infectious uveitis.
Preceding studies into human CD1c+ DCs revealed functionally distinct subsets termed ‘DC2’ and ‘DC3’, with the DC3 showing both transcriptomic features reminiscent of cDC2s and monocytes – such as elevated CD36 (Villani et al., 2017; Dutertre et al., 2019). DC3s also have distinct developmental pathways and transcriptional regulators compared to DC2 (Villani et al., 2017; Dutertre et al., 2019; Cytlak et al., 2020; Bourdely et al., 2020). Recently, Cytlak et al. revealed that lower expression of IRF8 is linked to DC3 (Cytlak et al., 2020), a transcription factor that was also decreased in non-infectious uveitis. According to Brown et al., 2019, CD1c+ DCs exhibit two subsets: cDC2A and cDC2B, whereas cDC2B exhibits higher expression of CX3CR1 and produces more TNF-alpha and IL-6 than cDC2A upon stimulation. Accordingly, uveitis-associated CX3CR1+ DC3s described in this study exhibit similar phenotypical and functional features.
Dutertre et al., 2019 showed that the phenotype of peripheral blood CD1c+ DCs can be further segregated according to the expression of CD163 and CD5, with ‘DC3’ cells being characterized as CD5−CD163− or CD5−CD163+ cells and ‘DC2’ as CD5+CD163 cells. Our flow-cytometry results confirm these findings, but we also show that CD5−CD163+ DC3s that express CD14 are composed of CX3CR1-positive and CX3CR1-negative cells, of which the CX3CR1+ population is implicated in non-infectious uveitis.
Single-cell analysis supported that CD1c+ DCs in eye fluid of patients with non-infectious uveitis contain also a population that has a gene profile reminiscent of CX3CR1+ DC3s, with relatively higher levels of CX3CR1, CD36, CCR2, and lower levels of RUNX3. Patients with SLE display accumulation of CD14+ DC3s in blood (Dutertre et al., 2019), while the population of CD14+ DC3 cells was decreased in non-infectious uveitis patients. The differences between non-infectious uveitis and SLE may be related to distinct (i.e., opposite) immunopathological mechanisms; Type I IFNs drive the maturation of cDC2s into ‘inflammatory cDC2s’ (infcDC2s) (Bosteels et al., 2020) and can induce CD1c+ DCs to express a distinct set of surface receptors (Girard et al., 2020). The type I IFN-α drives immunopathology of SLE and administration of type I IFN therapy can induce lupus-like disease (Rönnblom et al., 1991; Rönnblom et al., 1990). In favor of attributing the seemingly contrasting observations in blood CD1c+ subsets between SLE and non-infectious uveitis to distinct biology is the fact that, in contrast to elevated IFN-α in patients with SLE, in non-infectious uveitis patient’s disease exacerbations correlate with reduced blood type I IFN concentrations (Wang et al., 2019; Kuiper et al., 2022; Obermoser and Pascual, 2010). Despite the importance of type I IFN signaling on DC3s, our results suggest that DC3s are also dysregulated in conditions associated with decreased type I IFNs (Wang et al., 2019; Kuiper et al., 2022), supporting additional pathways involved in DC3 regulation during chronic inflammation.
We showed that the gene module of CD1c+ DCs showed overlap with the gene signature of disrupted NOTCH2 signaling in cDC2s. Notch2 signaling is mediated via the NF-κB family member Relb in murine cDC2s (Diener et al., 2021). NF-κB signaling via RelB suppresses type I IFN signaling in cDC2s (Saha et al., 2020) while selective deletion of RelB in dendritic cells protects against autoimmunity (Diener et al., 2021; Andreas et al., 2019). It is tempting to speculate that the enrichment for NOTCH gene signatures implies altered NF-κB-Relb signaling in CD1c+ DCs, with a mechanism that varies between diseases mediated by type I IFNs (e.g., SLE) and type I IFN-negative diseases (e.g., uveitis). Although this warrants further investigation, some circumstantial evidence for this is the presence of NF-κB family members in the black module, such as NFKB1 and NFKBIA, the former associated previously with a CX3CR1+ cDC2s, while the latter can regulate Relb function in dendritic cells (Shih et al., 2012; Sun, 2011). NF-κB-Relb signaling has been shown to suppress type I IFN via a histone demethylase encoded by KDM4A which was also in the black module (Jin et al., 2014). Regardless, the NF-κB pathways are regulated by the TNFR1, the main receptor for TNF-alpha (Hayden and Ghosh, 2014; Maney et al., 2014). TNFR1 is expressed at the cell surface of cDC2s and its ectodomain cleaved by the NOTCH2-pathway regulator ADAM10 (Maney et al., 2014; Yang et al., 2016; Iberg et al., 2022). We showed that both TNF-alpha and sTNFR1 were higher in the secretome of activated CX3CR1+ DC3s. This is in agreement with previous studies on CD36+ DC3s (Villani et al., 2017) or CX3CR1+ cDC2B (T-bet−) that also produced higher levels of TNF-alpha (Brown et al., 2019). Interestingly, altered NF-κB signaling specifically in cDC2 is associated with clinical response to anti-TNF-alpha therapy (Andres-Ejarque et al., 2021). Anti-TNF therapy is effective for treatment of non-infectious uveitis (Touhami et al., 2019), while anti-TNF therapy may also result in a dysregulated type I IFN response (Conrad et al., 2018) indicating potentially cross regulatory mechanisms via NF-κB signaling and type I IFN signaling affecting cDC2s. More research is needed to resolve the regulatory mechanisms driving CD1+ DC changes in type I IFN-negative inflammation, including non-infectious uveitis. It is also possible that the change in CD1c+ DCs observed in this study results from cytokine-induced precursor emigration or differentiation or that the affected peripheral blood DC3s marked by CX3CR1 are in a precursor or pre-activation state.
Other disease modifying factors possibly affect the CD1c+ DC pool in uveitis patients. In mice, antibiotic treatment to experimentally disturb the microbiota affects a cDC2 subset phenotypically similar to CD1c+ DCs and decreases their frequency in the intestine of mice, which suggests microbiota-dependent signals involved in the maintenance of cDC2 subsets (Brown et al., 2019). This is especially interesting in light of growing evidence that microbiota-dependent signals cause autoreactive T cells to trigger uveitis (Horai et al., 2015), which makes it tempting to speculate that gut-resident cDC2 subsets contribute to the activation of T cells in uveitis models. Dietary components can influence subsets of intestinal dendritic cells (Ko et al., 2020). Regardless, most likely, an ensemble of disease modulating factors is involved. For example, myeloid cytokines, such as GM-CSF, contribute to autoimmunity of the eye (Croxford et al., 2015) and GM-CSF has been shown to stimulate the differentiation of human CD1c+ subset from progenitors (Bourdely et al., 2020). However, GM-CSF signaling in conventional dendritic cells has a minor role in the inception of EAU (Bing et al., 2020). Our data support that stimulation of CD1c+subsets with GM-CSF or TLR ligands does not induce the transcriptional features of CD1c+ DCs during non-infectious uveitis, which is in line with previous observations that support that stimulated cDC2s do not convert from one into another subset (Bourdely et al., 2020).
Note that our results of a decreased subset of CD1c+ DCs in non-infectious uveitis are in contrast with previous flow-cytometry reports in non-infectious uveitis (Chen et al., 2016; Chen et al., 2015a; Chen et al., 2014). Chen et al., 2014 reported an increase of CD1c+ myeloid dendritic cells in non-infectious uveitis. It is important to note, however, that their study did not include the DC marker CD11c, thereby including CD1c+CD11c− populations that do not cluster phenotypically with CD1c+ (CD11c+) DCs (e.g., cluster 46, se Figure 3—figure supplement 1A), which may explain the differences compared to our study.
Better understanding of the changes in the CD1c+ DC pool during human non-infectious uveitis will help develop strategies to pharmacologically influence putative disease pathways involved at an early disease stage, which may lay the foundation for the design of effective strategies to halt progress toward severe visual complications or blindness. Perhaps targeting CD1c+ DCs may be achieved by dietary (microbiome) strategies and provide relatively safe preventive strategies for non-infectious uveitis.
To conclude, we have found that peripheral blood CD1c+ DCs have a gene module linked to a CX3CR1-positive CD1c+ DC subset implicated in non-infectious uveitis.
Appendix 1
Data availability
All raw data and data scripts are available via dataverseNL: https://doi.org/10.34894/9Q0FVO and deposited in NCBI's Gene Expression Omnibus accessible through GEO Series accession numbers GSE195501 and GSE194060.
-
NCBI Gene Expression OmnibusID GSE195501. Whole transcriptome-sequencing of CD1c+ conventional type 2 dendritic cells of human non-infectious uveitis patients [Replication cohort].
-
NCBI Gene Expression OmnibusID GSE194060. Whole transcriptome-sequencing of CD1c+ conventional type 2 dendritic cells of human non-infectious uveitis patients.
-
DataverseNLData and R Scripts for: "Transcriptome network analysis implicates CX3CR1-positive type 3 dendritic cells in non-infectious uveitis".https://doi.org/10.34894/9Q0FVO
-
NCBI Gene Expression OmnibusID GSE48590. The affect of specific ablation of Runx3 from Esam splenic dendritic cells.
-
NCBI Gene Expression OmnibusID GSE119242. Trancriptional profile of WT and Notch2 cDC2s after immunization with SRBC.
-
NCBI Gene Expression OmnibusID GSE149619. Inflammatory Type 2 cDCs Acquire Features of cDC1s and Macrophages to Orchestrate Immunity to Respiratory Virus Infection.
-
NCBI Gene Expression OmnibusID GSE110577. Notch signaling facilitates in vitro generation of cross-presenting classical dendritic cells.
-
NCBI Gene Expression OmnibusID GSE178833. Intraocular dendritic cells characterize HLA-B27-associated acute anterior uveitis.
References
-
Unraveling heterogeneity in pediatric atopic dermatitis: Identification of serum biomarker based patient clustersThe Journal of Allergy and Clinical Immunology 149:125–134.https://doi.org/10.1016/j.jaci.2021.06.029
-
A look at autoimmunity and inflammation in the eyeThe Journal of Clinical Investigation 120:3073–3083.https://doi.org/10.1172/JCI42440
-
Levels of blood cd1c+ mdc1 and cd1chi mdc1 subpopulation reflect disease activity in noninfectious uveitisInvestigative Ophthalmology & Visual Science 56:346–352.https://doi.org/10.1167/iovs.14-15416
-
Elevated cd1c+ myeloid dendritic cell proportions associate with clinical activity and predict disease reactivation in noninfectious uveitisInvestigative Ophthalmology & Visual Science 57:1765–1772.https://doi.org/10.1167/iovs.15-18357
-
STAR: Ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Specific fibroblastic niches in secondary lymphoid organs orchestrate distinct notch-regulated immune responsesThe Journal of Experimental Medicine 211:2265–2279.https://doi.org/10.1084/jem.20132528
-
CD83+ CCR7+ NK cells induced by interleukin 18 by dendritic cells promote experimental autoimmune uveitisJournal of Cellular and Molecular Medicine 1:1827–1839.https://doi.org/10.1111/jcmm.14081
-
Regulation of NF-κB by TNF family cytokinesSeminars in Immunology 26:253–266.https://doi.org/10.1016/j.smim.2014.05.004
-
Genomewide association study of acute anterior uveitis identifies new susceptibility lociInvestigative Ophthalmology & Visual Science 61:3.https://doi.org/10.1167/iovs.61.6.3
-
Standardization of uveitis nomenclature for reporting clinical data. Results of the first international workshopAmerican Journal of Ophthalmology 140:509–516.https://doi.org/10.1016/j.ajo.2005.03.057
-
Functionally distinct erap1 and erap2 are a hallmark of hla-a29-(birdshot) uveitisHuman Molecular Genetics 27:4333–4343.https://doi.org/10.1093/hmg/ddy319
-
HLA-A29 and birdshot uveitis: Further down the rabbit holeFrontiers in Immunology 11:599558.https://doi.org/10.3389/fimmu.2020.599558
-
Autoimmune and autoinflammatory mechanisms in uveitisSeminars in Immunopathology 36:581–594.https://doi.org/10.1007/s00281-014-0433-9
-
Dendritic cell maturation and survival are differentially regulated by TNFR1 and TNFR2Journal of Immunology 193:4914–4923.https://doi.org/10.4049/jimmunol.1302929
-
New insights into the genetic component of non-infectious uveitis through an immunochip strategyJournal of Medical Genetics 54:38–46.https://doi.org/10.1136/jmedgenet-2016-104144
-
The pathogenic role of dendritic cells in non-infectious anterior uveitisExperimental Eye Research 173:121–128.https://doi.org/10.1016/j.exer.2018.05.008
-
Autoimmunity after alpha-interferon therapy for malignant carcinoid tumorsAnnals of Internal Medicine 115:178–183.https://doi.org/10.7326/0003-4819-115-3-178
-
RelB suppresses type i interferon signaling in dendritic cellsCellular Immunology 349:104043.https://doi.org/10.1016/j.cellimm.2020.104043
-
The possible impact of uveitis in blindness: A literature surveyThe British Journal of Ophthalmology 80:844–848.https://doi.org/10.1136/bjo.80.9.844
-
Expert opinion on the use of biological therapy in non-infectious uveitisExpert Opinion on Biological Therapy 19:477–490.https://doi.org/10.1080/14712598.2019.1595578
-
Vitamin D-induced ectodomain shedding of TNF receptor 1 as a nongenomic action: D3 vs D2 derivativesThe Journal of Steroid Biochemistry and Molecular Biology 155:18–25.https://doi.org/10.1016/j.jsbmb.2015.09.019
Article and author information
Author details
Funding
UitZicht (#2014-4)
- Jonas JW Kuiper
UitZicht (#2019-10)
- Jonas JW Kuiper
UitZicht (#2021-4)
- Jonas JW Kuiper
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Ethics
This study was conducted in compliance with the Helsinki principles. Ethical approval was requested and obtained from the Medical Ethical Research Committee in Utrecht. All patients signed written informed consent before participation (METC protocol number #14-065/M).
Copyright
© 2023, Hiddingh, Pandit, Verhagen et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 805
- views
-
- 101
- downloads
-
- 5
- 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
Here, we sequenced rearranged TCRβ and TCRα chain sequences in CD4+CD8+ double positive (DP), CD4+CD8- single positive (SP4) and CD4-CD8+ (SP8) thymocyte populations from the foetus and young adult mouse. We found that life-stage had a greater impact on TCRβ and TCRα gene segment usage than cell-type. Foetal repertoires showed bias towards 3’TRAV and 5’TRAJ rearrangements in all populations, whereas adult repertoires used more 5’TRAV gene segments, suggesting that progressive TCRα rearrangements occur less frequently in foetal DP cells. When we synchronised young adult DP thymocyte differentiation by hydrocortisone treatment the new recovering DP thymocyte population showed more foetal-like 3’TRAV and 5’TRAJ gene segment usage. In foetus we identified less influence of MHC-restriction on α-chain and β-chain combinatorial VxJ usage and CDR1xCDR2 (V region) usage in SP compared to adult, indicating weaker impact of MHC-restriction on the foetal TCR repertoire. The foetal TCRβ repertoire was less diverse, less evenly distributed, with fewer non-template insertions, and all foetal populations contained more clonotypic expansions than adult. The differences between the foetal and adult thymus TCR repertoires are consistent with the foetal thymus producing αβT-cells with properties and functions that are distinct from adult T-cells: their repertoire is less governed by MHC-restriction, with preference for particular gene segment usage, less diverse with more clonotypic expansions, and more closely encoded by genomic sequence.