Tianhua Niu

The DosR Regulon Modulates Adaptive Immunity and Is Essential for Mycobacterium tuberculosis Persistence

RATIONALE Hypoxia promotes dormancy by causing physiologic changes to actively replicating Mycobacterium tuberculosis. DosR controls the response of M. tuberculosis to hypoxia. OBJECTIVES To understand DosRs contribution in the persistence of M. tuberculosis, we compared the phenotype of various DosR regulon mutants and a complemented strain to M. tuberculosis in macaques, which faithfully model M. tuberculosis infection. METHODS We measured clinical and microbiologic correlates of infection with M. tuberculosis relative to mutant/complemented strains in the DosR regulon, studied lung pathology and hypoxia, and compared immune responses in lung using transcriptomics and flow cytometry. MEASUREMENTS AND MAIN RESULTS Despite being able to replicate initially, mutants in DosR regulon failed to persist or cause disease. On the contrary, M. tuberculosis and a complemented strain were able to establish infection and tuberculosis. The attenuation of pathogenesis in animals infected with the mutants coincided with the appearance of a Th1 response and organization of hypoxic lesions wherein M. tuberculosis expressed dosR. The lungs of animals infected with the mutants (but not the complemented strain) exhibited early transcriptional signatures of T-cell recruitment, activation, and proliferation associated with an increase of T cells expressing homing and proliferation markers. CONCLUSIONS Delayed adaptive responses, a hallmark of M. tuberculosis infection, not only lead to persistence but also interfere with the development of effective antituberculosis vaccines. The DosR regulon therefore modulates both the magnitude and the timing of adaptive immune responses in response to hypoxia in vivo, resulting in persistent infection. Hence, DosR regulates key aspects of the M. tuberculosis life cycle and limits lung pathology.

DosS Is Required for the Complete Virulence of Mycobacterium tuberculosis in Mice with Classical Granulomatous Lesions

Mycobacterium tuberculosis (Mtb) must counter hypoxia within granulomas to persist. DosR, in concert with sensor kinases DosS and DosT, regulates the response to hypoxia. Yet Mtb lacking functional DosR colonize the lungs of C57Bl/6 mice, presumably owing to the lack of organized lesions with sufficient hypoxia in that model. We compared the phenotype of the Δ-dosR, Δ-dosS, and Δ-dosT mutants to Mtb using C3HeB/FeJ mice, an alternate mouse model where lesions develop hypoxia. C3HeB/FeJ mice were infected via aerosol. The progression of infection was analyzed by tissue bacterial burden and histopathology. A measure of the comparative global immune responses was also analyzed. Although Δ-dosR and Δ-dosT grew comparably to wild-type Mtb, Δ-dosS exhibited a significant defect in bacterial burden and pathology in vivo, accompanied by ablated proinflammatory response. Δ-dosS retained the ability to induce DosR. The Δ-dosS mutant was also attenuated in murine macrophages ex vivo, with evidence of reduced expression of the proinflammatory signature. Our results show that DosS, but not DosR and DosT, is required by Mtb to survive in C3HeB/FeJ mice. The attenuation of Δ-dosS is not due to its inability to induce the DosR regulon, nor is it a result of the accumulation of hypoxia. That the in vivo growth restriction of Δ-dosS could be mimicked ex vivo suggested sensitivity to macrophage oxidative burst. Anoxic caseous centers within tuberculosis lesions eventually progress to cavities. Our results provide greater insight into the molecular mechanisms of Mtb persistence within host lungs.

Role of TNF in the Altered Interaction of Dormant Mycobacterium tuberculosis with Host Macrophages

Mycobacterium tuberculosis (Mtb) persists within lung granulomas, despite being subjected to diverse stress conditions, including hypoxia. We hypothesized that the response of host phagocytes to Mtb experiencing hypoxia is radically altered and designed in vitro experiment to study this phenomenon. Hypoxia-stressed (Mtb-H) and aerobically grown Mtb (Mtb-A) were used to infect Rhesus Macaque Bone Marrow Derived Macrophages (Rh-BMDMs) and the comparative host response to Mtb infection studied. Mechanistic insights were gained by employing RNAi. Mtb-H accumulated significantly lower bacterial burden during growth in Rh-BMDMs, concomitantly generating a drastically different host transcriptional profile (with only <2% of all genes perturbed by either infection being shared between the two groups). A key component of this signature was significantly higher TNF and apopotosis in Mtb-H- compared to Mtb-A-infected Rh-BMDMs. Silencing of TNF by RNAi reversed the significant control of Mtb replication. These results indicate a potential mechanism for the rapid clearance of hypoxia-conditioned bacilli by phagocytes. In conclusion, hypoxia-conditioned Mtb undergo significantly different interactions with host macrophages compared to Mtb grown in normoxia. These interactions result in the induction of the TNF signaling pathway, activation of apoptosis, and DNA-damage stress response. Our results show that Mtb-H bacilli are particularly susceptible to killing governed by TNF.

Identification of IDUA and WNT16 Phosphorylation-Related Non-Synonymous Polymorphisms for Bone Mineral Density in Meta-Analyses of Genome-Wide Association Studies

Protein phosphorylation regulates a wide variety of cellular processes. Thus, we hypothesize that single‐nucleotide polymorphisms (SNPs) that may modulate protein phosphorylation could affect osteoporosis risk. Based on a previous conventional genome‐wide association (GWA) study, we conducted a three‐stage meta‐analysis targeting phosphorylation‐related SNPs (phosSNPs) for femoral neck (FN)‐bone mineral density (BMD), total hip (HIP)‐BMD, and lumbar spine (LS)‐BMD phenotypes. In stage 1, 9593 phosSNPs were meta‐analyzed in 11,140 individuals of various ancestries. Genome‐wide significance (GWS) and suggestive significance were defined by α = 5.21 × 10–6 (0.05/9593) and 1.00 × 10–4, respectively. In stage 2, nine stage 1–discovered phosSNPs (based on α = 1.00 × 10–4) were in silico meta‐analyzed in Dutch, Korean, and Australian cohorts. In stage 3, four phosSNPs that replicated in stage 2 (based on α = 5.56 × 10–3, 0.05/9) were de novo genotyped in two independent cohorts. IDUA rs3755955 and rs6831280, and WNT16 rs2707466 were associated with BMD phenotypes in each respective stage, and in three stages combined, achieving GWS for both FN‐BMD (p = 8.36 × 10–10, p = 5.26 × 10–10, and p = 3.01 × 10–10, respectively) and HIP‐BMD (p = 3.26 × 10–6, p = 1.97 × 10–6, and p = 1.63 × 10–12, respectively). Although in vitro studies demonstrated no differences in expressions of wild‐type and mutant forms of IDUA and WNT16B proteins, in silico analyses predicts that WNT16 rs2707466 directly abolishes a phosphorylation site, which could cause a deleterious effect on WNT16 protein, and that IDUA phosSNPs rs3755955 and rs6831280 could exert indirect effects on nearby phosphorylation sites. Further studies will be required to determine the detailed and specific molecular effects of these BMD‐associated non‐synonymous variants.

Androgen receptor-deficient islet β-cells exhibit alteration in genetic markers of insulin secretion and inflammation. A transcriptome analysis in the male mouse

AIMS Testosterone action is mediated via the androgen receptor (AR). We have reported that male mice lacking AR selectively in β-cells (βARKO-/y) develop decreased glucose-stimulated insulin secretion (GSIS), producing glucose intolerance. We showed that testosterone action on AR in β-cells amplifies the insulinotropic action of GLP-1 on its receptor via a cAMP-dependent protein kinase-A pathway. METHODS To investigate AR-dependent gene networks in β-cells, we performed a high throughput whole transcriptome sequencing (RNA-Seq) in islets from male βARKO-/y and control mice. RESULTS We identified 214 differentially expressed genes (DEGs) (158 up- and 56 down-regulated) with a false discovery rate (FDR) < 0.05 and a fold change (FC) > 2. Our analysis of individual transcripts revealed alterations in β-cell genes involved in cellular inflammation/stress and insulin secretion. Based on 312 DEGs with an FDR < 0.05, the pathway analysis revealed 23 significantly enriched pathways, including cytokine-cytokine receptor interaction, Jak-STAT signaling, insulin signaling, MAPK signaling, type 2 diabetes (T2D) and pancreatic secretion. The gene ontology analysis confirmed the results of the individual DEGs and the pathway analysis in showing enriched biological processes encompassing inflammation, ion transport, exocytosis and insulin secretion. CONCLUSIONS AR-deficient islets exhibit altered expression of genes involved in inflammation and insulin secretion demonstrating the importance of androgen action in β-cell health in the male with implications for T2D development in men.

Network-Based Meta-Analyses of Associations of Multiple Gene Expression Profiles with Bone Mineral Density Variations in Women

Background Existing microarray studies of bone mineral density (BMD) have been critical for understanding the pathophysiology of osteoporosis, and have identified a number of candidate genes. However, these studies were limited by their relatively small sample sizes and were usually analyzed individually. Here, we propose a novel network-based meta-analysis approach that combines data across six microarray studies to identify functional modules from human protein-protein interaction (PPI) data, and highlight several differentially expressed genes (DEGs) and a functional module that may play an important role in BMD regulation in women. Methods Expression profiling studies were identified by searching PubMed, Gene Expression Omnibus (GEO) and ArrayExpress. Two meta-analysis methods were applied across different gene expression profiling studies. The first, a nonparametric Fisher’s method, combined p-values from individual experiments to identify genes with large effect sizes. The second method combined effect sizes from individual datasets into a meta-effect size to gain a higher precision of effect size estimation across all datasets. Genes with Q test’s p-values < 0.05 or I2 values > 50% were assessed by a random effects model and the remainder by a fixed effects model. Using Fisher’s combined p-values, functional modules were identified through an integrated analysis of microarray data in the context of large protein–protein interaction (PPI) networks. Two previously published meta-analysis studies of genome-wide association (GWA) datasets were used to determine whether these module genes were genetically associated with BMD. Pathway enrichment analysis was performed with a hypergeometric test. Results Six gene expression datasets were identified, which included a total of 249 (129 high BMD and 120 low BMD) female subjects. Using a network-based meta-analysis, a consensus module containing 58 genes (nodes) and 83 edges was detected. Pathway enrichment analysis of the 58 module genes revealed that these genes were enriched in several important KEGG pathways including Osteoclast differentiation, B cell receptor signaling pathway, MAPK signaling pathway, Chemokine signaling pathway and Insulin signaling pathway. The importance of module genes was replicated by demonstrating that most module genes were genetically associated with BMD in the GWAS data sets. Meta-analyses were performed at the individual gene level by combining p-values and effect sizes. Five candidate genes (ESR1, MAP3K3, PYGM, RAC1 and SYK) were identified based on gene expression meta-analysis, and their associations with BMD were also replicated by two BMD meta-analysis studies. Conclusions In summary, our network-based meta-analysis not only identified important differentially expressed genes but also discovered biologically meaningful functional modules for BMD determination. Our study may provide novel therapeutic targets for osteoporosis in women.

KCNJ11 , ABCC8 and TCF7L2 polymorphisms and the response to sulfonylurea treatment in patients with type 2 diabetes: a bioinformatics assessment

BackgroundType 2 diabetes (T2D) is a worldwide epidemic with considerable health and economic consequences. Sulfonylureas are widely used drugs for the treatment of patients with T2D. KCNJ11 and ABCC8 encode the Kir6.2 (pore-forming subunit) and SUR1 (regulatory subunit that binds to sulfonylurea) of pancreatic β cell KATP channel respectively with a critical role in insulin secretion and glucose homeostasis. TCF7L2 encodes a transcription factor expressed in pancreatic β cells that regulates insulin production and processing. Because mutations of these genes could affect insulin secretion stimulated by sulfonylureas, the aim of this study is to assess associations between molecular variants of KCNJ11, ABCC8 and TCF7L2 genes and response to sulfonylurea treatment and to predict their potential functional effects.MethodsBased on a comprehensive literature search, we found 13 pharmacogenetic studies showing that single nucleotide polymorphisms (SNPs) located in KCNJ11: rs5219 (E23K), ABCC8: rs757110 (A1369S), rs1799854 (intron 15, exon 16 -3C/T), rs1799859 (R1273R), and TCF7L2: rs7903146 (intron 4) were significantly associated with responses to sulfonylureas. For in silico bioinformatics analysis, SIFT, PolyPhen-2, PANTHER, MutPred, and SNPs3D were applied for functional predictions of 36 coding (KCNJ11: 10, ABCC8: 24, and TCF7L2: 2; all are missense), and HaploReg v4.1, RegulomeDB, and Ensembl’s VEP were used to predict functions of 7 non-coding (KCNJ11: 1, ABCC8: 1, and TCF7L2: 5) SNPs, respectively.ResultsBased on various in silico tools, 8 KCNJ11 missense SNPs, 23 ABCC8 missense SNPs, and 2 TCF7L2 missense SNPs could affect protein functions. Of them, previous studies showed that mutant alleles of 4 KCNJ11 missense SNPs and 5 ABCC8 missense SNPs can be successfully rescued by sulfonylurea treatments. Further, 3 TCF7L2 non-coding SNPs (rs7903146, rs11196205 and rs12255372), can change motif(s) based on HaploReg v4.1 and are predicted as risk factors by Ensembl’s VEP.ConclusionsOur study indicates that a personalized medicine approach by tailoring sulfonylurea therapy of T2D patients according to their genotypes of KCNJ11, ABCC8, and TCF7L2 could attain an optimal treatment efficacy.

LAG-3 potentiates the survival of Mycobacterium tuberculosis in host phagocytes by modulating mitochondrial signaling in an in-vitro granuloma model

CD4+ T-cell mediated Th1 immune responses are critical for immunity to TB. The immunomodulatory protein, lymphocyte activation gene-3 (LAG-3) decreases Th1-type immune responses in T-cells. LAG-3 expression is significantly induced in the lungs of macaques with active TB and correlates with increased bacterial burden. Overproduction of LAG-3 can greatly diminish responses and could lead to uncontrolled Mtb replication. To assess the effect of LAG-3 on the progression of Mtb infection, we developed a co-culture system wherein blood-derived macrophages are infected with Mtb and supplemented with macaque blood or lung derived CD4+ T-cells. Silencing LAG-3 signaling in macaque lung CD4+ T-cells enhanced killing of Mtb in co-cultures, accompanied by reduced mitochondrial electron transport and increased IFN-γ expression. Thus, LAG-3 may modulate adaptive immunity to Mtb infection by interfering with the mitochondrial apoptosis pathway. Better understanding this pathway could allow us to circumvent immune features that promote disease.

Sequencing-relative to hybridization-based transcriptomics approaches better define Mycobacterium tuberculosis stress-response regulons

Mycobacterium tuberculosis (Mtb) infections cause tuberculosis (TB), an infectious disease which causes ∼1.5 million deaths annually. The ability of this pathogen to evade, escape and encounter immune surveillance is fueled by its adaptability. Thus, Mtb induces a transition in its transcriptome in response to environmental changes. Global transcriptome profiling has been key to our understanding of how Mtb responds to the different stress conditions it faces during its life cycle. While this was initially achieved using microarray technology, RNAseq is now widely employed. It is important to understand the correlation between the large amount of microarray based transcriptome data, which continues to shape our understanding of Mtb stress networks, and newer data being generated using RNAseq. We assessed how well the two platforms correlate using three well-defined stress conditions: diamide, hypoxia, and re-aeration. The data used here was generated by different individuals over time using distinct samples, providing a stringent test of platform correlation. While correlation between microarrays and sequencing was high upon diamide treatment, which causes a rapid reprogramming of the transcriptome, RNAseq allowed a better definition of the hypoxic response, characterized by subtle changes in the magnitude of gene-expression. RNAseq also allows for the best cross-platform reproducibility.