Analysis and interpretation of high-throughput transcriptional and chromatin accessibility data at single-cell (sc) resolution are still open challenges in the biomedical field. The existence of countless bioinformatics tools, for the different analytical steps, increases the complexity of data interpretation and the difficulty to derive biological insights. In this article, we present SCALA, a bioinformatics tool for analysis and visualization of single-cell RNA sequencing (scRNA-seq) and Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq) datasets, enabling either independent or integrative analysis of the two modalities. SCALA combines standard types of analysis by integrating multiple software packages varying from quality control to the identification of distinct cell populations and cell states. Additional analysis options enable functional enrichment, cellular trajectory inference, ligand-receptor analysis, and regulatory network reconstruction. SCALA is fully parameterizable, presenting data in tabular format and producing publication-ready visualizations. The different available analysis modules can aid biomedical researchers in exploring, analyzing, and visualizing their data without any prior experience in coding. We demonstrate the functionality of SCALA through two use-cases related to TNF-driven arthritic mice, handling both scRNA-seq and scATAC-seq datasets. SCALA is developed in R, Shiny and JavaScript and is mainly available as a standalone version, while an online service of more limited capacity can be found at http://scala.pavlopouloslab.info or https://scala.fleming.gr.

Tumor progression locus-2 (Tpl2) kinase is a major inflammatory mediator in immune cell types recently found to be genetically associated with inflammatory bowel diseases (IBDs). Here we show that Tpl2 may exert a dominant homeostatic rather than inflammatory function in the intestine mediated specifically by subepithelial intestinal myofibroblasts (IMFs). Mice with complete or IMF-specific Tpl2 ablation are highly susceptible to epithelial injury-induced colitis showing impaired compensatory proliferation in crypts and extensive ulcerations without significant changes in inflammatory responses. Following epithelial injury, IMFs sense innate or inflammatory signals and activate, via Tpl2, the cyclooxygenase-2 (Cox-2)-prostaglandin E2 (PGE2) pathway, which we show here to be essential for the epithelial homeostatic response. Exogenous PGE2 administration rescues mice with complete or IMF-specific Tpl2 ablation from defects in crypt function and susceptibility to colitis. We also show that Tpl2 expression is decreased in IMFs isolated from the inflamed ileum of IBD patients indicating that Tpl2 function in IMFs may be highly relevant to human disease. The IMF-mediated mechanism we propose also involves the IBD-associated genes IL1R1, MAPK1, and the PGE2 receptor-encoding PTGER4. Our results establish a previously unidentified myofibroblast-specific innate pathway that regulates intestinal homeostasis and may underlie IBD susceptibility in humans.

The reasons behind the clinical variability of SARS-CoV-2 infection, ranging from asymptomatic infection to lethal disease, are still unclear. We performed genome-wide transcriptional whole-blood RNA sequencing, bioinformatics analysis and PCR validation to test the hypothesis that immune response-related gene signatures reflecting baseline may differ between healthy individuals, with an equally robust antibody response, who experienced an entirely asymptomatic (n=17) versus clinical SARS-CoV-2 infection (n=15) in the past months (mean of 14 weeks). Among 12.789 protein-coding genes analysed, we identified six and nine genes with significantly decreased or increased expression, respectively, in those with prior asymptomatic infection relatively to those with clinical infection. All six genes with decreased expression (IFIT3, IFI44L, RSAD2, FOLR3, PI3, ALOX15), are involved in innate immune response while the first two are interferon-induced proteins. Among genes with increased expression six are involved in immune response (GZMH, CLEC1B, CLEC12A), viral mRNA translation (GCAT), energy metabolism (CACNA2D2) and oxidative stress response (ENC1). Notably, 8/15 differentially expressed genes are regulated by interferons. Our results suggest that subtle differences at baseline expression of innate immunity-related genes may be associated with an asymptomatic disease course in SARS-CoV-2 infection. Whether a certain gene signature predicts, or not, those who will develop a more efficient immune response upon exposure to SARS-CoV-2, with implications for prioritization for vaccination, warrant further study.

Paneth cells (PCs), a specialized secretory cell type in the small intestine, are increasingly recognized as having an essential role in host responses to microbiome and environmental stresses. Whether and how commensal and pathogenic microbes modify PC composition to modulate inflammation remain unclear. Using newly developed PC-reporter mice under conventional and gnotobiotic conditions, we determined PC transcriptomic heterogeneity in response to commensal and invasive microbes at single cell level. Infection expands the pool of CD74+ PCs, whose number correlates with auto or allogeneic inflammatory disease progressions in mice. Similar correlation was found in human inflammatory disease tissues. Infection-stimulated cytokines increase production of reactive oxygen species (ROS) and expression of a PC-specific mucosal pentraxin (Mptx2) in activated PCs. A PC-specific ablation of MyD88 reduced CD74+ PC population, thus ameliorating pathogen-induced systemic disease. A similar phenotype was also observed in mice lacking Mptx2. Thus, infection stimulates expansion of a PC subset that influences disease progression.