Cellular stress promotes NOD1/2-dependent inflammation via the endogenous metabolite sphingosine-1-phosphate.

Cellular stress has been associated with inflammation, yet precise underlying mechanisms remain elusive. In this study, various unrelated stress inducers were employed to screen for sensors linking altered cellular homeostasis and inflammation. We identified the intracellular pattern recognition receptors NOD1/2, which sense bacterial peptidoglycans, as general stress sensors detecting perturbations of cellular homeostasis. NOD1/2 activation upon such perturbations required generation of the endogenous metabolite sphingosine-1-phosphate (S1P). Unlike peptidoglycan sensing via the leucine-rich repeats domain, cytosolic S1P directly bound to the nucleotide binding domains of NOD1/2, triggering NF-κB activation and inflammatory responses. In sum, we unveiled a hitherto unknown role of NOD1/2 in surveillance of cellular homeostasis through sensing of the cytosolic metabolite S1P. We propose S1P, an endogenous metabolite, as a novel NOD1/2 activator and NOD1/2 as molecular hubs integrating bacterial and metabolic cues.

Pharmacological inhibition of bromodomain and extra-terminal proteins induces an NRF-2-mediated antiviral state that is subverted by SARS-CoV-2 infection.

Inhibitors of bromodomain and extra-terminal proteins (iBETs), including JQ-1, have been suggested as potential prophylactics against SARS-CoV-2 infection. However, molecular mechanisms underlying JQ-1-mediated antiviral activity and its susceptibility to viral subversion remain incompletely understood. Pretreatment of cells with iBETs inhibited infection by SARS-CoV-2 variants and SARS-CoV, but not MERS-CoV. The antiviral activity manifested itself by reduced reporter expression of recombinant viruses, and reduced viral RNA quantities and infectious titers in the culture supernatant. While we confirmed JQ-1-mediated downregulation of expression of angiotensin-converting enzyme 2 (ACE2) and interferon-stimulated genes (ISGs), multi-omics analysis addressing the chromatin accessibility, transcriptome and proteome uncovered induction of an antiviral nuclear factor erythroid 2-related factor 2 (NRF-2)-mediated cytoprotective response as an additional mechanism through which JQ-1 inhibits SARS-CoV-2 replication. Pharmacological inhibition of NRF-2, and knockdown of NRF-2 and its target genes reduced JQ-1-mediated inhibition of SARS-CoV-2 replication. Serial passaging of SARS-CoV-2 in the presence of JQ-1 resulted in predominance of ORF6-deficient variant, which exhibited resistance to JQ-1 and increased sensitivity to exogenously administered type I interferon (IFN-I), suggesting a minimised need for SARS-CoV-2 ORF6-mediated repression of IFN signalling in the presence of JQ-1. Importantly, JQ-1 exhibited a transient antiviral activity when administered prophylactically in human airway bronchial epithelial cells (hBAECs), which was gradually subverted by SARS-CoV-2, and no antiviral activity when administered therapeutically following an established infection. We propose that JQ-1 exerts pleiotropic effects that collectively induce an antiviral state in the host, which is ultimately nullified by SARS-CoV-2 infection, raising questions about the clinical suitability of the iBETs in the context of COVID-19.

lncRNA CDKN2B-AS1 regulates collagen expression.

The long noncoding RNA CDKN2B-AS1 harbors a major coronary artery disease risk haplotype, which is also associated with progressive forms of the oral inflammatory disease periodontitis as well as myocardial infarction (MI). Despite extensive research, there is currently no broad consensus on the function of CDKN2B-AS1 that would explain a common molecular role of this lncRNA in these diseases. Our aim was to investigate the role of CDKN2B-AS1 in gingival cells to better understand the molecular mechanisms underlying the increased risk of progressive periodontitis. We downregulated CDKN2B-AS1 transcript levels in primary gingival fibroblasts with LNA GapmeRs. Following RNA-sequencing, we performed differential expression, gene set enrichment analyses and Western Blotting. Putative causal alleles were searched by analyzing associated DNA sequence variants for changes of predicted transcription factor binding sites. We functionally characterized putative functional alleles using luciferase-reporter and antibody electrophoretic mobility shift assays in gingival fibroblasts and HeLa cells. Of all gene sets analysed, collagen biosynthesis was most significantly upregulated (Padj=9.7 × 10- 5 (AUC > 0.65) with the CAD and MI risk gene COL4A1 showing strongest upregulation of the enriched gene sets (Fold change = 12.13, Padj = 4.9 × 10- 25). The inflammatory "TNFA signaling via NFKB" gene set was downregulated the most (Padj=1 × 10- 5 (AUC = 0.60). On the single gene level, CAPNS2, involved in extracellular matrix organization, was the top upregulated protein coding gene (Fold change = 48.5, P < 9 × 10- 24). The risk variant rs10757278 altered a binding site of the pathogen responsive transcription factor STAT1 (P = 5.8 × 10- 6). rs10757278-G allele reduced STAT1 binding 14.4% and rs10757278-A decreased luciferase activity in gingival fibroblasts 41.2% (P = 0.0056), corresponding with GTEx data. CDKN2B-AS1 represses collagen gene expression in gingival fibroblasts. Dysregulated collagen biosynthesis through allele-specific CDKN2B-AS1 expression in response to inflammatory factors may affect collagen synthesis, and in consequence tissue barrier and atherosclerotic plaque stability.

Case-only design identifies interactions of genetic risk variants at SIGLEC5 and PLG with the lncRNA CTD-2353F22.1 implying the importance of periodontal wound healing for disease aetiology.

AIM: The basis of phenotypic variation of periodontitis is genetic variability. Disease relevant effects of individual risk alleles are considered to result from genetic interactions. We investigated gene × gene (G×G) interactions of suggestive periodontitis susceptibility alleles. MATERIALS AND METHODS: We used the case-only design and investigated single-nucleotide polymorphism (SNPs) that showed associations in our recent genome-wide association study (GWAS) and GWAS meta-analysis with p < 5 × 10-6 . CRISPR-dCas9 gene activation followed by RNA-sequencing and gene-set enrichment analyses elucidated differentially expressed genes and gene networks. With the databases of SNPInspector and Transfac professional, luciferase reporter gene assays and antibody electrophoretic mobility shift experiments, we analysed allele-specific effects on transcription factor binding. RESULTS: SNPs at the genes sialic acid binding Ig-like lectin 5 (SIGLEC5) and plasminogen (PLG) showed G×G interactions with rs1122900 at the long non-coding RNA (lncRNA) CTD-2353F22. Associated chromatin cis-activated CTD-2353F22.1 6.5-fold (p = .003), indicating CTD-2353F22.1 as target gene of this interaction. CTD-2353F22.1 regulated GADD45A (padj < 4.9 × 10-11 , log2 fold change (FC) = -0.55), THBS1, SERPINE1 and Tissue Factor F3 (padj < 5 × 10-7 , log2 FC ≥ -0.35) and the gene set "angiogenesis" (area under the curve = 0.71, padj  = 8.2 × 10-5 ). rs1122900 effect C-allele decreased reporter gene activity (5.5-fold, p = .0003) and PRDM14 binding (76%). CONCLUSIONS: CTD-2353F22.1 mediates interaction of SIGLEC5 and PLG, together with genes that function in periodontal wound healing.

RSPO4 is a potential risk gene of stages III-IV, grade C periodontitis through effects on innate immune response and oral barrier integrity.

AIM: R-spondin 4 (RSPO4) is a suggestive risk gene of stage III-IV, grade C periodontitis and upregulated in gingiva of mice resistant to bacteria-induced alveolar bone loss. We aimed to replicate the association, identify and characterize the putative causal variant(s) and molecular effects, and understand the downstream effects of RSPO4 upregulation. MATERIALS AND METHODS: We performed a two-step association study for RSPO4 with imputed genotypes of a German-Dutch (896 stage III-IV, grade C periodontitis cases, 7104 controls) and Spanish sample (441 cases and 1141 controls). We analysed the allelic effects on transcription factor binding sites with reporter gene and antibody electrophoretic mobility shift assays. We used CRISPR/dCas9 activation and RNA sequencing to pinpoint RSPO4 as the target gene and to analyse downstream effects. RESULTS: RSPO4 was associated with periodontitis (rs6056178, pmeta  = 4.6 × 10-5 ). rs6056178 contains a GATA-binding motif. The rs6056178 T-allele abolished reporter activity (p = .004) and reduced GATA binding (-14.5%). CRISPRa of the associated region increased RSPO4 expression (25.8 ± 6.5-fold, p = .003). RSPO4 activation showed strongest induction of Gliomedin (439-fold) and Mucin 21 (178-fold) and of the gene set "response to interferon-alpha" (area under the curve [AUC] = 0.8, p < 5 × 10-6 ). The most repressed gene set was "extracellular matrix interactions" (AUC = 0.8, padj  = .00016). CONCLUSION: RSPO4 is a potential periodontitis risk gene and modifies host defence and barrier integrity.

Innate-like Gene Expression of Lung-Resident Memory CD8+ T Cells during Experimental Human Influenza: A Clinical Study.

Rationale: Suboptimal vaccine immunogenicity and antigenic mismatch, compounded by poor uptake, means that influenza remains a major global disease. T cells recognizing peptides derived from conserved viral proteins could enhance vaccine-induced cross-strain protection. Objectives: To investigate the kinetics, phenotypes, and function of influenza virus-specific CD8+ resident memory T (Trm) cells in the lower airway and infer the molecular pathways associated with their response to infection in vivo. Methods: Healthy volunteers, aged 18-55, were inoculated intranasally with influenza A/California/4/09(H1N1). Blood, upper airway, and (in a subgroup) lower airway samples were obtained throughout infection. Symptoms were assessed by using self-reported diaries, and the nasal viral load was assessed by using quantitative PCR. T-cell responses were analyzed by using a three-color FluoroSpot assay, flow cytometry with MHC I-peptide tetramers, and RNA sequencing, with candidate markers being confirmed by using the immunohistochemistry results for endobronchial biopsy specimens. Measurements and Main Results: After challenge, 57% of participants became infected. Preexisting influenza-specific CD8+ T cells in blood correlated strongly with a reduced viral load, which peaked at Day 3. Influenza-specific CD8+ T cells in BAL fluid were highly enriched and predominantly expressed the Trm markers CD69 and CD103. Comparison between preinfection CD8+ T cells in BAL fluid and blood by using RNA sequencing revealed 3,928 differentially expressed genes, including all major Trm-cell markers. However, gene set enrichment analysis of BAL-fluid CD8+ T cells showed primarily innate cell-related pathways and, during infection, included upregulation of innate chemokines (Cxcl1, Cxcl10, and Cxcl16) that were also expressed by CD8+ cells in bronchial tissues. Conclusions: CD8+ Trm cells in the human lung display innate-like gene and protein expression that demonstrates blurred divisions between innate and adaptive immunity. Clinical study registered with www.clinicaltrials.gov (NCT02755948).

The potential of metabolic profiling for vaccine development.

Recent technological advances have provided deeper insights into the role of small molecules in biological processes. Metabolic profiling has thus entered the arena of -omics studies and rapidly proven its value both as stand-alone and as complement to other more advanced approaches, notably transcriptomics. Here we describe the potential of metabolic profiling for vaccinology embedded in the context of infection and immunity. This discussion is preceded by a description of the relevant technical and analytical tools for biological interpretation of metabolic data. Although not as widely applied as other -omics technologies, we believe that metabolic profiling can make important contributions to the better understanding of mechanisms underlying vaccine-induced responses and their effects on the prevention of infection or disease.

Heart-Specific Immune Responses in an Animal Model of Autoimmune-Related Myocarditis Mitigated by an Immunoproteasome Inhibitor and Genetic Ablation.

BACKGROUND: Immune checkpoint inhibitor (ICI) therapy is often accompanied by immune-related pathology, with an increasing occurrence of high-risk ICI-related myocarditis. Understanding the mechanisms involved in this side effect could enable the development of management strategies. In mouse models, immune checkpoints, such as PD-1 (programmed cell death protein 1), control the threshold of self-antigen responses directed against cardiac TnI (troponin I). We aimed to identify how the immunoproteasome, the main proteolytic machinery in immune cells harboring 3 distinct protease activities in the LMP2 (low-molecular-weight protein 2), LMP7 (low-molecular-weight protein 7), and MECL1 (multicatalytic endopeptidase complex subunit 1) subunit, affects TnI-directed autoimmune pathology of the heart. METHODS: TnI-directed autoimmune myocarditis (TnI-AM), a CD4+ T-cell-mediated disease, was induced in mice lacking all 3 immunoproteasome subunits (triple-ip-/-) or lacking either the gene encoding LMP2 and LMP7 by immunization with a cardiac TnI peptide. Alternatively, before induction of TnI-AM or after establishment of autoimmune myocarditis, mice were treated with the immunoproteasome inhibitor ONX 0914. Immune parameters defining heart-specific autoimmunity were investigated in experimental TnI-AM and in 2 cases of ICI-related myocarditis. RESULTS: All immunoproteasome-deficient strains showed mitigated autoimmune-related cardiac pathology with less inflammation, lower proinflammatory and chemotactic cytokines, less interleukin-17 production, and reduced fibrosis formation. Protection from TnI-directed autoimmune heart pathology with improved cardiac function in LMP7-/- mice involved a changed balance between effector and regulatory CD4+ T cells in the spleen, with CD4+ T cells from LMP7-/- mice showing a higher expression of inhibitory PD-1 molecules. Blocked immunoproteasome proteolysis, by treatment of TLR2 (Toll-like receptor 2)-engaged and TLR7 (Toll-like receptor 7)/TLR8 (Toll-like receptor 8)-engaged CD14+ monocytes with ONX 0914, diminished proinflammatory cytokine responses, thereby reducing the boost for the expansion of self-reactive CD4+ T cells. Correspondingly, in mice, ONX 0914 treatment reversed cardiac autoimmune pathology, preventing the induction and progression of TnI-AM when self-reactive CD4+ T cells were primed. The autoimmune signature during experimental TnI-AM, with high immunoproteasome expression, immunoglobulin G deposition, interleukin-17 production in heart tissue, and TnI-directed humoral autoimmune responses, was also present in 2 cases of ICI-related myocarditis, demonstrating the activation of heart-specific autoimmune reactions by ICI therapy. CONCLUSIONS: By reversing heart-specific autoimmune responses, immunoproteasome inhibitors applied to a mouse model demonstrate their potential to aid in the management of autoimmune myocarditis in humans, possibly including patients with ICI-related heart-specific autoimmunity.

Impact of Intermediate Hyperglycemia and Diabetes on Immune Dysfunction in Tuberculosis.

BACKGROUND: People with diabetes have an increased risk of developing active tuberculosis (TB) and are more likely to have poor TB-treatment outcomes, which may impact on control of TB as the prevalence of diabetes is increasing worldwide. Blood transcriptomes are altered in patients with active TB relative to healthy individuals. The effects of diabetes and intermediate hyperglycemia (IH) on this transcriptomic signature were investigated to enhance understanding of immunological susceptibility in diabetes-TB comorbidity. METHODS: Whole blood samples were collected from active TB patients with diabetes (glycated hemoglobin [HbA1c] ≥6.5%) or IH (HbA1c = 5.7% to <6.5%), TB-only patients, and healthy controls in 4 countries: South Africa, Romania, Indonesia, and Peru. Differential blood gene expression was determined by RNA-seq (n = 249). RESULTS: Diabetes increased the magnitude of gene expression change in the host transcriptome in TB, notably showing an increase in genes associated with innate inflammatory and decrease in adaptive immune responses. Strikingly, patients with IH and TB exhibited blood transcriptomes much more similar to patients with diabetes-TB than to patients with only TB. Both diabetes-TB and IH-TB patients had a decreased type I interferon response relative to TB-only patients. CONCLUSIONS: Comorbidity in individuals with both TB and diabetes is associated with altered transcriptomes, with an expected enhanced inflammation in the presence of both conditions, but also reduced type I interferon responses in comorbid patients, suggesting an unexpected uncoupling of the TB transcriptome phenotype. These immunological dysfunctions are also present in individuals with IH, showing that altered immunity to TB may also be present in this group. The TB disease outcomes in individuals with IH diagnosed with TB should be investigated further.

Prediction of anti-tuberculosis treatment duration based on a 22-gene transcriptomic model.

BACKGROUND: The World Health Organization recommends standardised treatment durations for patients with tuberculosis (TB). We identified and validated a host-RNA signature as a biomarker for individualised therapy durations for patients with drug-susceptible (DS)- and multidrug-resistant (MDR)-TB. METHODS: Adult patients with pulmonary TB were prospectively enrolled into five independent cohorts in Germany and Romania. Clinical and microbiological data and whole blood for RNA transcriptomic analysis were collected at pre-defined time points throughout therapy. Treatment outcomes were ascertained by TBnet criteria (6-month culture status/1-year follow-up). A whole-blood RNA therapy-end model was developed in a multistep process involving a machine-learning algorithm to identify hypothetical individual end-of-treatment time points. RESULTS: 50 patients with DS-TB and 30 patients with MDR-TB were recruited in the German identification cohorts (DS-GIC and MDR-GIC, respectively); 28 patients with DS-TB and 32 patients with MDR-TB in the German validation cohorts (DS-GVC and MDR-GVC, respectively); and 52 patients with MDR-TB in the Romanian validation cohort (MDR-RVC). A 22-gene RNA model (TB22) that defined cure-associated end-of-therapy time points was derived from the DS- and MDR-GIC data. The TB22 model was superior to other published signatures to accurately predict clinical outcomes for patients in the DS-GVC (area under the curve 0.94, 95% CI 0.9-0.98) and suggests that cure may be achieved with shorter treatment durations for TB patients in the MDR-GIC (mean reduction 218.0 days, 34.2%; p<0.001), the MDR-GVC (mean reduction 211.0 days, 32.9%; p<0.001) and the MDR-RVC (mean reduction of 161.0 days, 23.4%; p=0.001). CONCLUSION: Biomarker-guided management may substantially shorten the duration of therapy for many patients with MDR-TB.

Changes in Transcript, Metabolite, and Antibody Reactivity During the Early Protective Immune Response in Humans to Mycobacterium tuberculosis Infection.

BACKGROUND: Strategies to prevent Mycobacterium tuberculosis (Mtb) infection are urgently required. In this study, we aimed to identify correlates of protection against Mtb infection. METHODS: Two groups of Mtb-exposed contacts of tuberculosis (TB) patients were recruited and classified according to their Mtb infection status using the tuberculin skin test (TST; cohort 1) or QuantiFERON (QFT; cohort 2). A negative reading at baseline with a positive reading at follow-up classified TST or QFT converters and a negative reading at both time points classified TST or QFT nonconverters. Ribonucleic acid sequencing, Mtb proteome arrays, and metabolic profiling were performed. RESULTS: Several genes were found to be differentially expressed at baseline between converters and nonconverters. Gene set enrichment analysis revealed a distinct B-cell gene signature in TST nonconverters compared to converters. When infection status was defined by QFT, enrichment of type I interferon was observed. A remarkable area under the curve (AUC) of 1.0 was observed for IgA reactivity to Rv0134 and an AUC of 0.98 for IgA reactivity to both Rv0629c and Rv2188c. IgG reactivity to Rv3223c resulted in an AUC of 0.96 and was markedly higher compared to TST nonconverters. We also identified several differences in metabolite profiles, including changes in biomarkers of inflammation, fatty acid metabolism, and bile acids. Pantothenate (vitamin B5) was significantly increased in TST nonconverters compared to converters at baseline (q = 0.0060). CONCLUSIONS: These data provide new insights into the early protective response to Mtb infection and possible avenues to interfere with Mtb infection, including vitamin B5 supplementation.Analysis of blood from highly exposed household contacts from The Gambia who never develop latent Mycobacterium tuberculosis infection shows distinct transcriptomic, antibody, and metabolomic profiles compared to those who develop latent tuberculosis infection but prior to any signs of infection.

Gene set enrichment for reproducible science: comparison of CERNO and eight other algorithms.

MOTIVATION: Analysis of gene set (GS) enrichment is an essential part of functional omics studies. Here, we complement the established evaluation metrics of GS enrichment algorithms with a novel approach to assess the practical reproducibility of scientific results obtained from GS enrichment tests when applied to related data from different studies. RESULTS: We evaluated eight established and one novel algorithm for reproducibility, sensitivity, prioritization, false positive rate and computational time. In addition to eight established algorithms, we also included Coincident Extreme Ranks in Numerical Observations (CERNO), a flexible and fast algorithm based on modified Fisher P-value integration. Using real-world datasets, we demonstrate that CERNO is robust to ranking metrics, as well as sample and GS size. CERNO had the highest reproducibility while remaining sensitive, specific and fast. In the overall ranking Pathway Analysis with Down-weighting of Overlapping Genes, CERNO and over-representation analysis performed best, while CERNO and GeneSetTest scored high in terms of reproducibility. AVAILABILITY AND IMPLEMENTATION: tmod package implementing the CERNO algorithm is available from CRAN (cran.r-project.org/web/packages/tmod/index.html) and an online implementation can be found at http://tmod.online/. The datasets analyzed in this study are widely available in the KEGGdzPathwaysGEO, KEGGandMetacoreDzPathwaysGEO R package and GEO repository. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Mycobacterium tuberculosis infection modulates adipose tissue biology.

Mycobacterium tuberculosis (Mtb) primarily resides in the lung but can also persist in extrapulmonary sites. Macrophages are considered the prime cellular habitat in all tissues. Here we demonstrate that Mtb resides inside adipocytes of fat tissue where it expresses stress-related genes. Moreover, perigonadal fat of Mtb-infected mice disseminated the infection when transferred to uninfected animals. Adipose tissue harbors leukocytes in addition to adipocytes and other cell types and we observed that Mtb infection induces changes in adipose tissue biology depending on stage of infection. Mice infected via aerosol showed infiltration of inducible nitric oxide synthase (iNOS) or arginase 1 (Arg1)-negative F4/80+ cells, despite recruitment of CD3+, CD4+ and CD8+ T cells. Gene expression analysis of adipose tissue of aerosol Mtb-infected mice provided evidence for upregulated expression of genes associated with T cells and NK cells at 28 days post-infection. Strikingly, IFN-γ-producing NK cells and Mtb-specific CD8+ T cells were identified in perigonadal fat, specifically CD8+CD44-CD69+ and CD8+CD44-CD103+ subpopulations. Gene expression analysis of these cells revealed that they expressed IFN-γ and the lectin-like receptor Klrg1 and down-regulated CD27 and CD62L, consistent with an effector phenotype of Mtb-specific CD8+ T cells. Sorted NK cells expressed higher abundance of Klrg1 upon infection, as well. Our results reveal the ability of Mtb to persist in adipose tissue in a stressed state, and that NK cells and Mtb-specific CD8+ T cells infiltrate infected adipose tissue where they produce IFN-γ and assume an effector phenotype. We conclude that adipose tissue is a potential niche for Mtb and that due to infection CD8+ T cells and NK cells are attracted to this tissue.

Four-Gene Pan-African Blood Signature Predicts Progression to Tuberculosis.

Rationale: Contacts of patients with tuberculosis (TB) constitute an important target population for preventive measures because they are at high risk of infection with Mycobacterium tuberculosis and progression to disease.Objectives: We investigated biosignatures with predictive ability for incident TB.Methods: In a case-control study nested within the Grand Challenges 6-74 longitudinal HIV-negative African cohort of exposed household contacts, we employed RNA sequencing, PCR, and the pair ratio algorithm in a training/test set approach. Overall, 79 progressors who developed TB between 3 and 24 months after diagnosis of index case and 328 matched nonprogressors who remained healthy during 24 months of follow-up were investigated.Measurements and Main Results: A four-transcript signature derived from samples in a South African and Gambian training set predicted progression up to two years before onset of disease in blinded test set samples from South Africa, the Gambia, and Ethiopia with little population-associated variability, and it was also validated in an external cohort of South African adolescents with latent M. tuberculosis infection. By contrast, published diagnostic or prognostic TB signatures were predicted in samples from some but not all three countries, indicating site-specific variability. Post hoc meta-analysis identified a single gene pair, C1QC/TRAV27 (complement C1q C-chain / T-cell receptor-α variable gene 27) that would consistently predict TB progression in household contacts from multiple African sites but not in infected adolescents without known recent exposure events.Conclusions: Collectively, we developed a simple whole blood-based PCR test to predict TB in recently exposed household contacts from diverse African populations. This test has potential for implementation in national TB contact investigation programs.

Ranking metrics in gene set enrichment analysis: do they matter?

BACKGROUND: There exist many methods for describing the complex relation between changes of gene expression in molecular pathways or gene ontologies under different experimental conditions. Among them, Gene Set Enrichment Analysis seems to be one of the most commonly used (over 10,000 citations). An important parameter, which could affect the final result, is the choice of a metric for the ranking of genes. Applying a default ranking metric may lead to poor results. METHODS AND RESULTS: In this work 28 benchmark data sets were used to evaluate the sensitivity and false positive rate of gene set analysis for 16 different ranking metrics including new proposals. Furthermore, the robustness of the chosen methods to sample size was tested. Using k-means clustering algorithm a group of four metrics with the highest performance in terms of overall sensitivity, overall false positive rate and computational load was established i.e. absolute value of Moderated Welch Test statistic, Minimum Significant Difference, absolute value of Signal-To-Noise ratio and Baumgartner-Weiss-Schindler test statistic. In case of false positive rate estimation, all selected ranking metrics were robust with respect to sample size. In case of sensitivity, the absolute value of Moderated Welch Test statistic and absolute value of Signal-To-Noise ratio gave stable results, while Baumgartner-Weiss-Schindler and Minimum Significant Difference showed better results for larger sample size. Finally, the Gene Set Enrichment Analysis method with all tested ranking metrics was parallelised and implemented in MATLAB, and is available at https://github.com/ZAEDPolSl/MrGSEA . CONCLUSIONS: Choosing a ranking metric in Gene Set Enrichment Analysis has critical impact on results of pathway enrichment analysis. The absolute value of Moderated Welch Test has the best overall sensitivity and Minimum Significant Difference has the best overall specificity of gene set analysis. When the number of non-normally distributed genes is high, using Baumgartner-Weiss-Schindler test statistic gives better outcomes. Also, it finds more enriched pathways than other tested metrics, which may induce new biological discoveries.

Macrophage arginase-1 controls bacterial growth and pathology in hypoxic tuberculosis granulomas.

Lung granulomas develop upon Mycobacterium tuberculosis (Mtb) infection as a hallmark of human tuberculosis (TB). They are structured aggregates consisting mainly of Mtb-infected and -uninfected macrophages and Mtb-specific T cells. The production of NO by granuloma macrophages expressing nitric oxide synthase-2 (NOS2) via l-arginine and oxygen is a key protective mechanism against mycobacteria. Despite this protection, TB granulomas are often hypoxic, and bacterial killing via NOS2 in these conditions is likely suboptimal. Arginase-1 (Arg1) also metabolizes l-arginine but does not require oxygen as a substrate and has been shown to regulate NOS2 via substrate competition. However, in other infectious diseases in which granulomas occur, such as leishmaniasis and schistosomiasis, Arg1 plays additional roles such as T-cell regulation and tissue repair that are independent of NOS2 suppression. To address whether Arg1 could perform similar functions in hypoxic regions of TB granulomas, we used a TB murine granuloma model in which NOS2 is absent. Abrogation of Arg1 expression in macrophages in this setting resulted in exacerbated lung granuloma pathology and bacterial burden. Arg1 expression in hypoxic granuloma regions correlated with decreased T-cell proliferation, suggesting that Arg1 regulation of T-cell immunity is involved in disease control. Our data argue that Arg1 plays a central role in the control of TB when NOS2 is rendered ineffective by hypoxia.

Type I IFN signaling triggers immunopathology in tuberculosis-susceptible mice by modulating lung phagocyte dynamics.

General interest in the biological functions of IFN type I in Mycobacterium tuberculosis (Mtb) infection increased after the recent identification of a distinct IFN gene expression signature in tuberculosis (TB) patients. Here, we demonstrate that TB-susceptible mice lacking the receptor for IFN I (IFNAR1) were protected from death upon aerogenic infection with Mtb. Using this experimental model to mimic primary progressive pulmonary TB, we dissected the immune processes affected by IFN I. IFNAR1 signaling did not affect T-cell responses, but markedly altered migration of inflammatory monocytes and neutrophils to the lung. This process was orchestrated by IFNAR1 expressed on both immune and tissue-resident radioresistant cells. IFNAR1-driven TB susceptibility was initiated by augmented Mtb replication and in situ death events, along with CXCL5/CXCL1-driven accumulation of neutrophils in alveoli, followed by the discrete compartmentalization of Mtb in lung phagocytes. Early depletion of neutrophils rescued TB-susceptible mice to levels observed in mice lacking IFNAR1. We conclude that IFN I alters early innate events at the site of Mtb invasion leading to fatal immunopathology. These data furnish a mechanistic explanation for the detrimental role of IFN I in pulmonary TB and form a basis for understanding the complex roles of IFN I in chronic inflammation.

DoMosaics: software for domain arrangement visualization and domain-centric analysis of proteins.

UNLABELLED: DoMosaics is an application that unifies protein domain annotation, domain arrangement analysis and visualization in a single tool. It simplifies the analysis of protein families by consolidating disjunct procedures based on often inconvenient command-line applications and complex analysis tools. It provides a simple user interface with access to domain annotation services such as InterProScan or a local HMMER installation, and can be used to compare, analyze and visualize the evolution of domain architectures. AVAILABILITY AND IMPLEMENTATION: DoMosaics is licensed under theApache License, Version 2.0, and binaries can be freely obtained from www.domosaics.net.

Rapid similarity search of proteins using alignments of domain arrangements.

MOTIVATION: Homology search methods are dominated by the central paradigm that sequence similarity is a proxy for common ancestry and, by extension, functional similarity. For determining sequence similarity in proteins, most widely used methods use models of sequence evolution and compare amino-acid strings in search for conserved linear stretches. Probabilistic models or sequence profiles capture the position-specific variation in an alignment of homologous sequences and can identify conserved motifs or domains. While profile-based search methods are generally more accurate than simple sequence comparison methods, they tend to be computationally more demanding. In recent years, several methods have emerged that perform protein similarity searches based on domain composition. However, few methods have considered the linear arrangements of domains when conducting similarity searches, despite strong evidence that domain order can harbour considerable functional and evolutionary signal. RESULTS: Here, we introduce an alignment scheme that uses a classical dynamic programming approach to the global alignment of domains. We illustrate that representing proteins as strings of domains (domain arrangements) and comparing these strings globally allows for a both fast and sensitive homology search. Further, we demonstrate that the presented methods complement existing methods by finding similar proteins missed by popular amino-acid-based comparison methods. AVAILABILITY: An implementation of the presented algorithms, a web-based interface as well as a command-line program for batch searching against the UniProt database can be found at http://rads.uni-muenster.de. Furthermore, we provide a JAVA API for programmatic access to domain-string­based search methods.

Common patterns and disease-related signatures in tuberculosis and sarcoidosis.

In light of the marked global health impact of tuberculosis (TB), strong focus has been on identifying biosignatures. Gene expression profiles in blood cells identified so far are indicative of a persistent activation of the immune system and chronic inflammatory pathology in active TB. Definition of a biosignature with unique specificity for TB demands that identified profiles can differentiate diseases with similar pathology, like sarcoidosis (SARC). Here, we present a detailed comparison between pulmonary TB and SARC, including whole-blood gene expression profiling, microRNA expression, and multiplex serum analytes. Our analysis reveals that previously disclosed gene expression signatures in TB show highly similar patterns in SARC, with a common up-regulation of proinflammatory pathways and IFN signaling and close similarity to TB-related signatures. microRNA expression also presented a highly similar pattern in both diseases, whereas cytokines in the serum of TB patients revealed a slightly elevated proinflammatory pattern compared with SARC and controls. Our results indicate several differences in expression between the two diseases, with increased metabolic activity and significantly higher antimicrobial defense responses in TB. However, matrix metallopeptidase 14 was identified as the most distinctive marker of SARC. Described communalities as well as unique signatures in blood profiles of two distinct inflammatory pulmonary diseases not only have considerable implications for the design of TB biosignatures and future diagnosis, but they also provide insights into biological processes underlying chronic inflammatory disease entities of different etiology.