Identification of a systemic interferon-γ inducible antimicrobial gene signature in leprosy patients undergoing reversal reaction

Reversal reactions (RRs) in leprosy are characterized by a reduction in the number of bacilli in lesions associated with an increase in cell-mediated immunity against the intracellular bacterium Mycobacterium leprae, the causative pathogen of leprosy. To identify the mechanisms that contribute to cell-mediated immunity in leprosy, we measured changes in the whole blood-derived transcriptome of patients with leprosy before, during and after RR. We identified an ‘RR signature’ of 1017 genes that were upregulated at the time of the clinical diagnosis of RR. Using weighted gene correlated network analysis (WGCNA), we detected a module of 794 genes, bisque4, that was significantly correlated with RR, of which 434 genes were part of the RR signature. An enrichment for both IFN-γ and IFN-β downstream gene pathways was present in the RR signature as well as the RR upregulated genes in the bisque4 module, including those encoding proteins of the guanylate binding protein (GBP) family that contributes to antimicrobial responses against mycobacteria. Specifically, GBP1, GBP2, GBP3 and GBP5 mRNAs were upregulated in the RR peripheral blood transcriptome, with GBP1, GBP2 and GBP5 mRNAs also upregulated in the RR disease lesion transcriptome. These data indicate that RRs involve a systemic upregulation of IFN-γ downstream genes including GBP family members as part of the host antimicrobial response against mycobacteria.

Reversal reaction (RR) is a major cause of tissue injury and disabilities in leprosy, resulting from the rapid onset of cell-mediated immune responses to the intracellular bacterium Mycobacterium leprae. To identify the mechanisms related to the increase in cell-mediated immunity in RR, we measured changes in the whole blood-derived transcriptome of patients with leprosy before, during and after RR. We identified that RRs are associated with an IFN-γ induced inflammatory response including an antimicrobial gene network containing the GBP1, GBP2, GBP3 and GBP5 mRNAs in the peripheral blood. Furthermore, we show that GBP1, GBP2 and GBP5 mRNAs are also upregulated at the site of disease in RR patients. In summary, our study suggests that RR involves a systemic induction of IFN-γ regulated genes that contributes to an antimicrobial response against the pathogen, releasing ligands and antigens that can further amplify the inflammatory response.

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).

Genomic Characterization of Mycobacterium leprae to Explore Transmission Patterns Identifies New Subtype in Bangladesh

Mycobacterium leprae, the causative agent of leprosy, is an unculturable bacterium with a considerably reduced genome (3.27 Mb) compared to homologues mycobacteria from the same ancestry. In 2001, the genome of M. leprae was first described and subsequently four genotypes (1-4) and 16 subtypes (A-P) were identified providing means to study global transmission patterns for leprosy. In order to understand the role of asymptomatic carriers we investigated M. leprae carriage as well as infection in leprosy patients (n = 60) and healthy household contacts (HHC; n = 250) from Bangladesh using molecular detection of the bacterial element RLEP in nasal swabs (NS) and slit skin smears (SSS). In parallel, to study M. leprae genotype distribution in Bangladesh we explored strain diversity by whole genome sequencing (WGS) and Sanger sequencing. In the studied cohort in Bangladesh, M. leprae DNA was detected in 33.3% of NS and 22.2% of SSS of patients with bacillary index of 0 whilst in HHC 18.0% of NS and 12.3% of SSS were positive. The majority of the M. leprae strains detected in this study belonged to genotype 1D (55%), followed by 1A (31%). Importantly, WGS allowed the identification of a new M. leprae genotype, designated 1B-Bangladesh (14%), which clustered separately between the 1A and 1B strains. Moreover, we established that the genotype previously designated 1C, is not an independent subtype but clusters within the 1D genotype. Intraindividual differences were present between the M. leprae strains obtained including mutations in hypermutated genes, suggesting mixed colonization/infection or in-host evolution. In summary, we observed that M. leprae is present in asymptomatic contacts of leprosy patients fueling the concept that these individuals contribute to the current intensity of transmission. Our data therefore emphasize the importance of sensitive and specific tools allowing post-exposure prophylaxis targeted at M. leprae-infected or -colonized individuals.

Lack of evidence for the presence of leprosy bacilli in red squirrels from North-West Europe

Leprosy is a human infectious disease caused by Mycobacterium leprae or Mycobacterium lepromatosis that can also occur in animals and even manifest as zoonosis. Recently, both mycobacteria were detected in red squirrels (Sciurus vulgaris) from the British Isles. To further explore the presence of leprosy bacilli in North-West Europe, we screened Belgian and Dutch squirrels. Tissue samples from 115 animals tested by qPCR were negative for both pathogens. No molecular or pathological evidence was found of the presence of these zoonotic pathogens in North-West Europe.

Household Contacts of Leprosy Patients in Endemic Areas Display a Specific Innate Immunity Profile

Leprosy is a chronic infectious disease, caused by Mycobacterium leprae, that can lead to severe life-long disabilities. The transmission of M. leprae is continuously ongoing as witnessed by the stable new case detection rate. The majority of exposed individuals does, however, not develop leprosy and is protected from infection by innate immune mechanisms. In this study the relation between innate immune markers and M. leprae infection as well as the occurrence of leprosy was studied in household contacts (HCs) of leprosy patients with high bacillary loads. Serum proteins associated with innate immunity (ApoA1, CCL4, CRP, IL-1Ra, IL-6, IP-10, and S100A12) were determined by lateral flow assays (LFAs) in conjunction with the presence of M. leprae DNA in nasal swabs (NS) and/or slit-skin smears (SSS). The HCs displayed ApoA1 and S100A12 levels similar to paucibacillary patients and could be differentiated from endemic controls based on the levels of these markers. In the 31 households included the number (percentage) of HCs that were concomitantly diagnosed with leprosy, or tested positive for M. leprae DNA in NS and SSS, was not equally divided. Specifically, households where M. leprae infection and leprosy disease was not observed amongst members of the household were characterized by higher S100A12 and lower CCL4 levels in whole blood assays of HCs in response to M. leprae. Lateral flow assays provide a convenient diagnostic tool to quantitatively measure markers of the innate immune response and thereby detect individuals which are likely infected with M. leprae and at risk of developing disease or transmitting bacteria. Low complexity diagnostic tests measuring innate immunity markers can therefore be applied to help identify who should be targeted for prophylactic treatment.

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.

Archival, paleopathological and aDNA-based techniques in leprosy research and the case of Father Petrus Donders at the Leprosarium 'Batavia', Suriname

OBJECTIVE

We assessed whether Petrus Donders (died 1887), a Dutch priest who for 27 years cared for people with leprosy in the leprosarium Batavia, Suriname, had evidence of Mycobacterium (M.) leprae infection. A positive finding of M. leprae ancient (a)DNA would contribute to the origin of leprosy in Suriname.

MATERIALS

Skeletal remains of Father Petrus Donders; two additional skeletons excavated from the Batavia cemetery were used as controls.

METHODS

Archival research, paleopathological evaluation and aDNA-based testing of skeletal remains.

RESULTS

Neither archives nor inspection of Donders skeletal remains revealed evidence of leprosy, and aDNA-based testing for M. leprae was negative. We detected M. leprae aDNA by RLEP PCR in one control skeleton, which also displayed pathological lesions compatible with leprosy. The M. leprae aDNA was genotyped by Sanger sequencing as SNP type 4; the skeleton displayed mitochondrial haplogroup L3.

CONCLUSION

We found no evidence that Donders contracted leprosy despite years of intense leprosy contact, but we successfully isolated an archaeological M. leprae aDNA sample from a control skeleton from South America.

SIGNIFICANCE

We successfully genotyped recovered aDNA to a M. leprae strain that likely originated in West Africa. The detected human mitochondrial haplogroup L3 is also associated with this geographical region. This suggests that slave trade contributed to leprosy in Suriname.

LIMITATIONS

A limited number of skeletons was examined.

SUGGESTIONS FOR FURTHER RESEARCH

Broader review of skeletal collections is advised to expand on diversity of the M. leprae aDNA database.

Mycobacterium leprae and host immune transcriptomic signatures for reactional states in leprosy

Background

Mycobacterium leprae transcriptomic and human host immune gene expression signatures that demonstrate a plausible association with type I (T1R) and type II reactions (T2R) aid in early diagnosis, prevention of nerve damage and consequent demyelinating neuropathy in leprosy. The aim of the study is to identify M. leprae and host-associated gene-expression signatures that are associated with reactional states in leprosy.

Methods

The differentially expressed genes from the whole transcriptome of M. leprae were determined using genome-wide hybridization arrays with RNA extracted from skin biopsies of 20 T1R, 20 T2R and 20 non reactional controls (NR). Additionally, human immune gene-expressions were profiled using RT2-PCR profiler arrays and real-time qPCRs.

Results

The RNA quality was optimal in 16 NR, 18 T1R and 19 T2R samples. Whole transcriptome expression array of these samples revealed significant upregulation of the genes that encode integral and intrinsic membrane proteins, hydrolases and oxidoreductases. In T1R lesional skin biopsy specimens, the top 10 significantly upregulated genes are ML2064, ML1271, ML1960, ML1220, ML2498, ML1996, ML2388, ML0429, ML2030 and ML0224 in comparison to NR. In T2R, genes ML2498, ML1526, ML0394, ML1960, ML2388, ML0429, ML0281, ML1847, ML1618 and ML1271 were significantly upregulated. We noted ML2664 was significantly upregulated in T1R and repressed in T2R. Conversely, we have not noted any genes upregulated in T2R and repressed in T1R. In both T1R and T2R, ML2388 was significantly upregulated. This gene encodes a probable membrane protein and epitope prediction using Bepipred-2.0 revealed a distinct B-cell epitope. Overexpression of ML2388 was noted consistently across the reaction samples. From the host immune gene expression profiles, genes for CXCL9, CXCL10, CXCL2, CD40LG, IL17A and CXCL11 were upregulated in T1R when compared to the NR. In T2R, CXCL10, CXCL11, CXCL9, CXCL2 and CD40LG were upregulated when compared to the NR group.

Conclusion

A gene set signature involving bacterial genes ML2388, ML2664, and host immune genes CXCL10 and IL-17A can be transcriptomic markers for reactional states in leprosy.

Family-based genetics identifies association of CUBN, IL1RL1 and PRKN variants with leprosy in Bangladesh

Previous studies have demonstrated the association between the occurrence of leprosy and several host genetic variants that can vary across different ethnic groups and populations. Since leprosy susceptibility genes have not been investigated in detail in the Bangladeshi population that still includes highly endemic areas, we here studied known associations between leprosy and 13 genetic markers located in 11 genes in a leprosy endemic area in Bangladesh for which an association with leprosy was previously found in an isolated former leprosy colony in the Brazilian Amazon. For this study, a family-based analysis using 60 parent-affected child trios was performed, followed by a case-control study of 210 leprosy patients and 189 healthy controls from the same area (endemic controls). Genotypes for 11 markers were determined by TaqMan SNP genotyping. In the family-based study, a significant association was found between leprosy and three single nucleotide polymorphisms, rs1801224, rs13001714, and rs1801582, located in the CUBN, IL1RL1, and PRKN genes, respectively. These findings were not replicated in the case-control sample. Variants in the CUBN, IL1RL1, and PRKN genes were associated with leprosy in Bangladesh, validating the initial Brazilian finding in a population sample of distinct ethnic background.