Inferring ligand-receptor cellular networks from bulk and spatial transcriptomic datasets with BulkSignalR.

The study of cellular networks mediated by ligand-receptor interactions has attracted much attention recently owing to single-cell omics. However, rich collections of bulk data accompanied with clinical information exists and continue to be generated with no equivalent in single-cell so far. In parallel, spatial transcriptomic (ST) analyses represent a revolutionary tool in biology. A large number of ST projects rely on multicellular resolution, for instance the Visium™ platform, where several cells are analyzed at each location, thus producing localized bulk data. Here, we describe BulkSignalR, a R package to infer ligand-receptor networks from bulk data. BulkSignalR integrates ligand-receptor interactions with downstream pathways to estimate statistical significance. A range of visualization methods complement the statistics, including functions dedicated to spatial data. We demonstrate BulkSignalR relevance using different datasets, including new Visium liver metastasis ST data, with experimental validation of protein colocalization. A comparison with other ST packages shows the significantly higher quality of BulkSignalR inferences. BulkSignalR can be applied to any species thanks to its built-in generic ortholog mapping functionality.

Targeting the splicing isoforms of spleen tyrosine kinase affects the viability of colorectal cancer cells.

Spleen tyrosine kinase (Syk) expression have been both positively and negatively associated with tumorigenesis. Our goal was to evaluate the contribution of Syk and its two splice variants, full length Syk (L) and short isoform Syk (S), in the tumor biology of colorectal cancer cells (CRC). The analysis of Syk expression in primary human colorectal tumors, as well as the analysis of TCGA database, revealed a high Syk mRNA expression score in colorectal cancer tumors, suggesting a tumor promotor role of Syk in CRC. Our analysis showed that Syk (L) isoform is highly expressed in the majority of the tumor tissues and that it remains expressed in tumors in which global Syk expression is downregulated, suggesting the dependence of tumors to Syk (L) isoform. We also identified a small cluster of tumor tissues, which express a high proportion of Syk (S) isoform. This specific cluster is associated with overexpressed genes related to translation and mitochondria, and down regulated genes implicated in the progression of mitosis. For our functional studies, we used short hairpin RNA tools to target the expression of Syk in CRC cells bearing the activating K-Ras (G13D) mutation. Our results showed that while global Syk knock down increases cell proliferation and cell motility, Syk (L) expression silencing affects the viability and induces the apoptosis of the cells, confirming the dependence of cells on Syk (L) isoform for their survival. Finally, we report the promising potential of compound C-13, an original non-enzymatic inhibitor of Syk isolated in our group. In vitro studies showed that C-13 exerts cytotoxic effects on Syk-positive CRC cells by inhibiting their proliferation and their motility, and by inducing their apoptosis, while Syk-negative cell lines viability was not affected. Moreover, the oral and intraperitoneal administration of C-13 reduced the tumor growth of CRC DLD-1 cells xenografts in Nude mice in vivo.

Chemoradiation triggers antitumor Th1 and tissue resident memory-polarized immune responses to improve immune checkpoint inhibitors therapy.

BACKGROUND: Multiple synergistic combination approaches with cancer drugs are developed to overcome primary resistance to immunotherapy; however, the mechanistic rationale to combine chemoradiotherapy (CRT) with immune checkpoint inhibitors remains elusive. METHODS: This study described the immunological landscape of tumor microenvironment (TME) exposed to CRT. Tumor samples from patients with rectal cancer (n=43) treated with neoadjuvant CRT or radiotherapy were analyzed by nanostring and immunohistochemistry. Studies in mice were performed using three syngeneic tumors (TC1, CT26 and MC38). Tumor-bearing mice were treated either with platinum-based CRT, radiotherapy or chemotherapy. Anti-CTLA-4 and/or anti-Programmed Cell Death Receptor-1 (PD-1) therapy was used in combination with CRT. The therapy-exposed TME was screened by RNA sequencing and flow cytometry and tumor-infiltrating T lymphocyte functionality was evaluated by interferon (IFN)-γ ELIspot and intracellular cytokine staining. RESULTS: Front-to-front comparison analysis revealed the synergistic effect of CRT to establish a highly inflamed and Th1-polarized immune signature in the TME of patients and mice. In both settings, CRT-exposed TMEs were highly enriched in newly-infiltrated tumor-specific CD8+ T cells as well as tissue resident memory CD103+CD8+ T cells. In mice, CD8 T cells were involved in the antitumor response mediated by CRT and were primed by CRT-activated CD103+ dendritic cells. In the three tumor models, we showed that concurrent combination of CRT with a dual CTLA-4 and PD-1 blockade was required to achieve an optimal antitumor effect and to establish a broad and long-lasting protective antitumor T cell immunity. CONCLUSIONS: Our results highlight the ability of CRT to stimulate strong antitumor T-cell-mediated immunity and tissue resident memory T activation in TME, to foster immune checkpoint inhibitors action. These findings have implications in clinic for the design clinical trials combining chemoradiation with immunotherapy.

The immune contexture of primary central nervous system diffuse large B cell lymphoma associates with patient survival and specific cell signaling.

Rationale: Primary central nervous system diffuse large B-cell lymphoma (PCNSL) is a rare and aggressive entity that resides in an immune-privileged site. The tumor microenvironment (TME) and the disruption of the immune surveillance influence lymphoma pathogenesis and immunotherapy resistance. Despite growing knowledge on heterogeneous therapeutic responses, no comprehensive description of the PCNSL TME is available. We hence investigated the immune subtypes of PCNSL and their association with molecular signaling and survival. Methods: Analysis of PCNSL transcriptomes (sequencing, n = 20; microarrays, n = 34). Integrated correlation analysis and signaling pathway topology enabled us to infer intercellular interactions. Immunohistopathology and digital imaging were used to validate bioinformatic results. Results: Transcriptomics revealed three immune subtypes: immune-rich, poor, and intermediate. The immune-rich subtype was associated to better survival and characterized by hyper-activation of STAT3 signaling and inflammatory signaling, e.g., IFNγ and TNF-α, resembling the hot subtype described in primary testicular lymphoma and solid cancer. WNT/ß-catenin, HIPPO, and NOTCH signaling were hyper-activated in the immune-poor subtype. HLA down-modulation was clearly associated with a low or intermediate immune infiltration and the absence of T-cell activation. Moreover, HLA class I down-regulation was also correlated with worse survival with implications on immune-intermediate PCNSL that frequently feature reduced HLA expression. A ligand-receptor intercellular network revealed high expression of two immune checkpoints, i.e., CTLA-4/CD86 and TIM-3/LAGLS9. TIM-3 and galectin-9 proteins were clearly upregulated in PCNSL. Conclusion: Altogether, our study reveals that patient stratification according to immune subtypes, HLA status, and immune checkpoint molecule quantification should be considered prior to immune checkpoint inhibitor therapy. Moreover, TIM-3 protein should be considered an axis for future therapeutic development.

The molecular landscape and microenvironment of salivary duct carcinoma reveal new therapeutic opportunities.

Purpose: Salivary duct carcinoma (SDC) is a rare and aggressive salivary gland cancer subtype with poor prognosis. The mutational landscape of SDC has already been the object of several studies, however little is known regarding the functional genomics and the tumor microenvironment despite their importance in oncology. Our investigation aimed at describing both the functional genomics of SDC and the SDC microenvironment, along with their clinical relevance. Methods: RNA-sequencing (24 tumors), proteomics (17 tumors), immunohistochemistry (22 tumors), and multiplexed immunofluorescence (3 tumors) data were obtained from three different patient cohorts and analyzed by digital imaging and bioinformatics. Adjacent non-tumoral tissue from patients in two cohorts were used in transcriptomic and proteomic analyses. Results: Transcriptomic and proteomic data revealed the importance of Notch, TGF-ß, and interferon-γ signaling for all SDCs. We confirmed an overall strong desmoplastic reaction by measuring α-SMA abundance, the level of which was associated with recurrence-free survival (RFS). Two distinct immune phenotypes were observed: immune-poor SDCs (36%) and immune-infiltrated SDCs (64%). Advanced bioinformatics analysis of the transcriptomic data suggested 72 ligand-receptor interactions occurred in the microenvironment and correlated with the immune phenotype. Among these interactions, three immune checkpoints were validated by immunofluorescence, including CTLA-4/DC86 and TIM-3/galectin-9 interactions, previously unidentified in SDC. Immunofluorescence analysis also confirmed an important immunosuppressive role of macrophages and NK cells, also supported by the transcriptomic data. Conclusions: Together our data significantly increase the understanding of SDC biology and open new perspectives for SDC tumor treatment. Before applying immunotherapy, patient stratification according to the immune infiltrate should be taken into account. Immune-infiltrated SDC could benefit from immune checkpoint-targeting therapy, with novel options such as anti-CTLA-4. Macrophages or NK cells could also be targeted. The dense stroma, i.e., fibroblasts or hyaluronic acid, may also be the focus for immune-poor SDC therapies, e.g. in combination with Notch or TGF-ß inhibitors, or molecules targeting SDC mutations.

Babesia divergens glycosylphosphatidylinositols modulate blood coagulation and induce Th2-biased cytokine profiles in antigen presenting cells.

Glycosylphosphatidylinositols (GPIs) are glycolipids described as toxins of protozoan parasites due to their inflammatory properties in mammalian hosts characterized by the production of interleukin (IL)-1, IL-12 and tumor necrosis factor (TNF)-α. In the present work, we studied the cytokines produced by antigen presenting cells in response to ten different GPI species extracted from Babesia divergens, responsible for babesiosis. Interestingly, B. divergens GPIs induced the production of anti-inflammatory cytokines (IL-2, IL-5) and of the regulatory cytokine IL-10 by macrophages and dendritic cells. In contrast to all protozoan GPIs studied until now, GPIs from B. divergens did not stimulate the production of TNF-α and IL-12, leading to a unique Th1/Th2 profile. Analysis of the carbohydrate composition of the B. divergens GPIs indicated that the di-mannose structure was different from the evolutionary conserved tri-mannose structure, which might explain the particular cytokine profile they induce. Expression of major histocompatibility complex (MHC) molecules on dendritic cells and apoptosis of mouse peritoneal cells were also analysed. B. divergens GPIs did not change expression of MHC class I, but decreased expression of MHC class II at the cell surface, while GPIs slightly increased the percentages of apoptotic cells. During pathogenesis of babesiosis, the inflammation-coagulation auto-amplification loop can lead to thrombosis and the effect of GPIs on coagulation parameters was investigated. Incubation of B. divergens GPIs with rat plasma ex vivo led to increase of fibrinogen levels and to prolonged activated partial thromboplastin time, suggesting a direct modulation of the extrinsic coagulation pathway by GPIs.

Sequencing of the smallest Apicomplexan genome from the human pathogen Babesia microti.

We have sequenced the genome of the emerging human pathogen Babesia microti and compared it with that of other protozoa. B. microti has the smallest nuclear genome among all Apicomplexan parasites sequenced to date with three chromosomes encoding ∼3500 polypeptides, several of which are species specific. Genome-wide phylogenetic analyses indicate that B. microti is significantly distant from all species of Babesidae and Theileridae and defines a new clade in the phylum Apicomplexa. Furthermore, unlike all other Apicomplexa, its mitochondrial genome is circular. Genome-scale reconstruction of functional networks revealed that B. microti has the minimal metabolic requirement for intraerythrocytic protozoan parasitism. B. microti multigene families differ from those of other protozoa in both the copy number and organization. Two lateral transfer events with significant metabolic implications occurred during the evolution of this parasite. The genomic sequencing of B. microti identified several targets suitable for the development of diagnostic assays and novel therapies for human babesiosis.

Description of a xenoma-inducing microsporidian, Microgemma tincae n. sp., parasite of the teleost fish Symphodus tinca from Tunisian coasts.

A xenoma-inducing microsporidian species was found to infect the liver of the teleost fish, peacock wrasse Symphodus (Crenilabrus) tinca. Minimal estimates of the prevalence of the parasite in fishes caught along Tunisian coasts were as high as 43 % for Bizerte samples (over 2 yr) and 72% for Monastir samples (over 3 yr). Developmental stages were dispersed within a xenoma structure that was bounded only by the plasma membrane of the hypertrophic host cell. Ultrastructural features support allocation to the genus Microgemma Ralphs and Matthews, 1986. Meronts were multinucleate plasmodia and were surrounded by rough endoplasmic reticulum (RER) of the host cell. Merogonic plasmodia developed into sporogonic plasmodia, with loss of the RER interface. Sporogony was polysporoblastic. Ovocylindrical spores (3.6 x 1.2 microm) harbored a lamellar polaroplast and a polar tube that was coiled 9 times. Spore features and host specificity led us to propose a new species, Microgemma tincae. The conversion of M. tincae xenomas into well-visible cyst structures or granulomas reflected an efficient host response involving the infiltration of phagocytic cells, degradation of various parasite stages and formation of a thick fibrous wall. The small subunit rDNA gene of M. tincae was partially sequenced. Phylogenetic analysis confirms the placement within the family Tetramicriidae represented by the genera Tetramicra and Microgemma.