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.

In vitro efficacy of next-generation dihydrotriazines and biguanides against babesiosis and malaria parasites.

Babesia and Plasmodium pathogens, the causative agents of babesiosis and malaria, are vector-borne intraerythrocytic protozoan parasites, posing significant threats to both human and animal health. The widespread resistance exhibited by these pathogens to various classes of antiparasitic drugs underscores the need for the development of novel and more effective therapeutic strategies. Antifolates have long been recognized as attractive antiparasitic drugs as they target the folate pathway, which is essential for the biosynthesis of purines and pyrimidines, and thus is vital for the survival and proliferation of protozoan parasites. More efficacious and safer analogs within this class are needed to overcome challenges due to resistance to commonly used antifolates, such as pyrimethamine, and to address liabilities associated with the dihydrotriazines, WR99210 and JPC-2067. Here, we utilized an in vitro culture condition suitable for the continuous propagation of Babesia duncani, Babesia divergens, Babesia MO1, and Plasmodium falciparum in human erythrocytes to screen a library of 50 dihydrotriazines and 29 biguanides for their efficacy in vitro and compared their potency and therapeutic indices across different species and isolates. We identified nine analogs that inhibit the growth of all species, including the P. falciparum pyrimethamine-resistant strain HB3, with IC50 values below 10 nM, and display excellent in vitro therapeutic indices. These compounds hold substantial promise as lead antifolates for further development as broad-spectrum antiparasitic drugs.

FT-GPI, a highly sensitive and accurate predictor of GPI-anchored proteins, reveals the composition and evolution of the GPI proteome in Plasmodium species.

BACKGROUND: Protozoan parasites are known to attach specific and diverse group of proteins to their plasma membrane via a GPI anchor. In malaria parasites, GPI-anchored proteins (GPI-APs) have been shown to play an important role in host-pathogen interactions and a key function in host cell invasion and immune evasion. Because of their immunogenic properties, some of these proteins have been considered as malaria vaccine candidates. However, identification of all possible GPI-APs encoded by these parasites remains challenging due to their sequence diversity and limitations of the tools used for their characterization. METHODS: The FT-GPI software was developed to detect GPI-APs based on the presence of a hydrophobic helix at both ends of the premature peptide. FT-GPI was implemented in C ++and applied to study the GPI-proteome of 46 isolates of the order Haemosporida. Using the GPI proteome of Plasmodium falciparum strain 3D7 and Plasmodium vivax strain Sal-1, a heuristic method was defined to select the most sensitive and specific FT-GPI software parameters. RESULTS: FT-GPI enabled revision of the GPI-proteome of P. falciparum and P. vivax, including the identification of novel GPI-APs. Orthology- and synteny-based analyses showed that 19 of the 37 GPI-APs found in the order Haemosporida are conserved among Plasmodium species. Our analyses suggest that gene duplication and deletion events may have contributed significantly to the evolution of the GPI proteome, and its composition correlates with speciation. CONCLUSION: FT-GPI-based prediction is a useful tool for mining GPI-APs and gaining further insights into their evolution and sequence diversity. This resource may also help identify new protein candidates for the development of vaccines for malaria and other parasitic diseases.

In Vivo Large-Scale Mapping of Protein Turnover in Human Cerebrospinal Fluid.

The extraction of accurate physiological parameters from clinical samples provides a unique perspective to understand disease etiology and evolution, including under therapy. We introduce a new methodologic framework to map patient proteome dynamics in vivo, either proteome-wide or in large targeted panels. We applied it to ventricular cerebrospinal fluid (CSF) and could determine the turnover parameters of almost 200 proteins, whereas a handful were known previously. We covered a large number of neuron biology- and immune system-related proteins, including many biomarkers and drug targets. This first large data set unraveled a significant relationship between turnover and protein origin that relates to our ability to investigate organ physiology with protein-labeling strategy specifics. Our data constitute the first draft of CSF proteome dynamics as well as a repertoire of peptides for the community to design new analyses. The disclosed methods apply to other fluids or tissues provided sequential sample collection can be performed. We show that the proposed mathematical modeling applies to other analytical methods in the field.

Functional analysis of Plasmodium falciparum subpopulations associated with artemisinin resistance in Cambodia.

BACKGROUND: Plasmodium falciparum malaria is one of the most widespread parasitic infections in humans and remains a leading global health concern. Malaria elimination efforts are threatened by the emergence and spread of resistance to artemisinin-based combination therapy, the first-line treatment of malaria. Promising molecular markers and pathways associated with artemisinin drug resistance have been identified, but the underlying molecular mechanisms of resistance remains unknown. The genomic data from early period of emergence of artemisinin resistance (2008-2011) was evaluated, with aim to define k13 associated genetic background in Cambodia, the country identified as epicentre of anti-malarial drug resistance, through characterization of 167 parasite isolates using a panel of 21,257 SNPs. RESULTS: Eight subpopulations were identified suggesting a process of acquisition of artemisinin resistance consistent with an emergence-selection-diffusion model, supported by the shifting balance theory. Identification of population specific mutations facilitated the characterization of a core set of 57 background genes associated with artemisinin resistance and associated pathways. The analysis indicates that the background of artemisinin resistance was not acquired after drug pressure, rather is the result of fixation followed by selection on the daughter subpopulations derived from the ancestral population. CONCLUSIONS: Functional analysis of artemisinin resistance subpopulations illustrates the strong interplay between ubiquitination and cell division or differentiation in artemisinin resistant parasites. The relationship of these pathways with the P. falciparum resistant subpopulation and presence of drug resistance markers in addition to k13, highlights the major role of admixed parasite population in the diffusion of artemisinin resistant background. The diffusion of resistant genes in the Cambodian admixed population after selection resulted from mating of gametocytes of sensitive and resistant parasite populations.

Valine/isoleucine variants drive selective pressure in the VP1 sequence of EV-A71 enteroviruses.

BACKGROUND: In 2011-2012, Northern Vietnam experienced its first large scale hand foot and mouth disease (HFMD) epidemic. In 2011, a major HFMD epidemic was also reported in South Vietnam with fatal cases. This 2011-2012 outbreak was the first one to occur in North Vietnam providing grounds to study the etiology, origin and dynamic of the disease. We report here the analysis of the VP1 gene of strains isolated throughout North Vietnam during the 2011-2012 outbreak and before. METHODS: The VP1 gene of 106 EV-A71 isolates from North Vietnam and 2 from Central Vietnam were sequenced. Sequence alignments were analyzed at the nucleic acid and protein level. Gene polymorphism was also analyzed. A Factorial Correspondence Analysis was performed to correlate amino acid mutations with clinical parameters. RESULTS: The sequences were distributed into four phylogenetic clusters. Three clusters corresponded to the subgenogroup C4 and the last one corresponded to the subgenogroup C5. Each cluster displayed different polymorphism characteristics. Proteins were highly conserved but three sites bearing only Isoleucine (I) or Valine (V) were characterized. The isoleucine/valine variability matched the clusters. Spatiotemporal analysis of the I/V variants showed that all variants which emerged in 2011 and then in 2012 were not the same but were all present in the region prior to the 2011-2012 outbreak. Some correlation was found between certain I/V variants and ethnicity and severity. CONCLUSIONS: The 2011-2012 outbreak was not caused by an exogenous strain coming from South Vietnam or elsewhere but by strains already present and circulating at low level in North Vietnam. However, what triggered the outbreak remains unclear. A selective pressure is applied on I/V variants which matches the genetic clusters. I/V variants were shown on other viruses to correlate with pathogenicity. This should be investigated in EV-A71. I/V variants are an easy and efficient way to survey and identify circulating EV-A71 strains.

Babesia microti from humans and ticks hold a genomic signature of strong population structure in the United States.

BACKGROUND: Babesia microti is an emerging tick-borne apicomplexan parasite with increasing geographic range and incidence in the United States. The rapid expansion of B. microti into its current distribution in the northeastern USA has been due to the range expansion of the tick vector, Ixodes scapularis, upon which the causative agent is dependent for transmission to humans. RESULTS: To reconstruct the history of B. microti in the continental USA and clarify the evolutionary origin of human strains, we used multiplexed hybrid capture of 25 B. microti isolates obtained from I. scapularis and human blood. Despite low genomic variation compared with other Apicomplexa, B. microti was strongly structured into three highly differentiated genetic clusters in the northeastern USA. Bayesian analyses of the apicoplast genomes suggest that the origin of the current diversity of B. microti in northeastern USA dates back 46 thousand years with a signature of recent population expansion in the last 1000 years. Human-derived samples belonged to two rarely intermixing clusters, raising the possibility of highly divergent infectious phenotypes in humans. CONCLUSIONS: Our results validate the multiplexed hybrid capture strategy for characterizing genome-wide diversity and relatedness of B. microti from ticks and humans. We find strong population structure in B. microti samples from the Northeast indicating potential barriers to gene flow.

Subtelomere organization in the genome of the microsporidian Encephalitozoon cuniculi: patterns of repeated sequences and physicochemical signatures.

BACKGROUND: The microsporidian Encephalitozoon cuniculi is an obligate intracellular eukaryotic pathogen with a small nuclear genome (2.9 Mbp) consisting of 11 chromosomes. Although each chromosome end is known to contain a single rDNA unit, the incomplete assembly of subtelomeric regions following sequencing of the genome identified only 3 of the 22 expected rDNA units. While chromosome end assembly remains a difficult process in most eukaryotic genomes, it is of significant importance for pathogens because these regions encode factors important for virulence and host evasion. RESULTS: Here we report the first complete assembly of E. cuniculi chromosome ends, and describe a novel mosaic structure of segmental duplications (EXT repeats) in these regions. EXT repeats range in size between 3.5 and 23.8 kbp and contain four multigene families encoding membrane associated proteins. Twenty-one recombination sites were identified in the sub-terminal region of E. cuniculi chromosomes. Our analysis suggests that these sites contribute to the diversity of chromosome ends organization through Double Strand Break repair mechanisms. The region containing EXT repeats at chromosome extremities can be differentiated based on gene composition, GC content, recombination sites density and chromosome landscape. CONCLUSION: Together this study provides the complete structure of the chromosome ends of E. cuniculi GB-M1, and identifies important factors, which could play a major role in parasite diversity and host-parasite interactions. Comparison with other eukaryotic genomes suggests that terminal regions could be distinguished precisely based on gene content, genetic instability and base composition biais. The diversity of processes assciated with chromosome extremities and their biological consequences, as they are presented in the present study, emphasize the fact that great effort will be necessary in the future to characterize more carefully these regions during whole genome sequencing efforts.

A targeted immunomic approach identifies diagnostic antigens in the human pathogen Babesia microti.

BACKGROUND: Babesia microti is a protozoan parasite responsible for the majority of reported cases of human babesiosis and a major risk to the blood supply. Laboratory screening of blood donors may help prevent transfusion-transmitted babesiosis but there is no Food and Drug Administration-approved screening method yet available. Development of a sensitive, specific, and highly automated B. microti antibody assay for diagnosis of acute babesiosis and blood screening could have an important impact on decreasing the health burden of B. microti infection. STUDY DESIGN AND METHODS: Herein, we take advantage of recent advances in B. microti genomic analyses, field surveys of the reservoir host, and human studies in endemic areas to apply a targeted immunomic approach to the discovery of B. microti antigens that serve as signatures of active or past babesiosis infections. Of 19 glycosylphosphatidylinositol (GPI)-anchored protein candidates (BmGPI1-19) identified in the B. microti proteome, 17 were successfully expressed, printed on a microarray chip, and used to screen sera from uninfected and B. microti-infected mice and humans to determine immune responses that are associated with active and past infection. RESULTS: Antibody responses to various B. microti BmGPI antigens were detected and BmGPI12 was identified as the best biomarker of infection that provided high sensitivity and specificity when used in a microarray antibody assay. CONCLUSION: BmGPI12 alone or in combination with other BmGPI proteins is a promising candidate biomarker for detection of B. microti antibodies that might be useful in blood screening to prevent transfusion-transmitted babesiosis.

Plasmodium falciparum parasite population structure and gene flow associated to anti-malarial drugs resistance in Cambodia.

BACKGROUND: Western Cambodia is recognized as the epicentre of emergence of Plasmodium falciparum multi-drug resistance. The emergence of artemisinin resistance has been observed in this area since 2008-2009 and molecular signatures associated to artemisinin resistance have been characterized in k13 gene. At present, one of the major threats faced, is the possible spread of Asian artemisinin resistant parasites over the world threatening millions of people and jeopardizing malaria elimination programme efforts. To anticipate the diffusion of artemisinin resistance, the identification of the P. falciparum population structure and the gene flow among the parasite population in Cambodia are essential. METHODS: To this end, a mid-throughput PCR-LDR-FMA approach based on LUMINEX technology was developed to screen for genetic barcode in 533 blood samples collected in 2010-2011 from 16 health centres in malaria endemics areas in Cambodia. RESULTS: Based on successful typing of 282 samples, subpopulations were characterized along the borders of the country. Each 11-loci barcode provides evidence supporting allele distribution gradient related to subpopulations and gene flow. The 11-loci barcode successfully identifies recently emerging parasite subpopulations in western Cambodia that are associated with the C580Y dominant allele for artemisinin resistance in k13 gene. A subpopulation was identified in northern Cambodia that was associated to artemisinin (R539T resistant allele of k13 gene) and mefloquine resistance. CONCLUSIONS: The gene flow between these subpopulations might have driven the spread of artemisinin resistance over Cambodia.

Identification and functional analysis of the primary pantothenate transporter, PfPAT, of the human malaria parasite Plasmodium falciparum.

The human malaria parasite Plasmodium falciparum is absolutely dependent on the acquisition of host pantothenate for its development within human erythrocytes. Although the biochemical properties of this transport have been characterized, the molecular identity of the parasite-encoded pantothenate transporter remains unknown. Here we report the identification and functional characterization of the first protozoan pantothenate transporter, PfPAT, from P. falciparum. We show using cell biological, biochemical, and genetic analyses that this transporter is localized to the parasite plasma membrane and plays an essential role in parasite intraerythrocytic development. We have targeted PfPAT to the yeast plasma membrane and showed that the transporter complements the growth defect of the yeast fen2Δ pantothenate transporter-deficient mutant and mediates the entry of the fungicide drug, fenpropimorph. Our studies in P. falciparum revealed that fenpropimorph inhibits the intraerythrocytic development of both chloroquine- and pyrimethamine-resistant P. falciparum strains with potency equal or better than that of currently available pantothenate analogs. The essential function of PfPAT and its ability to deliver both pantothenate and fenpropimorph makes it an attractive target for the development and delivery of new classes of antimalarial drugs.

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.

Monitoring the Influence of Hand, Foot, and Mouth Disease: New Guidelines on Patient Care during the 2011-2012 Multiwaves and Multivariant Outbreak in Hai Phong City, Vietnam.

From 2011 to 2012, Northern Vietnam suffered its first large-scale hand, foot, and mouth disease (HFMD) epidemic. Two sets of official guidelines were issued during the outbreak to handle the HFMD crisis. The city of Hai Phong was used as a model to analyze the impact of the released guidelines. A total of 9621 HFMD cases were reported in Hai Phong city from April 2011 to December 2012. Three distinct waves of HFMD occurred. Enterovirus A71 and Coxsackievirus A16 were successively associated with the epidemics. Two periods, before and after the guidelines' release, could be distinguished and characterized by different patient patterns. The time to admission and severity changed notably. Guideline publications help the health system refocus on the 0.5-3 years age group with the highest incidence of the disease. The three waves showed different special distribution, but the main routes of infection were rivers and local secondary roads, most likely through local trade and occupational movements of people.

Modeling the Dynamic of Multiwave Diseases: The Model of Hand, Foot and Mouth Disease.

An HFMD outbreak spread over the city of Hai Phòng from summer 2011 to autumn 2012. This epidemic was chosen because it was the very first HFMD epidemic in North Vietnam, eliminating thus interferences with previous outbreaks. This epidemic displayed three separate waves. A complete dataset was collected for more than 9500 patients during this period, which enabled us to analyze this epidemic at different scales. Access to the healthcare system was crucial during this period, which was possible due to a reorganization of the system in February-March 2012. An analysis at the commune level enabled us to track the epidemic along certain communication routes. The three-waves structure reveals a wide disparity at the district level. We developed a mathematical model showing high accuracy at the adjustment of data for both the total number of cases and for the number of cases per week. As a consequence, the model was able to accurately determine the dates of the beginning and end of each wave and to show that they overlapped. Using mathematical functions associated with this model, it was possible to calculate the probability for a patient to belong to a specific wave.

Insights into the evolution, virulence and speciation of Babesia MO1 and Babesia divergens through multiomics analyses.

Babesiosis, caused by protozoan parasites of the genus Babesia, is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of new Babesia species underscores the ongoing risk of zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental changes. One such pathogen, Babesia MO1, previously implicated in severe cases of human babesiosis in the United States, was initially considered a subspecies of B. divergens, the predominant agent of human babesiosis in Europe. Here we report comparative multiomics analyses of B. divergens and B. MO1 that offer insight into their biology and evolution. Our analysis shows that despite their highly similar genomic sequences, substantial genetic and genomic divergence occurred throughout their evolution resulting in major differences in gene functions, expression and regulation, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens, and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

Multiomics analysis reveals B. MO1 as a distinct Babesia species and provides insights into its evolution and virulence.

Babesiosis, caused by protozoan parasites of the genus Babesia , is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of various Babesia species underscores the ongoing risk of new zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental shifts impacting the distribution and transmission dynamics of parasites, their vectors, and reservoir hosts. One such species, Babesia MO1, previously implicated in severe cases of human babesiosis in the midwestern United States, was initially considered closely related to B. divergens , the predominant agent of human babesiosis in Europe. Yet, uncertainties persist regarding whether these pathogens represent distinct variants of the same species or are entirely separate species. We show that although both B. MO1 and B. divergens share similar genome sizes, comprising three nuclear chromosomes, one linear mitochondrial chromosome, and one circular apicoplast chromosome, major differences exist in terms of genomic sequence divergence, gene functions, transcription profiles, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens , and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

Babesia duncani multi-omics identifies virulence factors and drug targets.

Babesiosis is a malaria-like disease in humans and animals that is caused by Babesia species, which are tick-transmitted apicomplexan pathogens. Babesia duncani causes severe to lethal infection in humans, but despite the risk that this parasite poses as an emerging pathogen, little is known about its biology, metabolic requirements or pathogenesis. Unlike other apicomplexan parasites that infect red blood cells, B. duncani can be continuously cultured in vitro in human erythrocytes and can infect mice resulting in fulminant babesiosis and death. We report comprehensive, detailed molecular, genomic, transcriptomic and epigenetic analyses to gain insights into the biology of B. duncani. We completed the assembly, 3D structure and annotation of its nuclear genome, and analysed its transcriptomic and epigenetics profiles during its asexual life cycle stages in human erythrocytes. We used RNA-seq data to produce an atlas of parasite metabolism during its intraerythrocytic life cycle. Characterization of the B. duncani genome, epigenome and transcriptome identified classes of candidate virulence factors, antigens for diagnosis of active infection and several attractive drug targets. Furthermore, metabolic reconstitutions from genome annotation and in vitro efficacy studies identified antifolates, pyrimethamine and WR-99210 as potent inhibitors of B. duncani to establish a pipeline of small molecules that could be developed as effective therapies for the treatment of human babesiosis.

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.

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.

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.