Whole blood transcriptomics reveals the enrichment of neutrophil activation pathways during erythema nodosum leprosum reaction.

Introduction: Patients with the multibacillary form of leprosy can develop reactional episodes of acute inflammation, known as erythema nodosum leprosum (ENL), which are characterized by the appearance of painful cutaneous nodules and systemic symptoms. Neutrophils have been recognized to play a role in the pathogenesis of ENL, and recent global transcriptomic analysis revealed neutrophil-related processes as a signature of ENL skin lesions. Methods: In this study, we expanded this analysis to the blood compartment, comparing whole blood transcriptomics of patients with non-reactional lepromatous leprosy at diagnosis (LL, n=7) and patients with ENL before administration of anti-reactional treatment (ENL, n=15). Furthermore, a follow-up study was performed with patients experiencing an ENL episode at the time of diagnosis and after 7 days of thalidomide treatment (THAL, n=10). Validation in an independent cohort (ENL=8; LL=7) was performed by RT-qPCR. Results: An enrichment of neutrophil activation and degranulation-related genes was observed in the ENL group, with the gene for the neutrophil activation marker CD177 being the most enriched gene of ENL episode when compared to its expression in the LL group. A more pro-inflammatory transcriptome was also observed, with increased expression of genes related to innate immunity. Validation in an independent cohort indicated that S100A8 expression could discriminate ENL from LL. Supernatants of blood cells stimulated in vitro with Mycobacterium leprae sonicate showed higher levels of CD177 compared to the level of untreated cells, indicating that the leprosy bacillus can activate neutrophils expressing CD177. Of note, suggestive higher CD177 protein levels were found in the sera of patients with severe/moderate ENL episodes when compared with patients with mild episodes and LL patients, highlighting CD177 as a potential systemic marker of ENL severity that deserves future confirmation. Furthermore, a follow-up study was performed with patients at the time of ENL diagnosis and after 7 days of thalidomide treatment (THAL, n=10). Enrichment of neutrophil pathways was sustained in the transcriptomic profile of patients undergoing treatment; however, important immune targets that might be relevant to the effect of thalidomide at a systemic level, particularly NLRP6 and IL5RA, were revealed. Discussion: In conclusion, our study reinforces the key role played by neutrophils in ENL pathogenesis and shed lights on potential diagnostic candidates and novel therapeutic targets that could benefit patients with leprosy.

MicroRNAs correlate with bacillary index and genes associated to cell death processes in leprosy.

Mycobacterium leprae infects skin and peripheral nerves causing a broad of clinical forms. MicroRNAs (miRNAs) control immune mechanisms such as apoptosis, autophagy as well as to target genes leading to abnormal proliferation, metastasis, and invasion of cells. Herein we evaluated miRNAs expression for leprosy phenotypes in biopsies obtained from patients with and without reactions. We also correlated those miRNAs with both, bacillary index (BI) and genes involved in the micobacteria elimination process. Our results show a significant increase in the miR-125a-3p expression in paucibacillary (PB) patients vs multibacillary (MB) subjects (p = 0.007) and vs reversal reactions (RR) (p = 0.005), respectively. Likewise, there was a higher expression of miR-125a-3p in patients with erythema nodosum leprosum (ENL) vs MB without reactions (p = 0.002). Furthermore, there was a positive correlation between miR-125a-3p, miR-146b-5p and miR-132-5p expression and BI in patients with RR and ENL. These miRNAS were also correlated with genes such as ATG12 (miR-125a-3p), TNFRSF10A (miR-146b-5p), PARK2, CFLAR and STX7 (miR-132-5p). All together we underpin a role for these miRNAs in leprosy pathogenesis, implicating mechanisms such as apoptosis and autophagy in skin. The miR-125a-3p might have a distinct role associated with PB phenotype and ENL in MB patients.

MicroRNAs regulating macrophages infected with Leishmania L. (V.) Braziliensis isolated from different clinical forms of American tegumentary leishmaniasis.

Background: Leishmaniasis is an infectious disease caused by protozoa of the genus Leishmania. There are still no vaccines, and therapeutic options are limited, indicating the constant need to understand the fine mechanisms of its pathophysiology. An approach that has been explored in leishmaniasis is the participation of microRNAs (miRNAs), a class of small non-coding RNAs that act, in most cases, to repress gene expression. miRNAs play a role in the complex and plastic interaction between the host and pathogens, either as part of the host's immune response to neutralize infection or as a molecular strategy employed by the pathogen to modulate host pathways to its own benefit. Methods: Monocyte-derived macrophages from healthy subjects were infected with isolates of three clinical forms of L. braziliensis: cutaneous (CL), mucosal (ML), and disseminated (DL) leishmaniasis. We compared the expression of miRNAs that take part in the TLR/NFkB pathways. Correlations with parasite load as well as immune parameters were analyzed. Results: miRNAs -103a-3p, -21-3p, 125a-3p -155-5p, -146a-5p, -132- 5p, and -147a were differentially expressed in the metastatic ML and DL forms, and there was a direct correlation between miRNAs -103a-3p, -21-3p, -155-5p, -146a-5p, -132-5p, and -9-3p and parasite load with ML and DL isolates. We also found a correlation between the expression of miR-21-3p and miR-146a-5p with the antiapoptotic gene BCL2 and the increase of viable cells, whereas miR-147a was indirectly correlated with CXCL-9 levels. Conclusion: The expression of miRNAs is strongly correlated with the parasite load and the inflammatory response, suggesting the participation of these molecules in the pathogenesis of the different clinical forms of L. braziliensis.

Reanalysis and validation of the transcriptional pleural fluid signature in pleural tuberculosis.

Introduction: Pleural tuberculosis (PlTB), the most common site of extrapulmonary TB, is characterized by a paucibacillary nature and a compartmentalized inflammatory response in the pleural cavity, both of which make diagnosis and management extremely challenging. Although transcriptional signatures for pulmonary TB have already been described, data obtained by using this approach for extrapulmonary tuberculosis and, specifically, for pleural tuberculosis are scarce and heterogeneous. In the present study, a set of candidate genes previously described in pulmonary TB was evaluated to identify and validate a transcriptional signature in clinical samples from a Brazilian cohort of PlTB patients and those with other exudative causes of pleural effusion. Methods: As a first step, target genes were selected by a random forest algorithm with recursive feature elimination (RFE) from public microarray datasets. Then, peripheral blood (PB) and pleural fluid (PF) samples from recruited patients presenting exudative pleural effusion were collected during the thoracentesis procedure. Transcriptional analysis of the selected top 10 genes was performed by quantitative RT-PCR (RT-qPCR). Results: Reanalysis of the public datasets identified a set of candidate genes (CARD17, BHLHE40, FCGR1A, BATF2, STAT1, BTN3A1, ANKRD22, C1QB, GBP2, and SEPTIN4) that demonstrated a global accuracy of 89.5% in discriminating pulmonary TB cases from other respiratory diseases. Our validation cohort consisted of PlTB (n = 35) patients and non-TB (n = 34) ones. The gene expressions of CARD17, GBP2, and C1QB in PF at diagnosis were significantly different between the two (PlTB and non-TB) groups (p < 0.0001). It was observed that the gene expressions of CARD17 and GBP2 were higher in PlTB PF than in non-TB patients. C1QB showed the opposite behavior, being higher in the non-TB PF. After anti-TB therapy, however, GBP2 gene expression was significantly reduced in PlTB patients (p < 0.001). Finally, the accuracy of the three above-cited highlighted genes in the PF was analyzed, showing AUCs of 91%, 90%, and 85%, respectively. GBP2 was above 80% (sensitivity = 0.89/specificity = 0.81), and CARD17 showed significant specificity (Se = 0.69/Sp = 0.95) in its capacity to discriminate the groups. Conclusion: CARD17, GBP2, and C1QB showed promise in discriminating PlTB from other causes of exudative pleural effusion by providing accurate diagnoses, thus accelerating the initiation of anti-TB therapy.

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.

The Type I Interferon Pathway Is Upregulated in the Cutaneous Lesions and Blood of Multibacillary Leprosy Patients With Erythema Nodosum Leprosum.

In leprosy patients, acute inflammatory episodes, known as erythema nodosum leprosum (ENL), are responsible for high morbidity and tissue damage that occur during the course of Mycobacterium leprae infection. In a previous study, we showed evidence implicating DNA-sensing via TLR9 as an important inflammatory pathway in ENL. A likely important consequence of TLR9 pathway activation is the production of type I interferons (IFN-I) by plasmacytoid dendritic cells (pDCs), also implicated in the pathogenesis of several chronic inflammatory diseases. In this study, we investigated whether the IFN-I pathway is activated during ENL. Blood samples and skin lesions from multibacillary patients diagnosed with ENL were collected and the expression of genes of the IFN-I pathway and interferon-stimulated genes were compared with samples collected from non-reactional multibacillary (NR) patients. Whole blood RNAseq analysis suggested higher activation of the IFN-I pathway in ENL patients, confirmed by RT-qPCR. Likewise, significantly higher mRNA levels of IFN-I-related genes were detected in ENL skin biopsies when compared to NR patient lesions. During thalidomide administration, the drug of choice for ENL treatment, a decrease in the mRNA and protein levels of some of these genes both in the skin and blood was observed. Indeed, in vitro assays showed that thalidomide was able to block the secretion of IFN-I by peripheral blood mononuclear cells in response to M. leprae sonicate or CpG-A, a TLR9 ligand. Finally, the decreased frequencies of peripheral pDCs in ENL patients, along with the higher TLR9 expression in ENL pDCs and the enrichment of CD123+ cells in ENL skin lesions, suggest the involvement of these cells as IFN-I producers in this type of reaction. Taken together, our data point to the involvement of the pDC/type I IFN pathway in the pathogenesis of ENL, opening new avenues in identifying biomarkers for early diagnosis and new therapeutic targets for the better management of this reactional episode.

Tissue Specific Dual RNA-Seq Defines Host-Parasite Interplay in Murine Visceral Leishmaniasis Caused by Leishmania donovani and Leishmania infantum.

Visceral leishmaniasis is associated with hepato-splenomegaly and altered immune and hematological parameters in both preclinical animal models and humans. We studied mouse experimental visceral leishmaniasis caused by Leishmania infantum and Leishmania donovani in BALB/c mice using dual RNA-seq to investigate the transcriptional response of host and parasite in liver and spleen. We identified only 4 species-specific parasite expressed genes (SSPEGs; log2FC >1, FDR <0.05) in the infected spleen, and none in the infected liver. For the host transcriptome, we found 789 differentially expressed genes (DEGs; log2FC >1, FDR <0.05) in the spleen that were common to both infections, with IFNγ signaling and complement and coagulation cascade pathways highly enriched, and an additional 286 and 186 DEGs that were selective to L. donovani and L. infantum infection, respectively. Among those, there were network interactions between genes of amino acid metabolism and PPAR signaling in L. donovani infection and increased IL1ß and positive regulation of fatty acid transport in L. infantum infection, although no pathway enrichment was observed. In the liver, there were 1,939 DEGs in mice infected with either L. infantum or L. donovani in comparison to uninfected mice, and the most enriched pathways were IFNγ signaling, neutrophil mediated immunity, complement and coagulation, cytokine-chemokine responses, and hemostasis. Additionally, 221 DEGs were selective in L. donovani and 429 DEGs in L. infantum infections. These data show that the host response for these two visceral leishmaniasis infection models is broadly similar, and ∼10% of host DEGs vary in infections with either parasite species. IMPORTANCE Visceral leishmaniasis (VL) is caused by two species of Leishmania parasites, L. donovani in the Old World and L. infantum in the New World and countries bordering the Mediterranean. Although cardinal features such as hepato-splenomegaly and alterations in blood and immune function are evident, clinical presentation may vary by geography, with for example severe bleeding often associated with VL in Brazil. Although animal models of both L. donovani and L. infantum have been widely used to study disease pathogenesis, a direct side-by-side comparison of how these parasites species impact the infected host and/or how they might respond to the stresses of mammalian infection has not been previously reported. Identifying common and distinct pathways to pathogenesis will be important to ensure that new therapeutic or prophylactic approaches will be applicable across all forms of VL.

Development and validation of a multiplex real-time qPCR assay using GMP-grade reagents for leprosy diagnosis.

Leprosy is a chronic dermato-neurological disease caused by Mycobacterium leprae, an obligate intracellular bacterium. Timely detection is a challenge in leprosy diagnosis, relying on clinical examination and trained health professionals. Furthermore, adequate care and transmission control depend on early and reliable pathogen detection. Here, we describe a qPCR test for routine diagnosis of leprosy-suspected patients. The reaction simultaneously amplifies two specific Mycobacterium leprae targets (16S rRNA and RLEP), and the human 18S rRNA gene as internal control. The limit of detection was estimated to be 2.29 copies of the M. leprae genome. Analytical specificity was evaluated using a panel of 20 other skin pathogenic microorganisms and Mycobacteria, showing no cross-reactivity. Intra- and inter-operator Cp variation was evaluated using dilution curves of M. leprae DNA or a synthetic gene, and no significant difference was observed between three operators in two different laboratories. The multiplex assay was evaluated using 97 patient samples with clinical and histopathological leprosy confirmation, displaying high diagnostic sensitivity (91%) and specificity (100%). Validation tests in an independent panel of 50 samples confirmed sensitivity and specificity of 97% and 98%, respectively. Importantly, assay performance remained stable for at least five months. Our results show that the newly developed multiplex qPCR effectively and specifically detects M. leprae DNA in skin samples, contributing to an efficient diagnosis that expedites the appropriate treatment.

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.

Autophagy-Associated IL-15 Production Is Involved in the Pathogenesis of Leprosy Type 1 Reaction.

Leprosy reactional episodes are acute inflammatory events that may occur during the clinical course of the disease. Type 1 reaction (T1R) is associated with an increase in neural damage, and the understanding of the molecular pathways related to T1R onset is pivotal for the development of strategies that may effectively control the reaction. Interferon-gamma (IFN-γ) is a key cytokine associated with T1R onset and is also associated with autophagy induction. Here, we evaluated the modulation of the autophagy pathway in Mycobacterium leprae-stimulated cells in the presence or absence of IFN-γ. We observed that IFN-γ treatment promoted autophagy activation and increased the expression of genes related to the formation of phagosomes, autophagy regulation and function, or lysosomal pathways in M. leprae-stimulated cells. IFN-γ increased interleukin (IL)-15 secretion in M. leprae-stimulated THP-1 cells in a process associated with autophagy activation. We also observed higher IL15 gene expression in multibacillary (MB) patients who later developed T1R during clinical follow-up when compared to MB patients who did not develop the episode. By overlapping gene expression patterns, we observed 13 common elements shared between T1R skin lesion cells and THP-1 cells stimulated with both M. leprae and IFN-γ. Among these genes, the autophagy regulator Translocated Promoter Region, Nuclear Basket Protein (TPR) was significantly increased in T1R cells when compared with non-reactional MB cells. Overall, our results indicate that IFN-γ may induce a TPR-mediated autophagy transcriptional program in M. leprae-stimulated cells similar to that observed in skin cells during T1R by a pathway that involves IL-15 production, suggesting the involvement of this cytokine in the pathogenesis of T1R.

Putative pathogen-selected polymorphisms in the PKLR gene are associated with mycobacterial susceptibility in Brazilian and African populations.

Pyruvate kinase (PK), encoded by the PKLR gene, is a key player in glycolysis controlling the integrity of erythrocytes. Due to Plasmodium selection, mutations for PK deficiency, which leads to hemolytic anemia, are associated with resistance to malaria in sub-Saharan Africa and with susceptibility to intracellular pathogens in experimental models. In this case-control study, we enrolled 4,555 individuals and investigated whether PKLR single nucleotide polymorphisms (SNPs) putatively selected for malaria resistance are associated with susceptibility to leprosy across Brazil (Manaus-North; Salvador-Northeast; Rondonópolis-Midwest and Rio de Janeiro-Southeast) and with tuberculosis in Mozambique. Haplotype T/G/G (rs1052176/rs4971072/rs11264359) was associated with leprosy susceptibility in Rio de Janeiro (OR = 2.46, p = 0.00001) and Salvador (OR = 1.57, p = 0.04), and with tuberculosis in Mozambique (OR = 1.52, p = 0.07). This haplotype downregulates PKLR expression in nerve and skin, accordingly to GTEx, and might subtly modulate ferritin and haptoglobin levels in serum. Furthermore, we observed genetic signatures of positive selection in the HCN3 gene (xpEHH>2 -recent selection) in Europe but not in Africa, involving 6 SNPs which are PKLR/HCN3 eQTLs. However, this evidence was not corroborated by the other tests (FST, Tajima's D and iHS). Altogether, we provide evidence that a common PKLR locus in Africans contribute to mycobacterial susceptibility in African descent populations and also highlight, for first, PKLR as a susceptibility gene for leprosy and TB.

Interferon-lambda 3 and 4 Polymorphisms Increase Sustained Virological Responses and Regulate Innate Immunity in Antiviral Therapy With Pegylated Interferon-Alpha.

Sustained virologic response (SVR) in chronic hepatitis C (CHC) treatment denotes that the host genetics controls the immune response and unequivocally contribute to viral clearance or disease severity. In this context, single nucleotide polymorphisms (SNPs) in the locus of interferon lambda 3 and 4 genes (IFNL3/4) have been important genetic markers of responsiveness to CHC as prognostic markers for the pegylated-Interferon-alpha/ribavirin (Peg-IFN-α/RBV). Here, we analyzed 12 SNPs at the IFNL3/4 region in 740 treatment-naïve patients with CHC infected with hepatitis C virus (HCV) genotypes 1, 2, or 3 treated with Peg-IFN-α/RBV. Individually, rs12979860-CC, rs8109886-CC, or rs8099917-TT were predictive markers of SVR, while rs12979860-CC demonstrated the stronger effect. Besides, the genotypic combination of these three predictors' genotypes, CC/CC/TT, increased the rate of SVR. Serum levels of cytokines and gene expression analysis on the genes IFNL3, IFNL4, IFNA1, and some of the IFN-stimulated genes (ISGs) were measured in a subgroup of 24 treated patients and 24 healthy volunteers. An antagonist effect was highlighted between the expression of IFNL3/4 and IFNA1 mRNA among patients. Besides, a prominent production of the pro-inflammatory chemokines CCL4 and CXCL10 was observed at a 12-week treatment follow-up. Lower serum levels of these chemokines were detected in patients with an rs12979860-CC genotype associated with the better treatment outcome. Also, lower expression levels of the IFI6, IFI16, IRF9 genes were observed among rs12979860-CC individuals. In conclusion, a combination of the genotypes at the IFNL3/4 locus can act as a better marker for the prognosis for virological responses in an admixed Brazilian population presenting the modulating effect over innate immunity and inflammation that are controlling the outcome of the viral infection, but also other infectious diseases. This study is registered on the ClinicalTrials.gov platform (accession number NCT01889849 and NCT01623336).

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.

Presence of Senescent and Memory CD8+ Leukocytes as Immunocenescence Markers in Skin Lesions of Elderly Leprosy Patients.

Leprosy is an infectious disease that remains endemic in approximately 100 developing countries, where about 200,000 new cases are diagnosed each year. Moreover, multibacillary leprosy, the most contagious form of the disease, has been detected at continuously higher rates among Brazilian elderly people. Due to the so-called immunosenescence, characterized by several alterations in the quality of the immune response during aging, this group is more susceptible to infectious diseases. In view of such data, the purpose of our work was to investigate if age-related alterations in the immune response could influence the pathogenesis of leprosy. As such, we studied 87 individuals, 62 newly diagnosed and untreated leprosy patients distributed according to the age range and to the clinical forms of the disease and 25 healthy volunteers, who were studied as controls. The frequency of senescent and memory CD8+ leukocytes was assessed by immunofluorescence of biopsies from cutaneous lesions, while the serum levels of IgG anti-CMV antibodies were analyzed by chemiluminescence and the gene expression of T cell receptors' inhibitors by RT-qPCR. We noted an accumulation of memory CD8+ T lymphocytes, as well as reduced CD8+CD28+ cell expression in skin lesions from elderly patients, when compared to younger people. Alterations in LAG3 and PDCD1 gene expression in cutaneous lesions of young MB patients were also observed, when compared to elderly patients. Such data suggest that the age-related alterations of T lymphocyte subsets can facilitate the onset of leprosy in elderly patients, not to mention other chronic inflammatory diseases.

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.

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

Due to multiple hypothesis testing with often limited sample size, microarrays and other-omics technologies can sometimes produce irreproducible findings. Complementary to better experimental design, reanalysis and integration of gene expression datasets may help overcome reproducibility issues by identifying consistent differentially expressed genes from independent studies. In this work, after a systematic search, nine microarray datasets evaluating host gene expression in leprosy were reanalyzed and the information was integrated to strengthen evidence of differential expression for several genes. Our results are relevant in prioritizing genes and pathways for further investigation, whether in functional studies or in biomarker discovery. Reanalysis of individual datasets revealed several differentially expressed genes (DEGs) in accordance with original reports. Then, five integration methods (P value and effect size based) were tested. In the end, random-effects model and ratio association were selected as the main methods to pinpoint DEGs. Overall, classic pathways were found corroborating previous findings and validating this approach. Also, we identified some novel DEG involved especially with skin development processes (AQP3, AKR1C3, CYP27B1, LTB, VDR) and keratinocyte biology (CSTA, DSG1, KRT14, KRT5, PKP1, IVL), both still poorly understood in leprosy context. In addition, here we provide aggregated evidence towards some gene candidates that should be prioritized in further leprosy research, as they are likely important in immunopathogenesis. Altogether, these data are useful in better understanding host responses to the disease and, at the same time, provide a list of potential host biomarkers that could be useful in complementing leprosy diagnosis based on transcriptional levels.

Population Genomics of Mycobacterium leprae Reveals a New Genotype in Madagascar and the Comoros.

Human settlement of Madagascar traces back to the beginning of the first millennium with the arrival of Austronesians from Southeast Asia, followed by migrations from Africa and the Middle East. Remains of these different cultural, genetic, and linguistic legacies are still present in Madagascar and other islands of the Indian Ocean. The close relationship between human migration and the introduction and spread of infectious diseases, a well-documented phenomenon, is particularly evident for the causative agent of leprosy, Mycobacterium leprae. In this study, we used whole-genome sequencing (WGS) and molecular dating to characterize the genetic background and retrace the origin of the M. leprae strains circulating in Madagascar (n = 30) and the Comoros (n = 3), two islands where leprosy is still considered a public health problem and monitored as part of a drug resistance surveillance program. Most M. leprae strains (97%) from Madagascar and Comoros belonged to a new genotype as part of branch 1, closely related to single nucleotide polymorphism (SNP) type 1D, named 1D-Malagasy. Other strains belonged to the genotype 1A (3%). We sequenced 39 strains from nine other countries, which, together with previously published genomes, amounted to 242 genomes that were used for molecular dating. Specific SNP markers for the new 1D-Malagasy genotype were used to screen samples from 11 countries and revealed this genotype to be restricted to Madagascar, with the sole exception being a strain from Malawi. The overall analysis thus ruled out a possible introduction of leprosy by the Austronesian settlers and suggests a later origin from East Africa, the Middle East, or South Asia.

Ultra-sensitive detection of Mycobacterium leprae: DNA extraction and PCR assays.

Leprosy urgently needs a precise and early diagnostic tool. The sensitivity of the direct (bacilli staining, Mycobacterium leprae DNA) and indirect (antibody levels, T cell assays) diagnostics methods vary based on the clinical form. Recently, PCR-based M. leprae DNA detection has been shown to differentially diagnose leprosy from other dermatological conditions. However, accuracy can still be improved, especially for use with less invasive clinical samples. We tested different commercial DNA extraction kits: DNeasy Blood & Tissue, QIAamp DNA Microbiome, Maxwell 16 DNA Purification, PowerSoil DNA Isolation; as well as in-house phenol-chloroform and Trizol/FastPrep methods. Extraction was performed on M. leprae-infected mouse footpads and different clinical samples of leprosy patients (skin biopsies and scrapings, lesion, oral and nasal swabs, body hair, blood on FTA cards, peripheral whole blood). We observed that the Microbiome kit was able to enrich for mycobacterial DNA, most likely due the enzymatic digestion cocktail along with mechanical disruption involved in this method. Consequently, we had a significant increase in sensitivity in skin biopsies from paucibacillary leprosy patients using a duplex qPCR targeting 16S rRNA (M. leprae) and 18S rRNA (mammal) in the StepOnePlus system. Our data showed that the presence of M. leprae DNA was best detected in skin biopsies and skin scrapings, independent of the extraction method or the clinical form. For multibacillary patients, detection of M. leprae DNA in nasal swabs indicates the possibility of having a much less invasive sample that can be used for the purposes of DNA sequencing for relapse analysis and drug resistance monitoring. Overall, DNA extracted with the Microbiome kit presented the best bacilli detection rate for paucibacillary cases, indicating that investments in extraction methods with mechanical and DNA digestion should be made.

Ultra-sensitive detection of Mycobacterium leprae: DNA extraction and PCR assays

Leprosy is difficult to diagnose since it is caused by a bacterium that does not grow in vitro. Bacilli direct detection or the presence of specific antibodies can vary greatly depending on the clinical form. M. leprae direct DNA detection can aid clinical diagnosis, although invasive skin biopsies are still necessary to detect the pathogen or histological features consistent with leprosy. Here we show that a kit combining mechanical and chemical lysis efficiently removes host DNA and enriches for M. leprae DNA, allowing better detection of paucibacillary cases. We believe our findings can contribute to improving disease diagnosis, as well as early detection and that could help monitoring strategies(AU).

Reprogramming of Small Noncoding RNA Populations in Peripheral Blood Reveals Host Biomarkers for Latent and Active Mycobacterium tuberculosis Infection.

In tuberculosis (TB), as in other infectious diseases, studies of small noncoding RNAs (sncRNA) in peripheral blood have focused on microRNAs (miRNAs) but have neglected the other major sncRNA classes in spite of their potential functions in host gene regulation. Using RNA sequencing of whole blood, we have therefore determined expression of miRNA, PIWI-interacting RNA (piRNA), small nucleolar RNA (snoRNA), and small nuclear RNA (snRNA) in patients with TB (n = 8), latent TB infection (LTBI; n = 21), and treated LTBI (LTBItt; n = 6) and in uninfected exposed controls (ExC; n = 14). As expected, sncRNA reprogramming was greater in TB than in LTBI, with the greatest changes seen in miRNA populations. However, substantial dynamics were also evident in piRNA and snoRNA populations. One miRNA and 2 piRNAs were identified as moderately accurate (area under the curve [AUC] = 0.70 to 0.74) biomarkers for LTBI, as were 1 miRNA, 1 piRNA, and 2 snoRNAs (AUC = 0.79 to 0.91) for accomplished LTBI treatment. Logistic regression identified the combination of 4 sncRNA (let-7a-5p, miR-589-5p, miR-196b-5p, and SNORD104) as a highly sensitive (100%) classifier to discriminate TB from all non-TB groups. Notably, it reclassified 8 presumed LTBI cases as TB cases, 5 of which turned out to have features of Mycobacterium tuberculosis infection on chest radiographs. SNORD104 expression decreased during M. tuberculosis infection of primary human peripheral blood mononuclear cells (PBMC) and M2-like (P = 0.03) but not M1-like (P = 0.31) macrophages, suggesting that its downregulation in peripheral blood in TB is biologically relevant. Taken together, the results demonstrate that snoRNA and piRNA should be considered in addition to miRNA as biomarkers and pathogenesis factors in the various stages of TB.IMPORTANCE Tuberculosis is the infectious disease with the worldwide largest disease burden and there remains a great need for better diagnostic biomarkers to detect latent and active M. tuberculosis infection. RNA molecules hold great promise in this regard, as their levels of expression may differ considerably between infected and uninfected subjects. We have measured expression changes in the four major classes of small noncoding RNAs in blood samples from patients with different stages of TB infection. We found that, in addition to miRNAs (which are known to be highly regulated in blood cells from TB patients), expression of piRNA and snoRNA is greatly altered in both latent and active TB, yielding promising biomarkers. Even though the functions of many sncRNA other than miRNA are still poorly understood, our results strongly suggest that at least piRNA and snoRNA populations may represent hitherto underappreciated players in the different stages of TB infection.