Abstract

One-fourth of the human population is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) and carries the infection in its latent form. This latent infection presents a lifelong risk of developing active tuberculosis (TB) disease, and persons with latent TB infection (LTBI) are significant contributors to the pool of active TB cases. Genetic polymorphisms among hosts have been shown to contribute to the outcome of Mtb infection. The SP110 gene, which encodes an interferon-induced nuclear protein, has been shown to control host innate immunity to Mtb infection. In this study, we provide experimental data demonstrating the ability of the gene to control genetic susceptibility to latent and active TB infection. Genetic variants of the SP110 gene were investigated in the Taiwanese population (including 301 pulmonary TB patients, 68 LTBI individuals, and 278 healthy household contacts of the TB patients), and their association with susceptibility to latent and active TB infection was examined by performing an association analysis in a case-control study. We identified several SNPs (rs7580900, rs7580912, rs9061, rs11556887, and rs2241525) in the SP110 gene that are associated with susceptibility to LTBI and/or TB disease. Our studies further showed that the same SNPs may have opposite effects on the control of susceptibility to LTBI versus TB. In addition, our analyses demonstrated that the SP110 rs9061 SNP was associated with tumor necrosis factor-α (TNFα) levels in plasma in LTBI subjects. The results suggest that the polymorphisms within SP110 have a role in controlling genetic susceptibility to latent and active TB infection in humans. To the best of our knowledge, this is the first report showing that the SP110 variants are associated with susceptibility to LTBI. Our study also demonstrated that the identified SP110 SNPs displayed the potential to predict the risk of LTBI and subsequent TB progression in Taiwan.

1. Introduction

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb) infection, remains one of the top ten causes of death in the world [1]. Approximately one-fourth of the world population has been infected with Mtb [2], but only 10% of these infected persons develop progressive disease during their lifetimes [3]. The majority of infected individuals remain healthy and noninfectious but carry Mtb in a latent form. Latent TB infection (LTBI) is a state during which a persistent host immune response to stimulation by Mtb antigens is sustained without evidence of clinically manifested active TB [4, 5]. As many as 10% of people with LTBI will go on to develop progressive disease in the near or remote future (a process named “TB reactivation”), and the risk is significantly higher in the presence of predisposing factors, such as coinfection with human immunodeficiency virus (HIV).

Persons with LTBI can be diagnosed by skin (tuberculin skin test (TST)) and/or blood (interferon-gamma (IFNγ) release assay (IGRA)) tests [6]. A positive result from these assays indicates an immune response to stimulation by Mtb antigens in the LTBI population, despite the fact that these individuals have a negative bacteriological test. However, neither TST nor IGRA can distinguish between LTBI and active TB disease or predict who will progress to active TB [7]. Given that the LTBI individuals represent a potential reservoir for future active TB cases, preventive therapy for LTBI is as important a goal as timely anti-TB treatment to reduce the burden of TB [8]. Therefore, identifying and treating cases with LTBI will contribute to TB elimination. Fundamental research for the development of diagnostic assays with improved performance and predictive assessment for TB reactivation will have practical applications and offer a substantial benefit for LTBI management.

In both humans and experimental animal models, genetic polymorphisms among hosts have been shown to contribute to the outcome of Mtb infection [913]. In mice, the Ipr1 (intracellular pathogen resistance 1) gene is located within the sst1 (supersusceptibility to tuberculosis 1) locus on chromosome 1 (49-54 cM) [14] and has been identified as a genetic determinant conferring host innate immunity to Mtb infection [15]. The previous studies indicate that the Ipr1 gene may function to integrate mechanisms on controlling cell death, innate immunity, and pathogenesis during intracellular pathogen infections [15, 16]. The gene orthologous to the mouse Ipr1 in humans is SP110, located on chromosome 2q37.1. Expression of both Ipr1 and SP110 genes is intensively regulated by IFNs, suggesting that the function of both genes is related to the IFN-mediated immune response [17].

Genetic defects in the SP110 gene have been found to be responsible for hepatic veno-occlusive disease and immunodeficiency [18, 19], indicating that the gene plays important roles in immunity [20]. The SP110 gene encodes the SP110 nuclear body protein, which has at least three isoforms, including the dominantly expressed SP110a, b, and c isoforms that are believed to be the result of alternative mRNA splicing. Our recent study demonstrated that SP110b, which is most similar to mouse Ipr1 and is expressed more abundantly than SP110a and SP110c, modulates nuclear factor-κB (NF-κB) activity resulting in the downregulation of tumor necrosis factor-α (TNFα) production and concomitant upregulation of NF-κB-induced antiapoptotic gene expression, thereby suppressing IFNγ-mediated monocyte/macrophage cell death [21]. This indicates that the protein is crucial in the control of the activation of macrophages, the reservoir for Mtb persistence.

Although a number of genetic variants of the SP110 gene have been reported to be associated with susceptibility to human TB, the results of studies regarding the relationship between SP110 polymorphisms and TB susceptibility are inconsistent [2229]. A family-based study in West Africa identified 3 SP110 polymorphisms that are associated with TB susceptibility [22]; however, no significant associations between SP110 and disease susceptibility were identified by other, larger case-control studies conducted on various populations [2426]. After screening Taiwanese populations for polymorphisms in SP110, we identified some single-nucleotide polymorphisms (SNPs) in SP110 that are significantly associated with susceptibility to LTBI as well as TB disease. These results suggest that the SP110 variants may provide novel predictive markers for TB infection status and disease outcome.

2. Materials and Methods

2.1. Human Subject Study

The study was conducted in accordance with the terms of the informed consent that was provided to, and received from, participants prior to inclusion in the study. This study was approved by the National Taiwan University Hospital Institutional Review Board (IRB No. 200612009M and IRB No. 201512169RINA). Human blood was ethically collected from patients with culture-confirmed pulmonary TB and their household contacts as described previously [30]. Briefly, the participating contacts received chest radiography, and mycobacteriology studies were conducted for 3 sputum samples (including acid-fast smear and mycobacterial culture) to exclude the possibility of active TB disease. Because routine BCG vaccination for newborns in Taiwan could affect the accuracy of tuberculin skin test, all enrolled contacts were then tested for LTBI using a T-SPOT.TB assay (Oxford Immunotec Ltd., Abingdon, UK) or QuantiFERON-TB Gold In-Tube assay (QFT) (Qiagen, Hilden, Germany), and the assays were interpreted according to the manufacturers’ criteria. Both patients and contacts were excluded if they were tested positive for HIV infection. In total, 301 pulmonary TB patients, 68 individuals with LTBI, and 278 healthy household contacts were included in the study, and genomic DNA was extracted from their peripheral blood (1-2 mL) using a kit from Qiagen according to the manufacturer’s protocol. SP110 polymorphisms were identified from the National Center for Biotechnology Information dbSNP database (https://www.ncbi.nlm.nih.gov/snp). The SNPs were genotyped using the MassARRAY System (Sequenom, San Diego, CA, US), and the primer extension products were analyzed by MALDI-TOF mass spectrometry as previously described [31, 32]. Details of the primers that were used are listed in Table S1 in the Supplementary Materials.

Plasma samples were prepared from blood samples by centrifugation and then stored at −80°C until analysis. The TNFα levels in plasma samples were determined by a MAGPIX® platform (Luminex Corp., Austin, TX. US) with a MILLIPLEX MAP Human Cytokine/Chemokine Magnetic Bead Panel I kit (Merck Millipore, Billerica, MA, US) according to the manufacturer’s instructions.

2.2. Statistical Analyses

The associations of gene polymorphisms with LTBI and TB disease were analyzed by SAS 9.4 software (SAS Institute, Cary, NC, US) [33]. Linkage disequilibrium and haplotype analyses were performed using Haploview (https://www.broadinstitute.org/haploview) [34]. Chi-square tests were used to compare frequencies. Odds ratios (ORs) were calculated with 95% confidence intervals (CIs). Bonferroni correction or false discovery rate (FDR) correction was applied for multiple comparison adjustments as indicated. The difference of the TNFα levels in plasma between two genotype groups of samples was calculated using a two-tailed unpaired -test with Prism software (GraphPad, San Diego, CA, US). values less than 0.05 were considered to indicate statistical significance, and the number of asterisks represents the degree of significance with regard to the values.

3. Results

3.1. Association of Polymorphisms in SP110 with LTBI and TB Susceptibility

To investigate the association between the SP110 gene and control of Mtb infection, we examined polymorphisms in the gene of members of the Taiwanese population for genetic association with TB disease status. In this study, 301 pulmonary TB patients (202 males and 99 females; mean age: 63.1 ± 19.9 years), 68 individuals with LTBI (35 males and 33 females; mean age: 46.9 ± 17.7 years), and 278 healthy household contacts of the patients (92 males and 186 females; mean age: 47.1 ± 17.3 years) were included (Table 1). In total, 20 SNPs in the SP110 gene were selected for analysis. Of these, 10 SNPs that were not polymorphic or had a minor allele frequency less than 1% were not further included in the analysis, and the remaining 10 SNPs were then analyzed. We found that 3 SNPs (rs7580912, , OR: 1.52, 95% CI: 1.01–2.39; rs7580900, , OR: 1.68, 95% CI: 1.14–2.48; and rs9061, , OR: 0.39, 95% CI: 0.21–0.73) showed differential allele frequency distributions in LTBI cases vs. healthy controls (Table 2). After Bonferroni correction, rs9061 remained significant in this analysis (). Although allele frequencies of none of the 10 SNPs differed significantly in TB patients vs. healthy controls (Table 2), 3 SNPs (rs7580900, , OR: 0.66, 95% CI: 0.45–0.97; rs9061, , OR: 2.21, 95% CI: 1.18–4.15; and rs2241525, , OR: 0.6, 95% CI: 0.38–0.96) exhibited differential allele frequency distributions in TB cases vs. LTBI individuals (Table 2). All of these SNPs were in accordance with the Hardy-Weinberg equilibrium (HWE).

The associations between SP110 genotypes and susceptibility to LTBI and TB were then analyzed. In LTBI cases vs. healthy controls, we found that genotypes “GG” in rs7580912 (, OR: 2.451, 95% CI: 1.12–5.364) and “GG” in rs7580900 (, OR: 2.584, 95% CI: 1.208–5.53) were associated with LTBI risk and that genotype “GA” in rs9061 exhibited a protective effect on LTBI (, OR: 0.494, 95% CI: 0.239–0.981) (Table 3). We also found that genotypes “GG” in rs7580912 (, OR: 0.392, 95% CI: 0.179–0.86) and “GG” in rs7580900 (, OR: 0.392, 95% CI: 0.182–0.848) in TB patients vs. LTBI cases (Table 3), as well as “CT” in rs11556887 (, OR: 0.626, 95% CI: 0.401–0.976) in TB patients vs. healthy controls had a protective effect on TB (Table 3). These results indicated that several SNPs (rs7580900, rs7580912, rs9061, and rs11556887) in SP110 were associated with susceptibility to LTBI and/or TB disease.

3.2. Association Analyses in Various Inheritance Models

The minor allele of each SNP was presumed as a risk factor compared to the major allele, and the associations between SNPs and susceptibility to LTBI and TB were analyzed in various inheritance models (Table 4). We found that in LTBI cases vs. healthy controls, rs9061 showed a protective effect on LTBI in both additive (, OR: 0.41, 95% CI: 0.22–0.78) and dominant (, OR: 0.42, 95% CI: 0.21–0.82) models, while rs7580900 was associated with LTBI risk in both additive (, OR: 1.62, 95% CI: 1.1–2.39) and recessive (, OR: 2.13, 95% CI: 1.13–4.03) models. In addition to this finding, rs7580912 was also associated with LTBI risk in a recessive model (, OR: 2.49, 95% CI: 1.2–5.18). In TB patients vs. LTBI cases, rs9061 was associated with TB risk in both additive (, OR: 2.17, 95% CI: 1.16–4.06) and dominant (, OR: 2.14, 95% CI: 1.09–4.22) models, while rs7580900 exhibited a protective effect on TB in both additive (, OR: 0.64, 95% CI: 0.43–0.95) and recessive (, OR: 0.45, 95% CI: 0.24–0.84) models and rs2241525 was associated with protection from TB in both additive (, OR: 0.62, 95% CI: 0.4–0.98) and dominant (, OR: 0.55, 95% CI: 0.31–0.97) models. Additionally, rs7580912 had a protective effect on TB in a recessive model (, OR: 0.35, 95% CI: 0.17–0.72). After false discovery rate (FDR) correction, rs7580900 and rs7580912 remained significant ( and 0.0301, respectively) in a recessive model in this analysis.

3.3. Linkage Disequilibrium and Haplotype Analyses

We then examined linkage disequilibrium (LD) with the SNP markers in the SP110 gene using a Haploview analysis. We found that the “ATATACGCGG” and “ATGTACGCGA” haplotypes met statistical significance (, OR: 2.11, 95% CI: 1.11–3.99 and , OR: 2.76, 95% CI: 1.12–6.84, respectively) for association with LTBI risk in LTBI cases vs. healthy controls (see Figure S1 in the Supplementary Materials) and that the “ATGTACGCGA” and “ATGAAAGCGA” haplotypes were statistically significant (, OR: 0.38, 95% CI: 0.16–0.95 and , OR: 0.24, 95% CI: 0.06–0.99, respectively) with a protective effect on TB in TB cases vs. LTBI individuals (see Figure S2 in the Supplementary Materials). Interestingly, the “GCGTACGCGG” haplotype was associated with TB risk (, OR: 3.81, 95% CI: 1.18–12.31) in TB cases vs. healthy controls (see Figure S3 in the Supplementary Materials), although none of the SNPs studied show a significant difference in frequency distributions in this comparison (Table 2). Noteworthily, when we compared “ATGTACGCAA” (the most frequent haplotype) with “ATGTACGCGA” (the haplotype that was associated with LTBI risk in LTBI cases vs. healthy controls and had a protective effect on TB in TB cases vs. LTBI individuals), the latter was found to be significantly affected by the rs7580900 SNP in both comparisons (, OR: 2.76, 95% CI: 1.12–6.84 for LTBI cases vs. healthy controls; , OR: 0.38, 95% CI: 0.16–0.95 for TB cases vs. LTBI individuals) (see Figures S1 and S2 in the Supplementary Materials).

We analyzed a block that includes 3 SNPs (rs9061, rs7580900, and rs7580912) in two comparisons (LTBI cases vs. healthy controls and TB cases vs. LTBI individuals) (Figure 1). We found that the “GAA” (, OR: 1.77, 95% CI: 1.20–2.60), “GGG” (, OR: 2.06, 95% CI: 1.34–3.17), “GGA” (, OR: 2.81, 95% CI: 1.49–5.27), and “GAG” (, OR: 10, 95% CI: 2.36–42.32) haplotypes were associated with disease risk and that the “AGG” haplotype had a protective effect on LTBI (, OR: 0.41, 95% CI: 0.18–0.92) for LTBI cases vs. healthy controls (Table 5). In addition, the “GAA” (, OR: 0.61, 95% CI: 0.42–0.90), “GGG” (, OR: 0.65, 95% CI: 0.43–0.98), “GGA” (, OR: 0.3, 95% CI: 0.16–0.57), and “GAG” (, OR: 0.1, 95% CI: 0.02–0.4) haplotypes showed a protective effect on TB disease for TB patients vs. LTBI individuals (Table 5). These genetic studies suggest that the SP110 gene plays a key role in modulating susceptibility to latent and active TB infection.

3.4. Association between the SP110 rs9061 SNP and the TNFα Production in LTBI Subjects

TNFα plays a crucial role in controlling Mtb infection and TB reactivation; however, overproduction of TNFα may cause pathology [35, 36]. Our previous studies demonstrated that the SP110b protein, which is encoded by the SP110 gene and whose expression is upregulated by IFNs, downregulates TNFα production in monocyte/macrophage cells activated by IFNγ, thereby alleviating cell death [21]. This finding indicates that the protein functions as a regulator of proinflammatory cytokines of host immunity contributing to a reduction in tissue damage caused by excessive inflammation [21]. To further investigate the potential association between the SP110 SNPs and disease status, we next analyzed clinical parameters in the studied subjects. We measured TNFα levels in plasma from LTBI individuals who carry the different genotypes of the SP110 SNPs and demonstrated that the “GA” genotype of rs9061 in LTBI individuals was associated with lower TNFα levels in plasma compared to “GG” LTBI subjects (Figure 2). The TNFα levels of healthy controls with both genotypes were undetectable (not shown). These data are in agreement with the protective role of the “GA” genotype of rs9061 in LTBI (Table 3) and further support our recent finding showing that the SP110b protein prevents cell death and tissue damage by downregulating TNFα production.

4. Discussion

SP110 is strongly regulated by IFNs [17], suggesting its possible role in microbial immunity. Although many groups have studied the associations between the gene and TB susceptibility in a variety of populations, these studies show inconclusive results [2229]. In our study, we recruited healthy household contacts of TB patients as controls, as these contacts are at a high risk of exposure to Mtb. It has been reported that approximately 80-100% of the contacts may have Mtb infection, and on average, 20% of them may develop disease [37], indicating that household contacts are at high risk of LTBI and active TB disease. Therefore, examination of this group carries a considerable importance for prevention and control of TB disease. To the best of our knowledge, this is the first study to demonstrate an association between SP110 and LTBI, and it is the first study of the gene in the Taiwanese population. This work will help clarify the relationship between genetic variation in SP110 with latent and active TB infection in an Asian population.

Several SNPs (rs7580900, rs7580912, rs9061, rs11556887, and rs2241525) in SP110 showed an association with susceptibility to LTBI and/or TB disease in our study. In Table 3, we found that genotypes “GG” in rs7580912 (, OR: 2.451, 95% CI: 1.12–5.364) and “GG” in rs7580900 (, OR: 2.584, 95% CI: 1.208–5.53) were associated with LTBI susceptibility in LTBI cases vs. healthy controls, while the same SNP genotypes exhibited a protective effect on LTBI (“GG” in rs7580912, , OR: 0.392, 95% CI: 0.179–0.86; “GG” in rs7580900, , OR: 0.392, 95% CI: 0.182–0.848) in TB patients vs. LTBI individuals. In the Haploview analysis, we also found that the haplotype of multiple SNPs “ATGTACGCGA” in the SP110 gene was significantly associated with LTBI risk (, OR: 2.76, 95% CI: 1.12–6.84) in LTBI cases vs. healthy controls, while the same haplotypes had a protective effect on TB disease (, OR: 0.38, 95% CI: 0.16–0.95) in LTBI individuals vs. TB patients (see Figures S1 and S2 in the Supplementary Materials). In addition, 3 SNPs (rs9061, rs7580900, and rs7580912) with the haplotype “GAA” (, OR: 1.77, 95% CI: 1.20–2.60), “GGG” (, OR: 2.06, 95% CI: 1.34–3.17), “GGA” (, OR: 2.81, 95% CI: 1.49–5.27), and “GAG” (, OR: 10, 95% CI: 2.36–42.32) were associated with disease risk in LTBI cases vs. healthy controls, while the same haplotypes showed a protective effect on TB disease for TB patients vs. LTBI individuals (, OR: 0.61, 95% CI: 0.42–0.90 for “GAA”; , OR: 0.65, 95% CI: 0.43–0.98 for “GGG”; , OR: 0.3, 95% CI: 0.16–0.57 for “GGA”; and , OR: 0.1, 95% CI: 0.02–0.4 for “GAG”) (Table 5). These results revealed that the same SNP genotypes or haplotypes in the SP110 gene had opposite effects on the control of susceptibility to LTBI and TB disease, suggesting that the gene may have differential roles in the control of susceptibility to LTBI and TB disease.

SNP rs9061 (G→A) introduces an amino acid change from glutamic acid to lysine at codon position 207 of the SP110 protein. Transforming an acidic amino acid to a basic amino acid may alter the protein structure or posttranslational modification of the SP110 protein leading to better downregulation of TNFα production. It has been suggested that the A allele may cause changes in α-helices and β-sheets in the secondary structure of the SP110 protein compared with the G allele [38]. In addition, we analyzed the SP110 protein using the ELM database (http://elm.eu.org/) and found that this amino acid change generates a potential binding motif for the C-terminal ubiquitin-like domain (CTD) of ubiquitin specific protease 7 (USP7), one of the most abundant deubiquitinases [39]. USP7 plays important roles in various biological activities, including cell survival, proliferation, apoptosis, and tumorigenesis [40, 41]. As shown in Figure 2, the data demonstrated an association between the “GA” genotype of rs9061 with lower TNFα levels in plasma from LTBI individuals compared to “GG” LTBI subjects. One possible explanation for this result is that SP110 with the “GA” genotype at rs9061 may interact with USP7 and thus is more stable than the SP110 with the “GG” genotype, resulting in a more efficient downregulation of TNFα production. However, further functional studies are needed to elucidate the exact effects of SNP rs9061 on the SP110 gene and disease risk.

Individuals with LTBI, based on a positive result in the TST or the IGRA, usually show no disease symptoms and acquire an effective adaptive immunity. However, a proportion of people with LTBI might reactivate and develop clinical disease. The positive results from these assays indicate an immune response to stimulation by Mtb antigens in the LTBI population; however, neither TST nor IGRA can distinguish between LTBI and active TB. In addition, these assays also cannot predict which LTBI cases will progress to active TB [7]. To control and eliminate TB, worldwide TB eradication endeavors have been focused on identifying and treating cases with LTBI. Therefore, our studies identified several SNPs in SP110 that are significantly associated with controlling susceptibility to LTBI and TB disease and thus may provide novel predictive markers for latent and active TB infection. This study may also yield an improved strategy that can identify persons at increased risk of the disease.

There were some limitations in our study. First, this study dealt with relatively small sample sizes (301, 68, and 278 for TB patients, LTBI cases, and healthy household contacts, respectively). Further large-scale studies are, therefore, needed to validate the predictive value of the SP110 SNPs identified in this study. Second, the LTBI number of our study group is small, and besides, the results on plasma TNFα levels were analyzed from groups with even smaller sample sizes; therefore, the difference in TNFα levels in plasma of LTBI cases with two genotypes may be underestimated and should be verified in a different population. In addition, the analysis of TNFα levels in plasma of active TB patients with the different genotypes is in progress. The result may help support our finding in LTBI individuals and clarify the potential association between the SP110 SNPs and disease status.

5. Conclusions

The results suggest that the SP110 variants have a role in controlling genetic susceptibility to latent and active TB infection in humans. To the best of our knowledge, this is the first report demonstrating associations of polymorphisms in SP110 with LTBI susceptibility. Additionally, we provide evidence that the SP110 rs9061 SNP is associated with plasma TNFα levels in LTBI individuals. These data suggest that the identified SP110 SNPs may serve as a biomarker for latent and active TB infection in Taiwan.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

This research was funded by the Ministry of Science and Technology, Taiwan (NSC 97-2320-B-010-006-MY3, NSC 99-2321-B-002-022-, NSC 100-2321-B-002-009-, NSC 100-2320-B-002-106-, NSC 101-2321-B-002-001-, NSC 101-2320-B-002-022-, NSC 103-2320-B-002-044-MY3, and 106-2320-B-002-020- to BSY; NSC 99-2320-B-002-078-MY3, NSC 103-2314-B-002-083-, and 104-2314-B-002-075- to MLC), and Institute for Biotechnology and Medicine Industry, Taiwan. We thank Hsiang-Ting Hsu and Bo-Wen Chen for technical assistance. We also thank the expert assistance provided by Microarray Core Facility of National Research Program for Genomic Medicine of National Science Council in Taiwan.

Supplementary Materials

Table S1: primers for genotype analysis for SP110. Figure S1: association of SP110 SNP haplotypes with LTBI risk in LTBI cases vs. healthy controls. (a) Haplotype block map for SP110 with 10 SNPs. (b) Association of haplotype frequencies with LTBI risk in LTBI cases and healthy controls. Figure S2: association of SP110 SNP haplotypes with TB risk in TB cases vs. LTBI individuals. (a) Haplotype block map for SP110 with 10 SNPs. (b) Association of haplotype frequencies with TB risk in TB cases vs. LTBI individuals. Figure S3: association of SP110 SNP haplotypes with TB risk in TB cases vs. healthy controls. (a) Haplotype block map for SP110 with 10 SNPs. (b) Association of haplotype frequencies with TB risk in TB cases and healthy controls. (Supplementary Materials)