Deciphering the genetic control of gene expression following Mycobacterium leprae antigen stimulation

Genotype × environment interactions in gene regulation and complex traits

Genotype × environment interactions (GxE) have long been recognized as a key mechanism underlying human phenotypic variation. Technological developments over the past 15 years have dramatically expanded our appreciation of the role of GxE in both gene regulation and complex traits. The richness and complexity of these datasets also required parallel efforts to develop robust and sensitive statistical and computational approaches. Although our understanding of the genetic architecture of molecular and complex traits has been maturing, a large proportion of complex trait heritability remains unexplained. Furthermore, there are increasing efforts to characterize the effect of environmental exposure on human health. We therefore review GxE in human gene regulation and complex traits, advocating for a comprehensive approach that jointly considers genetic and environmental factors in human health and disease.

Shaping immunity: The influence of natural selection on population immune diversity

Humans exhibit considerable variability in their immune responses to the same immune challenges. Such variation is widespread and affects individual and population-level susceptibility to infectious diseases and immune disorders. Although the factors influencing immune response diversity are partially understood, what mechanisms lead to the wide range of immune traits in healthy individuals remain largely unexplained. Here, we discuss the role that natural selection has played in driving phenotypic differences in immune responses across populations and present-day susceptibility to immune-related disorders. Further, we touch on future directions in the field of immunogenomics, highlighting the value of expanding this work to human populations globally, the utility of modeling the immune response as a dynamic process, and the importance of considering the potential polygenic nature of natural selection. Identifying loci acted upon by evolution may further pinpoint variants critically involved in disease etiology, and designing studies to capture these effects will enrich our understanding of the genetic contributions to immunity and immune dysregulation.

Characterization of caffeine response regulatory variants in vascular endothelial cells

Genetic variants in gene regulatory sequences can modify gene expression and mediate the molecular response to environmental stimuli. In addition, genotype-environment interactions (GxE) contribute to complex traits such as cardiovascular disease. Caffeine is the most widely consumed stimulant and is known to produce a vascular response. To investigate GxE for caffeine, we treated vascular endothelial cells with caffeine and used a massively parallel reporter assay to measure allelic effects on gene regulation for over 43,000 genetic variants. We identified 665 variants with allelic effects on gene regulation and 6 variants that regulate the gene expression response to caffeine (GxE, false discovery rate [FDR] < 5%). When overlapping our GxE results with expression quantitative trait loci colocalized with coronary artery disease and hypertension, we dissected their regulatory mechanisms and showed a modulatory role for caffeine. Our results demonstrate that massively parallel reporter assay is a powerful approach to identify and molecularly characterize GxE in the specific context of caffeine consumption.

Genome-wide meta-analysis and fine-mapping prioritize potential causal variants and genes related to leprosy

To date, genome-wide association studies (GWASs) have discovered 35 susceptible loci of leprosy; however, the cumulative effects of these loci can only partially explain the overall risk of leprosy, and the causal variants and genes within these loci remain unknown. Here, we conducted out new GWASs in two independent cohorts of 5007 cases and 4579 controls and then a meta-analysis in these newly generated and multiple previously published (2277 cases and 3159 controls) datasets were performed. Three novel and 15 previously reported risk loci were identified from these datasets, increasing the known leprosy risk loci of explained genetic heritability from 23.0 to 38.5%. A comprehensive fine-mapping analysis was conducted, and 19 causal variants and 14 causal genes were identified. Specifically, manual checking of epigenomic information from the Epimap database revealed that the causal variants were mainly located within the immune-relevant or immune-specific regulatory elements. Furthermore, by using gene-set, tissue, and cell-type enrichment analyses, we highlighted the key roles of immune-related tissues and cells and implicated the PD-1 signaling pathways in the pathogenetic mechanism of leprosy. Collectively, our study identified candidate causal variants and elucidated the potential regulatory and coding mechanisms for genes associated with leprosy.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single-cell resolution

Synthetic glucocorticoids, such as dexamethasone, have been used as a treatment for many immune conditions, such as asthma and, more recently, severe COVID-19. Single-cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single-cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated peripheral blood mononuclear cells from 96 African American children. We used novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we showed that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and showed widespread genetic regulation of the transcriptional dynamics of the gene expression response.

Stem Cell Models for Context-Specific Modeling in Psychiatric Disorders

Genome-wide association studies reveal the complex polygenic architecture underlying psychiatric disorder risk, but there is an unmet need to validate causal variants, resolve their target genes(s), and explore their functional impacts on disorder-related mechanisms. Disorder-associated loci regulate transcription of target genes in a cell type- and context-specific manner, which can be measured through expression quantitative trait loci. In this review, we discuss methods and insights from context-specific modeling of genetically and environmentally regulated expression. Human induced pluripotent stem cell-derived cell type and organoid models have uncovered context-specific psychiatric disorder associations by investigating tissue-, cell type-, sex-, age-, and stressor-specific genetic regulation of expression. Techniques such as massively parallel reporter assays and pooled CRISPR (clustered regularly interspaced short palindromic repeats) screens make it possible to functionally fine-map genome-wide association study loci and validate their target genes at scale. Integration of disorder-associated contexts with these patient-specific human induced pluripotent stem cell models makes it possible to uncover gene by environment interactions that mediate disorder risk, which will ultimately improve our ability to diagnose and treat psychiatric disorders.

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.

Analysis of transcriptional changes in the immune system associated with pubertal development in a longitudinal cohort of children with asthma

Puberty is an important developmental period marked by hormonal, metabolic and immune changes. Puberty also marks a shift in sex differences in susceptibility to asthma. Yet, little is known about the gene expression changes in immune cells that occur during pubertal development. Here we assess pubertal development and leukocyte gene expression in a longitudinal cohort of 251 children with asthma. We identify substantial gene expression changes associated with age and pubertal development. Gene expression changes between pre- and post-menarcheal females suggest a shift from predominantly innate to adaptive immunity. We show that genetic effects on gene expression change dynamically during pubertal development. Gene expression changes during puberty are correlated with gene expression changes associated with asthma and may explain sex differences in prevalence. Our results show that molecular data used to study the genetics of early onset diseases should consider pubertal development as an important factor that modifies the transcriptome.

Genome-wide association study of leprosy in Malawi and Mali

Leprosy is a chronic infection of the skin and peripheral nerves caused by Mycobacterium leprae. Despite recent improvements in disease control, leprosy remains an important cause of infectious disability globally. Large-scale genetic association studies in Chinese, Vietnamese and Indian populations have identified over 30 susceptibility loci for leprosy. There is a significant burden of leprosy in Africa, however it is uncertain whether the findings of published genetic association studies are generalizable to African populations. To address this, we conducted a genome-wide association study (GWAS) of leprosy in Malawian (327 cases, 436 controls) and Malian (247 cases, 368 controls) individuals. In that analysis, we replicated four risk loci previously reported in China, Vietnam and India; MHC Class I and II, LACC1 and SLC29A3. We further identified a novel leprosy susceptibility locus at 10q24 (rs2015583; combined p = 8.81 × 10-9; OR = 0.51 [95% CI 0.40 - 0.64]). Using publicly-available data we characterise regulatory activity at this locus, identifying ACTR1A as a candidate mediator of leprosy risk. This locus shows evidence of recent positive selection and demonstrates pleiotropy with established risk loci for inflammatory bowel disease and childhood-onset asthma. A shared genetic architecture for leprosy and inflammatory bowel disease has been previously described. We expand on this, strengthening the hypothesis that selection pressure driven by leprosy has shaped the evolution of autoimmune and atopic disease in modern populations. More broadly, our data highlights the importance of defining the genetic architecture of disease across genetically diverse populations, and that disease insights derived from GWAS in one population may not translate to all affected populations.

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.

Transcriptomic Analysis of Mycobacterium leprae-Stimulated Response in Peripheral Blood Mononuclear Cells Reveal Potential Biomarkers for Early Diagnosis of Leprosy

We aimed to identify an unique host transcriptional signature in peripheral blood mononuclear cells (PBMCs) in response to Mycobacterium leprae antigens to distinguish between patients with leprosy and non-leprosy controls for early diagnosis of the disease. Sixteen individuals were enrolled in the discovery cohort [eight patients with leprosy, comprising four multibacillary (MB) and four paucibacillary (PB); and eight non-leprosy controls, comprising four healthy house contacts (HHCs) and four endemic controls (ECs)]. The differences in the transcriptome response of PBMCs to M. leprae sonicate antigen were evaluated between leprosy patients and non-leprosy controls, and 12 differentially expressed genes (CCL2/MCP-1, IL-8, JAKM, ATP, ND1, SERP, FLJ10489, LINC00659, LOC34487, LOC101928143, MIR22, and NCF1C) were identified. The accuracy of the 12 differentially expressed genes was further validated for the diagnosis of leprosy using real-time quantitative PCR in 82 individuals (13 MB, 10 PB, 37 HHCs, and 22 ECs) in the validation cohort. We found that a 5 gene signature set IL-8, CCL2/MCP-1, SERP, LINC00659 and FLJ10489 had a suitable performance in discriminating leprosy from ECs. In addition, elevated expression of IL-8, CCL2/MCP-1, SERP and LINC00659 was associated with MB diagnosis compared with ECs, whereas increased expression of IL-8, CCL2/MCP-1, SERP and FLJ10489 was found to be useful biomarkers for PB diagnosis from ECs. Moreover, we found decreased expression of NCF1C among leprosy patients could distinguish leprosy from HHCs, whereas higher expression of CCL2 among MB than PB could distinguish different leprosy patients. In conclusion, among the 12 candidate host genes identified, a three gene signature IL-8, CCL2/MCP-1, and SERP showed the best performance in distinguishing leprosy patients from healthy controls. These findings may have implications for developing a rapid blood-based test for early diagnosis of leprosy.

Factors associated with the development of leprosy in Brazilian contacts: a systematic review

People who interact with leprosy patients in their environment, neighborhood, family, or social relationships are at risk to develop the disease. This systematic review investigated the risk and protective factors associated with the development of leprosy in Brazilian contacts. The studies were found in Cochrane Library, PubMed (MEDLINE), Embase, Virtual Health Library, grey literature and hand search until July 2021. The study selection, data extraction and quality assessment were independently performed by two investigators. The quality assessment was performed using the Newcastle-Ottawa Scale (NOS). This review was registered in PROSPERO (CRD42020160680). Seventeen articles fulfilled the inclusion criteria (n=544). The immunological and molecular factors, such as Anti-phenolic Glycolipid Antibodies (Anti-PGL-1) seropositivity, negative Mitsuda test, absence of Bacillus Calmette-Guérin (BCG) scar, positive Polymerase Chain Reaction (PCR) in blood; age and race; conviviality, education, contact time and type of contact, as well as elements related to the index case (bacilloscopic index; genetic conditions, family relationships), and some combined factors were shown to be relevant risk factors associated with the development of the disease in Brazilian leprosy contacts. The protective factors reported were the presence of one or more BCG scars, positive Mitsuda test, and education level. All selected studies were considered of high quality according to NOS. The knowledge of disease-related risk and protective factors provides the scientific basis for decision-making in the management of the disease in leprosy contacts.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single cell resolution.. [PMID: 35313584 PMCID: PMC8936121 DOI: 10.1101/2021.09.30.462672] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Synthetic glucocorticoids, such as dexamethasone, have been used as treatment for many immune conditions, such as asthma and more recently severe COVID-19. Single cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated Peripheral Blood Mononuclear Cells from 96 African American children. We employed novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we demonstrated that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and demonstrated widespread genetic regulation of the transcriptional dynamics of the gene expression response.

A signature of Neanderthal introgression on molecular mechanisms of environmental responses

Ancient human migrations led to the settlement of population groups in varied environmental contexts worldwide. The extent to which adaptation to local environments has shaped human genetic diversity is a longstanding question in human evolution. Recent studies have suggested that introgression of archaic alleles in the genome of modern humans may have contributed to adaptation to environmental pressures such as pathogen exposure. Functional genomic studies have demonstrated that variation in gene expression across individuals and in response to environmental perturbations is a main mechanism underlying complex trait variation. We considered gene expression response to in vitro treatments as a molecular phenotype to identify genes and regulatory variants that may have played an important role in adaptations to local environments. We investigated if Neanderthal introgression in the human genome may contribute to the transcriptional response to environmental perturbations. To this end we used eQTLs for genes differentially expressed in a panel of 52 cellular environments, resulting from 5 cell types and 26 treatments, including hormones, vitamins, drugs, and environmental contaminants. We found that SNPs with introgressed Neanderthal alleles (N-SNPs) disrupt binding of transcription factors important for environmental responses, including ionizing radiation and hypoxia, and for glucose metabolism. We identified an enrichment for N-SNPs among eQTLs for genes differentially expressed in response to 8 treatments, including glucocorticoids, caffeine, and vitamin D. Using Massively Parallel Reporter Assays (MPRA) data, we validated the regulatory function of 21 introgressed Neanderthal variants in the human genome, corresponding to 8 eQTLs regulating 15 genes that respond to environmental perturbations. These findings expand the set of environments where archaic introgression may have contributed to adaptations to local environments in modern humans and provide experimental validation for the regulatory function of introgressed variants.

Functional dynamic genetic effects on gene regulation are specific to particular cell types and environmental conditions

Genetic effects on gene expression and splicing can be modulated by cellular and environmental factors; yet interactions between genotypes, cell type, and treatment have not been comprehensively studied together. We used an induced pluripotent stem cell system to study multiple cell types derived from the same individuals and exposed them to a large panel of treatments. Cellular responses involved different genes and pathways for gene expression and splicing and were highly variable across contexts. For thousands of genes, we identified variable allelic expression across contexts and characterized different types of gene-environment interactions, many of which are associated with complex traits. Promoter functional and evolutionary features distinguished genes with elevated allelic imbalance mean and variance. On average, half of the genes with dynamic regulatory interactions were missed by large eQTL mapping studies, indicating the importance of exploring multiple treatments to reveal previously unrecognized regulatory loci that may be important for disease.

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.

Mapping leprosy-associated coding variants of interleukin genes by targeted sequencing

Previous genotyping-based assays have identified non-coding variants of several interleukins (ILs) being associated with genetic susceptibility to leprosy. However, understanding of the involvement of coding variants within all IL family genes in leprosy was still limited. To obtain the full mutation spectrum of all ILs in leprosy, we performed a targeted deep sequencing of coding regions of 58 ILs genes in 798 leprosy patients (age 56.2 ± 14.4; female 31.5%) and 990 healthy controls (age 38.1 ± 14.0; female 44.3%) from Yunnan, Southwest China. mRNA expression alterations of ILs in leprosy skin lesions or in response to M. leprae treatment were estimated by using publicly available expression datasets. Two coding variants in IL27 (rs17855750, p.S59A, p = 4.02 × 10^-8 , odds ratio [OR] = 1.748) and IL1RN (rs45507693, p.A106T, p = 1.45 × 10^-5 , OR = 3.629) were significantly associated with leprosy risk. mRNA levels of IL27 and IL1RN were upregulated in whole blood cells after M. leprae stimulation. These data showed that IL27 and IL1RN are leprosy risk genes. Further functional study is required for characterizing the exact role of ILs in leprosy.

Dynamic effects of genetic variation on gene expression revealed following hypoxic stress in cardiomyocytes

One life-threatening outcome of cardiovascular disease is myocardial infarction, where cardiomyocytes are deprived of oxygen. To study inter-individual differences in response to hypoxia, we established an in vitro model of induced pluripotent stem cell-derived cardiomyocytes from 15 individuals. We measured gene expression levels, chromatin accessibility, and methylation levels in four culturing conditions that correspond to normoxia, hypoxia, and short- or long-term re-oxygenation. We characterized thousands of gene regulatory changes as the cells transition between conditions. Using available genotypes, we identified 1,573 genes with a cis expression quantitative locus (eQTL) in at least one condition, as well as 367 dynamic eQTLs, which are classified as eQTLs in at least one, but not in all conditions. A subset of genes with dynamic eQTLs is associated with complex traits and disease. Our data demonstrate how dynamic genetic effects on gene expression, which are likely relevant for disease, can be uncovered under stress.

Polymorphisms in mitochondrial ribosomal protein S5 (MRPS5) are associated with leprosy risk in Chinese

Leprosy is an infectious disease caused by Mycobacterium leprae (M. leprae), with about 210,000 new cases per year worldwide. Although numerous risk loci have been uncovered by genome-wide association studies, the effects of common genetic variants are relatively modest. To identify possible new genetic locus involved in susceptibility to leprosy, whole exome sequencing was performed for 28 subjects including 14 patients and 12 unaffected members from 8 leprosy-affected families as well as another case and an unrelated control, and then the follow-up SNP genotyping of the candidate variants was studied in case-control sample sets. A rare missense variant in mitochondrial ribosomal protein S5 (MRPS5), rs200730619 (c. 95108402T>C [p. Tyr137Cys]) was identified and validated in 369 cases and 270 controls of Chinese descent (Padjusted = 0.006, odds ratio [OR] = 2.74) as a contributing factor to leprosy risk. Moreover, the mRNA level of MRPS5 was downregulated in M. leprae sonicate-stimulated peripheral blood mononuclear cells. Our results indicated that MRPS5 may be involved in leprosy pathogenesis. Further studies are needed to determine if defective MRPS5 could lead to impairment of energy metabolism of host immune cells, which could further cause defect in clearing M. leprae and increase susceptibility to infection.

Mycolactone toxin induces an inflammatory response by targeting the IL-1β pathway: Mechanistic insight into Buruli ulcer pathophysiology

Mycolactone, a lipid-like toxin, is the major virulence factor of Mycobacterium ulcerans, the etiological agent of Buruli ulcer. Its involvement in lesion development has been widely described in early stages of the disease, through its cytotoxic and immunosuppressive activities, but less is known about later stages. Here, we revisit the role of mycolactone in disease outcome and provide the first demonstration of the pro-inflammatory potential of this toxin. We found that the mycolactone-containing mycobacterial extracellular vesicles produced by M. ulcerans induced the production of IL-1β, a potent pro-inflammatory cytokine, in a TLR2-dependent manner, targeting NLRP3/1 inflammasomes. We show our data to be relevant in a physiological context. The in vivo injection of these mycolactone-containing vesicles induced a strong local inflammatory response and tissue damage, which were prevented by corticosteroids. Finally, several soluble pro-inflammatory factors, including IL-1β, were detected in infected tissues from mice and Buruli ulcer patients. Our results revisit Buruli ulcer pathophysiology by providing new insight, thus paving the way for the development of new therapeutic strategies taking the pro-inflammatory potential of mycolactone into account.

Buruli ulcer is a neglected tropical disease occurring mainly in poor rural areas of West and Central Africa. This cutaneous disease is caused by Mycobacterium ulcerans, a bacterium belonging to the same family as M. tuberculosis and M. leprae. The skin lesions are caused by a cytotoxic toxin named mycolactone, also known to act as an immunosuppressor and an anti-inflammatory molecule. However, Buruli ulcer lesions are characterized by a chronic cutaneous inflammation with a recruitment of cellular immune cells trying to counteract M. ulcerans. Our work allows for a reconcilitation of previous observations. We found by in vitro experiment on macrophages that the mycolactone-containing mycobacterial extracellular vesicles produced by M. ulcerans induced the production of IL-1β, a potent pro-inflammatory molecule, while other pro-inflammatory soluble factors are inhibited. We also detected IL-1β protein in a mouse model of M. ulcerans infection as well as in biopsies of Buruli ulcer patients. The pro-inflammatory potential of mycolacone has to be taken into account to understand the full pathophysiology of Buruli ulcer.

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.

MiCA: An extended tool for microarray gene expression analysis

Gene expression microarrays capture a complete image of all the transcriptional activity in a biological sample. Microarrays produce a large amount of data, which becomes a challenge when it comes to exploring and interpreting using modern computational and statistical tools. We propose the Microarray Analysis (MiCA) tool that outperforms other similar tools both in terms of ease of use and statistical features requiring minimal input to conduct an analysis. MiCA is an integrated, interactive, and streamlined desktop software for the analysis of microarray gene expression data. MiCA consists of a complete microarray analysis pipeline including but not limited to fetching data directly from GEO, normalization, interactive quality control, batch-effect correction, regression analysis, surrogate variable analysis and functional annotation methods such as GSVA using known existing R packages. We compare the features offered by MiCA and other similar tools while performing differential expression analysis using previously published datasets. MiCA offers additional statistical and visualization methods to conduct a microarray data analysis compared to other available microarray analysis tools. MiCA minimizes the need for technical knowledge by providing a very intuitive and versatile interface that integrates all necessary tasks and features required for basic microarray data analysis. We analyzed multiple published datasets and showed that the features offered by MiCA not only simplify the analysis pipeline but also provide additional interpretation to the data.

Genetics of leprosy: today and beyond

Leprosy is a chronic infectious disease of the skin and peripheral nerves that presents a strong link with the host genetic background. Different approaches in genetic studies have been applied to leprosy and today leprosy is among the infectious diseases with the greatest number of genetic risk variants identified. Several leprosy genes have been implicated in host immune response to pathogens and point to specific pathways that are relevant for host defense to infection. In addition, host genetic factors are also involved in the heterogeneity of leprosy clinical manifestations and in excessive inflammatory responses that occur in some leprosy patients. Finally, genetic studies in leprosy have provided strong evidence of pleiotropic effects between leprosy and other complex diseases, such as immune-mediated or neurodegenerative diseases. These findings not only impact on the field of leprosy and infectious diseases but also make leprosy a good model for the study of complex immune-mediated diseases. Here, we summarize recent genetic findings in leprosy susceptibility and discuss the overlap of the genetic control in leprosy with Parkinson's disease and inflammatory bowel disease. Moreover, some limitations, challenges, and potential new avenues for future genetics studies of leprosy are also discussed in this review.

Human genetics of mycobacterial disease

Mycobacterial diseases are caused by members of the genus Mycobacterium, acid-fast bacteria characterized by the presence of mycolic acids within their cell walls. Claiming almost 2 million lives every year, tuberculosis (TB) is the most common mycobacterial disease and is caused by infection with M. tuberculosis and, in rare cases, by M. bovis or M. africanum. The second and third most common mycobacterial diseases are leprosy and buruli ulcer (BU), respectively. Both diseases affect the skin and can lead to permanent sequelae and deformities. Leprosy is caused by the uncultivable M. leprae while the etiological agent of BU is the environmental bacterium M. ulcerans. After exposure to these mycobacterial species, a majority of individuals will not progress to clinical disease and, among those who do, inter-individual variability in disease manifestation and outcome can be observed. Susceptibility to mycobacterial diseases carries a human genetic component and intense efforts have been applied over the past decades to decipher the exact nature of the genetic factors controlling disease susceptibility. While for BU this search was mostly conducted on the basis of candidate genes association studies, genome-wide approaches have been widely applied for TB and leprosy. In this review, we summarize some of the findings achieved by genome-wide linkage, association and transcriptome analyses in TB disease and leprosy and the recent genetic findings for BU susceptibility.

AmpliSeq transcriptome analysis of human alveolar and monocyte-derived macrophages over time in response to Mycobacterium tuberculosis infection

Human alveolar macrophages (HAM) are primary bacterial niche and immune response cells during Mycobacterium tuberculosis (M.tb) infection, and human blood monocyte-derived macrophages (MDM) are a model for investigating M.tb-macrophage interactions. Here, we use a targeted RNA-Seq method to measure transcriptome-wide changes in RNA expression patterns of freshly obtained HAM (used within 6 h) and 6 day cultured MDM upon M.tb infection over time (2, 24 and 72 h), in both uninfected and infected cells from three donors each. The Ion AmpliSeq™ Transcriptome Human Gene Expression Kit (AmpliSeq) uses primers targeting 18,574 mRNAs and 2,228 non-coding RNAs (ncRNAs) for a total of 20,802 transcripts. AmpliSeqTM yields highly precise and reproducible gene expression profiles (R2 >0.99). Taking advantage of AmpliSeq's reproducibility, we establish well-defined quantitative RNA expression patterns of HAM versus MDM, including significant M.tb-inducible genes, in networks and pathways that differ in part between MDM and HAM. A similar number of expressed genes are detected at all time-points between uninfected MDM and HAM, in common pathways including inflammatory and immune functions, but canonical pathway differences also exist. In particular, at 2 h, multiple genes relevant to the immune response are preferentially expressed in either uninfected HAM or MDM, while the HAM RNA profiles approximate MDM profiles over time in culture, highlighting the unique RNA expression profile of freshly obtained HAM. MDM demonstrate a greater transcriptional response than HAM upon M.tb infection, with 2 to >10 times more genes up- or down-regulated. The results identify key genes involved in cellular responses to M.tb in two different human macrophage types. Follow-up bioinformatics analysis indicates that approximately 30% of response genes have expression quantitative trait loci (eQTLs in GTEx), common DNA variants that can influence host gene expression susceptibility or resistance to M.tb, illustrated with the TREM1 gene cluster and IL-10.

Missense Variants in HIF1A and LACC1 Contribute to Leprosy Risk in Han Chinese

Genome-wide association studies (GWASs) and genome-wide linkage studies (GWLSs) have identified numerous risk genes affecting the susceptibility to leprosy. However, most of the reported GWAS hits are noncoding variants and account for only part of the estimated heritability for this disease. In order to identify additional risk genes and map the potentially functional variants within the GWAS loci, we performed a three-stage study combining whole-exome sequencing (WES; discovery stage), targeted next-generation sequencing (NGS; screening stage), and refined validation of risk missense variants in 1,433 individuals with leprosy and 1,625 healthy control individuals from Yunnan Province, Southwest China. We identified and validated a rare damaging variant, rs142179458 (c.1045G>A [p.Asp349Asn]) in HIF1A, as contributing to leprosy risk (p = 4.95 × 10^-9, odds ratio [OR] = 2.266). We were able to show that affected individuals harboring the risk allele presented with multibacillary leprosy at an earlier age (p = 0.025). We also confirmed the association between missense variant rs3764147 (c.760A>G [p.Ile254Val]) in the GWAS hit LACC1 (formerly C13orf31) and leprosy (p = 6.11 × 10^-18, OR = 1.605). By using the population attributable fraction, we have shown that HIF1A and LACC1 are the major genes with missense variants contributing to leprosy risk in our study groups. Consistently, mRNA expression levels of both HIF1A and LACC1 were upregulated in the skin lesions of individuals with leprosy and in Mycobacterium leprae-stimulated cells, indicating an active role of HIF1A and LACC1 in leprosy pathogenesis.