Identification and Classification of Differentially Expressed Genes and Network Meta-Analysis Reveals Potential Molecular Signatures Associated With Tuberculosis

Tuberculosis (TB) is one of deadly transmissible disease that causes death worldwide; however, only 10% of people infected with Mycobacteriumtuberculosis develop disease, indicating that host genetic factors may play key role in determining susceptibility to TB disease. In this way, the analysis of gene expression profiling of TB infected individuals can give us a snapshot of actively expressed genes and transcripts under various conditions. In the present study, we have analyzed microarray data set and compared the gene expression profiles of patients with different datasets of healthy control, latent infection, and active TB. We observed the transition of genes from normal condition to different stages of the TB and identified and annotated those genes/pathways/processes that have important roles in TB disease during its cyclic interventions in the human body. We identified 488 genes that were differentially expressed at various stages of TB and allocated to pathways and gene set enrichment analysis. These pathways as well as GSEA’s importance were evaluated according to the number of DEGs presents in both. In addition, we studied the gene regulatory networks that may help to further understand the molecular mechanism of immune response against the TB infection and provide us a new angle for future biomarker and therapeutic targets. In this study, we identified 26 leading hubs which are deeply rooted from top to bottom in the gene regulatory network and work as the backbone of the network. These leading hubs contains 31 key regulator genes, of which 14 genes were up-regulated and 17 genes were down-regulated. The proposed approach is based on gene-expression profiling, and network analysis approaches predict some unknown TB-associated genes, which can be considered (or can be tested) as reliable candidates for further (in vivo/in vitro) studies.

Introduction

Tuberculosis (TB) is a communicable disease generally caused by the bacterium Mycobacterium tuberculosis (MTB). The bacteria typically affects the lungs (pulmonary TB), but other body parts can be also affected (extra pulmonary TB) (Ahmad, 2011). The disease is communicable and spread through the air by expelling out the active MTB while coughing and sneezing (Schnappinger et al., 2003). In 2016, 10.4 million individuals were infected with TB, and 1.7 million died from the disease (including 0.4 million among people with HIV i.e., 40% of HIV deaths were due to TB). TB kills more adults in India than any other infectious disease (In 2016, an estimated 28 lakh cases occurred, and 4.5 lakh people died due to TB). India has the highest burden of both TB and advanced TB (like MDR TB) and second highest of HIV-associated TB. In India, the major challenge to curb the TB is poor primary healthcare system in rural areas due to deregulation of private health care leading to indiscriminate use of first- and second-line TB drugs, poverty, spreading HIV infection, and lack of administrative coordination among government functionary bodies. In our current study, we used meta-analysis of individual raw microarray data (GSE series) deposited in the GEO database which are obtained from various blood cell types (macrophages, monocytes, and CD4+) and cell lines (THP1) of individuals with different datasets (e.g., controls vs. TB disease, control vs. latent TB, latent TB vs. TB disease). We have performed gene-transition study of differentially expressed genes (DEGs) data and text mining between different stages of TB and classified the DEGs in stage-specific manner like normal to latent TB infections, normal to active TB, and latent infection to active TB. It is widely thought that, to know that function of a gene, it must be analyzed in the context of gene interaction network, because gene networks are commonly interpreted as encoding functional information in their connections. So, our study usually focused on the “guilt-by-association (GBA)” presumption which state that physically and functionally linked genes are possibly participating in the same biological pathways having comparable effects on the phenotypes (Lee et al., 2011). The concept of network theory is an imperative method to know the topological properties and the complex-system dynamics that correspond to their functional modules. The complex networks may be classified into four types of networks: (a) scale-free network, (b) small-world network, (c) random network, and (d) hierarchical network. For the biologist, hierarchical network has special interest because it incorporates the mien of modules, and distributed hubs (sparsely) regulate the network. So, we constructed the gene regulatory network and then analyzed its topological properties, because it helps to understand the structure of a network which facilitates in understanding the hidden mechanisms. Further, we identified important network modules which contain fundamental key regulators that have a fundamental importance.

Material and Method

Inclusion Criteria for Differentially Expressed Genes

A set of 11 microarray data sets were selected from the NCBI GEO repository database (http://www.ncbi.nlm.nih.gov/geo/). The datasets include GSE54992 (Cai et al., 2014), GSE52819 (Verway et al., 2013), GSE64335 (Bouttier et al., 2016), GSE11199 (Thuong et al., 2008), GSE98750 (Lavalett et al., 2017), GSE78233 (Cheng et al., 2017), GSE57736 (Guerra-Laso et al., 2015), GSE16250 (Reyes et al., 2015), GSE27882 (Liu et al., 2012), GSE57028 (Salamon et al., 2014), and GSE83456 (Blankley et al., 2016). The background data correction and data normalization were done by the robust multiarray average (RMA) in R affy and lumi packages to ensure unbiased and dysregulated gene expression data. The RMA method, which performs quantile normalization, was used to minimize inconsistencies due to normalization of the each Affymetrix GSE series. The RMA method was chosen over others due to its fine differential change detection and stable variance on a log scale and minimizes the false positive results. The specificity and sensitivity of RMA method are good during fold-change calculation to identify DEGs. Similarly, we have used lumi pipeline (Bioconductor package) which is exclusively developed to analyze Illumina data (BeadChip). It checks the data quality and data normalization and stabilizes the data variance. The gene expression data (GSE57736 and GSE27882) which is generated by Agilent platforms were already present as normalized data. In the present study, we used linear model for microarray analysis (LIMMA) package which is a highly recommended method to measure the differential expression of genes. It not only calculates simple t-test but also calculates moderate t-test and f-test by applying the Empirical Bayes approach and reducing the standard errors and gives us steady and reproducible outcomes even with a less quantity of arrays. So, limma package (Smyth, 2004) was used to identify the DEGs among the healthy control, latent TB, and active TB. We set the criteria to select significant differentially expressed genes i.e., the adjusted p-value was <0.05, and the fold change was >1.5. We have used the BRCW [Bioinformatics & Research Computing website (http://jura.wi.mit.edu/bioc/tools/compare.php)] to select the DEGs which is common in at least two datasets of gene expression profile. Doing this, we became more specific in the selection of DEGs and the chances of biased data compilation became negligible. According to the corresponding correlation between the probe and gene from the data, the probe numbers of the expression profile were changed into the corresponding gene symbols using the Synergizer database (Huang et al., 2009).

Gene Classification, Ontology, and Pathway Analysis of Degs

To know the significance of the identified DEGs, we have categorized them by GO-molecular function, GO-biological process, and protein class using the Protein ANalysis THrough Evolutionary Relationships (PANTHER v.13.0) Classification System and analysis tools and DAVID database (https://david.ncifcrf.gov/) to enrich the given set of DEGs to possible GO terms (Mi et al., 2013) (Huang et al., 2009). The PANTHER overrepresentation study (Fisher’s exact with FDR multiple test correction) was used to search the data against the PANTHER, and GO databases and p-values were set according to Bonferroni correction.

The GO analysis (gene ontology) is the useful method for annotation of genes and its products and characterization of biological attributes for high-throughput genome or transcriptome data (Ashburner et al., 2000). The differentially expressed genes among “active-TB” and “latent TB” on compared with “normal condition” were over-represented in various GO classes. The gene ontology provides and visualizes us the basic terms subdivided into three important categories, namely, BP (Biological process), MF (molecular functions), and biological pathways among these groups.

Gene Transition Among Different Stages of TB

In order to get the behavior of normal gene expression perturbation, we tried a normal way for finding the associated genes while moving from one stage to another. In all the transitions that we took into consideration have been provided by a list of up- and down-regulated genes. These gene(s) in both the cases (i.e., stage from which is transferred to the targeted stage) have its own meaning. This meaning to a gene(s) regulation gives us an opportunity to say something about the ongoing mechanobiology inside the cell. So, to observe these transitions, we framed our study in such a way (based on the data we’ve got) which considers every possible transition in view. These transitions are discussed in brief as follows. In this study, we made comparison of the gene-expression profiles among individuals with normal conditions, latent infection, and active TB. Thus, we observed the expression of genes from normal to different stages of the TB and try to arrest those genes which play a key role in susceptibility to TB.

● Normal to Latent Infections: In this section, we have taken those DEGs which are involved in between normal and latent TB conditions. In LTBI, the bacteria remain in inactive form for many years (years to decades) before transforming into TB disease. In this study, we have studied several BPs and important pathways that lead to the latent infection for identification of latently infected individuals.

● Normal to Active TB: In this section, we have taken those DEGs which are differentially expressed in TB disease condition as compared to healthy control and recognized those immune process and pathways which are prominent in TB disease.

● Latent Infection to Active TB: In this section, we have taken those DEGs which are involved in between latent TB and active TB diseases. The individual with LTBI eventually reactivates and becomes infectious, seriously influencing epidemiological situation. Mechanisms of MTB transition to dormancy and TB reactivation are inadequately understood, and biomarkers of latency remain largely mysterious (Kondratieva et al., 2014).

Topological Characterization of Networks

The classified genes of various stages of TB were used to construct their regulatory networks. The networks were built by STRING database (https://string-db.org/) and then visualized in Cytoscape v3.4 (Shannon et al., 2003). As the network was built, the first and foremost basic analysis are its topological properties. Topological analysis helps to understand the structure of a network which facilitates in understanding the hidden mechanisms. The following networks properties were analyzed to seek the important behaviors of the network:

● Degree distribution: The degree (k) of a node (gene) in a biological network is the number of links with other nodes, the probability distribution of these degrees called degree distribution (P[k]).

$\begin{array}{ll}P\left(k\right)=\frac{n_k}{N}& (1)\end{array}$

wheren_k =no. of nodes with degree k.

N= the total number of nodes in the network.

● Neighborhood connectivity: A node (gene) has number of neighbors and is considered as a connectivity of node. So, neighborhood connectivity (C_N[K]) is:

$\begin{array}{ll}{C}_N\left(k\right)={\displaystyle {\sum }_{q}qP}\left(\frac{q}{k}\right)& (2)\end{array}$

where $$P\left(\frac{q}{k}\right)$$ = the conditional probability.

● Clustering co-efficient: The ratio of number of edges (e_i) between the node’s neighbor or maximum numbers of edges that could possibly exist between the nodes. So, the total network cluster coefficient is the average of cluster coefficient of all nodes (i^th ) in the network.

$\begin{array}{ll}C\left(k_i\right)=\frac{2e_i}{k_i\left(k_i-1\right)}& (3)\end{array}$

● Betweenness centrality: In the complex network, a node’s betweenness centrality (CB) represents the prominence of information flow through one node to another via shortest path (Brandes, 2001) (Mason and Verwoerd, 2007). From node (i) to node (j), the geodesic paths are shown by “dij(v)” passing via node “v” and “dij.”

$\begin{array}{ll}{C}_B\left(v\right)={\displaystyle {\sum }_{i,j,i\ne j\ne k}\frac{d_{ij\left(v\right)}}{d_{ij}}}& (4)\end{array}$

● Closeness centrality: In the network, how quick info is circulated from one node (i) to another (j) is measured by closeness centrality (C_C). Closeness centrality is given by:

$\begin{array}{ll}{C}_C\left(i\right)=\frac{n}{{\displaystyle {\sum }_{j}dij}}& (5)\end{array}$

● Eigenvector centrality: In the network, eigenvector centrality of a node “i(C_E[i])” is proportionate to the total of i’s neighbor centralities (Bonacich, 1987).

$\begin{array}{ll}{C}_E\left(i\right)=\frac{1}{\lambda }{\displaystyle {\sum }_{j=nn\left(i\right)}{v}_j}& (6)\end{array}$

where nn(i) = closest neighbors of nodes (“i”).

λ = eigenvalue of the eigenvector.

v_i = “Av_i=λv_i” where “A” (adjacency matrix).

Community Analysis: Leading Eigenvector Method

In hierarchical network, to distinguish the nature of modular and its properties is important to explain the activities of network at various levels of hierarchy. In our study, the leading eigenvector method (LEV) (Newman, 2006a; Newman, 2006b) was used to detect the communities in R from package “igraph” (Gabor and Nepusz, 2006) (the community detection script is given Supplementary Data 3).

● Modularity: Modularity determines to measure the strength of division of a network into clusters or communities (Newman and Girvan, 2004).

$\begin{array}{ll}Q=\frac{1}{2m}{\displaystyle {\sum }_{ij}\left(A_{ij}-\frac{k_i{k}_j}{2m}\right)}\delta \left(C_i,C_j\right)& (7)\end{array}$

where m = Total no. of edges.

A_ij = adjacency matrix of size “i × j.”

k = degrees.

δ = function yields 1 if nodes “i” and “j” are in the same community.

Genes Tracing

To access the regulation of network, we first tried to find out the most influential nodes (genes) within the network. The gene tracing (up to motif level) was done purely on the appearance of the respective genes in various submodules obtained from the clustering. Then, these genes were used to get the picture of changes in the network organization in their absence.

Hub Gene Knockout

To know the changes of organization within the complex network in the absence of most influencing gene (node), we must remove constructed leading hubs (rich clubs) in the networks. We consecutively eliminated all the important hubs (one by one) from each network and measured the topological properties of this reorganized network to observe the regulating abilities of these hubs by calculating the degree of structural change due to their absence.

LCP-DP Approach to Estimate the Network Compactness

The LCP-decomposition-plot (LCP-DP) provides a way to characterize the topological properties of the network in 2D space of common neighbors (CN) index of connecting nodes and local community links (LCL) of each pair of interacting nodes in the network (Cannistraci et al., 2013). The LCP correlation defined by:

$\begin{array}{ll}LCP-corr=\frac{cov\left(CN,LCL\right)}{{\sigma }_{CN}{\sigma }_{LCL}}withCN>1& (8)\end{array}$

where cov(CN, LCL) = the covariance among LCL and CN.

σ_LCL and σ_CN = standard deviations.

(Note: MATLAB script is given in Supplementary Data 2).

Hamiltonian Energy Estimation: Energy Distribution in the Network

The HE is used to organize a network at a certain level by following the formalism of constant Potts model (Traag et al., 2011; Traag, 2013). The energy distribution at the global and modular levels of the network is given by HE. HE of a network can be measured by:

$\begin{array}{ll}{H}^{\left[c\right]}=-{\displaystyle {\sum }_c[e_c-\gamma n_{c }^2]}& (9)\end{array}$

where e_c = no. of edges.

n_c = no. nodes.

“c” and “γ” = the resolution parameter

Results

Differentially Expressed Genes (DEGs) Classification and Overrepresentation Analysis

A total of 5, 680 differentially expressed genes (DEGs) were identified after the extensive analysis of all the 11 GSE series, of which 2,660 were up-regulated, and 3,020 were down-regulated genes. These differentially expressed genes were clustered according to “GO-MF (molecular function),” “GO-BP (biological process),” and “PANTHER protein class” shown in Figure 1. All these predictive DEGs showed a broad range of protein classes which involved in various processes. The helicases and nucleases were found within the “nucleic acid binding” protein class. The “enzyme modulator” category features kinase, G-protein, phosphatase, and protease modulators. The structural motif and nuclear hormone receptors are part of the “transcription factor” protein class. The “hydrolases” is a sub-category of proteases and phosphatases. The “receptor” protein class includes cytokine receptors, protein-kinase receptors, ligand-gated ion channels, nuclear-hormone receptors, and G-protein-coupled receptors. Besides of these above protein classes, signaling molecules, transferase, oxidoreductase, and transporter are also most abundant protein classes. The two most abundant GO-biological process groups—“metabolic process” and “cellular process” which are not astonishing as these processes carry those genes which involved in the most basic life processes. The “cellular process” includes cell cycle, cell–cell signaling, cell component movement, cytokinesis, and proliferation. The “metabolic process” section includes metabolism of lipid, protein, carbohydrate, cellular amino acid, and nucleobase-containing compound metabolism. The “biological regulation” includes metabolism, cell cycle, the regulation of apoptosis, catalytic activity, translation, and homeostasis. Further, these biological process responses to stimulus, localization, and developmental process also represent the significant number of proteins; the complete categories details are given in Supplementary Data 1. To know the probability of largely occupied protein classes and GO categories among the differentially expressed genes, we used PANTHER’s overrepresentation analysis. When we compared with the reference genome, we found that any of the richest classes are overrepresented in the data (Table 1). The most abundant protein classes “chemokine,” “cytokine,” “hydrolase,” “ribosomal protein,” and “RNA-binding protein” were enriched along with the classes “signaling molecule” and “cell adhesion molecule.” The highly populated GO-biological processes were enriched in “cytokine-mediated signaling pathway,” “response to external stimulus,” “immune response,” “locomotion,” “signal transduction,” “cell communication,” “developmental process,” “cellular process,” and “MAPK cascade.” Similarly, the most abundant GO-molecular functions were enriched in “chemokine activity,” “cytokine activity,” “cytokine receptor binding,” “oxidoreductase activity,” and “receptor binding.”

FIGURE 1

Figure 1 Functional Classification of Differentially expressed genes from various GSE series associated with TB disease according to (A) Molecular functions (Transporter activity [(U¹⁸⁸,D₂₈₆), Translation regulator activity(U⁰³,D₀₉), Catalytic activity(U⁹⁷⁸,D₈₀₇), Channel regulator activity(U⁰³,D₁₃) Receptor activity(U²¹⁹,D₂₉₅), Signal transducer activity(U¹²⁷,D₁₄₄), Antioxidant activity(U⁰⁸,D₁₂), Structural molecule activity(U⁹⁶,D₁₀₁), Binding(U⁹⁵⁷,D₁₁₅₀) and unclassified(U⁵³,D₂₅)], (B) Biological processes [(cellular component organization or biogenesis(U¹⁰⁸,D₁₃₄), cellular process(U⁷⁰⁹,D₇₉₀), localization(U¹⁹⁸,D₂₀₇), reproduction(U¹⁹,D₂₃), biological regulation(U³¹⁰,D₃₈₃), response to stimulus(U²⁶⁷,D₂₉₉), developmental process(U¹⁷²,D₂₁₃), rhythmic process(U⁰¹,D₀₀), multicellular organismal process (U¹²⁴,D₁₂₉), locomotion(U⁴⁸,D₅₇), biological adhesion(U⁴⁴,D₃₄), metabolic process(U⁴⁵⁶,D₄₇₁), growth(U¹⁰,D₀₄), immune system process(U⁴⁷,D₄₈), cell killing (U⁰⁰,D₀₇) and unclassified(U¹¹⁹,D₄₃)] and (C) Protein classes [(extracellular matrix protein(U⁵²,D₃₄), cytoskeletal protein(U¹⁰¹,D₁₃₉), transporter(U¹²⁶,D₁₄₂), transmembrane receptor regulatory/adaptor protein(U¹¹,D₀₉), transferase(U²²⁶,D₂₂₇), oxidoreductase(U¹³²,D₇₇), lyase(U²¹,D₂₈), cell adhesion molecule(U⁷⁶,D₇₄), ligase(U⁶³,D₄₈), nucleic acid binding(U²⁵⁰,D₃₁₈), signaling molecule(U²⁰⁵,D₃₅₂), enzyme modulator(U²⁶⁷,D₃₀₄),calcium-binding protein(U⁶⁹,D₆₃), defense/immunity protein(U⁵⁸,D₇₄), hydrolase(U³⁶³,D₃₀₄), transfer/carrier protein(U⁵⁰,D₃₁), membrane traffic protein(U⁵²,D₆₀), transcription factor(U²⁰⁰,D₂₂₂), chaperone(U¹⁹,D₂₃), cell junction protein(U²⁴,D₂₈), surfactant(U⁰³,D₀₀), structural protein(U²⁷,D₂₀), storage protein(U⁰³,D₀₃), isomerase(U¹⁸,D₁₁), receptor(U²¹⁶,D₂₄₁)]. (U, Up-regulated; D, Down-regulated).

TABLE 1

Table 1 Overrepresented PANTHER protein class and GO ontology categories of all differentially expressed genes.

Gene-Transition and GO Enrichment Analysis

All the 5,680 differentially expressed genes (DEGs) were divided into three different groups, normal vs. latent infection, normal vs. active TB, and latent infection vs. active TB. We filtered out a total number of 488 DEGs which are listed in (Table 2), and then, we also identified few genes which are commonly differentially expressed among various stages of TB (including NC vs. LI, NC vs. ATB, and LI vs. ATB) shown in Figure 2. Further, we isolated up- and down-regulated genes from each stage and performed GO-enrichment analysis using DAVID tool (v6.7). It is well preferred annotation database which used novel algorithms to extract biological information from large gene lists. The DEGs are enriched with a certain biological process or molecular function that are given in Table 3. Further, the most enriched pathways associated with these differentially expressed genes were identified, given in Table 4.

● Normal to Latent infections: In the early infection (latent infection), we have identified 28 genes which are differentially expressed in the host system, of which 12 genes (IER5L, MS4A6A, DOK2, FZD2, NCKI-ASI, SNHG12, NLRC4, XPO7, SMA4, CD36, AFFI, and NDUFS8) were up-regulated, and 16 genes (IL1A, IL6, ACOD1, IL1B, ELOVL7, PTGS2, EREG, F3, IFIT1, TNF, KANK1, CCL4, CXCL11, PTX3, IRAK2, and AREG) were down-regulated. The GO analysis of these DEGs are enriched in “immune response,” “inflammatory response,” “chemokine activity,” “cellular process,” “biological regulation,” “regulation of metabolic process,” “protein binding,” “catalytic activity,” “bindings” etc. On the pathway analysis, most of the DEGs were found to be down-regulated (beneficial to pathogen) and enriched in very important pathways like toll-like receptors, NF-kappaB signaling, cytokine–cytokine receptor interaction, MAPK signaling pathway, TB, and TNF signaling.

● Normal to Active TB diseases: A total of 266 differentially expressed genes (DEGs) were identified, of which 149 were up-regulated, and 117 were down-regulated among the various cases. In this condition, the DEGs were enriched in “inflammatory response,” “immune response,” “signal transduction,” “response to stimulus,” “apoptotic process,” “cellular process,” “metabolic process,” “binding,” “catalytic activity,” “receptor activity” etc. For the up-regulated genes, the most abundant pathways were “cytokine–cytokine receptor interaction,” “chemokine signaling,” “toll-like receptor signaling pathway,” “NF-kappa B signaling pathway,” “transcriptional misregulation in cancer,” and “pathways in cancer.” Similarly, for down-regulated genes, the most enriched pathways were “cell cycle,” “chemokine signaling pathway,” “NF-kappa B signaling pathway,” “TNF signaling pathway,” “PI3K-Akt signaling pathway,” “metabolic pathways,” and “MAPK signaling pathway.”

● Latent infection to Active TB Disease:: A total of 127 differentially expressed genes (DEGs) were found to be up-regulated, and 69 were down-regulated. These differentially expressed genes were enriched in “binding,” “catalytic activity,” “receptor activity,” “transporter activity,” “cellular and metabolic process,” “biological regulation,” “immune response,” and “oxidoreductase activity” including others. The most abundant enriched pathway for up-regulated genes belonged to “metabolic pathway,” “NOD-like receptor signaling pathway,” “regulation of actin cytoskeleton,” “biosynthesis of antibiotics,” “metabolism of xenobiotics by cytochrome p450,” “thyroid hormone synthesis,” “regulation of autophagy,” and “peroxisome,” and for down-regulated genes, the enriched pathways were “chemokine signaling pathway,” “TNF signaling pathway,” “RAP1 signaling pathway,” “PI3K_AKT signaling,” “metabolic process,” “apoptotic,” “focal adhesion,” “MAPK signaling pathway,” “microRNAs in cancer,” and “amoebiases.”

FIGURE 2

Figure 2 Venn diagram showing the number of common genes (in intersect region) which are differentially expressed among the normal vs. latent infection, normal vs. active TB, and latent infection vs. active TB.

TABLE 2

Table 2 The total number of DEGs from various GSE series associated with TB disease.

TABLE 3

Table 3 Gene set enrichment analysis of differentially expressed genes (DEGs) among active TB, LTBI, and normal condition.

TABLE 4

Table 4 Pathways enriched by differentially expressed genes (DEGs) among TB, LTBI, and healthy control (HC).

Gene Interaction: Hierarchical Scale-Free Network

The classified genes of various stages of TB (Table 2) were used to construct their regulatory network. We constructed six networks for Up and Down-regulated genes separately. The topological parameters of the networks obey power law distributions. The probability of clustering co-efficient C(k), degree distributions P(k), and neighborhood connectivity C_N(k) exhibits fractal nature as shown in Figure 3. The results for all the networks are summarized as follows:

$$\begin{array}{ll}\mathbf{Normal to latent Inf.}\hfill & \mathbf{Normal to active TB}\hfill \\ \left[\begin{array}{c}P\\ C\\ C_N\end{array}\right]\sim \left[\begin{array}{c}{k}^{-A}\\ k^{-B}\\ k^{+C}\end{array}\right]\to \left[\begin{array}{c}0.418\\ 0.218\\ 2.099\end{array}\right]\hfill & \left[\begin{array}{c}P\\ C\\ C_N\end{array}\right]\sim \left[\begin{array}{c}{k}^{-A}\\ k^{-B}\\ k^{+C}\end{array}\right]\to \left[\begin{array}{c}0.459\\ 0.153\\ 1.755\end{array}\right]\hfill \\ \left[\begin{array}{c}P\\ C\\ C_N\end{array}\right]\sim \left[\begin{array}{c}{k}^{-A}\\ k^{-B}\\ k^{+C}\end{array}\right]\to \left[\begin{array}{c}0.315\\ 0.413\\ 1.33\end{array}\right]\hfill & \left[\begin{array}{c}P\\ C\\ C_N\end{array}\right]\sim \left[\begin{array}{c}{k}^{-A}\\ k^{-B}\\ k^{+C}\end{array}\right]\to \left[\begin{array}{c}0.533\\ 0.354\\ 0.983\end{array}\right]\hfill \\ \mathbf{Latent Inf. to active TB}\hfill & \hfill \\ \left[\begin{array}{c}P\\ C\\ C_N\end{array}\right]\sim \left[\begin{array}{c}{k}^{-A}\\ k^{-B}\\ k^{+C}\end{array}\right]\to \left[\begin{array}{c}0.732\\ 0.158\\ 1.549\end{array}\right]\hfill & \to \text{\hspace{1em}For Up-regulated genes}\hfill \\ \left[\begin{array}{c}P\\ C\\ C_N\end{array}\right]\sim \left[\begin{array}{c}{k}^{-A}\\ k^{-B}\\ k^{+C}\end{array}\right]\to \left[\begin{array}{c}0.719\\ 0.563\\ 1.093\end{array}\right]\hfill & \to \text{\hspace{1em}For Down-regulated genes}\hfill \end{array}$$

FIGURE 3

Figure 3 Topological properties of all the six networks. The behaviors of degree distributions (P[k]), clustering co-efficient (C[k]), neighborhood connectivity (C_N[k]), betweenness (C_B[k]), closeness (C_C[k]), and eigenvector (C_E[k]) measurements as a function of degree k. The lines are fitted lines with power laws in the data sets.

Our network behavior indicates hierarchical scale free network (Barabási and Oltvai, 2004) (Pastor-Satorras et al., 2001) (Ravasz and Barabási, 2003). The power law fits on the data points of the network’s topological parameters were done and confirmed by following the standard statistical fitting method given by Clauset et al. (Clauset et al., 2009). The p values for all data sets were calculated (against 2,500 random samplings) and found to be greater than 0.1, and data fitting goodness was less than 0.1. The values of P(k) and C(k) were negative, which implies that the network follows a hierarchical pattern, and C_N(k) was positive, which means that the network follows the assortativity that identifies a huge cluster of degree-nodes (rich club), which regulates the network. The centrality parameters: betweenness (C_B), closeness (C_C), and eigenvector (C_E) centralities of the network also showed fractal behavior, and good connectivity of nodes in a network is distinguished by eigenvector or centrality C_E(k). It calculates the effectiveness of the spreading (receiving) power of data of nodes from the network. These properties follow the power law behaviors as follows: -

$$\begin{array}{ll}\mathbf{Normal to latent Inf.}\hfill & \mathbf{Normal to active TB}\hfill \\ \left[\begin{array}{c}{C}_B\\ C_C\\ C_E\end{array}\right]\sim \left[\begin{array}{c}{k}^E\\ k^F\\ k^D\end{array}\right]\to \left[\begin{array}{c}2.099\\ 0.115\\ 1.019\end{array}\right]\hfill & \left[\begin{array}{c}{C}_B\\ C_C\\ C_E\end{array}\right]\sim \left[\begin{array}{c}{k}^E\\ k^F\\ k^D\end{array}\right]\to \left[\begin{array}{c}1.755\\ 0.117\\ 1.056\end{array}\right]\hfill \\ \left[\begin{array}{c}{C}_B\\ C_C\\ C_E\end{array}\right]\sim \left[\begin{array}{c}{k}^E\\ k^F\\ k^D\end{array}\right]\to \left[\begin{array}{c}2.634\\ 0.269\\ 0.989\end{array}\right]\hfill & \left[\begin{array}{c}{C}_B\\ C_C\\ C_E\end{array}\right]\sim \left[\begin{array}{c}{k}^E\\ k^F\\ k^D\end{array}\right]\to \left[\begin{array}{c}2.728\\ 0.107\\ 1.048\end{array}\right]\hfill \\ \mathbf{Latent Inf. to active TB}\hfill & \hfill \\ \left[\begin{array}{c}{C}_B\\ C_C\\ C_E\end{array}\right]\sim \left[\begin{array}{c}{k}^E\\ k^F\\ k^D\end{array}\right]\to \left[\begin{array}{c}1.549\\ 0.111\\ 1.021\end{array}\right]\hfill & \to \text{\hspace{1em}For Up-regulated genes}\hfill \\ \left[\begin{array}{c}{C}_B\\ C_C\\ C_E\end{array}\right]\sim \left[\begin{array}{c}{k}^E\\ k^F\\ k^D\end{array}\right]\to \left[\begin{array}{c}1.805\\ 0.131\\ 0.987\end{array}\right]\hfill & \to \text{\hspace{1em}For Down-regulated genes}\hfill \end{array}$$

Identification of Key Regulators and Properties

Since the popularity of leading hubs gets changed according to the gene activities and its regulation, we can’t say that all the predictive leading hubs are key regulators for disease progression, but some of these hubs can play a significant role, which we called them as fundamental key regulators (FKR). The structure of modular and its arrangement have been performed through Newman and Girvan’s standard community finding algorithm at various levels of the organization. Using this community finding algorithm, we found that our six networks are hierarchically organized at various levels (Supplementary Data 2). The Hamiltonian energy (HE) and corresponding modularity (Q_N) are reduced as one goes from top to down organization of network (Figure 4).

FIGURE 4

Figure 4 Energy and modularity distribution in all networks quantified by Hamiltonian (HE) and modularity calculation as a function of network levels.

We have calculated the probability to know the regulating ability of our 31 key regulators:

$$P_x\left(y^{\left[l\right]}\right)=\frac{y^{\left[l\right]}}{E^{\left[l\right]}}$$

where,

x = no. of edges (y^[l])

(E^[l]) = Total no. of edges of the network (modules and sub-modules).

The measured probability P_x(y^l) of all the KR showed an increase in P_x as the level the increases from top to bottom direction. At deeper level of the organization, regulation of FKR increases, and their activities become more prominent. Thus, these FKR becomes backbone of the network organization, stabilization, and active workers at grassroot level. The FKR are shown in Figures 5 and 6.

FIGURE 5

Figure 5 This figure shown the fundamental key regulator from various stage of TB obtained from main networks to motif/hub level through various modules/sub-modules at various level of organization. The probability distribution of the 14 up-regulated key genes as a function of level of organization.

FIGURE 6

Figure 6 This figure shown the fundamental key regulator from various stage of TB obtained from main networks to motif/hub level through various modules/sub-modules at various level of organization. The probability distribution of the 15 down-regulated key genes as a function of level of organization.

As the results, we identified a total of 31 FKR genes from various stages of TB. In normal to latent infections, we identified eight genes, of which three genes (FZD2, NDUFS8, NLRC4) were up-regulated, and another five genes (CCL4, IL1B, IL1A, TNF, AREG) were down-regulated. While in normal to active TB, we found 15 genes out of which eight genes (EIF2AK2, SAMD9L, IFI44L, DDX58, IFI44, HERC5, NFE2L3, LRRK2) were up-regulated, and seven genes (CSF3, MAP3K8, IRAK2, TNFAIP3, TRAF1, PLAUR, and CD44) were down-regulated, and similarly, for latent to active TB diseases, three genes (TST, MTHFD1, CARS2) were up-regulated, and five genes (ETV6, NFKB1, MET, DOCK4, PARD6G) were down-regulated. The functional interpretations of these 31 key regulators were annotated by GO-TermFinder (LAGO) (Boyle et al., 2004) using biological process terms at P-value 0.01; the results are interpreted in Figure 7. The gene specific pathways were identified by KEGG pathway database (Kanehisa et al., 2016), and the results are shown in Figure 8.

FIGURE 7

Figure 7 This heatmap shows the 31 hub genes with their involvement in various biological processes.

FIGURE 8

Figure 8 The molecular pathways associated with our key regulators were identified by KEGG pathway database. Out of the 31 key regulators, we did not find any pathways information of five genes namely: SAMD9L, IFI44L, IFI44, HERC5, and NF2L3.

Local Perturbations of Network

The gene knockout experiment of all the hubs/motif from the parent networks may able to highlight the local perturbations driven by these individual hub or motif, and their effect on global network properties. It has been revealed that the network is tolerant to hub’s deletion which implies that the important network elements are still remains after elimination of hubs at level 0 (parent network). One of main causes is that the P-P interactions networks are too dense to be broken into fragments by only removing hubs. However, the elimination of these hubs/motif from the complete network does cause significant variations in the network properties, where P(k) and C(k) change significantly in complete network level, whereas C_N(k) changes slightly. Likewise, the variations in the exponents of centrality measurements have also showed significant changes, as shown in Figure 9. Since, it is clear from the differences in the exponents of topological parameters that network perturbation increases when it goes to deeper level (top to down direction). In our case, most of the perturbation increases after the 3rd level; at this level elimination of key regulator from network almost breaks down the submodules existing in the deeper levels; such type of behavior shows that local perturbation of network is highest at deeper levels.

FIGURE 9

Figure 9 The changes in the exponents of the three important topological parameters (P[k]), C(k), and C_N[k]) due to hub genes knock-out experiment from parent network.

The Local-Community-Paradigm: Evidence of Self-Organization

The LCP architecture supports the quick transfer of data across several network modules and through local processing too. We have analyzed all the six networks to check self-organization behavior at different levels using LCP method. For different level, the calculated LCP-corr of all the modules or sub-modules are shown in Figures 10 and 11. The average values of LCP-corr (we ignored modules which having zero LCP-corr) were greater than 0.853 at each level. This shows that the network maintains compactness and self-organization and has efficient data processing. It characterizes robust LCP networks that are dynamic in nature and heterogeneous, which help in network re-organization and evolution.

FIGURE 10

Figure 10 This figure shows the LCP-corr of all the modules/sub-modules of up-regulated genes at various levels. The compactness characterized by 2D plots between √(LCL) versus CA.

FIGURE 11

Figure 11 This figure shows the LCP-corr of all the modules/sub-modules of down-regulated genes at various levels. The compactness characterized by 2D plots between √(LCL) versus CA.

Discussion

TB disease remains to be one of the most significant infectious diseases and a leading cause of death worldwide. Currently, the main challenge is to develop a delicate and effective method to identify the latent TB infection (LTBI) because 90% of cases of latent infection with MTB do not show any symptoms/signs but have a 10% lifetime possibility of transforming into active TB. As we know, many events have happened in response to TB infection cycle. According to the cycle of infection by Young et al., 2008, the exposure of microorganism up to its development to latent stage, there are stringent amount of processes that have appeared in the host cell. The long fight of MTB to establish its generation by stabilizing its environment for metabolism has evolved with an ability to overcome immunity. In this scenario of transferring from initial infection stage to latent stage has a profound increase in the overall immune response outside the host cell and an increased intensity of basic cellular machine establishment process for the pathogen cell. In case, the elimination eliminates the eliminator even after the induction of T-cell but able to suppress the effect; pathogen cells (Pcells) have two choices, either develop army or attack the adaption of immunity. In order to decide to enter the resting phase (i.e., latent), MTB’s survival instincts allow it to develop its metabolic environment by downregulating antigen presentation and the release of cytokines too (like IFNγ). On the other hand, if adaptive immunity has taken its chance, then MTB takes the defensive mode by subverting various normal cell cycle functioning inside the macrophages thus making a halt on the maturation of phagosomes.

To stop the disease epidemic, early diagnostic method or techniques are required. In this way, the gene expression profiling has uncovered the differences in the transcriptome among normal condition, latent infection, and active TB. These results not only revealed significant genetic biomarkers indicative of LTBI, TB-disease conditions but also recognized transcriptionally regulated genes that vary in biological functions. In the current study, we have identified differentially expressed genes (DEGs) from various stages of disease and found that most of the DEGs relative to normal vs. TB disease condition (including latent infection) are enriched in important pathway like-toll like receptor, NOD-like receptor signaling pathway, MAPK signaling pathway, TNF signaling, chemokine signaling pathway, PI3K-Akt signaling pathway, apoptosis etc. In fact, pathogen is identified by the receptors on the surface of immune cells, and toll-like receptors (TLRs) are one of them. Different TLRs including TLR2, 4, 9, and 8 play an important roles in TB infection (Faridgohar and Nikoueinejad, 2017). These receptors (toll-like, NOD-like etc.) are expressed irrespective of whether they participate in immune signaling or immunity. The interaction of MTB with these receptors (like TLRs) initiates an intercellular signaling cascade that culminates in a pro-inflammatory response (cytokines and chemokines that serve as a signal for infection). As we know that MAPK signaling pathway is important for transcription and non-transcription responses of the immune system, most of the pathogens during infection hijack the immune system by targeting the MAPK signaling pathway (Soares-Silva et al., 2016). The TNF signaling pathway is involved in the regulation of immune cells, and it can regulate immune responses (Kim, 2018). However, pathogen can target the TNF signaling pathway for its survival in the host cell. The PI3K-Akt signaling pathway plays an important role in apoptosis, autophagy, metabolism, cell growth, and differentiation. The expression of FoxP3 by inhibiting the activation of transcription factor Forkhead-O3a (Foxo1-3a) is negatively regulated by this pathway (Zhang et al., 2017). The FoxP3+Treg cell activation which will assist to set up a new target for the involvement of TB immunotherapy molecules as part of the immune-escape mechanism to provide a theoretical basis is inhibited by M. tuberculosis (Scott-Browne et al., 2007). On comparing between latent infection and TB disease, few additional pathways are involved, like regulation of actin cytoskeleton during MTB infection; to maintain the stability of the cytoskeleton, macrophages cells themselves are also trying to regulate cytoskeletal-associated proteins (Wang et al., 2015b). Thyroid hormones (hormones, T4 and T3) are produced by thyroid gland, which is essential for the regulation of metabolic processes throughout the body. The regulation of autophagy is important for host in response to invading mycobacteria; the host defense mechanism identifies pathogen motifs through innate receptors but also releases appropriate cytokines. The autophagic pathways are regulated by these innate signals by regulation of genes of this pathway during infection (Jo, 2013). Recently, a study reports that the enzyme Msm_ACTase (from bacteria) helps in scavenging increased amount of H₂O₂ due to upregulation of genes involved in peroxisome pathway which gives an insight into a new idea as to how the MTB bacteria surpasses the host defense in MTB infection (Ganguli, 2015). Moreover, macrophages are the main effector cells responsible for killing pathogen (MTB) via different mechanisms, including apoptosis. But important genes in this pathway were down-regulated by MTB that slow down or stop the apoptotic process for its survival in host system (Behler et al., 2015). Most of the microRNAs that were found deregulated in cancer were affected by MTB in a similar way as its infestations, like deletion, mutation, and epigenetic silencing (Garzon et al., 2009). The apoptotic pathways are important to kill the infected macrophage. During MTB infection, the pathogen try to control the timing and mode of host cell death for its persistence and replication (Schaaf et al., 2017)

In addition, we have built the networks which emphasize on genes that were regulated by network. The constructed network of classified genes from various stages of TB showed hierarchical nature, which indicated that the networks have system level organization including modules/sub-modules which are interconnected. Since the network’s nature is hierarchical, its synchronization confirms several important functional regulations of the network, but individual gene activities are not so important. In our networks (including up- and down-regulated genes), a total of 31 key regulators were identified by affecting motifs and module regulation, showing their biological importance and serve as the foundation of network activities and their regulations, and could be the most probable target gene for disease control. We performed extensive analysis of all the key regulators (31 genes) to find evidence of their association with TB or which are directly or indirectly involved in host immune response and confirmed by manually searching for evidence in the literatures as presented:

● Down-regulated genes: We found that TNF and IL1 (A and B) are key mediators present in severe inflammatory diseases. However, both IL-1 and TNF receptor pathways are important for the control of MTB infection, and it is critical to assess the respective role of IL1A, IL1B, and TNF (Bourigault et al., 2013) (Cavalcanti et al., 2012). The expression of CCL4 is high in the late phase of the active TB disease, but there are low levels of expression during early infection (Rangel-Santiago et al., 2016). It has been currently reported that AREG plays an important role in orchestrating both host resistance and tolerance mechanisms. Although AREG is known as epithelial cell-derived factor, the recent studies showed that AREG can be expressed by multiple populations of activated immune cells in inflammatory conditions (Zaiss et al., 2015). The NFKB1 activates the immune cells by up-regulating the expression of various cytokines in the host immune system (Ghosh et al., 1998). The gene CD44 plays an important role in innate and adaptive immune responses, in the acute inflammatory response to both infectious and sterile stimuli. Also during infection, CD44 may influence host defense by affecting phagocytosis (van der Windt et al., 2010). In our present study, we found that NFKB1 and CD44 are down-regulated which means that MTB is burglarized and seizing the host immune system by slowing down the expression of NFKB1. It is well known that MET (hepatocyte growth factor receptor) regulates several functions of immune cells, including cytokine production, differentiation and maturation, cellular migration and adhesion, and T-cell effector function, and HGF exerts anti-inflammatory activities through MET signaling (Sagi and Hieronymus, 2018). A group of researchers has elucidated that PLAUR domain containing Lypd8 that inhibits bacterial invasion of colonic epithelia (Okumura et al., 2016). The gene TRFA1 has an active or passive interaction with many tumor necrosis factor receptor (TNFR), and it has been also shown that TRAF1 plays as an antiapoptotic role in lymphoma cells by activation of NFKB (Wan et al., 2016). Further, the gene TNFAIP3 (A20) is a cytoplasmic zinc-finger protein which works as a negative-feedback regulator of NFKB activation (Vereecke et al., 2011). It has been shown that IRAK2 is an important factor and potential novel biomarkers for human antiviral innate immunity (Wang et al., 2015a). Also, MAP3K8 is a serine-threonine kinase and has a critical function in integrating host immune responses to complex pathogens (Mielke et al., 2009). Similarly, S. Srijata, et al. has suggested that AgNP attenuate the NF-κB pathway as indicated by the downregulation of NF-κB target genes CSF3 and subsequently inhibit MTB-induced proinflammatory responses (Sarkar et al., 2015). The downregulation of these genes has allowed the process of host defense to cease and give time to bacteria for progression.

● Up-regulated genes: We found that NLRC4 gene is associated with inflammasome signaling. Its upregulation means inflammasome activation (play an important role in host defense against MTB) not only leads to cytokine secretion but may also cause pyroptosis (Ting et al., 2008) (Bergsbaken et al., 2009). The NDUFS8 genes also play an important role in host immunity by increasing the expression level (Lohman et al., 2017). The gene FZD2 (frizzled class receptor 2) is Wnt pathway component (Wnt signaling), and Wnt signaling pathway plays a key role at different stages of TB development (Villaseñor et al., 2017). It has been reported that LRRK2 is involved in the IFN-γ response and host response to pathogens (Gardet et al., 2010), and LRRK2 also inhibits the immune response transcription factor NFAT1 (Liu et al., 2011). The genes DDX58 interact with IRGM and promote its K63-linked polyubiquitination, indicating that IRGM is positioned at a nexus of various innate immunity (Chauhan et al., 2016) while SAMD9L and EIF2AK2 play key roles in the innate immune responses to multiple stimuli. The gene EIF2AK2 is also involved in the regulation of signal transduction, apoptosis, cell proliferation, and differentiation (Lemos de Matos et al., 2013) (Onomoto et al., 2012) (Reineke and Lloyd, 2015). It is well known that type-I IFN induces potent cellular defense against viral infection mediated by up-regulation of ISGs like IFI44, IFI44L, and HERC5 (Jeon et al., 2010; Power et al., 2015). In our study, the expression of IFI44 and HERC5 were found to be up-regulated that means host system tried to resist bacterial infection.

We didn’t find any literature evidence that supports the role of four genes (NFE2L3, TST, MTHFD1, and CARS2) and three genes (ETV6, DOCK4, and PARD6G) which are down-regulated and up-regulated respectively in host immune response against the pathogen, but most of these genes are involved in various types of cancers like ETV6 is involved in prostate cancer as well as colorectal cancer susceptibility (Fararjeh and Liu, 2019; Wang et al., 2016); mutation in DOCK4 can cause ovarian, prostate, glioma, and colorectal cancers (Sundaravel et al., 2015), and PARD6G suppresses cell proliferation and is targeted by loss-of-function mutations in several cancers (Marques et al., 2016). Besides of these NFE2L3 (transcription factor) which involved in several cellular processes like inflammation, stress response, differentiation, and carcinogenesis (Kannan et al., 2015), MTHFD1 are key players in folate metabolism, which is essential for de novo purine synthesis, and several defects in this pathway have been associated with immunodeficiency (Keller et al., 2013), and CARS2 is associated with a severe progressive myoclonic epilepsy most resembling MERRF syndrome (Hallmann et al., 2014). The TST gene is a mitochondrial matrix enzyme (Pallini et al., 1991). It may play roles in cyanide detoxification, the formation of iron–sulfur proteins, and the modification of sulfur-containing enzymes.

Conclusion

In this study, we have compared gene-expression profiling of active TB patient and latent TB patients, as well as healthy controls. We have combined different data sets from multiple microarray experiment, and this combined effect size method for the gene expression analysis gave us more biologically consistent and conservative results. We have identified differentially expressed genes and their associated particular biological processes and pathways among healthy controls, active TB and LTBI. The gene interaction networks employed to obtain functional modules and also uncovered novel key regulators for active TB and LTBI which might play an important role in regulating the expression of multiple TB-related host factors. In conclusion, our study provides new insights into host genes and pathways important for TB infection. Experimental validation of these hypotheses would confirm the credibility of the key regulators and make easier the development of a cost-effective and sensitive molecular diagnostic platform for TB disease.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/gds

Author Contributions

RI conceived the study design instructed on data analysis. AA curated data and performed statistical and network analyses. NI, MM and AF implemented key statistical tools. AA, NI, MA, ST, AF, and NT drafted the manuscript. All the authors read, edited and approved of the manuscript.

Funding

This work was supported by the Department of Health Research (DHR), Ministry of Health and Family Welfare, Government of India. Under the scheme of “Young Scientist” fellowship (R.12014/06/2019-HR).

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.

Acknowledgments

The authors are grateful to the Centre for Interdisciplinary Research in Basic Sciences (CIRBSc), Jamia Millia Islamia-110025 for providing the research infrastructure. AA acknowledges the Department of Health Research (DHR) New Delhi, India for the award of “Young Scientist” fellowship (R.12014/06/2019-HR). NT is thankful to Council of Scientific and Industrial Research (CSIR), New Delhi, India, NI, MA, ST and AF are thankful to Indian Council of Medical Research (ICMR), New Delhi, India. RI duly acknowledges the Indian Council of Medical Research (ICMR), Council of Scientific and Industrial Research (CSIR) and Department of Health and Research (DHR).

Note: “This manuscript (Old version) has been released as a Pre-Print at [https://www.biorxiv.org/content/10.1101/414110v2]”.

Abbreviations

MTB, Mycobacterium tuberculosis; DEGS, Differentially expressed genes; KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, Gene Ontology; LCP, Local–community–paradigm; MDR, Multidrug-resistance TB; GBA, Guilt by association; MF, Molecular function; BP, Biological process; PC, Protein class; NC, Normal control; LI, Latent infection; ATB, Active TB; FKR, Fundamental key regulators; LAGO, Logically accelerated GO; GEO, Gene Expression Omnibus; GSE, GEO series; RMA, Robust multiarray average; LIMMA, Linear model for microarray analysis; CN, Common neighbors; HE, Hamiltonian energy; CPM, Constant Potts model; LCP-DP, LCP-decomposition-plot.

Supplementary Material

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

Supplementary Data 1 | Functional Classification of Differentially expressed according to (A) Molecular functions, (B) Biological processes and (C) Protein classes.

Supplementary Data 2 | The modular path of important hub genes from complete network to motif level with the structures of modules/sub-modules at various levels.

Supplementary Data 3 | (A) Community Finding (Network cluster Analysis) script in R. (B) LCP-correlation and LCP-decomposition-plot (LCP-DP approach to estimate the network compactness) method and MATLAB commonds.

References

Ahmad, S. (2011). Pathogenesis, immunology, and diagnosis of latent Mycobacterium tuberculosis Infection. Clin. Dev. Immunol. 2011, 1–17. doi: 10.1155/2011/814943

CrossRef Full Text | Google Scholar

Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 25, 25–29. doi: 10.1038/75556

PubMed Abstract | CrossRef Full Text | Google Scholar

Barabási, A.-L., Oltvai, Z. N. (2004). Network biology: understanding the cell’s functional organization. Nat. Rev. Genet. 5, 101–113. doi: 10.1038/nrg1272

PubMed Abstract | CrossRef Full Text | Google Scholar

Behler, F., Maus, R., Bohling, J., Knippenberg, S., Kirchhof, G., Nagata, M., et al. (2015). Macrophage-inducible C-type lectin mincle-expressing dendritic cells contribute to control of splenic Mycobacterium bovis BCG infection in mice. Infect. Immun. 83, 184–196. doi: 10.1128/IAI.02500-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergsbaken, T., Fink, S. L., Cookson, B. T. (2009). Pyroptosis: host cell death and inflammation. Nat. Rev. Microbiol. 7, 99–109. doi: 10.1038/nrmicro2070

PubMed Abstract | CrossRef Full Text | Google Scholar

Blankley, S., Graham, C. M., Turner, J., Berry, M. P. R., Bloom, C. I., Xu, Z., et al. (2016). The transcriptional signature of active tuberculosis reflects symptom status in extra-pulmonary and pulmonary tuberculosis. PloS One 11, e0162220. doi: 10.1371/journal.pone.0162220

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonacich, P. (1987). Power and centrality: a family of measures. Am. J. Sociol. 92, 1170–1182. doi: 10.1086/228631

CrossRef Full Text | Google Scholar

Bourigault, M.-L., Segueni, N., Rose, S., Court, N., Vacher, R., Vasseur, V., et al. (2013). Relative contribution of IL-1α, IL-1β and TNF to the host response to Mycobacterium tuberculosis and attenuated M. bovis BCG: IL-1α, IL-1β versus TNF in host response to M. tuberculosis.Immun. Inflamm. Dis. 1, 47–62. doi: 10.1002/iid3.9

PubMed Abstract | CrossRef Full Text | Google Scholar

Bouttier, M., Laperriere, D., Memari, B., Mangiapane, J., Fiore, A., Mitchell, E., et al. (2016). Alu repeats as transcriptional regulatory platforms in macrophage responses to M. tuberculosis infection. Nucleic Acids Res. 44, 10571–10587. doi: 10.1093/nar/gkw782

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyle, E. I., Weng, S., Gollub, J., Jin, H., Botstein, D., Cherry, J. M., et al. (2004). GO::TermFinder–open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinforma. Oxf. Engl. 20, 3710–3715. doi: 10.1093/bioinformatics/bth456

CrossRef Full Text | Google Scholar

Brandes, U. (2001). A faster algorithm for betweenness centrality*. J. Math. Sociol. 25, 163–177. doi: 10.1080/0022250X.2001.9990249

CrossRef Full Text | Google Scholar

Cai, Y., Yang, Q., Tang, Y., Zhang, M., Liu, H., Zhang, G., et al. (2014). Increased complement C1q level marks active disease in human tuberculosis. PloS One 9, e92340. doi: 10.1371/journal.pone.0092340

PubMed Abstract | CrossRef Full Text | Google Scholar

Cannistraci, C. V., Alanis-Lobato, G., Ravasi, T. (2013). From link-prediction in brain connectomes and protein interactomes to the local-community-paradigm in complex networks. Sci. Rep. 3, 1613. doi: 10.1038/srep01613

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavalcanti, Y. V. N., Brelaz, M. C. A., Lemoine Neves, J. K., de, A., Ferraz, J. C., Pereira, V. R. A. (2012). Role of TNF-alpha, IFN-gamma, and IL-10 in the development of pulmonary tuberculosis. Pulm. Med. 2012, 1–10. doi: 10.1155/2012/745483

CrossRef Full Text | Google Scholar

Chauhan, S., Mandell, M. A., Deretic, V. (2016). Mechanism of action of the tuberculosis and Crohn disease risk factor IRGM in autophagy. Autophagy 12, 429–431. doi: 10.1080/15548627.2015.1084457

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, C. Y., Gutierrez, N. M., Marzuki, M. B., Lu, X., Foreman, T. W., Paleja, B., et al. (2017). Host sirtuin 1 regulates mycobacterial immunopathogenesis and represents a therapeutic target against tuberculosis. Sci. Immunol. 2. doi: 10.1126/sciimmunol.aaj1789

CrossRef Full Text | Google Scholar

Clauset, A., Shalizi, C. R., Newman, M. E. J. (2009). Power-law distributions in empirical data. SIAM Rev. 51, 661–703. doi: 10.1137/070710111

CrossRef Full Text | Google Scholar

Fararjeh, A. S., Liu, Y.-N. (2019). ZBTB46, SPDEF, and ETV6: Novel potential biomarkers and therapeutic targets in castration-resistant prostate cancer. Int. J. Mol. Sci. 20, 2802. doi: 10.3390/ijms20112802

CrossRef Full Text | Google Scholar

Faridgohar, M., Nikoueinejad, H. (2017). New findings of toll-like receptors involved in Mycobacterium tuberculosis infection. Pathog. Glob. Health 111, 256–264. doi: 10.1080/20477724.2017.1351080

PubMed Abstract | CrossRef Full Text | Google Scholar

Gabor, C., Nepusz, T. (2006). The igraph software package for complex network research. Inter J Comp Syst 1695, 1–9.

Google Scholar

Ganguli, G. (2015). Mycobacterium tuberculosis acetyltransferase reduces the oxidative stress response through expression of peroxisomal membrane transporter protein. Open Forum Infect. Dis. 2. doi: 10.1093/ofid/ofv131.70

CrossRef Full Text | Google Scholar

Gardet, A., Benita, Y., Li, C., Sands, B. E., Ballester, I., Stevens, C., et al. (2010). LRRK2 is involved in the IFN-response and host response to pathogens. J. Immunol. 185, 5577–5585. doi: 10.4049/jimmunol.1000548

PubMed Abstract | CrossRef Full Text | Google Scholar

Garzon, R., Calin, G. A., Croce, C. M. (2009). MicroRNAs in cancer. Annu. Rev. Med. 60, 167–179. doi: 10.1146/annurev.med.59.053006.104707

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghosh, S., May, M. J., Kopp, E. B. (1998). NF-κB and Rel proteins: evolutionarily conserved mediators of immune responses. Annu. Rev. Immunol. 16, 225–260. doi: 10.1146/annurev.immunol.16.1.225

PubMed Abstract | CrossRef Full Text | Google Scholar

Guerra-Laso, J. M., Raposo-García, S., García-García, S., Diez-Tascón, C., Rivero-Lezcano, O. M. (2015). Microarray analysis of Mycobacterium tuberculosis-infected monocytes reveals IL26 as a new candidate gene for tuberculosis susceptibility. Immunology 144, 291–301. doi: 10.1111/imm.12371

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallmann, K., Zsurka, G., Moskau-Hartmann, S., Kirschner, J., Korinthenberg, R., Ruppert, A.-K., et al. (2014). A homozygous splice-site mutation in CARS2 is associated with progressive myoclonic epilepsy. Neurology 83, 2183–2187. doi: 10.1212/WNL.0000000000001055

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, D. W., Sherman, B. T., Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeon, Y. J., Yoo, H. M., Chung, C. H. (2010). ISG15 and immune diseases. Biochim. Biophys. Acta BBA - Mol. Basis Dis. 1802, 485–496. doi: 10.1016/j.bbadis.2010.02.006

CrossRef Full Text | Google Scholar

Jo, E.-K. (2013). Autophagy as an innate defense against mycobacteria. Pathog. Dis. 67, 108–118. doi: 10.1111/2049-632X.12023

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., Tanabe, M. (2016). KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, D457–D462. doi: 10.1093/nar/gkv1070

CrossRef Full Text | Google Scholar

Kannan, M. B., Dodard-Friedman, I., Blank, V. (2015). Stringent control of NFE2L3 (nuclear factor, erythroid 2-like 3; NRF3) protein degradation by FBW7 (F-box/WD repeat-containing protein 7) and glycogen synthase kinase 3 (GSK3). J. Biol. Chem. 290, 26292–26302. doi: 10.1074/jbc.M115.666446

PubMed Abstract | CrossRef Full Text | Google Scholar

Keller, M. D., Ganesh, J., Heltzer, M., Paessler, M., Bergqvist, A. G. C., Baluarte, H. J., et al. (2013). Severe combined immunodeficiency resulting from mutations in MTHFD1. Pediatrics 131, e629–e634. doi: 10.1542/peds.2012-0899

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J. (2018). Regulation of immune cell functions by metabolic reprogramming. J. Immunol. Res. 2018, 1–12. doi: 10.1155/2018/8605471

CrossRef Full Text | Google Scholar

Kondratieva, T., Azhikina, T., Nikonenko, B., Kaprelyants, A., Apt, A. (2014). Latent tuberculosis infection: what we know about its genetic control? Tuberc. Edinb. Scotl. 94, 462–468. doi: 10.1016/j.tube.2014.06.009

CrossRef Full Text | Google Scholar

Lavalett, L., Rodriguez, H., Ortega, H., Sadee, W., Schlesinger, L. S., Barrera, L. F. (2017). Alveolar macrophages from tuberculosis patients display an altered inflammatory gene expression profile. Tuberc. Edinb. Scotl. 107, 156–167. doi: 10.1016/j.tube.2017.08.012

CrossRef Full Text | Google Scholar

Lee, I., Blom, U. M., Wang, P. I., Shim, J. E., Marcotte, E. M. (2011). Prioritizing candidate disease genes by network-based boosting of genome-wide association data. Genome Res. 21, 1109–1121. doi: 10.1101/gr.118992.110

PubMed Abstract | CrossRef Full Text | Google Scholar

Lemos de Matos, A., Liu, J., McFadden, G., Esteves, P. J. (2013). Evolution and divergence of the mammalian SAMD9/SAMD9L gene family. BMC Evol. Biol. 13, 121. doi: 10.1186/1471-2148-13-121

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, K., Zhang, Y., Hu, S., Yu, Y., Yang, Q., Jin, D., et al. (2012). Increased levels of BAFF and APRIL related to human active pulmonary tuberculosis. PloS One 7, e38429. doi: 10.1371/journal.pone.0038429

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Lee, J., Krummey, S., Lu, W., Cai, H., Lenardo, M. J. (2011). The kinase LRRK2 is a regulator of the transcription factor NFAT that modulates the severity of inflammatory bowel disease. Nat. Immunol. 12, 1063–1070. doi: 10.1038/ni.2113

PubMed Abstract | CrossRef Full Text | Google Scholar

Lohman, B. K., Steinel, N. C., Weber, J. N., Bolnick, D. I. (2017). Gene expression contributes to the recent evolution of host resistance in a model host parasite system. Front. Immunol. 8, 107. doi: 10.3389/fimmu.2017.01071

PubMed Abstract | CrossRef Full Text | Google Scholar

Marques, E., Englund, J. I., Tervonen, T. A., Virkunen, E., Laakso, M., Myllynen, M., et al. (2016). Par6G suppresses cell proliferation and is targeted by loss-of-function mutations in multiple cancers. Oncogene 35, 1386–1398. doi: 10.1038/onc.2015.196

PubMed Abstract | CrossRef Full Text | Google Scholar

Mason, O., Verwoerd, M. (2007). Graph theory and networks in biology. IET Syst. Biol. 1, 89–119. doi: 10.1049/iet-syb:20060038

PubMed Abstract | CrossRef Full Text | Google Scholar

Mi, H., Muruganujan, A., Thomas, P. D. (2013). PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res. 41, D377–D386. doi: 10.1093/nar/gks1118

PubMed Abstract | CrossRef Full Text | Google Scholar

Mielke, L. A., Elkins, K. L., Wei, L., Starr, R., Tsichlis, P. N., O’Shea, J. J., et al. (2009). Tumor progression locus 2 (Map3k8) is critical for host defense against Listeria monocytogenes and IL-1 production. J. Immunol. 183, 7984–7993. doi: 10.4049/jimmunol.0901336

PubMed Abstract | CrossRef Full Text | Google Scholar

Newman, M. E. J. (2006a). Finding community structure in networks using the eigenvectors of matrices. Phys. Rev. E 74, 036104. doi: 10.1103/PhysRevE.74.036104

CrossRef Full Text | Google Scholar

Newman, M. E. J. (2006b). Modularity and community structure in networks. Proc. Natl. Acad. Sci. 103, 8577–8582. doi: 10.1073/pnas.0601602103

CrossRef Full Text | Google Scholar

Newman, M. E. J., Girvan, M. (2004). Finding and evaluating community structure in networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 69, 026113. doi: 10.1103/PhysRevE.69.026113

PubMed Abstract | CrossRef Full Text | Google Scholar

Okumura, R., Kurakawa, T., Nakano, T., Kayama, H., Kinoshita, M., Motooka, D., et al. (2016). Lypd8 promotes the segregation of flagellated microbiota and colonic epithelia. Nature 532, 117–121. doi: 10.1038/nature17406

PubMed Abstract | CrossRef Full Text | Google Scholar

Onomoto, K., Jogi, M., Yoo, J.-S., Narita, R., Morimoto, S., Takemura, A., et al. (2012). Critical role of an antiviral stress granule containing RIG-I and PKR in viral detection and innate immunity. PLoS One 7, e43031. doi: 10.1371/journal.pone.0043031

PubMed Abstract | CrossRef Full Text | Google Scholar

Pallini, R., Guazzi, G. C., Cannella, C., Cacace, M. G. (1991). Cloning and sequence analysis of the human liver rhodanese: comparison with the bovine and chicken enzymes. Biochem. Biophys. Res. Commun. 180, 887–893. doi: 10.1016/S0006-291X(05)81148-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Pastor-Satorras, R., Vázquez, A., Vespignani, A. (2001). Dynamical and correlation properties of the internet. Phys. Rev. Lett. 87, 258701. doi: 10.1103/PhysRevLett.87.258701

PubMed Abstract | CrossRef Full Text | Google Scholar

Power, D., Santoso, N., Dieringer, M., Yu, J., Huang, H., Simpson, S., et al. (2015). IFI44 suppresses HIV-1 LTR promoter activity and facilitates its latency. Virology 481, 142–150. doi: 10.1016/j.virol.2015.02.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Rangel-Santiago, J. F., Baay-Guzman, G. J., Duran-Padilla, M. A., Lopez-Bochm, K. A., Garcia-Romero, B. L., Hernandez-Cueto, D. D., et al. (2016). A novel role of yin-yang-1 in pulmonary tuberculosis through the regulation of the chemokine CCL4. Tuberculosis 96, 87–95. doi: 10.1016/j.tube.2015.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ravasz, E., Barabási, A.-L. (2003). Hierarchical organization in complex networks. Phys. Rev. E 67, 026112. doi: 10.1103/PhysRevE.67.026112

CrossRef Full Text | Google Scholar

Reineke, L. C., Lloyd, R. E. (2015). The stress granule protein G3BP1 recruits protein kinase R to promote multiple innate immune antiviral responses. J. Virol. 89, 2575–2589. doi: 10.1128/JVI.02791-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Reyes, N., Bettin, A., Reyes, I., Geliebter, J. (2015). Microarray analysis of the in vitro granulomatous response to Mycobacterium tuberculosis H37Ra. Colomb. Medica Cali Colomb. 46, 26–32. doi: 10.25100/cm.v46i1.1570

CrossRef Full Text | Google Scholar

Sagi, Z., Hieronymus, T. (2018). The impact of the epithelial–mesenchymal transition regulator hepatocyte growth factor receptor/met on skin immunity by modulating Langerhans cell migration. Front. Immunol. 9, 517. doi: 10.3389/fimmu.2018.00517

PubMed Abstract | CrossRef Full Text | Google Scholar

Salamon, H., Bruiners, N., Lakehal, K., Shi, L., Ravi, J., Yamaguchi, K. D., et al. (2014). Cutting edge: vitamin D regulates lipid metabolism in Mycobacterium tuberculosis infection. J. Immunol. Baltim. Md 1950 193, 30–34. doi: 10.4049/jimmunol.1400736

CrossRef Full Text | Google Scholar

Sarkar, S., Leo, B. F., Carranza, C., Chen, S., Rivas-Santiago, C., Porter, A. E., et al. (2015). Modulation of human macrophage responses to Mycobacterium tuberculosis by silver nanoparticles of different size and surface modification. PLOS One 10, e0143077. doi: 10.1371/journal.pone.0143077

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaaf, K., Smith, S. R., Duverger, A., Wagner, F., Wolschendorf, F., Westfall, A. O., et al. (2017). Mycobacterium tuberculosis exploits the PPM1A signaling pathway to block host macrophage apoptosis. Sci. Rep. 7, 42101. doi: 10.1038/srep42101

PubMed Abstract | CrossRef Full Text | Google Scholar

Schnappinger, D., Ehrt, S., Voskuil, M. I., Liu, Y., Mangan, J. A., Monahan, I. M., et al. (2003). Transcriptional adaptation of Mycobacterium tuberculosis within macrophages: insights into the phagosomal environment. J. Exp. Med. 198, 693–704. doi: 10.1084/jem.20030846

PubMed Abstract | CrossRef Full Text | Google Scholar

Scott-Browne, J. P., Shafiani, S., Tucker-Heard, G., Ishida-Tsubota, K., Fontenot, J. D., Rudensky, A. Y., et al. (2007). Expansion and function of Foxp3-expressing T regulatory cells during tuberculosis. J. Exp. Med. 204, 2159–2169. doi: 10.1084/jem.20062105

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Smyth, G. K. (2004). Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 3, Article3. doi: 10.2202/1544-6115.1027

CrossRef Full Text | Google Scholar

Soares-Silva, M., Diniz, F. F., Gomes, G. N., Bahia, D. (2016). The mitogen-activated protein kinase (MAPK) pathway: role in immune evasion by trypanosomatids. Front. Microbiol. 7, 183. doi: 10.3389/fmicb.2016.00183

PubMed Abstract | CrossRef Full Text | Google Scholar

Sundaravel, S., Duggan, R., Bhagat, T., Ebenezer, D. L., Liu, H., Yu, Y., et al. (2015). Reduced DOCK4 expression leads to erythroid dysplasia in myelodysplastic syndromes. Proc. Natl. Acad. Sci. 112, E6359–E6368. doi: 10.1073/pnas.1516394112

CrossRef Full Text | Google Scholar

Thuong, N. T. T., Dunstan, S. J., Chau, T. T. H., Thorsson, V., Simmons, C. P., Quyen, N. T. H., et al. (2008). Identification of tuberculosis susceptibility genes with human macrophage gene expression profiles. PLoS Pathog. 4, e1000229. doi: 10.1371/journal.ppat.1000229

PubMed Abstract | CrossRef Full Text | Google Scholar

Ting, J. P.-Y., Willingham, S. B., Bergstralh, D. T. (2008). NLRs at the intersection of cell death and immunity. Nat. Rev. Immunol. 8, 372–379. doi: 10.1038/nri2296

PubMed Abstract | CrossRef Full Text | Google Scholar

Traag, V. A., Krings, G., Van Dooren, P. (2013). Significant scales in community structure. Sci. Rep. 3, 2930. doi: 10.1038/srep02930

PubMed Abstract | CrossRef Full Text | Google Scholar

Traag, V. A., Van Dooren, P., Nesterov, Y. (2011). Narrow scope for resolution-limit-free community detection. Phys. Rev. E 84, 016114. doi: 10.1103/PhysRevE.84.016114

CrossRef Full Text | Google Scholar

van der Windt, G. J. W., Florquin, S., de Vos, A. F., van’t Veer, C., Queiroz, K. C. S., Liang, J., et al. (2010). CD44 deficiency is associated with increased bacterial clearance but enhanced lung inflammation during gram-negative pneumonia. Am. J. Pathol. 177, 2483–2494. doi: 10.2353/ajpath.2010.100562

PubMed Abstract | CrossRef Full Text | Google Scholar

Vereecke, L., Beyaert, R., van Loo, G. (2011). Genetic relationships between A20/TNFAIP3, chronic inflammation and autoimmune disease. Biochem. Soc. Trans. 39, 1086–1091. doi: 10.1042/BST0391086

PubMed Abstract | CrossRef Full Text | Google Scholar

Verway, M., Bouttier, M., Wang, T.-T., Carrier, M., Calderon, M., An, B.-S., et al. (2013). Vitamin D induces interleukin-1β expression: paracrine macrophage epithelial signaling controls M. tuberculosis infection. PLoS Pathog. 9, e1003407. doi: 10.1371/journal.ppat.1003407

PubMed Abstract | CrossRef Full Text | Google Scholar

Villaseñor, T., Madrid-Paulino, E., Maldonado-Bravo, R., Urbán-Aragón, A., Pérez-Martínez, L., Pedraza-Alva, G. (2017). Activation of the Wnt pathway by Mycobacterium tuberculosis: a Wnt–Wnt Situation. Front. Immunol. 8, 50. doi: 10.3389/fimmu.2017.00050

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, X.-K., Yuan, S.-L., Wang, Y.-C., Tao, H.-X., Jiang, W., Guan, Z.-Y., et al. (2016). Helicobacter pylori inhibits the cleavage of TRAF1 via a CagA-dependent mechanism. World J. Gastroenterol. 22, 10566. doi: 10.3748/wjg.v22.i48.10566

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., El Maadidi, S., Fischer, J., Grabski, E., Dickhöfer, S., Klimosch, S., et al. (2015a). A frequent hypofunctional IRAK2 variant is associated with reduced spontaneous hepatitis C virus clearance. Hepatology 62, 1375–1387. doi: 10.1002/hep.28105

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Wang, Z., Yao, Y., Wu, J., Tang, X., Gu, T., et al. (2015b). The fibroblast growth factor-2 arrests Mycobacterium avium sp. paratuberculosis growth and immunomodulates host response in macrophages. Tuberc. Edinb. Scotl. 95, 505–514. doi: 10.1016/j.tube.2015.04.006

CrossRef Full Text | Google Scholar

Wang, M., Gu, D., Du, M., Xu, Z., Zhang, S., Zhu, L., et al. (2016). Common genetic variation in ETV6 is associated with colorectal cancer susceptibility. Nat. Commun. 7, 11478. doi: 10.1038/ncomms11478

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaiss, D. M. W., Gause, W. C., Osborne, L. C., Artis, D. (2015). Emerging functions of amphiregulin in orchestrating immunity, inflammation, and tissue repair. Immunity 42, 216–226. doi: 10.1016/j.immuni.2015.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Huang, T., Wu, Y., Peng, W., Xie, H., Pan, M., et al. (2017). Inhibition of the PI3K-Akt-mTOR signaling pathway in T lymphocytes in patients with active tuberculosis. Int. J. Infect. Dis. 59, 110–117. doi: 10.1016/j.ijid.2017.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar