Skip to main content

Gene expression pattern analysis using dual-color RT-MLPA and integrative genome-wide association studies of eQTL for tuberculosis suscepitibility



When infected with Mycobacterium tuberculosis, only a small proportion of the population will develop active TB, and the role of host genetic factors in different TB infection status was not fully understood.


Forty-three patients with active tuberculosis and 49 with latent tuberculosis were enrolled in the prospective cohort. Expressing levels of 27 candidate mRNAs, which were previously demonstrated to differentially expressed in latent and active TB, were measured by dual color reverse transcription multiplex ligation dependent probe amplification assay (dcRT-MLPA). Using expression levels of these mRNAs as quantitative traits, associations between expression abundance and genome-wild single nucleotide polymorphisms (SNPs) were calculated. Finally, identified candidate SNPs were further assessed for their associations with TB infection status in a validation cohort with 313 Chinese Han cases.


We identified 9 differentially expressed mRNAs including il7r, il4, il8, tnfrsf1b, pgm5, ccl19, il2ra, marco and fpr1 in the prospective cohort. Through expression quantitative trait loci mapping, we screened out 8 SNPs associated with these mRNAs. Then, CG genotype of the SNP rs62292160 was finally verified to be significantly associated with higher transcription levels of IL4 in LTBI than in TB patients.


We reported that the SNP rs62292160 in Chinese Han population may link to higher expression of il4 in latent tuberculosis. Our findings provided a new genetic variation locus for further exploration of the mechanisms of TB and a possible target for TB genetic susceptibility studies, which might aid the clinical decision to precision treatment of TB.


Tuberculosis (TB) remains one of the most important threats to human health for thousands of years. According to a report of the World Health Organization (WHO) in 2017, about 10 million individuals newly infected with TB and 1.3 million died [1]. Most people infected with Mycobacterium tuberculosis (Mtb) can clear Mtb. However, approximately a quarter to three quarters of the population would keep a latent tuberculosis infection (LTBI) status and among them, about 5–10% persons will develop active tuberculosis (ATB) infections within their lifetimes [2].

Although living conditions, age, immune status could be predictors for ATB, with the development of molecule biology, the host genetic factors were thought to be substantial risk factors [3,4,5]. Kramnik et al. found that susceptibility to Mtb was a complex multigenic trait [6]. The study enrolling 202 co-twins suggested that the tuberculosis incidence among identical twins was higher than that among fraternal twins, which provided more evidence that the host genetic factors might affect the susceptibility to active tuberculosis infections [7].

In 2007, researchers firstly analyzed the gene-expression profiles of mRNA in whole blood samples from patients with ATB and LTBI [8]. Studies have shown that the expression levels of the host mRNAs vary at different post-infection status. It has been confirmed in several studies that expression levels of genes such as FcγRIB, FCGR1A, RAB33A, RAS and RIN3 could discriminate between ATB and LTBI patients [8]. Besides, some studies tried to analyze the association between human genetic variations and TB infection phenotype. In 2010, Thy et al. identified a single nucleotide polymorphism (SNP) rs4331426 located in a gene-poor region on chromosome 18q11.2 in African population by genome-wide association study (GWAS), which might make sense in infection of Mtb [9]. Then in 2012, they identified another SNP rs2057178 in Ghanaian, Gambian, Indonesian and Russian individuals using the same method [10].

Although GWAS can link SNPs with different disease phenotypes, some identified SNPs might locate on non-coding regions and thus provided limited information for the exploration of the function of the genetic variations. Mapping of expression quantitative trait loci (eQTL) can help to go a further step to associate differential mRNA expression with SNPs, which provides in-depth understanding in studying disease predisposition loci with potential function [11]. Previous studies have already identified SNPs associated with the discriminatory expression levels between Mtb infected and uninfected dendritic cells which were derived from peripheral blood mononuclear cells. Barreiro et al. performed eQTL mapping to identify DUSP14 and rs712039 as tuberculosis susceptibility genes [12]. Siddle et al. analyzed the expression levels of miRNA using eQTL and identified two SNPs, rs532751 and rs11159250 that could alter the expression of miR-326 and miR-1260[13]. These two studies analyzed the different expression profiles in healthy and infected cells, but they didn’t explore the differentiation of these profiles between ATB and LTBI. Besides, most of the previous studies were in Africa and Europe. Because the genetic backgrounds of different populations vary widely, conclusions from studies of other ethnic groups can not be applied in Chinese.

In our study, we directly analyzed the expression levels of candidate mRNAs in peripheral blood from ATB patients and LTBI patients using dcRT-MLPA, and conducted expression quantitative trait loci. Furthermore, we tried to find the associations between expression abundance and genome-wild single nucleotide polymorphisms (SNPs). Overall, this study aimed to comprehensively analyze the susceptible genes associated with tuberculosis infection and the possible mechanism of TB genetic susceptibility in Chinese Han population.


Ethics approval

This study was approved by ethical committees in Hua Shan Hospital (Shanghai, China) and we had already obtained written informed consent from all study subjects.

Prospective cohort enrollment

Patients and healthy donors

In this study, 92 adults aged 18 to 65 were enrolled in the prospective cohort.

ATB patients and LTBI patients were recruited at Zhuji People's Hospital of Zhejiang Province(Zhejiang, China) and 905th Hospital of PLA Navy (Shanghai, China). Individuals satisfying at least one of the following criteria were enrolled to ATB groups: (1) positive results of Mtb culture of respiratory specimens; (2) positive results of acid-fast bacilli smear microscopy; (3) positive histopathological findings together with TB clinical assessment: a history of close contact with ATB patients, positive clinical signs and chest radiography (e.g. cavities, hydrothorax, or changes on serial chest X-rays). All patients were diagnosed as active pulmonary tuberculosis and treated according to the guideline. LTBI patients were individuals who were closely household contact with the diagnosed ATB patients. They all had positive T-SPOT TB test findings and were without any clinical evidence of ATB.

Participants meeting any one of the following criteria were excluded: human immunodeficiency virus (HIV) testing positive, receiving immunosuppressive, chemical or biochemical therapy, treated for ATB for more than 8 weeks, or patients with abnormal liver function (AST, ALT ≥ 2ULN, TB ≥ 2.5 mg/dl).

Preparation of RNA

Two and half milliliters of whole blood were drawn from each individual in the prospective cohort with PAXgene blood RNA tube (Qiagen, German). According to the manufacturer’s instructions of kit, then we isolated and purified RNA using the PAXgene Blood RNA extraction kit(Qiagen, German). The concentration and the quality of RNA from whole blood were assessed by NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, USA).

Whole blood dual color reverse transcription multiplex ligation dependent probe amplification assay (dcRT-MLPA)

After reviewing and summarizing previous related literature and selected out 27 mRNAs as target mRNAs, all of which were reported to be differentially expressed between latent and active tuberculosis patients (il4, sec14l1, timp2, il7r, ltf, casp8, il8, mmp9, tnfrsf1a, tnfrsf1b, ccr7, ifng, pgm5, bcl2, tgfb1, ccl19, cd4, cd8a, il2ra, il6, foxp3, marco, tnfsf10c, blr1, fpr, fcgr1a and cxcl10) (Additional file 1: Table S1).

We designed immediately half probes for each target-specific sequence and also designed specific RT primers for the probe target sequence. All half probes and RT primers used in this study were described in details in Additional file 1: Table S2.

We conducted reverse transcription using the PrimeScriptTM RT reagent kit (Takara, Japan) and SLASA MLPA-RT RNA reagent kit (MRC-Holland, USA). The reaction conditions were described in Additional file 1: Table S3. Using HiDi formamide-containing 400 HD ROX size standard, the product of PCR amplification was diluted 10 times and then analyzed on 3730 capillary sequence in GeneScan mode (Applied Biosystems, USA). As for trace data analysis, we used the GeneMapper software package (Applied Biosystems, USA). We used B2M as reference gene and calculated the relative expression level of selected mRNAs to compare with B2M expressions.

DNA extraction and genome-wide SNP assay

For each sample, the genomic DNA was extracted from 200 µl of peripheral venous blood using PAXgene blood DNA tube (Qiagen, German) then extract genomic DNA according to manufacturer’s instructions. DNA content, purity, and integrity were determined using NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, USA). The genome-wide SNP assay was performed by the Infinium OmniZhongHua-8 BeadChip(Illumina, USA), which could detect 887,270 SNP genotype across the whole genome.

Mapping of expression quantitative trait loci (eQTLs)

EQTL mapping can associate the transcript abundance with specific genomic. Expression levels are viewed as quantitative traits, and associations between SNPs and expression abundance were calculated using a mixed-linear regression model:

$${\text{gene expression}} = \beta_{1} \times {\text{Genotype}} + \beta_{2} \times {\text{age}} + \beta_{3} \times {\text{gender}} + \beta_{4} \times {\text{kinship}} + \varepsilon$$

where patient age, gender and kinship information were considered as covariances in the above model. Gene expression phenotypes were mapped to particular genomic loci by combining studies of variation in gene expression patterns with genome-wide genotyping. The genome-wide association analysis was executed under Plink (v1.90b4.4 64-bit).

Statistical analysis

Microsoft excel spreadsheet version 2013 (Microsoft Corp., Redmond, WA, USA) and Microsoft access version 2010 (Microsoft Corp., Redmond, WA, USA) were used for data entry and management. As for the comparison between two groups, if the continuous variable conforms to the normal distribution, the t-test was used, while Mann Whitney U rank sum test was used for the non-normal distribution; the chi-square test or Fisher's exact test was used for the comparison of categorical variables. As for comparison between more than two clinical groups, a non-parametric Kruskal–Wallis test was performed. The statistical significance cut-off level was p < 0.05. All data was analyzed using SPSS 20.0 (IBM Corp., Armonk, NY, USA). Graphs was drawn with GraphPad 6.0(GraphPad Software, La Jolla, CA, USA). The regulatory network was drawn base on STRING database.

Validation of the identified SNPs

The validation cohort enrolled 204 ATB patients and 109 LTBI patients. Inclusion and exclusion criteria were the same as the prospective cohort. Patients were recruited from hospitals in Yangtze River Delta region (China) and pulmonary Hospital of Chonqing (Chongqing, China).

We genotyped the identified SNPs which were screened out by eQTLs mapping and further analyzed via multiplex SNaPshot® assay (Life Technologies, USA) in the validation cohort. Data was analyzed using SPSS 20.0 software (IBM Corp, USA) and histogram was performed to show the difference of selected SNPs between ATB and LTBI patients using GraphPad 6.0(La Jolla, USA).


We screened out 27 mRNAs for further analysis after reviewing the recent studies on biomarkers differently expressed in two TB infection status. Ninety-two individuals including 43 ATB patients and 49 LTBI patients were enrolled in the prospective cohort (Additional file 1: Figure S1). Using dcRT-MLPA, we identified 9 mRNAs with significant gene expression differences between the two groups. Furthermore, through eQTL mapping, we screened out SNPs significantly associated with the identified 9 mRNAs. And according to the function and interaction analysis, 8 SNPs were finally selected to be assessed in the validation cohort. In the validation cohort, we evaluated the identified SNPs and finally verified a new SNP rs62292160, which was shown to be associated with tuberculosis infection.

Identifying differentially expressed mRNAs between ATB and LTBI cohorts

The baseline information of 43 ATB and 49 LTBI patients in the prospective cohort was assessed. The mean age of ATB patients was 44 and the mean age of LTBI patients was 49.7. Significant gender difference was found between two groups (P = 0.012). The spot count of T-SPOT.TB and percentage of coughing and chest tightness or pain were significantly higher in the ATB groups than that in the LTBI group, while the percentage of enrolled subjects who admitted of BCG vaccination were significantly higher in the LTBI groups (P < 0.05). The overall baseline data was shown in Table 1.

Table 1 Baseline information of prospective and validation cohort

We analyzed the peripheral blood expression levels of 27 selected mRNAs using dcRT-MLPA and finally identified 9 differentially expressed (p < 0.05) mRNAs including il7r, il4, il8, tnfrsf1b, pgm5, ccl19, il2ra, marco, and fpr1. The peripheral blood gene expression levels of il7r, pgm5, ccl19, marco, and fpr1 were statistically up-regulated in ATB patients (p < 0.05), whereas the expression levels of il4, il8, tnfrsf1b, and il2ra were significantly up-regulated in LTBI groups (p < 0.05). There were no significant differences in the expression levels of the rest 18 genes (Table 2 and Additional file 1: Figure S2).

Table 2 Analysis of differential expression of whole blood mRNA in ATB and LTBI patients by dcRT-MLPA method

We further drew the receiver operating characteristic (ROC) curve and the area under the curves (AUCs) of identified 9 genes were 0.7635, 0.7633, 0.763, 0.7253, 0.6888, 0.6753, 0.6489, 0.643 and 0.6393, respectively (il7r, il4, il8, tnfrsf1b, pgm5, ccl19, il2ra, marco and fpr1) (Additional file 1: Figure S3).

Screening out specific SNPs associated with differentially expressed mRNAs using eQTLs analysis

We extracted DNA from peripheral blood mononuclear cells and obtained genome-wide SNPs information. Then we associated the SNPs with the 9 identified mRNAs. After eQTL mapping, we selected out the significant SNPs with the following strategy (Table 3). Firstly, 107 SNPs with p values less than E−05 were selected out (Fig. 1 and Additional file 1: Table S4). In the 107 SNPs, rs7203246, rs62292160, rs17170200, rs2189164 and rs3786247 were shown to have smallest p values (related genes were EEF2K, ADRA2C, TPK1, TACR3 and LIPG respectively). We then analyzed the function of 107 SNPs-related genes. Some were previously reported to participate in the process of human immune regulation or TB infection and these were also included for the future validation study. Among those, rs1667908-related gene was MBP and it translated protein which was found to express ubiquitously, especially in the immune system. And the protein might regulate proliferation and activation of T cells negatively [14]. Rs7203246-related gene was EEF2K, the gene coded the protein which was found to have effect in the differentiation of lung fibroblast through p38 mitogen activated protein kinases (MAPK) signal pathway [15]. Rs12712986 was related to gene TTC7A. TTC7A may have some effect in immunity and mutation of this gene was reported to be associated with inflammatory bowel disease [16]. What’s more, we analyzed the interaction relationship of SNPs-related genes and differentially expressed mRNAs previously identified (Fig. 2). Seven SNPs, rs4654400, rs7203246, rs62292160, rs2189164, rs3786247, rs1667908, and rs17170200 related genes were in the regulatory network. After comprehensive consideration of the p value, function of related genes and interaction among genes and identified mRNAs, we finally screened out 8 SNPs, rs7203246, rs62292160, rs17170200, rs2189164, rs3786247, rs12712986, rs4654400 and rs1667908 for the validation.

Table 3 Selected SNPs, related genes and corresponding selection strategies
Fig. 1

Manhattan plot. The expression levels of mRNAs were treated as quantitative traits and were associated with single nucleotide polymorphisms. The p value of green plots was less than E−05. The blue line represents E−05 and the red lines for E−07

Fig. 2

Regulatory network of significantly differentially expressed mRNAs and SNPs related genes. Regulatory network of mRNAs and SNPs related genes. The genes were surrounded by a red rectangle and the RNAs were surrounded by a blue rectangle. Seven genes and 9 mRNAs were in this network

Validation of the identified SNPs

Two-hundred and nine ATB and 104 LTBI patients were enrolled in the validation cohort. There were no significant differences in the mean age in ATB and LTBI groups. All patients were Han Chinese and were recruited from multi-centers. These patients were from Yangtze River Delta region and southwest area of China (Table 1). Eight selected SNPs were analyzed and successfully, we found the C/G genotype of rs62292160 showed significantly higher percentage in the LTBI patients and CC genotype of rs62292160 showed significantly higher in ATB patients. However, no variations of the rest 7 SNPs were found in the validation cohort (Fig. 3).

Fig. 3

Genotypes of 8 Sigle Nucleoside polymorphisms. We validated the identified SNPs and the histogram was performed to show the variations of selected SNPs in ATB and LTBI patients. The rs62292160 whose related gene was ADRA2C had significant differences in ATB and LTBI patient


Up to now, TB remains to be one of the top 10 causes of morbidity and mortality around the world [1]. The challenge from TB is due to many reasons, one of which is the existence of a large population infected with TB latently, and 10% of which will develop ATB during their lifetime [2]. The treatment for ATB and LTBI were different and high treatment coverage of LTBI is equally important to that of ATB. It can help to reduce the incidence of ATB and the burden of TB eventually. Currently, tuberculin skin test (TST) and interferon-gamma release assays (IGRAs) cannot discriminate between ATB and LTBI groups. And it is still unclear why some LTBI patients will develop ATB but others will not. The WHO guideline described at-risk groups who might progress to ATB, including patients with HIV, diabetes, chronic renal failure and so on [17]. Our study explored the characteristics of susceptible population at the molecular level and tried to explain the molecular mechanisms of susceptible population to TB in Chinese Han population.

Thye et al. have already identified rs4331426 and rs2057178 in TB patients by GWAS [9, 10]. The gene mutations of different ethnics were different and the same gene mutation might have opposite effects in different ethnics [18]. In recent years, researchers successively found the SNPs specific to Chinese Han population including rs1946518, rs5744247, rs1800796 and so on [18,19,20,21]. Although using GAWS or SNaPshot assay could associate gene variations with disease phenotypes, that could not explain the mechanisms underlying actions of SNPs. Chen et al. tried to explain the effects of SNPs in a Chinese cohort by correlating gene expression levels with previously identified rs4331426 and rs2057178, and finally screened out rs2057178 which might influence the expression levels of MAFB and SOCS2 in Chinese Han population [22]. Recently, researchers used eQTL mapping to identify the relationship between gene expression levels and genotypic variations in several diseases including chronic kidney disease, tumor and insomnia. And it helped to find plenty of significant SNPs [23,24,25]. In our study, we used eQTL mapping to analyze the association between SNPs, mRNA from peripheral blood and tuberculosis infection phenotypes in Chinese Han population.

In the prospective cohort, we selected 27 mRNAs after reviewing the recent studies on biomarkers which can discriminate between ATB and LTBI. The expression levels of these mRNAs were mainly analyzed in groups from Europe, Africa and America, but our study further tested the expression levels of theses mRNAs in Chinese Han population in different status of TB infection. Finally, we identified 9 mRNAs differentially expressed in ATB patients and LTBI patients, including il7r, il4, il8, tnfrsf1b, pgm5, ccl19, il2ra, marco and fpr1.

All these mRNAs have been found to play key roles in human immune responses or the progressing to ATB. Human il4 was the archetypal T-helper-2 (Th2) cytokine, serving as the critical element in the induction of Th2 response and blocked some of the effects of il2 [26]. il4 could subvert mycobacterial containment in human macrophages, likely via perturbations in Treg and Th1-linked pathways [27]. As for il7r, both Jenum et al. and Mihret et al. have found the different gene expression level of il7r in TB patients [28, 29]. Jenum et al. found that the tnfrsf1b separated ATB from LTBI patients in children [30]. Tnfr2 was an important TNF-α receptor and it was crucial to control pleural tuberculosis [31, 32]. Pgm5 made sense in the formation of carbohydrates and researchers have already found it expressed differentially in several diseases [33]. Ccl19 was involved in B cell-related pathways and played key role in immune response [34]. Marco was essential in the inflammatory responses to bacterial pathogens and several SNPs have been proved to be associated with susceptibility to infectious diseases [35]. And Fpr1 was detected to be important in neutrophils function-related pathways [36].

These selected mRNAs were previously proved to be critical elements in the development of infectious diseases including TB, and our cohort included a larger sample size of cases to further validate these mRNAs in Han population. Therefore, using these gene expression levels for further eQTLs mapping was more likely to find significant and meaningful genetic variations.

By the eQTL mapping analysis, we found out 8 SNPs in the prospective cohort. In the final validation step, we successfully identified a new SNP rs62292160, which showed a significant difference in ATB and LTBI patients in Chinese Han population. and that had not been reported the association with TB previously in any ethnics.

Rs62292160 was located at 4q12 and 62 Kb away from ADRA2C, which was totally a new locus. And the variation of rs62292160 was significantly associated with the higher expression level of IL4 in LTBI patients. The rs62292160 related to gene ADRA2C encodes adrenergic receptor alpha-2C (α2C-AR), which is one of the three subtypes of the alpha-2-adrenergic receptor(α2-AR). The protein α2-AR is one of the members of the G protein-coupled receptor superfamily and it has the general characteristics of G-protein coupled receptors. It can couple with G proteins of the Go/Gi-type, and inhibit the function of adenylate cyclase enzyme and the voltage-gated calcium channels. Furthermore, it plays a role in the activation of potassium channels and could also stimulate the mitogen-activated protein kinase (MAPK) [37]. The previous study suggested the decrease of alpha-adrenergic receptors in human lung tissue when patients suffered from destructive tuberculosis [38]. This suggested that the α2-ARs might be associated with TB. Bao et al. found that the α2-ARs participated in the down-regulation of cytokine production including IFN-γ and IL4 in the experiment using T lymphocytes from rats [39]. The cytokines played an essential role in combating Mtb. After the host is infected, Mtb interacts with the host to induce changes in cytokines release such as TNF-α, IFN-γ, IL8 and IL4 by affecting JAK/STAT, MAKP, calcium-dependent signaling pathways and so on [40,41,42]. The rs62292160-related gene ADRA2C might influence the expression of these cytokines through MAKP or calcium-dependent signaling pathways. Up to now, most studies focused on the adrenergic receptor beta(β2-AR) in lung tissue, more studies were needed to explore the function of α2-ARs in TB.

To our knowledge, this is the first application of dc-MLPA combined with eqtl mapping to analyze gene susceptibility of latent and active tuberculosis infection status. Using the technique, we have successfully found rs62292160 in the Chinese Han population, which might effect human immune responses in TB infection which leads to different infection status. On one hand, our result provides a new direction for further study of the mechanism of TB infection; and more importantly, our study raises the possible concept that receiving preventive anti-tuberculosis treatment might be more necessary for some LTBI patients with certain susceptible genes. However, such concept would still require a large amount of further clinical research and verification.


In this study, we validated the differential expression levels of 9 mRNAs in latent and active TB patients in the Chinese Han population. and we identified rs62292160, which might affect human immune responses in TB infection through regulating the expression of IL4. Not only did we provide a new variation locus for further exploration of mechanisms in TB infection, but also it might be helpful for the clinical researches. Because it holds the possibility to identify susceptible population who may progress to ATB, and guide preventive treatments for latent tuberculosis in the future.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.





World Health Organization


Mycobacterium tuberculosis


Latent tuberculosis infection


Active tuberculosis


Single nucleotide polymorphism


Genome-wide association study


Expression quantitative trait loci


Human immunodeficiency virus


Area under the curves


Mitogen activated protein kinases


Tuberculin skin test




Adrenergic receptor alpha-2C


Alpha-2-adrenergic receptor


Mitogen-activated protein kinase


  1. 1.

    Global tuberculosis report. 2018.

  2. 2.

    Lillebaek T, Dirksen A, Baess I, Strunge B, Thomsen VO, Andersen AB. Molecular evidence of endogenous reactivation of Mycobacterium tuberculosis after 33 years of latent infection. J Infect Dis. 2002;185(3):401–4.

    CAS  Article  Google Scholar 

  3. 3.

    Abel L, Fellay J, Haas DW, Schurr E, Srikrishna G, Urbanowski M, et al. Genetics of human susceptibility to active and latent tuberculosis: presentknowledge and future perspectives. Lancet Infect Dis. 2018;18(3):e64-75.

    CAS  Article  Google Scholar 

  4. 4.

    Basu Roy R, Whittaker E, Seddon JA, Kampmann B. Tuberculosis susceptibility and protection in children. Lancet Infect Dis. 2019;19(3):e96-108.

    Article  Google Scholar 

  5. 5.

    Furin J, Cox H, Pai M. Tuberculosis. Lancet (London, England). 2019;393(10181):1642–56.

    Article  Google Scholar 

  6. 6.

    Kramnik I, Demant P, Bloom BB. Susceptibility to tuberculosis as a complex genetic trait: analysis usingrecombinant congenic strains of mice. Novartis Found Symp. 1998;217(120–31):132–7.

    Google Scholar 

  7. 7.

    Comstock GW. Tuberculosis in twins: a re-analysis of the Prophit survey. Am Rev Respir Dis. 1978;117(4):621–4.

    CAS  PubMed  Google Scholar 

  8. 8.

    Mistry R, Cliff JM, Clayton CL, Beyers N, Mohamed YS, Wilson PA, et al. Gene-expression patterns in whole blood identify subjects at risk for recurrenttuberculosis. J Infect Dis. 2007;195(3):357–65.

    CAS  Article  Google Scholar 

  9. 9.

    Thye T, Vannberg FO, Wong SH, Owusu-Dabo E, Osei I, Gyapong J, et al. Genome-wide association analyses identifies a susceptibility locus for tuberculosis on chromosome 18q11.2. Nat Genet. 2010;42(9):739–41.

    CAS  Article  Google Scholar 

  10. 10.

    Thye T, Owusu-Dabo E, Vannberg FO, van Crevel R, Curtis J, Sahiratmadja E, et al. Common variants at 11p13 are associated with susceptibility to tuberculosis. Nat Genet. 2012;44(3):257–9.

    CAS  Article  Google Scholar 

  11. 11.

    Zou W, Aylor DL, Zeng Z. eQTL Viewer: visualizing how sequence variation affects genome-widetranscription. BMC Bioinf. 2007;8:7.

    Article  Google Scholar 

  12. 12.

    Barreiro LB, Tailleux L, Pai AA, Gicquel B, Marioni JC, Gilad Y. Deciphering the genetic architecture of variation in the immune response to Mycobacterium tuberculosis infection. Proc Natl Acad Sci. 2012;109:1204–9.

    CAS  Article  Google Scholar 

  13. 13.

    Siddle KJ, Deschamps M, Tailleux L, Nedelec Y, Pothlichet J, Lugo-Villarino G, et al. A genomic portrait of the genetic architecture and regulatory impact of microRNA expression in response to infection. Genome Res. 2014;24(5):850–9.

    CAS  Article  Google Scholar 

  14. 14.

    Xu J , Zhu D , Shan J , Fan Y. Changes of gene expression in T lymphocytes following GolliMBP gene RNA interference. Molecular Medicine Reports, 2016.

  15. 15.

    Wang Y, Huang G, Wang Z, Qin H, Mo B, Wang C. Elongation factor-2 kinase acts downstream of p38 MAPK to regulate proliferation, apoptosis and autophagy in human lung fibroblasts. Exp Cell Res. 2018;363(2):291–8.

    CAS  Article  Google Scholar 

  16. 16.

    Lawless D, Mistry A, Wood PM, Stahlschmidt J, Arumugakani G, Hull M, et al. Bialellic mutations in tetratricopeptide repeat domain 7A (TTC7A) cause common variable immunodeficiency-like phenotype with enteropathy. J Clin Immunol. 2017;37:617–22.

    CAS  Article  Google Scholar 

  17. 17.

    Global tuberculosis report. 2018.

  18. 18.

    Fornage M, Jacobs DR, Steffes MW, Gross MD, Bray MS, Schreiner PJ. Inverse effects of the PPARγ2 Pro12Ala polymorphism on measures of adiposity over 15 years in African Americans and whites. The CARDIA study. Metabolism. 2005;54(7):910–7.

    CAS  Article  Google Scholar 

  19. 19.

    Yu X, Zhao F, Zhang J, Pan X. Retracted article: IL-18 genetic polymorphisms may contribute to the pathogenesis of tuberculosis among Asians: a meta-analysis of case–control studies. Mol Biol Rep. 2014;41(9):6013–23.

    CAS  Article  Google Scholar 

  20. 20.

    Han M, Yue J, Lian Y, Zhao Y, Wang H, Liu L. Relationship between single nucleotide polymorphism of interleukin-18 and susceptibility to pulmonary tuberculosis in the Chinese Han population. Microbiol Immunol. 2011;55(6):388–93.

    CAS  Article  Google Scholar 

  21. 21.

    Zhang G, Zhou B, Wang W, Zhang M, Zhao Y, Wang Z, et al. A functional single-nucleotide polymorphism in the promoter of the gene encoding interleukin 6 is associated with susceptibility to tuberculosis. J Infect Dis. 2012;205(11):1697–704.

    CAS  Article  Google Scholar 

  22. 22.

    Chen C, Zhao Q, Hu Y, Shao Y, Li G, Zhu L, et al. A rare variant at 11p13 is associated with tuberculosis susceptibility in the Han Chinese population. Sci Rep-UK. 2016.

    Article  Google Scholar 

  23. 23.

    Jansen PR, Watanabe K, Stringer S, Skene N, Bryois J, Hammerschlag AR, et al. Genome-wide analysis of insomnia in 1,331,010 individuals identifies new riskloci and functional pathways. Nat Genet. 2019;51(3):394–403.

    CAS  Article  Google Scholar 

  24. 24.

    Qiu C, Huang S, Park J, Park Y, Ko Y, Seasock MJ, et al. Renal compartment-specific genetic variation analyses identify new pathways in chronic kidney disease. Nat Med. 2018;24(11):1721–31.

    CAS  Article  Google Scholar 

  25. 25.

    Zhang W, Bojorquez-Gomez A, Velez DO, Xu G, Sanchez KS, Shen JP, et al. A global transcriptional network connecting noncoding mutations to changes intumor gene expression. Nat Genet. 2018;50(4):613–20.

    CAS  Article  Google Scholar 

  26. 26.

    Darboe F, Adetifa JU, Reynolds J, Hossin S, Plebanski M, Netea MG, et al. Minimal sex-differential modulation of reactivity to pathogens and toll-likereceptor ligands following infant Bacillus Calmette-Guerin Russia Vaccination. Front Immunol. 2017;8:1092.

    Article  Google Scholar 

  27. 27.

    Malik A, Gupta M, Mani R, Bhatnagar R. Single-dose Ag85B-ESAT6-loaded poly(lactic-co-glycolic acid) nanoparticles conferprotective immunity against tuberculosis. Int J Nanomed. 2019;14:3129–43.

    CAS  Article  Google Scholar 

  28. 28.

    Jenum S, Dhanasekaran S, Lodha R, Mukherjee A, Kumar Saini D, Singh S, et al. Approaching a diagnostic point-of-care test for pediatric tuberculosis throughevaluation of immune biomarkers across the clinical disease spectrum. Sci Rep-UK. 2016a;6:18520.

    CAS  Article  Google Scholar 

  29. 29.

    Mihret A, Loxton AG, Bekele Y, Kaufmann SHE, Kidd M, Haks MC, et al. Combination of gene expression patterns in whole blood discriminate betweentuberculosis infection states. BMC Infect Dis. 2014;14:257.

    Article  Google Scholar 

  30. 30.

    Jenum S, Dhanasekaran S, Lodha R, Mukherjee A, Kumar Saini D, Singh S, et al. Approaching a diagnostic point-of-care test for pediatric tuberculosis through evaluation of immune biomarkers across the clinical disease spectrum. Sci Rep-UK. 2016b.

    Article  Google Scholar 

  31. 31.

    Chavez-Galan L, Vesin D, Uysal H, Blaser G, Benkhoucha M, Ryffel B, et al. Transmembrane tumor necrosis factor controls myeloid-derived suppressor cell activity via TNF receptor 2 and protects from excessive inflammation during BCG-induced pleurisy. Front Immunol. 2017;8:999.

    Article  Google Scholar 

  32. 32.

    Chavez-Galan L, Vesin D, Segueni N, Prasad P, Buser-Llinares R, Blaser G, et al. Tumor necrosis factor and its receptors are crucial to control Mycobacteriumbovis Bacillus Calmette-Guerin pleural infection in a murine model. Am J Pathol. 2016;186(9):2364–77.

    CAS  Article  Google Scholar 

  33. 33.

    Jiao Y, Li Y, Jiang P, Han W, Liu Y. PGM5: a novel diagnostic and prognostic biomarker for liver cancer. Peer J. 2019;7:e7070.

    Article  Google Scholar 

  34. 34.

    Tokunaga R, Naseem M, Lo JH, Battaglin F, Soni S, Puccini A, et al. B cell and B cell-related pathways for novel cancer treatments. Cancer Treat Rev. 2019;73:10–9.

    CAS  Article  Google Scholar 

  35. 35.

    Novakowski KE, Yap NVL, Yin C, Sakamoto K, Heit B, Golding GB, et al. Human-specific mutations and positively selected sites in MARCO Confer functional changes. Mol Biol Evol. 2018;35(2):440–50.

    CAS  Article  Google Scholar 

  36. 36.

    Schepetkin IA, Khlebnikov AI, Kirpotina LN, Quinn MT. Antagonism of human formyl peptide receptor 1 with natural compounds and their synthetic derivatives. Int Immunopharmacol. 2016;37:43–58.

    CAS  Article  Google Scholar 

  37. 37.

    Weber B, Schaper C, Wang Y, Scholz J, Bein B. Interaction of the ubiquitin carboxyl terminal esterase L1 with alpha(2)-adrenergic receptors inhibits agonist-mediated p44/42 MAP kinase activation. Cell Signal. 2009;21(10):1513–21.

    CAS  Article  Google Scholar 

  38. 38.

    Kondratenko TI, Kuzina NV, Severin ES, Kornilova ZK, Perel’man MI. Alpha-adrenergic receptors of the human lung parenchyma in destructive tuberculosis. Vopr Med Khim. 1991;37:69–70.

    PubMed  Google Scholar 

  39. 39.

    Bao J, Huang Y, Wang F, Peng Y, Qiu Y. Expression of alpha-AR subtypes in T lymphocytes and role of the alpha-ARs inmediating modulation of T cell function. NeuroImmunoModulation. 2007;14(6):344–53.

    CAS  Article  Google Scholar 

  40. 40.

    Rojas M, Olivier M, Garcia LF. Activation of JAK2/STAT1-alpha-dependent signaling events during Mycobacteriumtuberculosis-induced macrophage apoptosis. Cell Immunol. 2002;217(1–2):58–66.

    CAS  Article  Google Scholar 

  41. 41.

    Hasan Z, Shah BH, Mahmood A, Young DB, Hussain R. The effect of mycobacterial virulence and viability on MAP kinase signalling and TNF alpha production by human monocytes. Tuberculosis (Edinburgh, Scotland). 2003;83(5):299–309.

    CAS  Article  Google Scholar 

  42. 42.

    Rojas M, Garcia LF, Nigou J, Puzo G, Olivier M. Mannosylated lipoarabinomannan antagonizes Mycobacterium tuberculosis-induced macrophage apoptosis by altering Ca+2-dependent cell signaling. J Infect Dis. 2000;182(1):240–51.

    CAS  Article  Google Scholar 

Download references


Not applicable.


This work was supported by National Natural Science Foundation of China (82002141) Shanghai Youth Science and Technology Talents Sailing Project (20YF1404300).

Author information




J-WA: designed the study, collected the samples, and conducted the MLPA and eQTLs studies. Also, she is in charge of co-writing of this manuscript. HZ: statistics analysis and manuscript writing. ZZ: assisting in study design and samples collections. SW: conducting the eQTLs studies and statistics analysis. HH: in charge of clinical samples collections. SW: assisting the clinical samples collections. LS: MLPA experiment guidance. YG: Assisting the clinical samples collections. JW: MLPA experiment and statistics analysis guidance. QR: assisting the clinical samples collections in the validation cohort. FW: assisting the clinical samples collections in the validation cohort. NJ: guidance in the eQTLs analysis and manuscript writing. JC: guidance in the study design and MLPA experiment guidance. WZ: designed the study and in charge of the overall cohort management.

Corresponding authors

Correspondence to Ning Jiang or Jiazhen Chen or Wenhong Zhang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by ethical committees in Hua Shan Hospital (Shanghai, China) and we had already obtained written informed consent from all study subjects.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

The list of selected mRNAs that were reported to distinguish between ATB and LTBI in previous studies Table S2 The probes of 27 selected mRNAs for dcRT-MLPA. Table S3 The reaction steps in dcRT-MLPA. Table S4 The location and symbol of selected 107 SNPs with p value less than E-05. Figure S1. Flow chart. Figure S2. 9 mRNAs with statistically differential expression in ATB and LTBI patients. Figure S3. The receiver operating characteristic curve (ROC) of the selected mRNA in distinguishing ATB and LTBI patients.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ai, JW., Zhang, H., Zhou, Z. et al. Gene expression pattern analysis using dual-color RT-MLPA and integrative genome-wide association studies of eQTL for tuberculosis suscepitibility. Respir Res 22, 23 (2021).

Download citation


  • Latent tuberculosis
  • Active tuberculosis infection
  • Genetic susceptibility
  • Single nucleotide polymorphisms
  • Expression quantitative trait loci