Single-Cell Analysis Reveals Malignant Cells Reshape the Cellular Landscape and Foster an Immunosuppressive Microenvironment of Extranodal NK/T-Cell Lymphoma.

Extranodal natural killer/T-cell lymphoma (NKTCL) is an aggressive type of lymphoma associated with Epstein-Barr virus (EBV) and characterized by heterogeneous tumor behaviors. To better understand the origins of the heterogeneity, this study utilizes single-cell RNA sequencing (scRNA-seq) analysis to profile the tumor microenvironment (TME) of NKTCL at the single-cell level. Together with in vitro and in vivo models, the study identifies a subset of LMP1+ malignant NK cells contributing to the tumorigenesis and development of heterogeneous malignant cells in NKTCL. Furthermore, malignant NK cells interact with various immunocytes via chemokines and their receptors, secrete substantial DPP4 that impairs the chemotaxis of immunocytes and regulates their infiltration. They also exhibit an immunosuppressive effect on T cells, which is further boosted by LMP1. Moreover, high transcription of EBV-encoded genes and low infiltration of tumor-associated macrophages (TAMs) are favorable prognostic indicators for NKTCL in multiple patient cohorts. This study for the first time deciphers the heterogeneous composition of NKTCL TME at single-cell resolution, highlighting the crucial role of malignant NK cells with EBV-encoded LMP1 in reshaping the cellular landscape and fostering an immunosuppressive microenvironment. These findings provide insights into understanding the pathogenic mechanisms of NKTCL and developing novel therapeutic strategies against NKTCL.

A genomic-augmented multivariate prognostic model for the survival of natural-killer/T-cell lymphoma patients from an international cohort.

With lowering costs of sequencing and genetic profiling techniques, genetic drivers can now be detected readily in tumors but current prognostic models for Natural-killer/T cell lymphoma (NKTCL) have yet to fully leverage on them for prognosticating patients. Here, we used next-generation sequencing to sequence 260 NKTCL tumors, and trained a genomic prognostic model (GPM) with the genomic mutations and survival data from this retrospective cohort of patients using LASSO Cox regression. The GPM is defined by the mutational status of 13 prognostic genes and is weakly correlated with the risk-features in International Prognostic Index (IPI), Prognostic Index for Natural-Killer cell lymphoma (PINK), and PINK-Epstein-Barr virus (PINK-E). Cox-proportional hazard multivariate regression also showed that the new GPM is independent and significant for both progression-free survival (PFS, HR: 3.73, 95% CI 2.07-6.73; p < .001) and overall survival (OS, HR: 5.23, 95% CI 2.57-10.65; p = .001) with known risk-features of these indices. When we assign an additional risk-score to samples, which are mutant for the GPM, the Harrell's C-indices of GPM-augmented IPI, PINK, and PINK-E improved significantly (p < .001, χ2 test) for both PFS and OS. Thus, we report on how genomic mutational information could steer toward better prognostication of NKTCL patients.

HLA class I genes modulate disease risk and age at onset together with DR-DQ in Chinese patients with insulin-requiring type 1 diabetes.

AIMS/HYPOTHESIS: The study aimed to investigate the effects of HLA class I genes on susceptibility to type 1 diabetes with different onset ages, in addition to the well-established effects of HLA class II genes. METHODS: A total of 361 patients with type 1 diabetes (192 patients with onset <18 years and 169 patients with onset ≥18 years) and 500 healthy control participants from China were enrolled and genotyped for the HLA-A, -B, -C, -DQA1, -DQB1 and -DRB1 genes using next-generation sequencing. RESULTS: The susceptible DR3 (ß = -0.09, p = 0.0009) and DR4-DQ8 (ß = -0.13, p = 0.0059) haplotypes were negatively associated with onset age, while the protective DR11 (ß = 0.21, p = 0.0314) and DR12 (ß = 0.27, p < 0.0001) haplotypes were positively associated with onset age. After adjustment for linkage disequilibrium with DR-DQ haplotypes, A*11:01:01 was positively associated with onset age (ß = 0.06, p = 0.0370), while the susceptible C*15:02:01 was negatively associated with onset age (ß = -0.21, p = 0.0050). The unit for ß was double square-root (fourth root) transformed years of change in onset age associated with per copy of the HLA haplotype/allele. In addition, B*46:01:01 was protective (OR 0.41, 0.46; pc [corrected for multiple comparisons] = 0.0044, 0.0040), whereas A*24:02:01 (OR 2.71, 2.25; pc = 0.0003, 0.0002) and B*54:01:01 (OR 3.96, 3.79; pc = 0.0018, 0.0004) were predisposing in both the <18 group and the ≥18 group compared with healthy control participants. In the context of DR4-DQ4, A*11:01:01 (61.29% vs 28.26%, pc = 0.0144) was increased while the predisposing A*24:02:01 (19.35% vs 47.83%, pc = 0.0403) was decreased in patients with onset ≥18 years when compared with patients with onset <18 years. CONCLUSIONS/INTERPRETATION: In addition to DR-DQ haplotypes, novel HLA class I alleles were detected to play a role in susceptibility to type 1 diabetes with different onset ages, which could improve the understanding of disease heterogeneity and has implications for the design of future studies.

Tumour heterogeneity and intercellular networks of nasopharyngeal carcinoma at single cell resolution.

The heterogeneous nature of tumour microenvironment (TME) underlying diverse treatment responses remains unclear in nasopharyngeal carcinoma (NPC). Here, we profile 176,447 cells from 10 NPC tumour-blood pairs, using single-cell transcriptome coupled with T cell receptor sequencing. Our analyses reveal 53 cell subtypes, including tumour-infiltrating CD8+ T, regulatory T (Treg), and dendritic cells (DCs), as well as malignant cells with different Epstein-Barr virus infection status. Trajectory analyses reveal exhausted CD8+ T and immune-suppressive TNFRSF4+ Treg cells in tumours might derive from peripheral CX3CR1+CD8+ T and naïve Treg cells, respectively. Moreover, we identify immune-regulatory and tolerogenic LAMP3+ DCs. Noteworthily, we observe intensive inter-cell interactions among LAMP3+ DCs, Treg, exhausted CD8+ T, and malignant cells, suggesting potential cross-talks to foster an immune-suppressive niche for the TME. Collectively, our study uncovers the heterogeneity and interacting molecules of the TME in NPC at single-cell resolution, which provide insights into the mechanisms underlying NPC progression and the development of precise therapies for NPC.

Single-cell transcriptome profiling of an adult human cell atlas of 15 major organs.

BACKGROUND: As core units of organ tissues, cells of various types play their harmonious rhythms to maintain the homeostasis of the human body. It is essential to identify the characteristics of cells in human organs and their regulatory networks for understanding the biological mechanisms related to health and disease. However, a systematic and comprehensive single-cell transcriptional profile across multiple organs of a normal human adult is missing. RESULTS: We perform single-cell transcriptomes of 84,363 cells derived from 15 tissue organs of one adult donor and generate an adult human cell atlas. The adult human cell atlas depicts 252 subtypes of cells, including major cell types such as T, B, myeloid, epithelial, and stromal cells, as well as novel COCH+ fibroblasts and FibSmo cells, each of which is distinguished by multiple marker genes and transcriptional profiles. These collectively contribute to the heterogeneity of major human organs. Moreover, T cell and B cell receptor repertoire comparisons and trajectory analyses reveal direct clonal sharing of T and B cells with various developmental states among different tissues. Furthermore, novel cell markers, transcription factors, and ligand-receptor pairs are identified with potential functional regulations in maintaining the homeostasis of human cells among tissues. CONCLUSIONS: The adult human cell atlas reveals the inter- and intra-organ heterogeneity of cell characteristics and provides a useful resource in uncovering key events during the development of human diseases in the context of the heterogeneity of cells and organs.

Germline Polymorphisms and Length of Survival of Nasopharyngeal Carcinoma: An Exome-Wide Association Study in Multiple Cohorts.

Germline polymorphisms are linked with differential survival outcomes in cancers but are not well studied in nasopharyngeal carcinoma (NPC). Here, a two-phase association study is conducted to discover germline polymorphisms that are associated with the prognosis of NPC. The discovery phase includes two consecutive hospital cohorts of patients with NPC from Southern China. Exome-wide genotypes at 246 173 single nucleotide polymorphisms (SNPs) are determined, followed by survival analysis for each SNP under Cox proportional hazard regression model. Candidate SNP is replicated in another two independent cohorts from Southern China and Singapore. Meta-analysis of all samples (n = 5553) confirms that the presence of rs1131636-T, located in the 3'-UTR of RPA1, confers an inferior overall survival (HR = 1.33, 95% CI = 1.20-1.47, P = 6.31 × 10-8). Bioinformatics and biological assays show that rs1131636 has regulatory effects on upstream RPA1. Functional studies further demonstrate that RPA1 promotes the growth, invasion, migration, and radioresistance of NPC cells. Additionally, miR-1253 is identified as a suppressor for RPA1 expression, likely through regulation of its binding affinity to rs1131636 locus. Collectively, these findings provide a promising biomarker aiding in stratifying patients with poor survival, as well as a potential drug target for NPC.

Genetic risk of extranodal natural killer T-cell lymphoma: a genome-wide association study in multiple populations.

BACKGROUND: Extranodal natural killer T-cell lymphoma (NKTCL; nasal type) is an aggressive malignancy with a particularly high prevalence in Asian and Latin American populations. Epstein-Barr virus infection has a role in the pathogenesis of NKTCL, and HLA-DPB1 variants are risk factors for the disease. We aimed to identify additional novel genetic variants affecting risk of NKTCL. METHODS: We did a genome-wide association study of NKTCL in multiple populations from east Asia. We recruited a discovery cohort of 700 cases with NKTCL and 7752 controls without NKTCL of Han Chinese ancestry from 19 centres in southern, central, and northern regions of China, and four independent replication samples including 717 cases and 12 650 controls. Three of these independent samples (451 cases and 5301 controls) were from eight centres in the same regions of southern, central, and northern China, and the fourth (266 cases and 7349 controls) was from 11 centres in Hong Kong, Taiwan, Singapore, and South Korea. All cases had primary NKTCL that was confirmed histopathologically, and matching with controls was based on geographical region and self-reported ancestry. Logistic regression analysis was done independently by geographical regions, followed by fixed-effect meta-analyses, to identify susceptibility loci. Bioinformatic approaches, including expression quantitative trait loci, binding motif and transcriptome analyses, and biological experiments were done to fine-map and explore the functional relevance of genome-wide association loci to the development of NKTCL. FINDINGS: Genetic data were gathered between Jan 1, 2008, and Jan 23, 2019. Meta-analysis of all samples (a total of 1417 cases and 20 402 controls) identified two novel loci significantly associated with NKTCL: IL18RAP on 2q12.1 (rs13015714; p=2·83 × 10-16; odds ratio 1·39 [95% CI 1·28-1·50]) and HLA-DRB1 on 6p21.3 (rs9271588; 9·35 × 10-26 1·53 [1·41-1·65]). Fine-mapping and experimental analyses showed that rs1420106 at the promoter of IL18RAP was highly correlated with rs13015714, and the rs1420106-A risk variant had an upregulatory effect on IL18RAP expression. Cell growth assays in two NKTCL cell lines (YT and SNK-6 cells) showed that knockdown of IL18RAP inhibited cell proliferation by cell cycle arrest in NKTCL cells. Haplotype association analysis showed that haplotype 47F-67I was associated with reduced risk of NKTCL, whereas 47Y-67L was associated with increased risk of NKTCL. These two positions are component parts of the peptide-binding pocket 7 (P7) of the HLA-DR heterodimer, suggesting that these alterations might account for the association at HLA-DRB1, independent of the previously reported HLA-DPB1 variants. INTERPRETATION: Our findings provide new insights into the development of NKTCL by showing the importance of inflammation and immune regulation through the IL18-IL18RAP axis and antigen presentation involving HLA-DRB1, which might help to identify potential therapeutic targets. Taken in combination with additional genetic and other risk factors, our results could potentially be used to stratify people at high risk of NKTCL for targeted prevention. FUNDING: Guangdong Innovative and Entrepreneurial Research Team Program, National Natural Science Foundation of China, National Program for Support of Top-Notch Young Professionals, Chang Jiang Scholars Program, Singapore Ministry of Health's National Medical Research Council, Tanoto Foundation, National Research Foundation Singapore, Chang Gung Memorial Hospital, Recruitment Program for Young Professionals of China, First Affiliated Hospital and Army Medical University, US National Institutes of Health, and US National Cancer Institute.

MDH2 is an RNA binding protein involved in downregulation of sodium channel Scn1a expression under seizure condition.

Voltage-gated sodium channel α-subunit type I (NaV1.1, encoded by SCN1A gene) plays a critical role in the excitability of brain. Downregulation of SCN1A expression is associated with epilepsy, a common neurological disorder characterized by recurrent seizures. Here we reveal a novel role of malate dehydrogenase 2 (MDH2) in the posttranscriptional regulation of SCN1A expression under seizure condition. We identified that MDH2 was an RNA binding protein that could bind two of the four conserved regions in the 3' UTRs of SCN1A. We further showed that knockdown of MDH2 or inactivation of MDH2 activity in HEK-293 cells increased the reporter gene expression through the 3' UTR of SCN1A, and MDH2 overexpression decreased gene expression by affecting mRNA stability. In the hippocampus of seizure mice, the upregulation of MDH2 expression contributed to the decrease of the NaV1.1 levels at posttranscriptional level. In addition, we showed that the H2O2 levels increased in the hippocampus of the seizure mice, and H2O2 could promote the binding of MDH2 to the binding sites of Scn1a gene, whereas ß-mercaptoethanol decreased the binding capability, indicating an important effect of the seizure-induced oxidation on the MDH2-mediated downregulation of Scn1a expression. Taken together, these data suggest that MDH2, functioning as an RNA-binding protein, is involved in the posttranscriptional downregulation of SCN1A expression under seizure condition.

Epigenetic Downregulation of Scn3a Expression by Valproate: a Possible Role in Its Anticonvulsant Activity.

Upregulation of sodium channel SCN3A expression in epileptic tissues is known to contribute to enhancing neuronal excitability and the development of epilepsy. Therefore, certain strategies to reduce SCN3A expression may be helpful for seizure control. Here, we reveal a novel role of valproate (VPA) in the epigenetic downregulation of Scn3a expression. We found that VPA, instead of carbamazepine (CBZ) and lamotrigine (LTG), could significantly downregulate Scn3a expression in mouse Neuro-2a cells. Luciferase assays and CpG methylation analyses showed that VPA induced the methylation at the -39C site in Scn3a promoter which decreased the promoter activity. We further showed that VPA downregulated the expression of methyl-CpG-binding domain protein 2 (MBD2) at the posttranscriptional level and knockdown of MBD2 increased Scn3a expression. In addition, we found that VPA induced the expression of fat mass and obesity-associated (FTO) protein and FTO knockdown abolished the repressive effects of VPA on MBD2 and Nav1.3 expressions. Furthermore, VPA, instead of other two anticonvulsant drugs, induced the expressions of Scn3a and Mbd2 and reduced Fto expression in the hippocampus of VPA-treated seizure mice. Taken together, this study suggests an epigenetic pathway for the VPA-induced downregulation of Scn3a expression, which provides a possible role of this pathway in the anticonvulsant action of VPA.

GAPDH-mediated posttranscriptional regulations of sodium channel Scn1a and Scn3a genes under seizure and ketogenic diet conditions.

Abnormal expressions of sodium channel SCN1A and SCN3A genes alter neural excitability that are believed to contribute to the pathogenesis of epilepsy, a long-term risk of recurrent seizures. Ketogenic diet (KD), a high-fat and low-carbohydrate treatment for difficult-to-control (refractory) epilepsy in children, has been suggested to reverse gene expression patterns. Here, we reveal a novel role of GAPDH on the posttranscriptional regulation of mouse Scn1a and Scn3a expressions under seizure and KD conditions. We show that GAPDH binds to a conserved region in the 3' UTRs of human and mouse SCN1A and SCN3A genes, which decreases and increases genes' expressions by affecting mRNA stability through SCN1A 3' UTR and SCN3A 3' UTR, respectively. In seizure mice, the upregulation and phosphorylation of GAPDH enhance its binding to the 3' UTR, which lead to downregulation of Scn1a and upregulation of Scn3a. Furthermore, administration of KD generates ß-hydroxybutyric acid which rescues the abnormal expressions of Scn1a and Scn3a by weakening the GAPDH's binding to the element. Taken together, these data suggest that GAPDH-mediated expression regulation of sodium channel genes may be associated with epilepsy and the anticonvulsant action of KD.

The SCN1A mutation database: updating information and analysis of the relationships among genotype, functional alteration, and phenotype.

Mutations in the SCN1A gene have been identified in epilepsy patients with widely variable phenotypes and modes of inheritance and in asymptomatic carriers. This raises challenges in evaluating the pathogenicity of SCN1A mutations. We systematically reviewed all SCN1A mutations and established a database containing information on functional alterations. In total, 1,257 mutations have been identified, of which 81.8% were not recurrent. There was a negative correlation between phenotype severity and missense mutation frequency. Further analyses suggested close relationships among genotype, functional alteration, and phenotype. Missense mutations located in different sodium channel regions were associated with distinct functional changes. Missense mutations in the pore region were characterized by the complete loss of function, similar to haploinsufficiency. Mutations with severe phenotypes were more frequently located in the pore region, suggesting that functional alterations are critical in evaluating pathogenicity and can be applied to patient management. A negative correlation was found between phenotype severity and familial incidence, and incomplete penetrance was associated with missense and splice site mutations, but not truncations or genomic rearrangements, suggesting clinical genetic counseling applications. Mosaic mutations with a load of 12.5-25.0% were potentially pathogenic with low penetrance, suggesting the need for future studies on less pathogenic genomic variations.