Transcriptome-Based Molecular Networks Uncovered Interplay Between Druggable Genes of CD8+ T Cells and Changes in Immune Cell Landscape in Patients With Pulmonary Tuberculosis

Background: Tuberculosis (TB) is a major infectious disease, where incomplete information about host genetics and immune responses is hindering the development of transformative therapies. This study characterized the immune cell landscape and blood transcriptomic profile of patients with pulmonary TB (PTB) to identify the potential therapeutic biomarkers.

Methods: The blood transcriptome profile of patients with PTB and controls were used for fractionating immune cell populations with the CIBERSORT algorithm and then to identify differentially expressed genes (DEGs) with R/Bioconductor packages. Later, systems biology investigations (such as semantic similarity, gene correlation, and graph theory parameters) were implemented to prioritize druggable genes contributing to the immune cell alterations in patients with TB. Finally, real time-PCR (RT-PCR) was used to confirm gene expression levels.

Results: Patients with PTB had higher levels of four immune subpopulations like CD8⁺ T cells (P = 1.9 × 10⁻⁸), natural killer (NK) cells resting (P = 6.3 × 10⁻⁵), monocytes (P = 6.4 × 10⁻⁶), and neutrophils (P = 1.6 × 10⁻⁷). The functional enrichment of 624 DEGs identified in the blood transcriptome of patients with PTB revealed major dysregulation of T cell-related ontologies and pathways (q ≤ 0.05). Of the 96 DEGs shared between transcriptome and immune cell types, 39 overlapped with TB meta-profiling genetic signatures, and their semantic similarity analysis with the remaining 57 genes, yielded 45 new candidate TB markers. This study identified 9 CD8⁺ T cell-associated genes (ITK, CD2, CD6, CD247, ZAP70, CD3D, SH2D1A, CD3E, and IL7R) as potential therapeutic targets of PTB by combining computational druggability and co-expression (r² ≥ |0.7|) approaches.

Conclusion: The changes in immune cell proportion and the downregulation of T cell-related genes may provide new insights in developing therapeutic compounds against chronic TB.

Introduction

Tuberculosis (TB) is a chronic infectious lung disease caused by pathogenic Mycobacterium tuberculosis (MTB) belonging to the Mycobacteriaceae family. Despite the widespread use of antibiotics and live attenuated vaccine, TB remains to be the major cause of morbidity and death among all bacterial diseases (1). This is primarily due to the rapid emergence of drug-resistant MTB strains and the incomplete knowledge of complex host-pathogen interactions (2). In the initial stages of infection, MTB invades and replicates in the macrophages after reaching the alveolar air sacs of the lungs (3). Granulomas, hallmark of TB, are formed around the infected macrophages by the organized aggregation of immune cells (like T and B lymphocytes), multinucleated giant cells, dendritic cells, and fibroblasts. Granulomas also suppress the host immune responses, as dendritic cells and macrophages were unable to present antigen to lymphocytes (4). It is noteworthy to mention that mycobacteria can induce distinct host responses from asymptomatic conditions to severe pulmonary illness (5). However, underlying immune cell types and their association with the differentially expressed genes in TB and how they contribute to severe infection are not yet fully explored.

Over the past few decades, microarray-based genome-wide RNA profiling has evolved as a powerful approach to investigate the host transcriptional response (of ~19,000 genes) in infectious diseases (6). However, differences in the type of clinical samples, array platforms, and statistical approaches used, created a discordance in interpreting massive transcriptomics data. Advances in statistical modeling and bioinformatics approaches have accelerated the identification of disease-centric genes by employing gene networking methods based on graph topological parameters for many infectious diseases (7, 8). Moreover, the new bioinformatic methods like estimating relative subsets of RNA transcripts (CIBERSORT), Tumor Immune Estimation Resource (TIMER), and Estimating the Proportions of Immune and Cancer cells (EPIC) are developed to characterize immune cell composition using large-scale gene expression data (9, 10). These bioinformatic methods implement functional enrichment scores based on the presence of the query genes over reference gene sets. They perform variety of biological analyses including immune responses based on the defined gene sets. Exploring abnormal immune cell infiltration is critical for developing novel transformative therapies to combat diseases such as cancer, myocarditis, and TB (11, 12). Therefore, in order to characterize alterations in immune cell proportion landscape and transcriptomic profile, and to identify new molecular therapeutic targets, this study applied statistical and knowledge-based systemic investigations (such as semantic similarity, gene correlation, and graph theory parameters) to the blood transcription data of patients with TB.

Materials and Methods

Study Design and Global Expression Data

The genome-wide gene expression dataset (GSE83456) (13) was imported in raw format from the Gene Expression Omnibus (GEO) database (www.ncbi.nlm.nih.gov/geo). This dataset has expression profiles of 45 pulmonary TB (PTB) and 61 control blood samples generated on the Illumina Human HT-12 V4.0 expression bead chip (Illumina, Inc, USA). The detailed sample information is given in the Supplementary Material Table 1. Figure 1 depicts the overall work design employed in the current research analysis.

FIGURE 1

Figure 1. Study workflow. The gene expression profiles of patients with tuberculosis (TB) and controls were used to deconvolute immune cell fractions. Differentially expressed genes (DEGs) were mapped to functional pathways and then correlated with immune cell and TB meta-analysis gene signatures. The overlapping genes showing semantic similarity were explored by druggability and protein interaction analysis to identify novel candidate therapeutic targets/biomarkers for combating TB infection.

Global Data Preprocessing and Screening of Differentially Expressed Genes

R/Bioconductor packages were used to analyze microarray gene expression data. Raw data was fed into the R package limma (14) for the standardization and noise reduction of the probe data, and the raw signal levels for each probe set were standardized. The Quantile method was used to normalize the microarray datasets. The t-statistic was used to detect statistically significant differentially expressed genes (DEGs) between the PTB and control samples. To eliminate false positives, the Benjamini and Hochberg false discovery rate (FDR) with p < 0.05 was used as a cut-off point for gene data. Thereafter, probes were matched to Entrez Gene IDs, and duplicates (with the highest fold change difference) and unmatched transcripts were filtered out. In the final stage, all the DEGs were classified as up- and downregulated genes based on the fold change threshold (FC ≥ |1.5|).

Identification of Immune Cell Composition From Gene Expression Profiles

The fractions of 22 immune cell types in the PTB transcriptome profile were estimated using the CIBERSORT algorithm (15). This program employs linear support vector (SVR) regression to perform feature selection and to deconvolve the cell mixture from the gene expression profile. In this study, gene expression profiles of PTB and control samples were fed into the CIBERSORT algorithm where the algorithm converts the gene expression matrix into the immune cell-matrix and applies the filtering criteria of 1,000 permutations and significant p value set at ≤ 0.05.

MTB Meta-Profiling Genetic Signatures

We obtained 380 genetic signatures identified from the modular analysis and meta-profiling of 16 publicly available gene expression datasets (16) (Supplementary Material Table 2). We compared the overlapping genes between these 380 TB genetic signatures, DEGs, and genes associated with immune cell populations for downstream analysis.

Identification of Immune Pathways From Gene Expression Profiles

The functional enrichment analysis of the DEGs was performed using g:Profiler (17), a webserver to interpret the function of gene lists (https://biit.cs.ut.ee/gprofiler/gost). This server matches a queried gene list to established functional data sources and uncovers gene ontologies as well as pathway terms that are significantly enriched at q ≤ 0.05. Immune-related pathways were screened from the functional enrichment list. The DEGs which are contributing to immune-associated pathways were mapped to the known signature of TB and immune signature genes in the CIBERSORT to identify unreported genes in TB.

Identification of Semantic Similarity

Using encoded evidence in the Gene Ontology (GO) hierarchy, the functional similarity between unreported genes and known TB signatures is assessed. In this study, we used Wang's similarity metric to compare the biological process (BP) hierarchy. To quantify the semantic similarity between gene pairs, we used the R tool GoSemSim (18).

We employed Resnik's measure of Best-Match Average (BMA) method, which combines the semantic relationship scores of numerous GO terms and produces the average of all maximal similarities in each row and column because a gene can be annotated by many GO terms (19). Following that, gene pairs were selected based on a semantic score of ≤ 0.5, with a larger score indicating a stronger relationship. The following formula was used to calculate the semantic similarity among gene pairs:

$\begin{array}{ll}\begin{array}{c}{S}_{GO}(A,B)=\frac{{\displaystyle {\sum }_{t\in T_A{{\displaystyle \cap }}^{}{T}_B}(S_A(t)+S_B(t))}}{{\displaystyle {\sum }_{t\in T_A}{S}_A(t)+{\displaystyle {\sum }_{t\in T_B}{S}_B(t)}}}\end{array}& (1)\end{array}$

where T_A designates the contribution of t ∈ T_A term to the semantics of A based on the relative positions of t and A in the graph, and S_A(t) implies the role of t ∈ T_A term to the semantics of A.

Druggability Analysis

The Drug–Gene Interaction Database (DGIdb) (20) was used to assess the druggability of the genes. DGIdb is a central resource for drug-gene interaction data and the potential druggability of each query gene based on different databases. We included approved drugs, antineoplastic drugs, and Immunotherapeutic drug interactions filters and in advance filters, we selected 9 Disease-Agnositic sources databases, 43 gene categories, and 31 interaction types. We used drug target interaction with interactions score ≥0.03 to search the DEGs, which could act as potential drug target genes for MTB.

Correlation Among the Druggable Genes

The correlation between the druggable genes in PTB was investigated using Pearson's correlation method. The correlation (r) between each pair of gene matrices was ranked using Pearson's correlation coefficient (PCC). The formula used for computing the PCC existing between two genes is given below.

$\begin{array}{ll}PCC(r)=\frac{{\displaystyle {\sum }_{i=1}^n(x_i-\overline{x})(y_i-\overline{y})}}{\sqrt{{\displaystyle {\sum }_{i=1}^n{(x_i-\overline{x})}^2}}\sqrt{{\displaystyle {\sum }_{i=1}^n{(y_i-\overline{y})}^2}}}& (2)\end{array}$

where $\bar{x}$ and ȳ are the average of sample's gene expression signal in PTB of the two genes, respectively. The gene co-expression was confirmed using the Search Tool for the Retrieval of Interacting Genes (STRING) (21), an online protein interaction database, with high confidence interaction score of ≥0.7.

Real Time-PCR (RT-PCR) Validation of Druggable Genes

In order to verify our bioinformatics findings, we validated the expression of 9 druggable genes by the RT-PCR method. In brief, the RNA collected from THP-1 cell lines infected with the MTB strain (H37Rv) was collected after the post-incubation period as previously described (22). In brief, total RNA was reverse transcribed as complementary DNA (cDNA) and then amplified by the RT-PCR method using gene-specific oligonucleotide primers. The relative expression level of potentially druggable genes between the control and test cell lines was estimated by the 2^−ΔΔCT formula after normalizing their expression levels with the GAPDH internal reference gene. A p < 0.05 under the standard two-tailed t-test was considered a significant value.

Results

Immune Cell Proportion Analysis of PTB Gene Expression Profile

The immune cell proportion landscape of PTB is not yet fully revealed, particularly in low abundant cell subpopulations. In this study, the CIBERSORT algorithm has identified the enrichment of genes associated with 10 types of adaptive immune cells like B cells naive, plasma cells, T follicular helper cells, CD8⁺ T cells, resting memory CD4⁺ T cells, T cells, CD4⁺ memory T cells activated, memory B cells, naive CD4⁺ T cells, regulatory T cells (Tregs), and Gamma-delta (γδ) T cells. On the other hand, DEGs associated with 12 innate immune cell type categories were NK cells resting, macrophages M2, monocytes, macrophages M1, macrophages M0, resting dendritic cells, eosinophils, dendritic cells activated, mast cells resting, NK cells activated, mast cells activated, and neutrophils were also found be enriched. The immune cell proportions of adaptive immune cells and innate immune cells are represented in Supplementary Material Figure 1.

The genes associated with four immune cells; NK cells activated, T follicular helper cells, dendritic cells resting, and eosinophils, were not significantly enriched in both groups. The proportion plot of the enriched immune cell types is represented in Figure 2A. We observed higher relative proportion of genes enriched for cell types like CD8⁺ T cells (P = 1.9 × 10⁻⁸), NK cells resting (P = 6.3 × 10⁻⁵), monocytes (P = 6.4 × 10⁻⁶), and neutrophils (P = 1.6 × 10⁻⁷) in the PTB samples compared to the control samples (Figure 2B, Supplementary Material Figure 1). Among the 4 cell types with a higher relative proportion of enriched genes from DEGs, the genes of CD8⁺ T cells and NK resting cells were found to be downregulated in PTB when compared to the control samples. On other hand, monocytes and neutrophil-associated genes were highly active in PTB when compared to control samples.

FIGURE 2

Figure 2. The immune cell proportion landscape between pulmonary TB (PTB) and controls. (A) The relative proportion of immune cell subpopulations in GSE83456 dataset. (B) The difference of immune infiltration between PTB and normal controls (the control group was marked in blue color and the PTB group was marked in red color. P < 0.05 were considered as statistically significant).

Identification of DEGs From Gene Expression Profile

The standardized gene expression data of “PTB vs. controls” was used to identify the differentially expressed genes. The volcano plot representing the distribution of fold change and the significant p-value is given in Figure 3A. The PTB vs. control group analysis revealed 624 DEGs (FC |1.5|, adj p-value of 0.05), with 393 upregulated and 231 downregulated genes. The top 10 DEGs obtained from PTB vs controls are given in Table 1. The mean distribution of intensity of differentially expressed genes in PTB and control samples is represented in Figure 3B.

FIGURE 3

Figure 3. Graphical distribution of differentially expressed genes. (A) Volcano plot representing the distribution of fold change and p-value significance. (B) The distribution mean intensity of differentially expressed genes in the PTB and control samples. (C) Red and green nodes represent up and downregulation of genes and black nodes are the immune cell types. (D) The Venn diagram represents the overlap of DEGs with immune and TB signatures. (E) Semantic similarity of pairs of genes between differentially expressed immune signatures and TB signatures. The selected gene pairs with higher functional similarity (≥0.5) are highlighted in green color.

TABLE 1

Table 1. The top 10 differentially expressed gene list in pulmonary tuberculosis (PTB).

Functional Enrichment Analysis of DEGs

The differentially expressed genes enriched using g:Profiler with the statistical significance of q value ≤ 0.05, generated 309 ontologies of Biological Process (BP), 17 ontologies of Molecular Function (MF), 42 ontologies of Cellular Component (CC), and 85 terms in pathways (Supplementary Material Table 3). Overall, the enrichment analysis has shown the overlap with immune-related ontologies and pathways. We pooled immune-related pathways from enrichment terms to check how DEGs affect the immune system pathways. We observed the upregulation of pathways such as interferon signaling (q = 9.16 × 10⁻²⁷), cytokine signaling in the immune system (q = 2.34 × 10⁻²¹), neutrophil degranulation (q = 3.13 × 10⁻¹¹), viral genome replication (q = 2.47 × 10⁻⁷), and response to biotic stimulus, etc. (Supplementary Material Figure S1). On the other hand, pathways like T-cell antigen receptor signaling, antigen receptor-mediated signaling, NF-kappa B signaling, T cell activation, T cell receptor signaling, leukocyte differentiation, leukocyte activation, alpha-beta T cell activation, and T cell differentiation, were downregulated (Supplementary Material Figure S1). Overall, our functional enrichment analysis points to a major downregulation of T cell-related ontologies and pathways.

Mapping DEGs to Immune Cell Proportions in PTB

Here we investigated the genes overlapping between the CIBERSORT signature and DEGs. There are about 96 DEGs (Figure 3C) contributing to different immune cell types (Supplementary Material Table 4). Interestingly, we found that 31.25% of DEGs were contributing to the immune cell type “CD8⁺ T cells.” We also observed that all those genes contributing to the “CD8⁺ T cells” were downregulated in PTB samples, as shown in Figure 3C. The findings from the mapping of DEGs to immune cell proportions are consistent with functional enrichment analysis, where T cell-related pathways have shown major dysregulation.

Comparison of DEGs, Immune Cell Signatures With TB Meta-Analysis Signatures

The differentially expressed immune signatures in the sample of patients with TB were compared with the known signatures of TB (Supplementary Material Table 4). Here, we observed 39 (40.6%) differentially expressed immune signatures overlapping with TB meta-analysis gene signatures and 57 novel genes (59.3%) contributing to the immune cell proportion (Figure 3D). Further, semantic similarity (functional association) of these 57 novel genes with 39 overlapping with TB signatures was performed to identify the most predominant genes. The semantic similarity score of ≥0.5 among the gene pairs was considered as a highly significant score implying a stronger association. The semantic similarity of 45/57 genes has shown a stronger functional association with overlapping with TB signatures (Figure 3E). Again, it is important to pinpoint that 20 out of 45 genes (44%) were contributing to the immune cell type “CD8+ T cells.”

Druggability and Co-expression Analysis

Druggability analysis was performed on the 45 genes that had shown higher functional similarity to the known TB signature. We found that 21 druggable genes (Supplementary Material Figure S2) (46%) with an interaction score ≥ 0.03, were enriched against terms like an antibody, binder, inhibitor, antagonist, agonist, modulator, and activator. Of all the druggable genes, ITK and FCGR3B genes were observed to have the highest number of drug interactions (19 drugs) followed by PTGDR (13 drugs), TLR7 (12 drugs), and CD3E (10 drugs). The drug-target interaction network is represented in Figure 4A. Next, we checked the association of druggable targets to the immune cell types. Among the 21 druggable targets, 48% were contributing to the immune cell type “CD8⁺ T cells” followed by monocytes (9%), naïve B cells (9%), and neutrophils (9%) (Figure 4B). Interestingly, the expression patterns of these 21 druggable genes have shown a clear distinction in PTB when compared to control samples (Figure 4C).

FIGURE 4

Figure 4. The druggable targets with their expression, interaction and co-expression (A) The network depicts the drug-target interaction where black and pink nodes represent target and drugs respectively. (B) The co-expression (a similar pattern of gene expression) among the druggable genes where red and blue represent the positive and negative correlation. (C) The pattern of gene expression of druggable genes in control and PTB patient samples clearly depicts a distinction between the control and PTB groups.

To check the correlation of these 21 druggable targets in patients with PTB, we performed Pearson's correlation analysis (Figure 4B). The correlation analysis performed between druggable targets in PTB samples resulted in 15 genes (ITK, CD2, CD6, CD247, CD27, CD3D, ZAP70, SH2D1A, IL2RB, CEACAM3, IL7R, TLR2, CD3E, PTGDR, FCGR3B) with higher Pearson correlation coefficient (r² ≥ |0.7|). The co-expressed genes are shown in Table 2. Among the 15 co-expressed genes, 12 and 3 were down and upregulated respectively (Figure 5A). Upregulated genes were contributing to the immune cell types “Monocytes” (TLR2) and Neutrophils (FCGR3B and CEACAM3). Again, 10 downregulated genes were contributing to “T cells CD8” and 1 each for “T cells regulatory” and “B cells memory” (Figure 5B).

TABLE 2

Table 2. List of co-expressed genes and the druggable targets.

FIGURE 5

Figure 5. The interaction networks of druggable targets. (A) co-expressed network of druggable targets where red and green colored nodes represent up and downregulated genes. The red edge represents the association of targets to immune cells and the blue edges depict the interaction among the target genes. (B) The pie chart represents the composition of immune cell types with the co-expressed druggable targets. (C) The interaction of the co-expressed druggable targets in the STRING database.

Interestingly, we noticed a cluster of 12 co-expressed genes ITK, CD2, CD6, ZAP70, CD247, CD3D, SH2D1A, CD27, CD3E, IL2RB, IL7R, and PTGDR. To validate the co-expression among the 12 downregulated genes we queried them in the STRING database with a high confidence score of ≥ 0.7. The STRING database identified strong interaction among the 11 genes except for PTGDR (Figure 5C). Co-expression and protein interaction network from the STRING database has shown the mutual influence of the 11 genes in the expression and functional activities. Nine genes (expect CD27, PTGDR, and IL2RB) were contributing to “CD8⁺ T cells.” All the genes were predicted to be targeted by different drug molecules, most of which are monoclonal antibodies (Table 3). Hence, the integrated analysis has depicted predominant deregulation of “CD8⁺ T cells” as the key genetic signatures for active PTB.

TABLE 3

Table 3. The list of drugs shows direct interaction with 9 genes associated with CD+T cell functioning.

RT-PCR Validation of Druggable Genes

The real-time PCR gene expression results showed that the relative expression levels of 9 potential druggable genes (ITK, CD2, CD6, CD247, ZAP70, CD3D, SH2D1A, CD3E, and IL7R) were consistent with the findings of microarray hybridization. All the genes were differentially expressed between treated and untreated cell lines (p ≤ 0.01). These results confirm the dysregulated “CD8⁺ T cell signaling” plays important role in establishing TB infection (Supplementary Material Figure S3).

Discussion

Host genetic factors are known to play an important role in regulating the initial TB infection and determining the disease progression in the lungs (37). Genome-wide association studies have underlined the relevance of numerous polymorphisms in immune response-related genes in contributing to susceptibility or resistance to TB (38). However, polymorphism studies were unable to provide full insight into the complex molecular crosstalk between thousands of host genes involved in innate and adaptive immune responses. In this context, high throughput transcriptome approaches have shown great promise in dissecting the host-pathogen interactions thereby helping to develop a novel vaccine and therapeutic targets for several infectious diseases (39–41). Therefore, we explored host immune system response through integrated systems biology approach based on immune cell subtyping and differential gene expression profiles of patients with PTB to normal controls.

The cellular and molecular background of TB-induced systemic immunological dysregulation is poorly understood. Therefore, we screened the DEGs in PTB and deciphered their contribution to immune cell proportion alterations. Traditionally, host transcriptomics studies have relied on whole blood to characterize TB gene signatures by aggregating transcriptomic signals from many different cell types but were unable to identify specific immune cell type signatures (42). To overcome these constraints, we used a powerful computational technique called CIBERSORT to define the range of immune cell states in the blood of patients with TB. This method relies on linear support vector regression (SVR), a machine learning approach to deconvolute the gene expression signatures, known as “signature matrix” for determining the relative fraction of immune cell proportions in blood or tissues (15). The CIBERSORT method has been widely used to infer immune cell types from transcriptomics data to predict outcomes of different cancers (9, 43, 44) and infectious diseases (45–47). In this study, the CIBERSORT output identified the downregulation of “CD8⁺ T cells” in patients with PTB.

The GO functional enrichment analysis of gene expression profile revealed the upregulation of interferon signaling, cytokine signaling in the immune system, neutrophil degranulation, and response to biotic stimulus pathways. The MTB infection of primary human macrophages is shown to induce type I IFN signaling and limit the expression of IL-1β, which imparts immunity against the infection (48). Additionally, the downregulation of major pathways associated with T cells function like T-cell antigen receptor signaling, leukocyte differentiation, leukocyte activation, T cell activation, T cell differentiation, T cell receptor signaling, alpha-beta T cell activation, NF-kappa B signaling, and antigen receptor-mediated signaling pathway were noted in PTB samples.

Majority of the DEGs contributing to the immune cell type “CD8⁺ T cells” were clearly downregulated in the PTB samples indicating their potential roles in defense against TB. These cells are also known as killer or cytotoxic T lymphocytes, as they potentially destroy the infected cells by recruiting cytokines and other immune cells to the site of infection. The low abundance of blood CD8⁺ T lymphocytes may impair the effective immunity against pathogens, as they lack a sufficient cytotoxic T cells to recognize the MHC class I-restricted epitopes of MTB antigens, in the site of infection (49). A recent RNA transcriptome study used the positron emission tomography (PET) data collected from recovered patients with TB, at 4th and 24th weeks has also reported that genes associated with the overexpression of B cells and down expression of T cells and platelets confirms our findings (50). Thus, the downregulation that contribute to immune cell type is concordant with the pathway enrichment analysis findings of lower expression of T cell-related ontologies and pathways.

The druggability potential of any protein is attributed to its binding specificity with small compounds following Lipinski's rule-of-five for drug likeliness (51). Numerous bioinformatic and empirical methods which consume less time and provide faster prescreening of druggability of candidate proteins than conventional methods have been developed (52, 53). There are a variety of computational methods available, which can predict druggability and protein-binding sites by using energy dynamics to geometrical topological estimations, and from flexible to rigid proteins (54, 55).

By applying druggability and co-expression features we identified 9 CD8⁺ T cells associated genes (ITK, CD2, CD6, CD247, ZAP70, CD3D, SH2D1A, CD3E, and IL7R) as potential therapeutic targets of PTB. However, it is pivotal to carefully prioritize the drug molecules based on their mode of action, whether activator or inhibitor based on the gene expression status. For example, over-expressed genes can be targeted by inhibitory molecules, and downregulated genes can be targeted by activator molecules (56). From the above 9 genes, the therapeutic potential of ITK and IL7R has been characterized by experimental methods. ITK is a tyrosine kinase expressed on T-cells, which regulates its T-cell development and function. Human lungs with ITK deficiency impair early protection against MTB in vivo (57). Improving ITK signaling pathways could become an alternative approach for combating MTB infection. One study reveals the role of IL7R on T-cell immunity in human TB (58). The authors reported that patients with TB had lower IL7R concentrations and lower IL7R expression in T cells than healthy controls, indicating that patients with TB have impaired T-cell sensitivity. In addition, due to post-transcriptional processes, patients with TB had reduced amounts of IL7R in T cells. In vitro experiments revealed that MTB-specific T lymphocytes from patients with TB have reduced IL-7-induced STAT5 phosphorylation and IL-7-promoted cytokine production (59). The role of the remaining 7 genes (CD2, CD6, CD247, ZAP70, CD3D, SH2D1A, and CD3E) in T-cell signaling and modulation of host immune responses in mycobacterium infections is also supported (60–64).

Our results highlight the dysregulation of CD8⁺ T cells and the associated genes in PTB patients. These findings are exciting not just from the fact that CD8⁺ T cell-associated genes have the potential to act as potential therapeutic targets but prove that their role is not less important than CD4⁺ T cells in controlling MTB infection. We acknowledge that our strategy has some technical constraints. CIBERSORT was a convenient computational tool for determining infiltrating immune cell fractions, but it was still less precise than immunohistochemistry or flow cytometry, which could lead to inaccuracies in immune cell fractions. However, to overcome this limitation to some extent, we have linked gene expression profiles of immune signatures followed by functional enrichment, semantic similarity, druggability, and co-expression among the identified key signatures.

Conclusion

In this study, by coupling computational deconvolution algorithms and high throughput blood transcriptomics data, we identified the difference in T-cell-related immune cell populations among patients with PTB. The functional enrichment of 624 DEGs (393 over-expressed and 231 under-expressed) identified in the blood transcriptome of PTB patients revealed the major dysregulation of T cell-related ontologies and pathways. By linking DEGs against immune cell populations and TB gene signature, this study identified 9 CD8⁺ T cells associated genes (ITK, CD2, CD6, CD247, ZAP70, CD3D, SH2D1A, CD3E, and IL7R) as potential therapeutic targets of PTB. The expression levels of these 9 genes in MTB infection in cell lines were assessed by RT-PCR-based expression assay, confirming the experimental validation. However, further in vitro and in vivo studies are needed to establish the role of these genes in PTB infection, progression, and treatment.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: http://www.ncbi.nlm.nih.gov/geo; GSE83456.

Author Contributions

FA, BB, VV, and NS: conceptualization. FA, HA, SP, and BB: data curation. FA and BB: formal analysis. FA: funding acquisition and project administration. FA, NZ, SP, and BB: methodology. BB and SP: software and visualization. FA, NS, and BB: supervision. BB: validation. FA, NZ, HA, SP, PS, VV, RE, NS, and BB: writing original draft and editing. All authors contributed to the article and approved the submitted version.

Funding

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number IFPRC-063-247-2020 and King Abdulaziz University, DSR, Jeddah, Saudi Arabia.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.812857/full#supplementary-material

References

1. Makam P, Matsa R. “Big Three” infectious diseases: tuberculosis, malaria and HIV/AIDS. Curr Top Med Chem. (2021) 21:2779–99. doi: 10.2174/1568026621666210916170417

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Chakaya J, Khan M, Ntoumi F, Aklillu E, Fatima R, Mwaba P, et al. Global Tuberculosis Report 2020—reflections on the Global TB burden, treatment and prevention efforts. Int J Infect Dis. (2021). doi: 10.1016/j.ijid.2021.02.107

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Maertzdorf J, Tonnies M, Lozza L, Schommer-Leitner S, Mollenkopf H, Bauer TT, et al. Mycobacterium tuberculosis Invasion of the Human Lung: First Contact. Front Immunol. (2018) 9:1346. doi: 10.3389/fimmu.2018.01346

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Reece ST, Kaufmann SH. Floating between the poles of pathology and protection: can we pin down the granuloma in tuberculosis? Curr Opin Microbiol. (2012) 15:63–70. doi: 10.1016/j.mib.2011.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Morris SK, Giroux RJP, Consunji-Araneta R, Stewart K, Baikie M, Kakkar F, et al. Epidemiology, clinical features and outcomes of incident tuberculosis in children in Canada in 2013-2016: results of a national surveillance study. Arch Dis Child. (2021). doi: 10.1136/archdischild-2021-322092

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Von Both U, Kaforou M, Levin M, Newton SM. Understanding immune protection against tuberculosis using RNA expression profiling. Vaccine. (2015) 33:5289–93. doi: 10.1016/j.vaccine.2015.05.025

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Banaganapalli B, Al-Rayes N, Awan ZA, Alsulaimany FA, Alamri AS, Elango R, et al. Multilevel systems biology analysis of lung transcriptomics data identifies key miRNAs and potential miRNA target genes for SARS-CoV-2 infection. Comput Biol Med. (2021) 135:104570. doi: 10.1016/j.compbiomed.2021.104570

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Banaganapalli B, Mansour H, Mohammed A, Alharthi AM, Aljuaid NM, Nasser KK, et al. Exploring celiac disease candidate pathways by global gene expression profiling and gene network cluster analysis. Sci Rep. (2020) 10:16290. doi: 10.1038/s41598-020-73288-6

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Craven KE, Gokmen-Polar Y, Badve SS. CIBERSORT analysis of TCGA and METABRIC identifies subgroups with better outcomes in triple negative breast cancer. Sci Rep. (2021) 11:4691. doi: 10.1038/s41598-021-83913-7

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Xie Z, Bailey A, Kuleshov MV, Clarke DJB, Evangelista JE, Jenkins SL, et al. Gene set knowledge discovery with enrichr. Curr Protocol. (2021) 1:e90. doi: 10.1002/cpz1.90

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kawada JI, Takeuchi S, Imai H, Okumura T, Horiba K, Suzuki T, et al. Immune cell infiltration landscapes in pediatric acute myocarditis analyzed by CIBERSORT. J Cardiol. (2021) 77:174–8. doi: 10.1016/j.jjcc.2020.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Feng Z, Qu J, Liu X, Liang J, Li Y, Jiang J, et al. Integrated bioinformatics analysis of differentially expressed genes and immune cell infiltration characteristics in Esophageal Squamous cell carcinoma. Sci Rep. (2021) 11:16696. doi: 10.1038/s41598-021-96274-y

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Blankley S, Graham CM, Turner J, Berry MP, Bloom CI, Xu Z, et al. The transcriptional signature of active tuberculosis reflects symptom status in extra-pulmonary and pulmonary tuberculosis. PLoS ONE. (2016) 11:e0162220. doi: 10.1371/journal.pone.0162220

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Blankley S, Graham CM, Levin J, Turner J, Berry MP, Bloom CI, et al. A 380-gene meta-signature of active tuberculosis compared with healthy controls. Eur Respir J. (2016) 47:1873–6. doi: 10.1183/13993003.02121-2015

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H, et al. g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. (2019) 47:W191–8. doi: 10.1093/nar/gkz369

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Yu G, Li F, Qin Y, Bo X, Wu Y, Wang S. GOSemSim: an R package for measuring semantic similarity among GO terms and gene products. Bioinformatics. (2010) 26:976–8. doi: 10.1093/bioinformatics/btq064

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Pesquita C, Faria D, Bastos H, Ferreira AE, Falcao AO, Couto FM. Metrics for GO based protein semantic similarity: a systematic evaluation. BMC Bioinformatics. (2008) 9 Suppl 5:S4. doi: 10.1186/1471-2105-9-S5-S4

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Cotto KC, Wagner AH, Feng YY, Kiwala S, Coffman AC, Spies G, et al. (2018). DGIdb 3.0: a redesign and expansion of the drug-gene interaction database. Nucleic Acids Res 46:D1068–73. doi: 10.1093/nar/gkx1143

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. (2021) 49:D605–12. doi: 10.1093/nar/gkab835

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Zhang W, Niu C, Fu RY, Peng ZY. Mycobacterium tuberculosis H37Rv infection regulates alternative splicing in Macrophages. Bioengineered. (2018) 9:203–8. doi: 10.1080/21655979.2017.1387692

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zapf CW, Gerstenberger BS, Xing L, Limburg DC, Anderson DR, Caspers N, et al. Covalent inhibitors of interleukin-2 inducible T cell kinase (itk) with nanomolar potency in a whole-blood assay. J Med Chem. (2012) 55:10047–63. doi: 10.1021/jm301190s

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ortonne JP. Clinical response to alefacept: results of a phase 3 study of intramuscular administration of alefacept in patients with chronic plaque psoriasis. J Eur Acad Dermatol Venereol. (2003) 17 Suppl 2:12–6. doi: 10.1046/j.1468-3083.17.s2.3.x

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Goedkoop AY, De Rie MA, Picavet DI, Kraan MC, Dinant HJ, Van Kuijk AW, et al. Alefacept therapy reduces the effector T-cell population in lesional psoriatic epidermis. Arch Dermatol Res. (2004) 295:465–73. doi: 10.1007/s00403-004-0450-y

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Chamian F, Lin SL, Lee E, Kikuchi T, Gilleaudeau P, Sullivan-Whalen M, et al. Alefacept (anti-CD2) causes a selective reduction in circulating effector memory T cells (Tem) and relative preservation of central memory T cells (Tcm) in psoriasis. J Transl Med. (2007) 5:27. doi: 10.1186/1479-5876-5-27

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Chen X, Ji ZL, Chen YZ. TTD: therapeutic target database. Nucleic Acids Res. (2002) 30:412–5. doi: 10.1093/nar/30.1.412

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Da Silva AJ, Brickelmaier M, Majeau GR, Li Z, Su L, Hsu YM, et al. Alefacept, an immunomodulatory recombinant LFA-3/IgG1 fusion protein, induces CD16 signaling and CD2/CD16-dependent apoptosis of CD2(+) cells. J Immunol. (2002) 168:4462–71. doi: 10.4049/jimmunol.168.9.4462

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Larsen R, Ryder LP, Svejgaard A, Gniadecki R. Changes in circulating lymphocyte subpopulations following administration of the leucocyte function-associated antigen-3 (LFA-3)/IgG1 fusion protein alefacept. Clin Exp Immunol. (2007) 149:23–30. doi: 10.1111/j.1365-2249.2007.03380.x

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Todd PA, Brogden RN. Muromonab CD3. A review of its pharmacology and therapeutic potential. Drugs. (1989) 37:871–99. doi: 10.2165/00003495-198937060-00004

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Wilde MI, Goa KL. Muromonab CD3: a reappraisal of its pharmacology and use as prophylaxis of solid organ transplant rejection. Drugs. (1996) 51:865–94. doi: 10.2165/00003495-199651050-00010

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Chen IJ, Chen HL, Demetriou M. Lateral compartmentalization of T cell receptor versus CD45 by galectin-N-glycan binding and microfilaments coordinate basal and activation signaling. J Biol Chem. (2007) 282:35361–72. doi: 10.1074/jbc.M706923200

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Imming P, Sinning C, Meyer A. Drugs, their targets and the nature and number of drug targets. Nat Rev Drug Discov. (2006) 5:821–34. doi: 10.1038/nrd2132

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Overington JP, Al-Lazikani B, Hopkins AL. How many drug targets are there? Nat Rev Drug Discov. (2006) 5:993–6. doi: 10.1038/nrd2199

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Roberts KG, Morin RD, Zhang J, Hirst M, Zhao Y, Su X, et al. Genetic alterations activating kinase and cytokine receptor signaling in high-risk acute lymphoblastic leukemia. Cancer Cell. (2012) 22:153–66. doi: 10.1016/j.ccr.2012.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Maude SL, Tasian SK, Vincent T, Hall JW, Sheen C, Roberts KG, et al. Targeting JAK1/2 and mTOR in murine xenograft models of Ph-like acute lymphoblastic leukemia. Blood. (2012) 120:3510–8. doi: 10.1182/blood-2012-03-415448

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Chai Q, Lu Z, Liu CH. Host defense mechanisms against Mycobacterium tuberculosis. Cell Mol Life Sci. (2020) 77:1859–78. doi: 10.1007/s00018-019-03353-5

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Van Tong H, Velavan TP, Thye T, Meyer CG. Human genetic factors in tuberculosis: an update. Trop Med Int Health. (2017) 22:1063–71. doi: 10.1111/tmi.12923

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Domínguez Á, Muñoz E, López MC, Cordero M, Martínez JP, Viñas M. Transcriptomics as a tool to discover new antibacterial targets. Biotechnol Lett. (2017) 39:819–28. doi: 10.1007/s10529-017-2319-0

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Guleria V, Jaiswal V. Comparative transcriptome analysis of different stages of Plasmodium falciparum to explore vaccine and drug candidates. Genomics. (2020) 112:796–804. doi: 10.1016/j.ygeno.2019.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Le BL, Andreoletti G, Oskotsky T, Vallejo-Gracia A, Rosales R, Yu K, et al. Transcriptomics-based drug repositioning pipeline identifies therapeutic candidates for COVID-19. Sci Rep. (2021) 11:12310. doi: 10.1038/s41598-021-91625-1

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Blankley S, Berry MP, Graham CM, Bloom CI, Lipman M, O'garra A. The application of transcriptional blood signatures to enhance our understanding of the host response to infection: the example of tuberculosis. Philos Trans R Soc Lond B Biol Sci. (2014) 369:20130427. doi: 10.1098/rstb.2013.0427

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Huang R, Mao M, Lu Y, Yu Q, Liao L. A novel immune-related genes prognosis biomarker for melanoma: associated with tumor microenvironment. Aging. (2020) 12:6966–80. doi: 10.18632/aging.103054

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Wu Z, Chen H, Luo W, Zhang H, Li G, Zeng F, et al. The landscape of immune cells infiltrating in prostate cancer. Front Oncol. (2020) 10:517637. doi: 10.3389/fonc.2020.517637

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Reyes M, Filbin MR, Bhattacharyya RP, Billman K, Eisenhaure T, Hung DT, et al. An immune-cell signature of bacterial sepsis. Nat Med. (2020) 26:333–40. doi: 10.1038/s41591-020-0752-4

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Bossel Ben-Moshe N, Hen-Avivi S, Levitin N, Yehezkel D, Oosting M, Joosten L, et al. Predicting bacterial infection outcomes using single cell RNA-sequencing analysis of human immune cells. Nat Commun. (2019) 10:3266. doi: 10.1038/s41467-019-11257-y

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Wu YY, Wang SH, Wu CH, Yen LC, Lai HF, Ho CL, et al. In silico immune infiltration profiling combined with functional enrichment analysis reveals a potential role for naïve B cells as a trigger for severe immune responses in the lungs of COVID-19 patients. PLoS ONE. (2020) 15:e0242900. doi: 10.1371/journal.pone.0242900

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Novikov A, Cardone M, Thompson R, Shenderov K, Kirschman KD, Mayer-Barber KD, et al. Mycobacterium tuberculosis triggers host type I IFN signaling to regulate IL-1β production in human macrophages. J Immunol. (2011) 187:2540–7. doi: 10.4049/jimmunol.1100926

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Cho S, Mehra V, Thoma-Uszynski S, Stenger S, Serbina N, Mazzaccaro RJ, et al. Antimicrobial activity of MHC class I-restricted CD8⁺ T cells in human tuberculosis. Proc Natl Acad Sci USA. (2000) 97:12210–5. doi: 10.1073/pnas.210391497

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Odia T, Malherbe ST, Meier S, Maasdorp E, Kleynhans L, Du Plessis N, et al. The peripheral blood transcriptome is correlated with pet measures of lung inflammation during successful tuberculosis treatment. Front Immunol. (2020) 11:596173. doi: 10.3389/fimmu.2020.596173

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Abi Hussein H, Geneix C, Petitjean M, Borrel A, Flatters D, Camproux AC. Global vision of druggability issues: applications and perspectives. Drug Discov Today. (2017) 22:404–15. doi: 10.1016/j.drudis.2016.11.021

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Moumbock AFA, Li J, Mishra P, Gao M, Günther S. Current computational methods for predicting protein interactions of natural products. Comput Struct Biotechnol J. (2019) 17:1367–76. doi: 10.1016/j.csbj.2019.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Kozakov D, Grove LE, Hall DR, Bohnuud T, Mottarella SE, Luo L, et al. The FTMap family of web servers for determining and characterizing ligand-binding hot spots of proteins. Nat Protoc. (2015) 10:733–55. doi: 10.1038/nprot.2015.043

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Ozdemir ES, Halakou F, Nussinov R, Gursoy A, Keskin O. Methods for discovering and targeting druggable protein-protein interfaces and their application to repurposing. Methods Mol Biol. (2019) 1903:1–21. doi: 10.1007/978-1-4939-8955-3_1

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Fauman EB, Rai BK, Huang ES. Structure-based druggability assessment–identifying suitable targets for small molecule therapeutics. Curr Opin Chem Biol. (2011) 15:463–8. doi: 10.1016/j.cbpa.2011.05.020

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Gottesfeld JM, Neely L, Trauger JW, Baird EE, Dervan PB. Regulation of gene expression by small molecules. Nature. (1997) 387:202–5. doi: 10.1038/387202a0

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Huang L, Ye K, Mcgee MC, Nidetz NF, Elmore JP, Limper CB, et al. Interleukin-2-inducible T-cell kinase deficiency impairs early pulmonary protection against mycobacterium tuberculosis infection. Front Immunol. (2019) 10:3103. doi: 10.3389/fimmu.2019.03103

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Lundtoft C, Afum-Adjei Awuah A, Rimpler J, Harling K, Nausch N, Kohns M, et al. Aberrant plasma IL-7 and soluble IL-7 receptor levels indicate impaired T-cell response to IL-7 in human tuberculosis. PLoS Pathog. (2017) 13:e1006425. doi: 10.1371/journal.ppat.1006425

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Rochman Y, Spolski R, Leonard WJ. New insights into the regulation of T cells by gamma(c) family cytokines. Nat Rev Immunol. (2009) 9:480–90. doi: 10.1038/nri2580

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Mahon RN, Sande OJ, Rojas RE, Levine AD, Harding CV, Boom WH. Mycobacterium tuberculosis ManLAM inhibits T-cell-receptor signaling by interference with ZAP-70, Lck and LAT phosphorylation. Cell Immunol. (2012) 275:98–105. doi: 10.1016/j.cellimm.2012.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Bughani U, Saha A, Kuriakose A, Nair R, Sadashivarao RB, Venkataraman R, et al. T cell activation and differentiation is modulated by a CD6 domain 1 antibody Itolizumab. PLoS ONE. (2017) 12:e0180088. doi: 10.1371/journal.pone.0180088

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Cai Y, Dai Y, Wang Y, Yang Q, Guo J, Wei C, et al. Single-cell transcriptomics of blood reveals a natural killer cell subset depletion in tuberculosis. EBioMedicine. (2020) 53:102686. doi: 10.1016/j.ebiom.2020.102686

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Ju Y, Kim S-S, Lee K, Park S, Jin H. Association of the Genetic Polymorphisms for CD247 Gene and Tuberculosis Case. Biomed Sci Lett. (2020) 26:22–7. doi: 10.15616/BSL.2020.26.1.22

CrossRef Full Text | Google Scholar

64. Gebremicael G, Kassa D, Quinten E, Alemayehu Y, Gebreegziaxier A, Belay Y, et al. Host gene expression kinetics during treatment of tuberculosis in HIV-coinfected individuals is independent of highly active antiretroviral therapy. J Infect Dis. (2018) 218:1833–46. doi: 10.1093/infdis/jiy404

PubMed Abstract | CrossRef Full Text | Google Scholar