CHEX-seq detects single-cell genomic single-stranded DNA with catalytical potential.

Genomic DNA (gDNA) undergoes structural interconversion between single- and double-stranded states during transcription, DNA repair and replication, which is critical for cellular homeostasis. We describe "CHEX-seq" which identifies the single-stranded DNA (ssDNA) in situ in individual cells. CHEX-seq uses 3'-terminal blocked, light-activatable probes to prime the copying of ssDNA into complementary DNA that is sequenced, thereby reporting the genome-wide single-stranded chromatin landscape. CHEX-seq is benchmarked in human K562 cells, and its utilities are demonstrated in cultures of mouse and human brain cells as well as immunostained spatially localized neurons in brain sections. The amount of ssDNA is dynamically regulated in response to perturbation. CHEX-seq also identifies single-stranded regions of mitochondrial DNA in single cells. Surprisingly, CHEX-seq identifies single-stranded loci in mouse and human gDNA that catalyze porphyrin metalation in vitro, suggesting a catalytic activity for genomic ssDNA. We posit that endogenous DNA enzymatic activity is a function of genomic ssDNA.

HRS mediates tumor immune evasion by regulating proteostasis-associated interferon pathway activation.

By sorting receptor tyrosine kinases into endolysosomes, the endosomal sorting complexes required for transport (ESCRTs) are thought to attenuate oncogenic signaling in tumor cells. Paradoxically, ESCRT members are upregulated in tumors. Here, we show that disruption of hepatocyte growth factor-regulated tyrosine kinase substrate (HRS), a pivotal ESCRT component, inhibited tumor growth by promoting CD8+ T cell infiltration in melanoma and colon cancer mouse models. HRS ablation led to misfolded protein accumulation and triggered endoplasmic reticulum (ER) stress, resulting in the activation of the type I interferon pathway in an inositol-requiring enzyme-1α (IRE1α)/X-box binding protein 1 (XBP1)-dependent manner. HRS was upregulated in tumor cells with high tumor mutational burden (TMB). HRS expression associates with the response to PD-L1/PD-1 blockade therapy in melanoma patients with high TMB tumors. HRS ablation sensitized anti-PD-1 treatment in mouse melanoma models. Our study shows a mechanism by which tumor cells with high TMB evade immune surveillance and suggests HRS as a promising target to improve immunotherapy.

Upregulation of exosome secretion from tumor-associated macrophages plays a key role in the suppression of anti-tumor immunity.

Macrophages play a pivotal role in tumor immunity. We report that reprogramming of macrophages to tumor-associated macrophages (TAMs) promotes the secretion of exosomes. Mechanistically, increased exosome secretion is driven by MADD, which is phosphorylated by Akt upon TAM induction and activates Rab27a. TAM exosomes carry high levels of programmed death-ligand 1 (PD-L1) and potently suppress the proliferation and function of CD8+ T cells. Analysis of patient melanoma tissues indicates that TAM exosomes contribute significantly to CD8+ T cell suppression. Single-cell RNA sequencing analysis showed that exosome-related genes are highly expressed in macrophages in melanoma; TAM-specific RAB27A expression inversely correlates with CD8+ T cell infiltration. In a murine melanoma model, lipid nanoparticle delivery of small interfering RNAs (siRNAs) targeting macrophage RAB27A led to better T cell activation and sensitized tumors to anti-programmed cell death protein 1 (PD-1) treatment. Our study demonstrates tumors use TAM exosomes to combat CD8 T cells and suggests targeting TAM exosomes as a potential strategy to improve immunotherapies.

Astrocyte ethanol exposure reveals persistent and defined calcium response subtypes and associated gene signatures.

Astrocytes play a critical role in brain function, but their contribution during ethanol (EtOH) consumption remains largely understudied. In light of recent findings on the heterogeneity of astrocyte physiology and gene expression, an approach with the ability to identify subtypes and capture this heterogeneity is necessary. Here, we combined measurements of calcium signaling and gene expression to define EtOH-induced astrocyte subtypes. In the absence of a demonstrated EtOH receptor, EtOH is believed to have effects on the function of many receptors and downstream biological cascades that underlie calcium responsiveness. This mechanism of EtOH-induced calcium signaling is unknown and this study provides the first step in understanding the characteristics of cells displaying these observed responses. To characterize underlying astrocyte subtypes, we assessed the correlation between calcium signaling and astrocyte gene expression signature in response to EtOH. We found that various EtOH doses increased intracellular calcium levels in a subset of astrocytes, distinguishing three cellular response types and one nonresponsive subtype as categorized by response waveform properties. Furthermore, single-cell RNA-seq analysis of astrocytes from the different response types identified type-enriched discriminatory gene expression signatures. Combining single-cell calcium responses and gene expression analysis identified specific astrocyte subgroups among astrocyte populations defined by their response to EtOH. This result provides a basis for identifying the relationship between astrocyte susceptibility to EtOH and corresponding measurable markers of calcium signaling and gene expression, which will be useful to investigate potential subgroup-specific influences of astrocytes on the physiology and pathology of EtOH exposure in the brain.

Exosomal PD-L1 contributes to immunosuppression and is associated with anti-PD-1 response.

Tumour cells evade immune surveillance by upregulating the surface expression of programmed death-ligand 1 (PD-L1), which interacts with programmed death-1 (PD-1) receptor on T cells to elicit the immune checkpoint response1,2. Anti-PD-1 antibodies have shown remarkable promise in treating tumours, including metastatic melanoma2-4. However, the patient response rate is low4,5. A better understanding of PD-L1-mediated immune evasion is needed to predict patient response and improve treatment efficacy. Here we report that metastatic melanomas release extracellular vesicles, mostly in the form of exosomes, that carry PD-L1 on their surface. Stimulation with interferon-γ (IFN-γ) increases the amount of PD-L1 on these vesicles, which suppresses the function of CD8 T cells and facilitates tumour growth. In patients with metastatic melanoma, the level of circulating exosomal PD-L1 positively correlates with that of IFN-γ, and varies during the course of anti-PD-1 therapy. The magnitudes of the increase in circulating exosomal PD-L1 during early stages of treatment, as an indicator of the adaptive response of the tumour cells to T cell reinvigoration, stratifies clinical responders from non-responders. Our study unveils a mechanism by which tumour cells systemically suppress the immune system, and provides a rationale for the application of exosomal PD-L1 as a predictor for anti-PD-1 therapy.

Dissection of the module network implementation "LemonTree": enhancements towards applications in metagenomics and translation in autoimmune maladies.

Under the current deluge of omics, module networks distinctively emerge as methods capable of not only identifying inherently coherent groups (modules), thus reducing dimensionality, but also hypothesizing cause-effect relationships between modules and their regulators. Module networks were first designed in the transcriptomic era and further exploited in the multi-omic context to assess (for example) miRNA regulation of gene expression. Despite a number of available implementations, expansion of module networks to other omics is constrained by a limited characterization of the solutions' (modules plus regulators) accuracy and stability - an immediate need for the better characterization of molecular biology complexity in silico. We hence carefully assessed for LemonTree - a popular and open source module network implementation - the dependency of the software performances (sensitivity, specificity, false discovery rate, solutions' stability) on the input parameters and on the data quality (sample size, expression noise) based on synthetic and real data. In the process, we uncovered and fixed an issue in the code for the regulator assignment procedure. We concluded this evaluation with a table of recommended parameter settings. Finally, we applied these recommended settings to gut-intestinal metagenomic data from rheumatoid arthritis patients, to characterize the evolution of the gut-intestinal microbiome under different pharmaceutical regimens (methotrexate and prednisone) and we inferred innovative clinical recommendations with therapeutic potential, based on the computed module network.

Systemic Wound Healing Associated with local sub-Cutaneous Mechanical Stimulation.

Degeneration is a hallmark of autoimmune diseases, whose incidence grows worldwide. Current therapies attempt to control the immune response to limit degeneration, commonly promoting immunodepression. Differently, mechanical stimulation is known to trigger healing (regeneration) and it has recently been proposed locally for its therapeutic potential on severely injured areas. As the early stages of healing consist of altered intra- and inter-cellular fluxes of soluble molecules, we explored the potential of this early signal to spread, over time, beyond the stimulation district and become systemic, to impact on distributed or otherwise unreachable injured areas. We report in a model of arthritis in rats how stimulations delivered in the subcutaneous dorsal tissue result, over time, in the control and healing of the degeneration of the paws' joints, concomitantly with the systemic activation of wound healing phenomena in blood and in correlation with a more eubiotic microbiome in the gut intestinal district.

Mechanotransduction map: simulation model, molecular pathway, gene set.

MOTIVATION: Mechanotransduction--the ability to output a biochemical signal from a mechanical input--is related to the initiation and progression of a broad spectrum of molecular events. Yet, the characterization of mechanotransduction lacks some of the most basic tools as, for instance, it can hardly be recognized by enrichment analysis tools, nor could we find any pathway representation. This greatly limits computational testing and hypothesis generation on mechanotransduction biological relevance and involvement in disease or physiological mechanisms. RESULTS: We here present a molecular map of mechanotransduction, built in CellDesigner to warrant that maximum information is embedded in a compact network format. To validate the map's necessity we tested its redundancy in comparison with existing pathways, and to estimate its sufficiency, we quantified its ability to reproduce biological events with dynamic simulations, using Signaling Petri Networks. AVAILABILITY AND IMPLEMENTATION: SMBL language map is available in the Supplementary Data: core_map.xml, basic_map.xml. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Exploring the molecular causes of hepatitis B virus vaccination response: an approach with epigenomic and transcriptomic data.

BACKGROUND: Variable responses to the Hepatitis B Virus (HBV) vaccine have recently been reported as strongly dependent on genetic causes. Yet, the details on such mechanisms of action are still unknown. In parallel, altered DNA methylation states have been uncovered as important contributors to a variety of health conditions. However, methodologies for the analysis of such high-throughput data (epigenomic), especially from the computational point of view, still lack of a gold standard, mostly due to the intrinsic statistical distribution of methylomic data i.e. binomial rather than (pseudo-) normal, which characterizes the better known transcriptomic data.We present in this article our contribution to the challenge of epigenomic data analysis with application to the variable response to the Hepatitis B virus (HBV) vaccine and its most lethal degeneration: hepatocellular carcinoma (HCC). METHODS: Twenty-five infants were recruited and classified as good and non-/low- responders according to serological test results. Whole genome DNA methylation states were profiled by Illumina HumanMethylation 450 K beadchips. Data were processed through quality and dispersion filtering and with differential methylation analysis based on a combination of average methylation differences and non-parametric statistical tests. Results were finally associated to already published transcriptomics and post-transcriptomics to gain further insight. RESULTS: We highlight 2 relevant variations in poor-responders to HBV vaccination: the hypomethylation of RNF39 (Ring Finger Protein 39) and the complex biochemical alteration on SULF2 via hypermethylation, down-regulation and post-transcriptional control. CONCLUSIONS: Our approach appears to cope with the new challenges implied by methylomic data distribution to warrant a robust ranking of candidates. In particular, being RNF39 within the Major Histocompatibility Complex (MHC) class I region, its altered methylation state fits with an altered immune reaction compatible with poor responsiveness to vaccination. Additionally, despite SULF2 having been indicated as a potential target for HCC therapy, we can recommend that non-responders to HBV vaccine who develop HCC are quickly directed to other therapies, as SULF2 appears to be already under multiple molecular controls in such patients. Future research in this direction is warranted.