Reanalysis and integration of public microarray datasets reveals novel host genes modulated in leprosy

DExplore: An Online Tool for Detecting Differentially Expressed Genes from mRNA Microarray Experiments

Microarray experiments, a mainstay in gene expression analysis for nearly two decades, pose challenges due to their complexity. To address this, we introduce DExplore, a user-friendly web application enabling researchers to detect differentially expressed genes using data from NCBI's GEO. Developed with R, Shiny, and Bioconductor, DExplore integrates WebGestalt for functional enrichment analysis. It also provides visualization plots for enhanced result interpretation. With a Docker image for local execution, DExplore accommodates unpublished data. To illustrate its utility, we showcase two case studies on cancer cells treated with chemotherapeutic drugs. DExplore streamlines microarray data analysis, empowering molecular biologists to focus on genes of biological significance.

Impact of dietary fiber on gut microbiota composition, function and gut-brain-modules in healthy adults - a systematic review protocol

Background

The gut microbiota has been extensively implicated in health and disease. The functional outputs of the gut microbiota, such as microbial metabolites, are considered particularly important in this regard. Significant associations exist between alterations in the relative abundance of specific microbial taxa and mental health disorders. Dietary fiber has the potential to alter gut microbiota composition and function, modifying bacterial enzymatic function and the production of metabolites. As many taxa of microorganisms have enzymes capable of producing or degrading neurochemicals i.e. neuroactive gut brain modules, new predictive tools can be applied to existing datasets such as those harvested from dietary fiber interventions. We endeavor to perform a systematic review in order to identify studies reporting compositional gut microbiota alterations after interventions with dietary fiber in healthy individuals. We aim to also extract from the selected studies publicly available microbial genomic sequence datasets for reanalysis with a consistent bioinformatics pipeline, with the ultimate intention of identifying altered gut brain modules following dietary fiber interventions.

Methods

Interventional trials and randomized controlled studies that are originally published, including cross-over and non-crossover design and involving healthy adult humans will be included. A systematic search of PubMed/MEDLINE and EMBASE, two electronic databases, will be completed.

Discussion

Various types of dietary fiber have an impact on the gut microbiota composition, with some promoting the growth of particular taxa while others are reduced in relative abundance. Our search focuses on the impact of this food component on the microbiota of healthy individuals. Compositional gut microbial changes have been reported and our review will compile and update these observations after reanalysis of their datasets with a consistent bioinformatic pipeline. From this it may be possible to predict more detailed functional consequences in terms of neuroactive gut brain modules, of the compositional alterations in gut microbial taxa.

Leprosy-specific subsets of macrophages and Schwann cells identified by single-cell RNA-sequencing

In Mycobacterium leprae (M. leprae)-infection, inflammatory cells' subsets and dynamics as well as the interactions with Schwann cells have remained elusive. We investigated individual cells in M. leprae-inoculated nude mice by single-cell RNA-sequencing (scRNA-seq). For macrophages, we dissected two M1-like subsets and five M2-like subsets, where lipid-associated signatures were pervasive in both M1-like and M2-like subsets. There were four macrophage trajectories showing: (i) pro-inflammatory (M1), (ii) lipid metabolism-related (M2), (iii) anti-inflammatory (M2), and (iv) interferon-stimulated gene-related (M2) fates. They displayed early divergence without ever rejoining along the paths, suggesting simultaneous or continuous stimuli for macrophage activation in leprosy. The scRNA-seq predicted Schwann cell-macrophage interactions (Notch1-Jag1, Plxnb1-Sema4d interactions). An immature Schwann cell subset showing Tfap2a expression was identified, indicating Schwann cell dedifferentiation in leprosy tissues. Expressions of Notch1, Jag1, Plxnb1, Sema4d, and Tfap2a were validated in mouse or human leprosy tissues by immunohistochemistry. We identified both pro-inflammatory and inflammation-resolution signatures, where lipid-associated signatures were pervasive to the macrophages, representing leprosy-specific macrophage states for prolonged and repeated episodes of inflammation and resolution. Our study identified refined molecular states and interactions of macrophages and Schwann cells, suggesting novel insights into the pathogenesis of unhealed inflammation with neuropathy and potential therapeutic targets for leprosy.

BioTEA: Containerized Methods of Analysis for Microarray-Based Transcriptomics Data

Tens of thousands of gene expression data sets describing a variety of model organisms in many different pathophysiological conditions are currently stored in publicly available databases such as the Gene Expression Omnibus (GEO) and ArrayExpress (AE). As microarray technology is giving way to RNA-seq, it becomes strategic to develop high-level tools of analysis to preserve access to this huge amount of information through the most sophisticated methods of data preparation and processing developed over the years, while ensuring, at the same time, the reproducibility of the results. To meet this need, here we present bioTEA (biological Transcript Expression Analyzer), a novel software tool that combines ease of use with the versatility and power of an R/Bioconductor-based differential expression analysis, starting from raw data retrieval and preparation to gene annotation. BioTEA is an R-coded pipeline, wrapped in a Python-based command line interface and containerized with Docker technology. The user can choose among multiple options—including gene filtering, batch effect handling, sample pairing, statistical test type—to adapt the algorithm flow to the structure of the particular data set. All these options are saved in a single text file, which can be easily shared between different laboratories to deterministically reproduce the results. In addition, a detailed log file provides accurate information about each step of the analysis. Overall, these features make bioTEA an invaluable tool for both bioinformaticians and wet-lab biologists interested in transcriptomics. BioTEA is free and open-source.

Gene expression patterns associated with multidrug therapy in multibacillary leprosy

Multidrug therapy (MDT) has been successfully used in the treatment of leprosy. However, although patients are cured after the completion of MDT, leprosy reactions, permanent disability, and occasional relapse/reinfection are frequently observed in patients. The immune system of multibacillary patients (MB) is not able to mount an effective cellular immune response against M. leprae. Consequently, clearance of bacilli from the body is a slow process and after 12 doses of MDT not all MB patients reduce bacillary index (BI). In this context, we recruited MB patients at the uptake and after 12-month of MDT. Patients were stratified according to the level of reduction of the BI after 12 doses MDT. A reduction of at least one log in BI was necessary to be considered a responder patient. We evaluated the pattern of host gene expression in skin samples with RNA sequencing before and after MDT and between samples from patients with or without one log reduction in BI. Our results demonstrated that after 12 doses of MDT there was a reduction in genes associated with lipid metabolism, inflammatory response, and cellular immune response among responders (APOBEC3A, LGALS17A, CXCL13, CXCL9, CALHM6, and IFNG). Also, by comparing MB patients with lower BI reduction versus responder patients, we identified high expression of CDH19, TMPRSS4, PAX3, FA2H, HLA-V, FABP7, and SERPINA11 before MDT. From the most differentially expressed genes, we observed that MDT modulates pathways related to immune response and lipid metabolism in skin cells from MB patients after MDT, with higher expression of genes like CYP11A1, that are associated with cholesterol metabolism in the group with the worst response to treatment. Altogether, the data presented contribute to elucidate gene signatures and identify differentially expressed genes associated with MDT outcomes in MB patients.

Potential biomarkers and signaling pathways associated with the pathogenesis of primary salivary gland carcinoma: a bioinformatics study

Salivary gland carcinoma (SGC) is rare cancer, constituting 6% of neoplasms in the head and neck area. The most responsible genes and pathways involved in the pathology of this disorder have not been fully understood. We aimed to identify differentially expressed genes (DEGs), the most critical hub genes, transcription factors, signaling pathways, and biological processes (BPs) associated with the pathogenesis of primary SGC. The mRNA dataset GSE153283 in the Gene Expression Omnibus database was re-analyzed for determining DEGs in cancer tissue of patients with primary SGC compared to the adjacent normal tissue (adjusted p-value < 0.001; |Log2 fold change| > 1). A protein interaction map (PIM) was built, and the main modules within the network were identified and focused on the different pathways and BP analyses. The hub genes of PIM were discovered, and their associated gene regulatory network was built to determine the master regulators involved in the pathogenesis of primary SGC. A total of 137 genes were found to be differentially expressed in primary SGC. The most significant pathways and BPs that were deregulated in the primary disease condition were associated with the cell cycle and fibroblast proliferation procedures. TP53, EGF, FN1, NOTCH1, EZH2, COL1A1, SPP1, CDKN2A, WNT5A, PDGFRB, CCNB1, and H2AFX were demonstrated to be the most critical genes linked with the primary SGC. SPIB, FOXM1, and POLR2A significantly regulate all the hub genes. This study illustrated several hub genes and their master regulators that might be appropriate targets for the therapeutic aims of primary SGC.

Mycobacterium leprae Induces Neutrophilic Degranulation and Low-Density Neutrophil Generation During Erythema Nodosum Leprosum

Erythema Nodosum Leprosum (ENL) is a recurrent acute inflammatory complication of leprosy affecting up to 50% of all Borderline Lepromatous and Lepromatous Leprosy (BL/LL) patients. Although ENL is described as an immune reaction mediated by neutrophils, studies demonstrating the direct role of neutrophils in ENL are still rare. One subpopulation of low-density neutrophils (LDNs), present within the fraction of peripheral blood mononuclear cells (PBMC), has been associated with the pathogenesis and severity of diseases like sepsis, lupus, and tuberculosis. We herein analyzed LDNs and high-density neutrophils (HDNs) in terms of frequency, phenotype, and morphology. Serum levels of MMP-9 (a neutrophilic degranulation marker) were evaluated by ELISA; and LDNs were generated in vitro by stimulating healthy-donor, whole-blood cultures. PBMC layers of ENL patients presented segmented/hypersegmented cells that were morphologically compatible with neutrophils. Immunofluorescence analyses identified LDNs in ENL. Flow cytometry confirmed the elevated frequency of circulating LDNs (CD14⁻CD15⁺) in ENL patients compared to healthy donors and nonreactional Borderline Tuberculoid (BT) patients. Moreover, flow cytometry analyses revealed that ENL LDNs had a neutrophilic-activated phenotype. ENL patients under thalidomide treatment presented similar frequency of LDNs as observed before treatment but its activation status was lower. In addition, Mycobacterium leprae induced in vitro generation of LDNs in whole blood in a dose-dependent fashion; and TGF-β, an inhibitor of neutrophilic degranulation, prevented LDNs generation. MMP-9 serum levels of BL/LL patients with or without ENL correlated with LDNs frequency at the same time that ultrastructural observations of ENL LDNs showed suggestive signs of degranulation. Together, our data provide new insights into the knowledge and understanding of the pathogenesis of ENL while enriching the role of neutrophils in leprosy.

A new paradigm for leprosy diagnosis based on host gene expression

Transcriptional profiling is a powerful tool to investigate and detect human diseases. In this study, we used bulk RNA-sequencing (RNA-Seq) to compare the transcriptomes in skin lesions of leprosy patients or controls affected by other dermal conditions such as granuloma annulare, a confounder for paucibacillary leprosy. We identified five genes capable of accurately distinguishing multibacillary and paucibacillary leprosy from other skin conditions. Indoleamine 2,3-dioxygenase 1 (IDO1) expression alone was highly discriminatory, followed by TLR10, BLK, CD38, and SLAMF7, whereas the HS3ST2 and CD40LG mRNA separated multi- and paucibacillary leprosy. Finally, from the main differentially expressed genes (DEG) and enriched pathways, we conclude that paucibacillary disease is characterized by epithelioid transformation and granuloma formation, with an exacerbated cellular immune response, while multibacillary leprosy features epithelial-mesenchymal transition with phagocytic and lipid biogenesis patterns in the skin. These findings will help catalyze the development of better diagnostic tools and potential host-based therapeutic interventions. Finally, our data may help elucidate host-pathogen interplay driving disease clinical manifestations.

Neutrophil NETworking in ENL: Potential as a Putative Biomarker: Future Insights

Erythema nodosum leprosum (ENL), also known as type 2 reaction (T2R) is an immune complex mediated (type III hypersensitivity) reactional state encountered in patients with borderline lepromatous and lepromatous leprosy (BL and LL) either before, during, or after the institution of anti-leprosy treatment (ALT). The consequences of ENL may be serious, leading to permanent nerve damage and deformities, constituting a major cause of leprosy-related morbidity. The incidence of ENL is increasing with the increasing number of multibacillary cases. Although the diagnosis of ENL is not difficult to make for physicians involved in the care of leprosy patients, its management continues to be a most challenging aspect of the leprosy eradication program: the chronic and recurrent painful skin lesions, neuritis, and organ involvement necessitates prolonged treatment with prednisolone, thalidomide, and anti-inflammatory and immunosuppressive drugs, which further adds to the existing morbidity. In addition, the use of immunosuppressants like methotrexate, azathioprine, cyclosporine, or biologics carries a risk of reactivation of persisters (Mycobacterium leprae), apart from their own end-organ toxicities. Most ENL therapeutic guidelines are primarily designed for acute episodes and there is scarcity of literature on management of patients with chronic and recurrent ENL. It is difficult to predict which patients will develop chronic or recurrent ENL and plan the treatment accordingly. We need simple point-of-care or ELISA-based tests from blood or skin biopsy samples, which can help us in identifying patients who are likely to require prolonged treatment and also inform us about the prognosis of reactions so that appropriate therapy may be started and continued for better ENL control in such patients. There is a significant unmet need for research for better understanding the immunopathogenesis of, and biomarkers for, ENL to improve clinical stratification and therapeutics. In this review we will discuss the potential of neutrophils (polymorphonuclear granulocytes) as putative diagnostic and prognostic biomarkers by virtue of their universal abundance in human blood, functional versatility, phenotypic heterogeneity, metabolic plasticity, differential hierarchical cytoplasmic granule mobilization, and their ability to form NETs (neutrophil extracellular traps). We will touch upon the various aspects of neutrophil biology relevant to ENL pathophysiology in a step-wise manner. We also hypothesize about an element of metabolic reprogramming of neutrophils by M. leprae that could be investigated and exploited for biomarker discovery. In the end, a potential role for neutrophil derived exosomes as a novel biomarker for ENL will also be explored.

Blood RNA signature RISK4LEP predicts leprosy years before clinical onset

BACKGROUND

Leprosy, a chronic infectious disease caused by Mycobacterium leprae, is often late- or misdiagnosed leading to irreversible disabilities. Blood transcriptomic biomarkers that prospectively predict those who progress to leprosy (progressors) would allow early diagnosis, better treatment outcomes and facilitate interventions aimed at stopping bacterial transmission. To identify potential risk signatures of leprosy, we collected whole blood of household contacts (HC, n=5,352) of leprosy patients, including individuals who were diagnosed with leprosy 4-61 months after sample collection.

METHODS

We investigated differential gene expression (DGE) by RNA-Seq between progressors before presence of symptoms (n=40) and HC (n=40), as well as longitudinal DGE within each progressor. A prospective leprosy signature was identified using a machine learning approach (Random Forest) and validated using reverse transcription quantitative PCR (RT-qPCR).

FINDINGS

Although no significant intra-individual longitudinal variation within leprosy progressors was identified, 1,613 genes were differentially expressed in progressors before diagnosis compared to HC. We identified a 13-gene prospective risk signature with an Area Under the Curve (AUC) of 95.2%. Validation of this RNA-Seq signature in an additional set of progressors (n=43) and HC (n=43) by RT-qPCR, resulted in a final 4-gene signature, designated RISK4LEP (MT-ND2, REX1BD, TPGS1, UBC) (AUC=86.4%).

INTERPRETATION

This study identifies for the first time a prospective transcriptional risk signature in blood predicting development of leprosy 4 to 61 months before clinical diagnosis. Assessment of this signature in contacts of leprosy patients can function as an adjunct diagnostic tool to target implementation of interventions to restrain leprosy development.

FUNDING

This study was supported by R2STOP Research grant, the Order of Malta-Grants-for-Leprosy-Research, the Q.M. Gastmann-Wichers Foundation and the Leprosy Research Initiative (LRI) together with the Turing Foundation (ILEP# 702.02.73 and # 703.15.07).

Reduction of host cell mitochondrial activity as Mycobacterium leprae's strategy to evade host innate immunity

Leprosy is a much-feared incapacitating infectious disease caused by Mycobacterium leprae or M lepromatosis, annually affecting roughly 200,000 people worldwide. During host-pathogen interaction, M leprae subverts the immune response, leading to development of disease. Throughout the last few decades, the impact of energy metabolism on the control of intracellular pathogens and leukocytic differentiation has become more evident. Mitochondria play a key role in regulating newly-discovered immune signaling pathways by controlling redox metabolism and the flow of energy besides activating inflammasome, xenophagy, and apoptosis. Likewise, this organelle, whose origin is probably an alphaproteobacterium, directly controls the intracellular pathogens attempting to invade its niche, a feature conquered at the expense of billions of years of coevolution. In the present review, we discuss the role of reduced host cell mitochondrial activity during M leprae infection and the consequential fates of M leprae and host innate immunity. Conceivably, inhibition of mitochondrial energy metabolism emerges as an overlooked and novel mechanism developed by M leprae to evade xenophagy and the host immune response.

Large-Scale Gene Expression Signatures Reveal a Microbicidal Pattern of Activation in Mycobacterium leprae-Infected Monocyte-Derived Macrophages With Low Multiplicity of Infection

Leprosy is a disease with a clinical spectrum of presentations that is also manifested in diverse histological features. At one pole, lepromatous lesions (L-pole) have phagocytic foamy macrophages heavily parasitized with freely multiplying intracellular Mycobacterium leprae. At the other pole, the presence of epithelioid giant cells and granulomatous formation in tuberculoid lesions (T-pole) lead to the control of M. leprae replication and the containment of its spread. The mechanism that triggers this polarization is unknown, but macrophages are central in this process. Over the past few years, leprosy has been studied using large scale techniques to shed light on the basic pathways that, upon infection, rewire the host cellular metabolism and gene expression. M. leprae is particularly peculiar as it invades Schwann cells in the nerves, reprogramming their gene expression leading to a stem-like cell phenotype. This modulatory behavior exerted by M. leprae is also observed in skin macrophages. Here, we used live M. leprae to infect (10:1 multiplicity of infection) monocyte-derived macrophages (MDMs) for 48 h and analyzed the whole gene expression profile using microarrays. In this model, we observe an intense upregulation of genes consistent with a cellular immune response, with enriched pathways including peptide and protein secretion, leukocyte activation, inflammation, and cellular divalent inorganic cation homeostasis. Among the most differentially expressed genes (DEGs) are CCL5/RANTES and CYP27B1, and several members of the metallothionein and metalloproteinase families. This is consistent with a proinflammatory state that would resemble macrophage rewiring toward granulomatous formation observed at the T-pole. Furthermore, a comparison with a dataset retrieved from the Gene Expression Omnibus of M. leprae-infected Schwann cells (MOI 100:1) showed that the patterns among the DEGs are highly distinct, as the Schwann cells under these conditions had a scavenging and phagocytic gene profile similar to M2-like macrophages, with enriched pathways rearrangements in the cytoskeleton, lipid and cholesterol metabolism and upregulated genes including MVK, MSMO1, and LACC1/FAMIN. In summary, macrophages may have a central role in defining the paradigmatic cellular (T-pole) vs. humoral (L-pole) responses and it is likely that the multiplicity of infection and genetic polymorphisms in key genes are gearing this polarization.

In search of biomarkers for leprosy by unraveling the host immune response to Mycobacterium leprae

Mycobacterium leprae, the causative agent of leprosy, is still actively transmitted in endemic areas reflected by the fairly stable number of new cases detected each year. Recognizing the signs and symptoms of leprosy is challenging, especially at an early stage. Improved diagnostic tools, based on sensitive and specific biomarkers, that facilitate diagnosis of leprosy are therefore urgently needed. In this review, we address the challenges that leprosy biomarker research is facing by reviewing cell types reported to be involved in host immunity to M leprae. These cell types can be associated with different possible fates of M leprae infection being either protective immunity, or pathogenic immune responses inducing nerve damage. Unraveling these responses will facilitate the search for biomarkers. Implications for further studies to disentangle the complex interplay between host responses that lead to leprosy disease are discussed, providing leads for the identification of new biomarkers to improve leprosy diagnostics.

COVID-19, leprosy, and neutrophils

Coronavirus Disease 2019 (COVID-19), a disease caused by the betacoronavirus Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), has only recently emerged, while Mycobacterium leprae, the etiological agent of leprosy, has endured for more than 2,000 years. As soon as the initial reports of COVID-19 became public, several entities, including the Brazilian Leprosy Society, warned about the possible impact of COVID-19 on leprosy patients. It has been verified that COVID-19 carriers can be either asymptomatic or present varying degrees of severe respiratory failure in association with cytokine storm and death, among other diseases. Severe COVID-19 patients show increased numbers of neutrophils and serum neutrophil extracellular trap (NET) markers, in addition to alterations in the neutrophil-to-lymphocyte ratio (NLR). The absence of antiviral drugs and the speed of COVID-19 transmission have had a major impact on public health systems worldwide, leading to the almost total collapse of many national and local healthcare services. Leprosy, an infectious neurological and dermatological illness, is widely considered to be the most frequent cause of physical disabilities globally. The chronic clinical course of the disease may be interrupted by acute inflammatory episodes, named leprosy reactions. These serious immunological complications, characterized by cytokine storms, are responsible for amplifying peripheral nerve damage. From 30% to 40% of all multibacillary leprosy (MB) patients experience erythema nodosum leprosum (ENL), a neutrophilic immune-mediated condition. ENL patients often present these same COVID-19-like symptoms, including high levels of serum NET markers, altered NLR, and neutrophilia. Moreover, the consequences of a M. leprae-SARS-CoV-2 coinfection have yet to be fully investigated. The goal of the present viewpoint is to describe some of the similarities that may be found between COVID-19 and leprosy disease in the context of neutrophilic biology.