Impact of Genetic Polymorphisms on Human Immune Cell Gene Expression.

While many genetic variants have been associated with risk for human diseases, how these variants affect gene expression in various cell types remains largely unknown. To address this gap, the DICE (database of immune cell expression, expression quantitative trait loci [eQTLs], and epigenomics) project was established. Considering all human immune cell types and conditions studied, we identified cis-eQTLs for a total of 12,254 unique genes, which represent 61% of all protein-coding genes expressed in these cell types. Strikingly, a large fraction (41%) of these genes showed a strong cis-association with genotype only in a single cell type. We also found that biological sex is associated with major differences in immune cell gene expression in a highly cell-specific manner. These datasets will help reveal the effects of disease risk-associated genetic polymorphisms on specific immune cell types, providing mechanistic insights into how they might influence pathogenesis (https://dice-database.org).

Multi-cell type gene coexpression network analysis reveals coordinated interferon response and cross-cell type correlations in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is an incurable autoimmune disease disproportionately affecting women. A major obstacle in finding targeted therapies for SLE is its remarkable heterogeneity in clinical manifestations as well as in the involvement of distinct cell types. To identify cell-specific targets as well as cross-correlation relationships among expression programs of different cell types, we here analyze six major circulating immune cell types from SLE patient blood. Our results show that presence of an interferon response signature stratifies patients into two distinct groups (IFNneg vs. IFNpos). Comparing these two groups using differential gene expression and differential gene coexpression analysis, we prioritize a relatively small list of genes from classical monocytes including two known immune modulators: TNFSF13B/BAFF (target of belimumab, an approved therapeutic for SLE) and IL1RN (the basis of anakinra, a therapeutic for rheumatoid arthritis). We then develop a multi-cell type extension of the weighted gene coexpression network analysis (WGCNA) framework, termed mWGCNA. Applying mWGCNA to RNA-seq data from six sorted immune cell populations (15 SLE, 10 healthy donors), we identify a coexpression module with interferon-stimulated genes (ISGs) among all cell types and a cross-cell type correlation linking expression of specific T helper cell markers to B cell response as well as to TNFSF13B expression from myeloid cells, all of which in turn correlates with disease severity of IFNpos patients. Our results demonstrate the power of a hypothesis-free and data-driven approach to discover drug targets and to reveal novel cross-correlation across cell types in SLE with implications for other autoimmune diseases.

Fc-optimized NKG2D-Fc constructs induce NK cell antibody-dependent cellular cytotoxicity against breast cancer cells independently of HER2/neu expression status.

The ability of NK cells to mediate Ab-dependent cellular cytotoxicity (ADCC) largely contributes to the clinical success of antitumor Abs, including trastuzumab, which is approved for the treatment of breast cancer with HER2/neu overexpression. Notably, only ∼25% of breast cancer patients overexpress HER2/neu. Moreover, HER2/neu is expressed on healthy cells, and trastuzumab application is associated with side effects. In contrast, the ligands of the activating immunoreceptor NKG2D (NKG2DL) are selectively expressed on malignant cells. In this study, we took advantage of the tumor-associated expression of NKG2DL by using them as target Ags for NKG2D-IgG1 fusion proteins optimized by amino acid exchange S239D/I332E in their Fc part. Compared to constructs with wild-type Fc parts, fusion proteins carrying the S239D/I332E modification (NKG2D-Fc-ADCC) mediated highly enhanced degranulation, ADCC, and IFN-γ production of NK cells in response to breast cancer cells. NKG2D-Fc-ADCC substantially enhanced NK reactivity also against HER2/neu-low targets that were unaffected by trastuzumab, as both compounds mediated their immunostimulatory effects in strict dependence of target Ag expression levels. Thus, in line with the hierarchically organized potential of the various activating receptors governing NK reactivity and due to its highly increased affinity to CD16, NKG2D-Fc-ADCC potently enhances NK cell reactivity despite the inevitable reduction of activating signals upon binding to NKG2DL. Due to the tumor-restricted expression of NKG2DL, NKG2D-Fc-ADCC may constitute an attractive means for immunotherapy especially of HER2/neu-low or -negative breast cancer.

Intra-tumoral T cells in pediatric brain tumors display clonal expansion and effector properties.

Brain tumors in children are a devastating disease in a high proportion of patients. Owing to inconsistent results in clinical trials in unstratified patients, the role of immunotherapy remains unclear. We performed an in-depth survey of the single-cell transcriptomes and clonal relationship of intra-tumoral T cells from children with brain tumors. Our results demonstrate that a large fraction of T cells in the tumor tissue are clonally expanded with the potential to recognize tumor antigens. Such clonally expanded T cells display enrichment of transcripts linked to effector function, tissue residency, immune checkpoints and signatures of neoantigen-specific T cells and immunotherapy response. We identify neoantigens in pediatric brain tumors and show that neoantigen-specific T cell gene signatures are linked to better survival outcomes. Notably, among the patients in our cohort, we observe substantial heterogeneity in the degree of clonal expansion and magnitude of T cell response. Our findings suggest that characterization of intra-tumoral T cell responses may enable selection of patients for immunotherapy, an approach that requires prospective validation in clinical trials.

4-1BB ligand modulates direct and Rituximab-induced NK-cell reactivity in chronic lymphocytic leukemia.

NK cells play an important role in tumor immunosurveillance and largely contribute to the therapeutic success of anti-tumor antibodies like Rituximab. Here, we studied the role of the TNF family member 4-1BB ligand (4-1BBL) during the interaction of NK cells with chronic lymphocytic leukemia (CLL) cells. 4-1BBL was highly expressed on patient B-CLL cells in all 56 investigated cases. Signaling via 4-1BBL following interaction with 4-1BB, which was detected on NK cells of CLL patients but not healthy individuals, led to the release of immunoregulatory cytokines including TNF by CLL cells. CLL patient sera contained elevated levels of TNF and induced 4-1BB upregulation on NK cells, which in turn impaired direct and Rituximab-induced NK-cell reactivity against 4-1BBL-expressing targets. NK-cell reactivity was not only enhanced by blocking the interaction of NK cell-expressed 4-1BB with 4-1BBL expressed by CLL cells, but also by preventing 4-1BB upregulation on NK cells via neutralization of TNF in patient serum with Infliximab. Our data indicate that 4-1BBL mediates NK-cell immunosubversion in CLL, and thus might contribute to the reportedly compromised efficacy of Rituximab to induce NK-cell reactivity in the disease, and that TNF neutralization may serve to enhance the efficacy of Rituximab treatment in CLL.

Germline modifiers of the tumor immune microenvironment implicate drivers of cancer risk and immunotherapy response.

With the continued promise of immunotherapy for treating cancer, understanding how host genetics contributes to the tumor immune microenvironment (TIME) is essential to tailoring cancer screening and treatment strategies. Here, we study 1084 eQTLs affecting the TIME found through analysis of The Cancer Genome Atlas and literature curation. These TIME eQTLs are enriched in areas of active transcription, and associate with gene expression in specific immune cell subsets, such as macrophages and dendritic cells. Polygenic score models built with TIME eQTLs reproducibly stratify cancer risk, survival and immune checkpoint blockade (ICB) response across independent cohorts. To assess whether an eQTL-informed approach could reveal potential cancer immunotherapy targets, we inhibit CTSS, a gene implicated by cancer risk and ICB response-associated polygenic models; CTSS inhibition results in slowed tumor growth and extended survival in vivo. These results validate the potential of integrating germline variation and TIME characteristics for uncovering potential targets for immunotherapy.

Sequence and expression of the chicken membrane-associated phospholipases A1 alpha (LIPH) and beta (LIPI).

Cancer/testis antigens (CTA) are a heterogeneous group of antigens that are expressed preferentially in tumor cells and testis. Based on this definition the human membrane-associated phospholipase A1 beta (lipase family member I, LIPI) has been identified as CTA. The high homology of LIPI and the membrane-associated phospholipase A1 alpha (lipase family member H, LIPH) suggests that both genes are derived from a common ancestor by gene duplication. In contrast to human LIPI, human LIPH is expressed in several tissues. LIPI sequences have only been identified in mammals. Here, we describe the identification of LIPI in non-mammalian vertebrates. Based on the conserved genomic organization of LIPI and LIPH we identified sequences for both lipases in birds and fishes. In all vertebrates the LIPI locus is neighbored by a member of the RNA binding motif (RBM) family, RBM11. By sequencing of reverse transcriptase-polymerase chain reaction products we determined the sequences of LIPI and LIPH messenger RNA from broilers. We found that the sequence homology between LIPI and LIPH is much higher in non-mammalian species than in mammals. In addition, we found broad expression of LIPI in broilers, resembling the expression profile of LIPH. Our data suggest that LIPI is a CTA only in mammalian species and that the unique sequence features of the mammalian LIPI/RBM11 locus have evolved together with the CTA-like expression pattern of LIPI.

Single-cell eQTL analysis of activated T cell subsets reveals activation and cell type-dependent effects of disease-risk variants.

The impact of genetic variants on cells challenged in biologically relevant contexts has not been fully explored. Here, we activated CD4+ T cells from 89 healthy donors and performed a single-cell RNA sequencing assay with >1 million cells to examine cell type-specific and activation-dependent effects of genetic variants. Single-cell expression quantitative trait loci (sc-eQTL) analysis of 19 distinct CD4+ T cell subsets showed that the expression of over 4000 genes is significantly associated with common genetic polymorphisms and that most of these genes show their most prominent effects in specific cell types. These genes included many that encode for molecules important for activation, differentiation, and effector functions of T cells. We also found new gene associations for disease-risk variants identified from genome-wide association studies and highlighted the cell types in which their effects are most prominent. We found that biological sex has a major influence on activation-dependent gene expression in CD4+ T cell subsets. Sex-biased transcripts were significantly enriched in several pathways that are essential for the initiation and execution of effector functions by CD4+ T cells like TCR signaling, cytokines, cytokine receptors, costimulatory, apoptosis, and cell-cell adhesion pathways. Overall, this DICE (Database of Immune Cell Expression, eQTLs, and Epigenomics) subproject highlights the power of sc-eQTL studies for simultaneously exploring the activation and cell type-dependent effects of common genetic variants on gene expression (https://dice-database.org).

Azacytidine impairs NK cell reactivity while decitabine augments NK cell responsiveness toward stimulation.

Azacytidine and decitabine are approved for treatment of acute myeloid leukemias and myelodysplastic syndromes. While clinical responses are attributed to epigenetic effects and induction of apoptosis in malignant cells, these azanucleosides also affect antitumor immune responses. NK cells as components of innate immunity may confine development and progression of cancer. Numerous therapeutic strategies presently aim to reinforce NK reactivity against hematopoietic malignancies. We here comparatively analyzed the effect of the two clinically available azanucleosides and report that NK cytotoxicity and IFN-γ production are significantly impaired by pharmacological concentrations of azacytidine but enhanced by decitabine. This was not due to alterations in the target cells but caused by direct effects on NK cells depending on the chemical modifications by which azanucleosides differ from their physiological analogues. Although azacytidine impaired mRNA synthesis and induced apoptosis in NK cells, decitabine did not per se alter NK cell viability or reactivity but enhanced responsiveness to activating stimuli by inducing transcription of genes involved in NK reactivity. Tantalizingly, these effects were independent of incorporation of the azanucleosides into DNA during cell division. While azacytidine impairs NK antitumor immunity, decitabine augments NK reactivity by yet unidentified mechanisms and may thus serve well in therapeutic strategies combining its effects on malignant cells with its ability to enhance NK functions.

COVID-19 genetic risk variants are associated with expression of multiple genes in diverse immune cell types.

Common genetic polymorphisms associated with COVID-19 illness can be utilized for discovering molecular pathways and cell types driving disease pathogenesis. Given the importance of immune cells in the pathogenesis of COVID-19 illness, here we assessed the effects of COVID-19-risk variants on gene expression in a wide range of immune cell types. Transcriptome-wide association study and colocalization analysis revealed putative causal genes and the specific immune cell types where gene expression is most influenced by COVID-19-risk variants. Notable examples include OAS1 in non-classical monocytes, DTX1 in B cells, IL10RB in NK cells, CXCR6 in follicular helper T cells, CCR9 in regulatory T cells and ARL17A in TH2 cells. By analysis of transposase accessible chromatin and H3K27ac-based chromatin-interaction maps of immune cell types, we prioritized potentially functional COVID-19-risk variants. Our study highlights the potential of COVID-19 genetic risk variants to impact the function of diverse immune cell types and influence severe disease manifestations.

Promoter-interacting expression quantitative trait loci are enriched for functional genetic variants.

Expression quantitative trait loci (eQTLs) studies provide associations of genetic variants with gene expression but fall short of pinpointing functionally important eQTLs. Here, using H3K27ac HiChIP assays, we mapped eQTLs overlapping active cis-regulatory elements that interact with their target gene promoters (promoter-interacting eQTLs, pieQTLs) in five common immune cell types (Database of Immune Cell Expression, Expression quantitative trait loci and Epigenomics (DICE) cis-interactome project). This approach allowed us to identify functionally important eQTLs and show mechanisms that explain their cell-type restriction. We also devised an approach to eQTL discovery that relies on HiChIP-based promoter interaction maps as a structural framework for deciding which SNPs to test for association with gene expression, and observe ultra-long-distance pieQTLs (>1 megabase away), including several disease-risk variants. We validated the functional role of pieQTLs using reporter assays, CRISPRi, dCas9-tiling guides and Cas9-mediated base-pair editing. In this article we present a method for functional eQTL discovery and provide insights into relevance of noncoding variants for cell-specific gene regulation and for disease association beyond conventional eQTL mapping.

Neutralization of tumor-derived soluble glucocorticoid-induced TNFR-related protein ligand increases NK cell anti-tumor reactivity.

NK cell anti-tumor reactivity is governed by a balance of activating and inhibitory receptors including various TNF receptor (TNFR) family members. Here we report that human tumor cells release a soluble form of the TNF family member Glucocorticoid-Induced TNFR-Related Protein (GITR) ligand (sGITRL), which can be detected in cell culture supernatants. Tumor-derived sGITRL concentration-dependently reduced NK cell cytotoxicity and IFN-gamma production, which could be overcome by neutralization of sGITRL using a GITR-Ig fusion protein. Although sGITRL did not induce apoptosis in NK cells, it diminished nuclear localized RelB, indicating that sGITRL negatively modulates NK cell NF-kappaB activity. Furthermore, we detected substantial levels of sGITRL in sera of patients with various malignancies, but not in healthy controls. Presence of sGITRL-containing patient serum in cocultures with tumor cells significantly reduced NK cell cytotoxicity and IFN-gamma production, which could again be restored by neutralization of sGITRL. The strong correlation of tumor incidence and elevated sGITRL levels indicates that sGITRL is released from cancers in vivo, leading to impaired NK cell immunosurveillance of human tumors. Our data suggest that determination of sGITRL levels might be implemented as a tumor marker in patients, and GITRL neutralization may be used to improve immunotherapeutic strategies relying on NK cell reactivity.

COVID-19 genetic risk variants are associated with expression of multiple genes in diverse immune cell types.

Common genetic polymorphisms associated with severity of COVID-19 illness can be utilized for discovering molecular pathways and cell types driving disease pathogenesis. Here, we assessed the effects of 679 COVID-19-risk variants on gene expression in a wide-range of immune cell types. Severe COVID-19-risk variants were significantly associated with the expression of 11 protein-coding genes, and overlapped with either target gene promoter or cis -regulatory regions that interact with target promoters in the cell types where their effects are most prominent. For example, we identified that the association between variants in the 3p21.31 risk locus and the expression of CCR2 in classical monocytes is likely mediated through an active cis-regulatory region that interacted with CCR2 promoter specifically in monocytes. The expression of several other genes showed prominent genotype-dependent effects in non-classical monocytes, NK cells, B cells, or specific T cell subtypes, highlighting the potential of COVID-19 genetic risk variants to impact the function of diverse immune cell types and influence severe disease manifestations.

Single-cell transcriptomic analysis of allergen-specific T cells in allergy and asthma.

CD4+ T helper (TH) cells and regulatory T (Treg) cells that respond to common allergens play an important role in driving and dampening airway inflammation in patients with asthma. Until recently, direct, unbiased molecular analysis of allergen-reactive TH and Treg cells has not been possible. To better understand the diversity of these T cell subsets in allergy and asthma, we analyzed the single-cell transcriptome of ~50,000 house dust mite (HDM) allergen-reactive TH cells and Treg cells from asthmatics with HDM allergy and from three control groups: asthmatics without HDM allergy and nonasthmatics with and without HDM allergy. Our analyses show that HDM allergen-reactive TH and Treg cells are highly heterogeneous and certain subsets are quantitatively and qualitatively different in individuals with HDM-reactive asthma. The number of interleukin-9 (IL-9)-expressing HDM-reactive TH cells is greater in asthmatics with HDM allergy compared with nonasthmatics with HDM allergy, and this IL-9-expressing TH subset displays enhanced pathogenic properties. More HDM-reactive TH and Treg cells expressing the interferon response signature (THIFNR and TregIFNR) are present in asthmatics without HDM allergy compared with those with HDM allergy. In cells from these subsets (THIFNR and TregIFNR), expression of TNFSF10 was enriched; its product, tumor necrosis factor-related apoptosis-inducing ligand, dampens activation of TH cells. These findings suggest that the THIFNR and TregIFNR subsets may dampen allergic responses, which may help explain why only some people develop TH2 responses to nearly ubiquitous allergens.

Reduced expression of phosphatase PTPN2 promotes pathogenic conversion of Tregs in autoimmunity.

Genetic variants at the PTPN2 locus, which encodes the tyrosine phosphatase PTPN2, cause reduced gene expression and are linked to rheumatoid arthritis (RA) and other autoimmune diseases. PTPN2 inhibits signaling through the T cell and cytokine receptors, and loss of PTPN2 promotes T cell expansion and CD4- and CD8-driven autoimmunity. However, it remains unknown whether loss of PTPN2 in FoxP3+ regulatory T cells (Tregs) plays a role in autoimmunity. Here we aimed to model human autoimmune-predisposing PTPN2 variants, the presence of which results in a partial loss of PTPN2 expression, in mouse models of RA. We identified that reduced expression of Ptpn2 enhanced the severity of autoimmune arthritis in the T cell-dependent SKG mouse model and demonstrated that this phenotype was mediated through a Treg-intrinsic mechanism. Mechanistically, we found that through dephosphorylation of STAT3, PTPN2 inhibits IL-6-driven pathogenic loss of FoxP3 after Tregs have acquired RORγt expression, at a stage when chromatin accessibility for STAT3-targeted IL-17-associated transcription factors is maximized. We conclude that PTPN2 promotes FoxP3 stability in mouse RORγt+ Tregs and that loss of function of PTPN2 in Tregs contributes to the association between PTPN2 and autoimmunity.

Membrane-associated phospholipase A1 beta (LIPI) Is an Ewing tumour-associated cancer/testis antigen.

BACKGROUND: Cancer/testis antigens (CTA) represent a heterogeneous group of antigens expressed nearly exclusively in tumour cells and testis. Recently, we identified phospholipase A1 beta (a CTA also known as lipase member I, LIPI) as a gene with high expression in Ewing family tumours (EFT). In the present paper we analyzed expression of LIPI in a panel of normal tissues and tumour samples. PROCEDURE: The expression of CTA in EFT and normal tissues was analyzed by using DNA microarray datasets. Expression of LIPI in EFT, a panel of other tumour samples, and normal tissues was analyzed by using RT-PCR and quantitative RT-PCR. RESULTS: LIPI was expressed in EFT samples but not in other investigated tumour samples. Expression of LIPI in normal tissues was restricted to testis and thyroid. However, expression in these tissues was low compared with EFT. Interestingly testis as well as thyroid expressed all analyzed EFT-associated transcripts, suggesting that these tissues harbour a small cell population with molecular features of EFT. The sensitivity of the LIPI RT-PCR was similar to the sensitivity of the conventional EWSR1-FLI1 RT-PCR, suggesting that LIPI might be useful as additional diagnostic target structure. CONCLUSIONS: The human cancer/testis antigen LIPI is highly expressed in Ewing family tumours and can be easily detected by RT-PCR or quantitative RT-PCR. LIPI might be an interesting target for the development of future diagnostic tools or treatment strategies.

Precursors of human CD4+ cytotoxic T lymphocytes identified by single-cell transcriptome analysis.

CD4+ cytotoxic T lymphocytes (CD4-CTLs) have been reported to play a protective role in several viral infections. However, little is known in humans about the biology of CD4-CTL generation, their functional properties, and heterogeneity, especially in relation to other well-described CD4+ memory T cell subsets. We performed single-cell RNA sequencing in more than 9000 cells to unravel CD4-CTL heterogeneity, transcriptional profile, and clonality in humans. Single-cell differential gene expression analysis revealed a spectrum of known transcripts, including several linked to cytotoxic and costimulatory function that are expressed at higher levels in the TEMRA (effector memory T cells expressing CD45RA) subset, which is highly enriched for CD4-CTLs, compared with CD4+ T cells in the central memory (TCM) and effector memory (TEM) subsets. Simultaneous T cell antigen receptor (TCR) analysis in single cells and bulk subsets revealed that CD4-TEMRA cells show marked clonal expansion compared with TCM and TEM cells and that most of CD4-TEMRA were dengue virus (DENV)-specific in donors with previous DENV infection. The profile of CD4-TEMRA was highly heterogeneous across donors, with four distinct clusters identified by the single-cell analysis. We identified distinct clusters of CD4-CTL effector and precursor cells in the TEMRA subset; the precursor cells shared TCR clonotypes with CD4-CTL effectors and were distinguished by high expression of the interleukin-7 receptor. Our identification of a CD4-CTL precursor population may allow further investigation of how CD4-CTLs arise in humans and, thus, could provide insights into the mechanisms that may be used to generate durable and effective CD4-CTL immunity.

Prognostic relevance of HER2/neu in acute lymphoblastic leukemia and induction of NK cell reactivity against primary ALL blasts by trastuzumab.

The epidermal growth factor receptor HER2/neu is expressed on various cancers and represents a negative prognostic marker, but is also a target for the therapeutic monoclonal antibody Trastuzumab. In about 30% of cases, HER2/neu is expressed on acute lymphoblastic leukemia (ALL) cells and was proposed to be associated with a deleterious prognosis. Here we evaluated clinical data from 65 ALL patients (HER2/neu+, n = 17; HER2/neu-, n = 48) with a median follow-up of 19.4 months (range 0.6-176.5 months) and observed no association of HER2/neu expression with response to chemotherapy, disease free or overall survival. In vitro, treatment of primary ALL cells (CD20+HER2/neu+, CD20+HER2/neu- and CD20-HER2/neu-) with Rituximab and Trastuzumab led to activation of NK cells in strict dependence of the expression of the respective antigen. NK reactivity was more pronounced with Rituximab as compared to Trastuzumab, and combined application could lead to additive effects in cases where both antigens were expressed. Besides providing evidence that HER2/neu expression is no risk factor in ALL patients, our data demonstrates that HER2/neu can be a promising target for Trastuzumab therapy in the subset of ALL patients with the potential to improve disease outcome.