The immune-suppressive landscape in lepromatous leprosy revealed by single-cell RNA sequencing

Lepromatous leprosy (L-LEP), caused by the massive proliferation of Mycobacterium leprae primarily in macrophages, is an ideal disease model for investigating the molecular mechanism of intracellular bacteria evading or modulating host immune response. Here, we performed single-cell RNA sequencing of both skin biopsies and peripheral blood mononuclear cells (PBMCs) of L-LEP patients and healthy controls. In L-LEP lesions, we revealed remarkable upregulation of APOE expression that showed a negative correlation with the major histocompatibility complex II gene HLA-DQB2 and MIF, which encodes a pro-inflammatory and anti-microbial cytokine, in the subset of macrophages exhibiting a high expression level of LIPA. The exhaustion of CD8+ T cells featured by the high expression of TIGIT and LAG3 in L-LEP lesions was demonstrated. Moreover, remarkable enhancement of inhibitory immune receptors mediated crosstalk between skin immune cells was observed in L-LEP lesions. For PBMCs, a high expression level of APOE in the HLA-DRhighFBP1high monocyte subset and the expansion of regulatory T cells were found to be associated with L-LEP. These findings revealed the primary suppressive landscape in the L-LEP patients, providing potential targets for the intervention of intracellular bacteria caused persistent infections.

A multivariate modeling framework to quantify immune checkpoint context-dependent stimulation on T cells

Cells receive, and adjust to, various stimuli, which function as part of complex microenvironments forming their “context”. The possibility that a given context impacts the response to a given stimulus defines “context-dependency” and it explains large parts of the functional variability of physiopathological and pharmacological stimuli. Currently, there is no framework to analyze and quantify context-dependency over multiple contexts and cellular response outputs. We established an experimental system including a stimulus of interest, applied to an immune cell type in several contexts. We studied the function of OX40 ligand (OX40L) on T helper (Th) cell differentiation, in 4 molecular (Th0, Th1, Th2, and Th17) and 11 dendritic cell (DC) contexts (monocyte-derived DC and cDC2 conditions). We measured 17 Th output cytokines in 302 observations, and developed a statistical modeling strategy to quantify OX40L context-dependency. This revealed highly variable context-dependency, depending on the output cytokine and context type itself. Among molecular contexts, Th2 was the most influential on OX40L function. Among DC contexts, the DC type rather than the activating stimuli was dominant in controlling OX40L context-dependency. This work mathematically formalizes the complex determinants of OX40L functionality, and provides a unique framework to decipher and quantify the context-dependent variability of any biomolecule or drug function.

SARS-CoV-2 uses metabotropic glutamate receptor subtype 2 as an internalization factor to infect cells

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses angiotensin-converting enzyme 2 (ACE2) as a binding receptor to enter cells via clathrin-mediated endocytosis (CME). However, receptors involved in other steps of SARS-CoV-2 infection remain largely unknown. Here, we found that metabotropic glutamate receptor subtype 2 (mGluR2) is an internalization factor for SARS-CoV-2. Our results show that mGluR2 directly interacts with the SARS-CoV-2 spike protein and that knockdown of mGluR2 decreases internalization of SARS-CoV-2 but not cell binding. Further, mGluR2 is uncovered to cooperate with ACE2 to facilitate SARS-CoV-2 internalization through CME and mGluR2 knockout in mice abolished SARS-CoV-2 infection in the nasal turbinates and significantly reduced viral infection in the lungs. Notably, mGluR2 is also important for SARS-CoV spike protein- and Middle East respiratory syndrome coronavirus spike protein-mediated internalization. Thus, our study identifies a novel internalization factor used by SARS-CoV-2 and opens a new door for antiviral development against coronavirus infection.

LncRNA RP11-295G20.2 regulates hepatocellular carcinoma cell growth and autophagy by targeting PTEN to lysosomal degradation

PTEN is a crucial tumor suppressor and loss of PTEN protein is involved in various cancers. However, the detailed molecular mechanisms of PTEN loss in cancers remain elusive, especially the involvement of lncRNAs. Here, lncRNA RP11-295G20.2 is found to be significantly upregulated in hepatocellular carcinoma (HCC) and promotes the growth of liver cancer cells both in vitro and in vivo. Furthermore, RP11-295G20.2 inhibits autophagy in liver cancer cells. Interestingly, RP11-295G20.2 directly binds to the PTEN protein and leads to its degradation. RP11-295G20.2 expression is inversely correlated with PTEN protein expression in 82 TCGA/TCPA-LIHC samples. Surprisingly, RP11-295G20.2-induced PTEN degradation occurs through the lysosomal pathway instead of the proteasome pathway. RP11-295G20.2 binds to the N terminus of PTEN and facilitates the interaction of p62 with PTEN. Thus, PTEN is translocated into lysosomes and degraded. RP11-295G20.2 also influences AKT phosphorylation and forkhead box O 3a (FOXO3a) translocation into the nucleus, in turn regulating the transcription of autophagy-related genes. Collectively, RP11-295G20.2 directly binds to PTEN and enables its lysosomal degradation. This newly identified RP11-295G20.2/PTEN axis reveals an unexplored molecular mechanism regarding PTEN loss in liver cancer and might provide new therapeutic benefits for liver cancer patients.

Metagenome-genome-wide association studies reveal human genetic impact on the oral microbiome

The oral microbiota contains billions of microbial cells, which could contribute to diseases in many body sites. Challenged by eating, drinking, and dental hygiene on a daily basis, the oral microbiota is regarded as highly dynamic. Here, we report significant human genomic associations with the oral metagenome from more than 1915 individuals, for both the tongue dorsum (n = 2017) and saliva (n = 1915). We identified five genetic loci associated with oral microbiota at study-wide significance (p < 3.16 × 10−11). Four of the five associations were well replicated in an independent cohort of 1439 individuals: rs1196764 at APPL2 with Prevotella jejuni, Oribacterium uSGB 3339 and Solobacterium uSGB 315; rs3775944 at the serum uric acid transporter SLC2A9 with Oribacterium uSGB 1215, Oribacterium uSGB 489 and Lachnoanaerobaculum umeaense; rs4911713 near OR11H1 with species F0422 uSGB 392; and rs36186689 at LOC105371703 with Eggerthia. Further analyses confirmed 84% (386/455 for tongue dorsum) and 85% (391/466 for saliva) of host genome-microbiome associations including six genome-wide significant associations mutually validated between the two niches. As many of the oral microbiome-associated genetic variants lie near miRNA genes, we tentatively validated the potential of host miRNAs to modulate the growth of specific oral bacteria. Human genetics accounted for at least 10% of oral microbiome compositions between individuals. Machine learning models showed that polygenetic risk scores dominated over oral microbiome in predicting risk of dental diseases such as dental calculus and gingival bleeding. These findings indicate that human genetic differences are one explanation for a stable or recurrent oral microbiome in each individual.

Genome-wide association and functional interrogation identified a variant at 3p26.1 modulating ovarian cancer survival among Chinese women

Ovarian cancer survival varies considerably among patients, to which germline variation may also contribute in addition to mutational signatures. To identify genetic markers modulating ovarian cancer outcome, we performed a genome-wide association study in 2130 Chinese ovarian cancer patients and found a hitherto unrecognized locus at 3p26.1 to be associated with the overall survival (Pcombined = 8.90 × 10−10). Subsequent statistical fine-mapping, functional annotation, and eQTL mapping prioritized a likely casual SNP rs9311399 in the non-coding regulatory region. Mechanistically, rs9311399 altered its enhancer activity through an allele-specific transcription factor binding and a long-range interaction with the promoter of a lncRNA BHLHE40-AS1. Deletion of the rs9311399-associated enhancer resulted in expression changes in several oncogenic signaling pathway genes and a decrease in tumor growth. Thus, we have identified a novel genetic locus that is associated with ovarian cancer survival possibly through a long-range gene regulation of oncogenic pathways.

Histone lysine methacrylation is a dynamic post-translational modification regulated by HAT1 and SIRT2

Histone lysine crotonylation is a posttranslational modification with demonstrated functions in transcriptional regulation. Here we report the discovery of a new type of histone posttranslational modification, lysine methacrylation (Kmea), corresponding to a structural isomer of crotonyllysine. We validate the identity of this modification using diverse chemical approaches and further confirm the occurrence of this type of histone mark by pan specific and site-specific anti-methacryllysine antibodies. In total, we identify 27 Kmea modified histone sites in HeLa cells using affinity enrichment with a pan Kmea antibody and mass spectrometry. Subsequent biochemical studies show that histone Kmea is a dynamic mark, which is controlled by HAT1 as a methacryltransferase and SIRT2 as a de-methacrylase. Altogether, these investigations uncover a new type of enzyme-catalyzed histone modification and suggest that methacrylyl-CoA generating metabolism is part of a growing number of epigenome-associated metabolic pathways.