Research article | Open | Open Peer Review | Published:
2′-5′-Oligoadenylate synthetase 1 polymorphisms are associated with tuberculosis: a case-control study
BMC Pulmonary Medicinevolume 18, Article number: 180 (2018)
2′-5′-Oligoadenylate synthetase 1 (OAS1) plays an important role in inflammatory immune reactions. OAS1 polymorphisms have been associated with increased susceptibility to various diseases. We investigated the association of polymorphisms in OAS1 with tuberculosis (TB).
A total of 1215 TB cases and 1114 healthy controls were enrolled from two independent studies. Genotyping was conducted using the improved multiplex ligase detection reaction (iMLDR) method. Associations between OAS1 polymorphisms (rs2240190, rs1131454, 10,774,671 and 11,066,453) and TB risk were established based on distributions of allelic frequencies using different genetic models.
Significant association was observed between rs10774671, rs1131454 and TB. In the initial study, the G allele of rs10774671 was a significantly protective factor against TB (P = 0.006) and the genotype of GG differed significantly between TB patients and controls under the codominant model (P = 0.008) after Bonferroni correction. In the validation study, we also observed that the rs10774671 G allele (P = 0.001) and GG genotype (P = 0.001) were associated with TB. In addition, we found that the rs1131454 G allele (P = 0.004) and GG genotype (P = 0.001) were protective against TB in the Chinese Han population.
We report novel associations of polymorphisms in OAS1 with TB in the Chinese Tibetan and Han populations. Similar studies in different populations and functional studies are warranted to confirm our results.
Tuberculosis (TB) is an infectious disease that constitutes a major global health problem. It is a major cause of morbidity and mortality globally, particularly in Asia and Africa, and ranks alongside human immunodeficiency virus (HIV) as a leading cause of death worldwide . In 2017, an estimated 1.7 billion individuals were newly infected with the causative agent of TB . However, only 5–15% of them will develop TB during their lifetime . The outcome of TB infection is affected by many factors, such as malnutrition, co-infection with other pathogens, exposure to microbes and previous vaccination . It was reported that host genetic factors play a crucial role in determining an individual’s susceptibility to TB . Recently, a number of genes potentially associated with TB susceptibility have been analysed in case-control studies.
Studies performed in mouse models and in patients indicate an important role of interferon (IFN) during TB infection [5, 6]. IFN-α, IFN-β and IFN-γ play significant immunomodulatory roles in the protection of the host against infections. As the main mediator of the type II immune response, IFN-γ is critical for controlling infections by intracellular pathogens, such as Mycobacteria . Specifically, IFN-γ plays a vital role in resistance to Mycobacterium tuberculosis (MTB) by activating macrophages, monocytes and Th1 cells . 2′-5′-oligoadenylate synthetases (OASs) are IFN-inducible enzymes that play vital roles in the innate immune response against viruses. The OAS family includes OAS1, OAS2, OAS3, and OAS-like (OASL) protein. The OAS1 gene encodes OAS1, which is an extensively characterized enzyme induced by IFNs . Binding of IFNs to their specific receptors results in induction of OAS1 gene expression . OAS1 is activated by the presence of double-stranded RNA and stimulates the oligomerisation of ATP into 2′,5′-linked oligoadenylates (2-5A) . 2-5A can bind to latent ribonuclease L (RNase L), which then dimerizes into the active form. This active form facilitates apoptosis , attenuates proliferation  and inhibits protein synthesis . It is possible that polymorphisms in OAS1 influence the expression of IFN-γ, the elimination of MTB, and ultimately affect the development of TB [14, 15].
Other OAS family members including OAS2, OAS3, and OASL also have been the focus of much research. OAS2 encodes the p69 and p71 IFN-induced isoforms . OAS2 play an important role in the control of bacterial infection through direct interaction with pattern-recognition receptors. It also has immunomodulatory function through immune cell receptor interaction . In terms of the ability to synthesize 2-5A, OAS3 shows higher activity than OAS1 . OAS3 is specialized for binding long dsRNA. The major function of OAS3 during infection is to produce 2-5A activators of RNaseL . In contrast to other OAS family members, OASL lacks any synthetic activity. OASL has the ability to regulate type I IFN responses when infected with pathogens. It has been suggested that OASL could act as a regulator in the control of antiviral innate immunity through IFN signalling .
Single-nucleotide polymorphisms (SNPs) can regulate OAS1 function at multiple levels, including expression, alternative splicing and enzyme activity. So far, among the OAS gene cluster, OAS1 polymorphisms have been the most studied and were reported to influence susceptibility to various diseases. For example, SNPs in OAS1 have been identified as candidates for susceptibility to viral infections, such as West Nile virus , hepatitis C , Chikungunya , dengue  and measles . Previous data also indicated that OAS1 polymorphisms and haplotypes were potential risk factors for autoimmune conditions, including type 1 diabetes  and multiple sclerosis . In addition, OAS1 polymorphisms were reported to be associated with respiratory infection . However, no previous studies have examined the association between OAS1 and TB.
Collectively, these observations prompted us to propose the hypothesis that OAS1 polymorphisms confer susceptibility to TB. The purpose of the study was to evaluate the prevalence of polymorphisms within the OAS1 gene in TB cases and healthy controls from the Chinese Tibetan and Han populations.
The initial study comprised 613 patients with TB and 602 healthy controls recruited from the People’s Hospital of the Aba Tibetan Autonomous Prefecture. In the independent validation study, we selected 571 TB cases and 543 healthy controls based on the same inclusion and exclusion criteria from the Chinese Han population. All subjects were unrelated ethnic Chinese.
The patients’ selection criteria were: 1) written consent and agreement to participate the study; 2) having a diagnosis of TB; 3) age ≥ 18 years. The diagnosis of TB was based on sputum smear tests, sputum culture, clinical symptoms and radiological and histological pathologic examinations. Patients with HIV, hepatitis B or C, diabetes, immune-mediated disorders or other lung diseases were excluded from our study. Controls were individuals without history of TB and active TB infection that had undergone routine medical examination.
Informed consents were obtained from all individual participants included in the study. The study was approved by the ethical committee of the West China Hospital of Sichuan University.
SNP selection and genotyping
Tag SNPs in the study were downloaded from the Chinese Han in Beijing database of HapMap (http://hapmap.ncbi.nlm.nih.gov/index.html.en, HapMap Data Rel 27 Phase II + III, on NCBI B36 assembly, dbSNP b126) in a region stretching from 1000 base pairs upstream and 1000 base pairs downstream of OAS1; other SNPs of OAS1 were also selected from previous studies [27, 28]. Tag SNPs were filtered according to the following criteria: (i) Hardy-Weinberg equilibrium (HWE) test P ≥ 0.05; (ii) r2 of pairwise linkage disequilibrium (LD) ≤ 0.8 and (iii) minor allele frequency (MAF) ≥ 0.05.
A sample of 5 ml of peripheral blood was obtained from each participant and placed in a test tube containing EDTA. We then isolated genomic DNA from the collected whole blood samples using a genomic DNA purification kit, in accordance with the manufacturer’s instructions (Axygen Scientific Inc., Union City, CA, USA). Genotyping was conducted using PCR and the improved multiple ligase detection reaction (iMLDR). Five percent of samples were re-genotyped to verify the initial results.
HWE and the differences in sex between the case and control were evaluated by the Chi-squared test. Unpaired t-tests were used to evaluate the significance of differences in age between the case and control groups. Genotype and allele frequencies of SNPs in the TB group vs. those in the control group were calculated by logistic regression. As the genetic determinants of TB were age- and gender-dependent [29, 30], logistic regression was used to adjust for gender and age. Odds ratio (OR) and 95% confidence interval (95% CI) were examined to estimate the magnitude of the risk. LD and haplotype analyses were conducted using the SHEsis online software platform (http://analysis.bio-x.cn). A 2-sided P value of less than 0.05 was considered statistically significant and P values were conservatively Bonferroni-corrected for multiple comparisons. Power calculation was performed (α = 0.05) by Quanto software (version 1.2.4, May 2009). The power of the two cohorts to detect relative risks (RR) of 2.0, 3.0 and 4.0 using the four SNPs studied under an allele genetic model were estimated by the Quanto software and are listed in Table S1. The two cohorts have reasonable power (> 80%) to detect genetic factors with RR 2.0 and above under the allele model in all four SNPs (Table S1). Gene-gene interaction was detected using the multifactor dimensionality reduction (MDR) software and MDR permutation testing module (http://www.epistasis.org and http://sourceforge.net/projects/mdr/files/mdrpt). We performed an eQTL analysis by using the GTEx Portal database (https://gtexportal.org/home/). All analyses were conducted using the Statistical Package for the Social Sciences (SPSS; SPSS Inc., Chicago, IL, USA), version 19.0.
Genotype data of rs10774671 in the Han Chinese in Beijing (CHB) were obtained from the HapMap project and whole-genome Expression data were obtained from the Gene Expression Omnibus of PubMed. rs10774671 genotype data were downloaded from the 1000 Genomes Browser (https://www.ncbi.nlm.nih.gov/variation/tools/1000genomes/). DNA samples were collected from participants living in the residential community at Beijing Normal University. As a result, a total of 40 individuals provided genotypic data of rs10774671. We then used the mRNA data obtained from lymphoblastoid cell lines (extracted from Gene Expression Omnibus of Pubmed, accession number GSE6536) to analyze the associations between rs10774671 and OAS1 mRNA expression levels .
Demographic characteristics of the study groups
In the initial study, the number of cases was 613 (male = 392; female = 221) and the controls was 602 (male = 392; female = 210). The validation cohort comprised 584 cases (male = 299; female = 285) and 543 controls (male = 266; female = 277). The mean age was 34.53 ± 14.54 years for the case group and 34.63 ± 13.85 years for the control group in the initial study. There were no significant differences between the case and control groups regarding age and sex (P > 0.05) in this cohort. However, the mean age was significantly different between cases (mean age = 36.50 ± 15.49) and controls (mean age = 39.35 ± 15.13) in the validation cohort (P = 0.002).
Analysis of HWE and association of OAS1 SNPs with TB
A total of four SNPs (rs2240190, rs1131454, 10,774,671 and 11,066,453) were selected for analysis. Two of these four were tag SNPs (rs1131454 and rs10774671). All four SNPs did not deviate from HWE. The MAFs of SNPs were more than 0.05. Table 1 provides basic information on the four SNPs in our study.
The allele frequencies in the two groups are shown in Table 2. The genotype frequencies of the OAS1 polymorphisms under different genetic models are summarized in Table 3. In the initial cohort, a significant association was only detected in the allele frequency of rs10774671 and the G allele acted as a protective factor against TB (P = 0.006, OR = 0.74, 95% CI: 0.59–0.92). Compared with the AA genotype, the frequency of the GG genotype of rs10774671 polymorphism in the case group significantly differed from the controls under the codominant model (GG vs. AA: P = 0.008, OR = 0.41, 95% CI: 0.21–0.80), suggesting that rs10774671 decreased the risk of TB.
To verify the results of the Tibetan population, we analyzed an independent population of Chinese Han individuals. We found that rs10774671 was also significantly associated with TB in the Chinese Han population. Compared with allele A, the G allele was more common in the control group (P = 0.001, OR = 0.72, 95% CI: 0.59–0.87). We also found that the GG genotype was a protective factor against TB under codominant (GG VS. AA: P = 0.001, OR = 0.40, 95% CI: 0.23–0.69) and recessive (GG VS. AA+GA: P = 0.005, OR = 0.46, 95% CI: 0.27–0.79) models (Table 4). In addition, rs1131454 was associated with TB susceptibility in the validation study. Compared with the control group, the TB group had lower proportions of the G allele (P = 0.004, OR = 0.78, 95% CI: 0.66–0.93) and GG genotype (codominant model: GG vs. AA, P = 0.001, OR = 0.52, 95% CI: 0.36–0.76; recessive model: GG vs. AA+GA, P = 0.002, OR = 0.59, 95% CI: 0.42–0.82).
Linkage disequilibrium and haplotype analyses
Information on LD for these four SNPs in the Chinese Tibetan population is shown in Fig. 1. LD among these four SNPs in the OAS1 gene was evaluated by the LD test using r2, which confirmed that these SNPs were not in LD (r2 < 0.5). No haplotypes were associated with TB in the two studies (Table 4).
rs2240190, rs1131454 and rs10774671 formed the best gene-gene interaction model with a testing balanced accuracy of 57.61% and cross-validation consistency of 10/10. It was confirmed statistically significant using 1000-fold permutation testing (P = 0.001).
eQTL and gene expression analysis for rs10774671
eQTL analysis using the GTEx Portal database showed significant eQTL effects of rs10774671 (P < 0.001) on OAS1 expression in whole blood from 369 individuals. The 40 CHB were divided into three groups according to rs10774671genotypic data. In brief, 15 CHB had the AA genotype, 20 had the GA genotype and 5 had the GG genotype. As shown in Fig. 2, GEO2R analysis revealed that the rs10774671 polymorphism was associated with OAS1 expression levels after adjustment for multiple comparisons (P < 0.001).
In China, the ethnic Tibetan population is more likely to be infected with TB than the Han population [32, 33]. It is possible that host genetic difference is one of the bases of TB susceptibility. In this study, four SNPs in OAS1 (rs2240190, rs1131454, 10,774,671 and 11,066,453) were investigated in Chinese Han and Chinese Tibetan patients with TB. Our data showed, for the first time, that the G allele and the GG genotype of rs10774671 in OAS1 were associated with TB in the two populations.
OAS1 is an enzyme that plays an important role in innate antiviral defence . This enzyme is activated by the presence of double-stranded or single-stranded RNA with a secondary structure and catalyses the oligomerisation of ATP into 2-5A . These processes activate latent RNase L, which degrades viral and cellular RNA and blocks protein synthesis. In accordance with this function, human OAS1 mainly affects susceptibility to viral infections [35, 36]. Given the importance of the OAS1 protein, SNPs in the OAS1 gene may affect susceptibility to infectious diseases. In the present study, we demonstrated that polymorphisms in OAS1 were associated with TB. It is unclear why the heightened antiviral enzyme activity in individuals was associated with TB, but this association may be explained as follows. One explanation is that the major role of the OAS proteins is as immune regulators in innate immunity, although previous studies indicated that they also play important roles in other cellular functions. The role of OAS1 in immune-related diseases has been shown by a number of association studies [25, 26]. It is well known that the immune response against TB plays a critical role in the progression of MTB infection. Another explanation is that the OAS1 gene affects susceptibility to TB through the type II IFN pathway [14, 15]. Type II IFN (IFN-γ) is critical for host defence against certain bacterial and parasitic pathogens . IFN-γ can activate infected host macrophages to inhibit the replication of MTB directly . In addition, OAS gene expression was regulated by both type I and type II IFNs . A previous study suggested that IFN-γ could increase the levels of OAS mRNA . Furthermore, OAS1 gene up-regulation was observed in several gene expression signatures that differentiated active TB from latent TB infection [41, 42]. It was also established that there was a prominent correlation between OAS expression and TB . Moreover, data from Noguchi et al. revealed the potential role of OAS1 polymorphisms in respiratory infection . Therefore, we assume that the polymorphism in the OAS gene may be associated with the risk of TB infection through the type II IFN pathway.
As a functional OAS1 polymorphic marker, rs10774671 was reported to be associated with disease/viral infection [36, 43]. rs10774671 is located at the last nucleotide of intron 6 in OAS1, which acts as a splice-acceptor site for exon 7. The G allele was demonstrated to allow splicing to occur leading to the production of a p46 form with high enzymatic OAS activity . The rs10774671 A allele usually prevents splicing at this site and thus splicing happens further downstream, resulting in the p48 and p52 forms associated with low OAS enzymatic activity. rs10774671 could also control other splice variants of OAS including p42 and p44 . In addition, the HapMap genotypic data (Fig. 2) of rs10774671 suggested that individuals with GG genotype had lower expression levels than those with AA and GA genotypes. However, rs10774671 GG genotype was reported to be associated with the highest enzyme activity of OAS1 in unstimulated lymphocytes . The inconsistent results between mRNA expression and enzyme activity may attribute to the reason that gene expressions at the transcriptional and translational levels might be different. Our data were partly consistent with these observations, in that the G allele and GG genotype had significant protective effects against TB, whereas the AG genotype was not associated with TB. Furthermore, this result was validated in the Chinese Han population. Considering the aforementioned functional effects of rs10774671 on phenotypes, we speculate that rs10774671 could alter OAS1 expression levels precluding the type I IFN pathway resulting in a protective response.
In order to detect whether two or more SNPs interact either directly or indirectly to change TB risk separate from their independent effects, we conducted a gene-gene interaction analysis. We did not find any significant differences in single SNP analysis between cases and controls groups for rs2240190, rs1131454 and 11,066,453. However, rs2240190, rs1131454 and rs10774671 formed the best gene-gene interaction model. This result was consistent with the phenomenon that most diseases including TB are caused by multiple factors and gene-gene or gene-environment interactions are characteristic of the genetic factors. In addition, we also conducted a haplotype analysis which was based on the association between a polymorphism and the ancestral haplotype . However, inconsistent with the gene-gene interaction result, there was no significant association between haplotypes and TB. Since the haplotype analysis approach was used for analyzing haplotype effects at numerous closely linked loci, and the gene-gene interaction analysis could be used for analyzing the unlinked loci [46, 47]. We speculate the discrepant results are likely attributable to different methods.
rs2240190 is located in intron 1 and it may have an alternative splicing regulatory effect, based on the Functional Single Nucleotide Polymorphism database . To the best of our knowledge, no previous studies investigated the association between rs2240190 and disease. Only one study suggested that this genetic polymorphism was not in LD with rs10774671 . rs1131454 (formerly rs3741981), in the evolutionarily conserved exon 3, is nonsynonymous (the G → A transition results in a G → S amino acid change). This polymorphic locus is close to the dsRNA binding domain of all OAS1 isoforms. In general, the AA, GA and GG genotypes of rs1131454 lead to low, intermediate and high OAS1 enzyme activity, respectively . This polymorphic locus was demonstrated to be functional in the progression of the severe acute respiratory syndrome . In our study, we demonstrated that the G allele and GG genotype were associated with TB in the Han population.
Genome-wide association studies with longitudinal data found that rs11066453 was a significant disease susceptibility locus, as confirmed by the Health Examination cohort . Furthermore, the rs11066453 polymorphism was strongly related to OAS1 mRNA expression levels in Epstein–Barr virus infection , and the monoallelically expressed OAS1 gene was also significantly associated with the gamma-glutamyltransferase level . This strong association of the OAS1 gene that contains the rs11066453 polymorphism in the development and progression of human diseases may be due to it affecting OAS1 expression. However, inconsistent with these findings, there was no significant association with TB susceptibility detected in the present study.
There are several limitations in our study. First, lack of data regarding the association between other genes in the OAS cluster and TB limited our understanding of genetic mechanisms regulating the pathogenesis of TB. Second, our study lacked clinical information for the Tibetans, which limited the analysis of the association with clinical characteristics. Finally, functional validation of the included SNPs was not carried out. As a result, the true causal allele underlying the genetic association result is still unknown.
In conclusion, we were the first to investigate the association between OAS1 polymorphisms and TB in the Chinese Tibetan population and validated the results in an independent cohort. We found that rs10774671, located in the OAS1 gene region is associated with susceptibility to TB in Tibetans and may act as a protective factor against TB, which was validated in the Chinese Han population. In addition, we also found that rs1131454 was a protective factor against TB in the Chinese Han population. rs2240190, rs1131454 and rs10774671 formed the best gene-gene interaction model in MDR analysis and may be considered as one of the multiple contributors to the progress of this complicated disease. However, in order to get more detailed mechanistic insight into the pathway of how OAS1 may act on TB, there is a need for further research in this field. Similar studies in different populations, gene-environment interactions analysis and functional studies are warranted to confirm and reinforce our results.
Human immunodeficiency virus
Improved multiple ligase detection reaction
Pairwise linkage disequilibrium
Minor allele frequency
Multifactor dimensionality reduction
2′-5′-oligoadenylate synthetase 1
- RNase L:
World Health Organization (2018). Global tuberculosis report 2018. http://www.who.int/tb/publications/global_report/en/ (accessed Sep 26, 2018).
Comstock GW, et al. The prognosis of a positive tuberculin reaction in childhood and adolescence. Am J Epidemiol. 1974;99(2):131–8.
van de Vosse E, et al. Human genetics of intracellular infectious diseases: molecular and cellular immunity against mycobacteria and salmonellae. Lancet Infect Dis. 2004;4(12):739–49.
Newport MJ, et al. A mutation in the interferon-gamma-receptor gene and susceptibility to mycobacterial infection. N Engl J Med. 1996;335(26):1941–9.
Manca C, et al. Virulence of a Mycobacterium tuberculosis clinical isolate in mice is determined by failure to induce Th1 type immunity and is associated with induction of IFN-alpha /beta. Proc Natl Acad Sci U S A. 2001;98(10):5752–7.
de Paus RA, et al. Inhibition of the type I immune responses of human monocytes by IFN-alpha and IFN-beta. Cytokine. 2013;61(2):645–55.
van de Vosse E, et al. Genetic deficiencies of innate immune signalling in human infectious disease. Lancet Infect Dis. 2009;9(11):688–98.
Flynn JL, et al. An essential role for interferon gamma in resistance to Mycobacterium tuberculosis infection. J Exp Med. 1993;178(6):2249–54.
Hovnanian A, et al. The human 2′,5′-oligoadenylate synthetase locus is composed of three distinct genes clustered on chromosome 12q24.2 encoding the 100-, 69-, and 40-kDa forms. Genomics. 1998;52(3):267–77.
Sadler AJ, Williams BR. Interferon-inducible antiviral effectors. Nat Rev Immunol. 2008;8(7):559–68.
Hovanessian AG, Justesen J. The human 2′-5’oligoadenylate synthetase family: unique interferon-inducible enzymes catalyzing 2′-5′ instead of 3′-5′ phosphodiester bond formation. Biochimie. 2007;89(6–7):779–88.
Domingo-Gil E, Esteban M. Role of mitochondria in apoptosis induced by the 2-5A system and mechanisms involved. Apoptosis. 2006;11(5):725–38.
Yaffe A, et al. Inhibition of 2-5A synthetase expression by antisense RNA interferes with interferon-mediated antiviral and antiproliferative effects and induces anchorage-independent cell growth. Cell Growth Differ. 1996;7(8):969–78.
Leisching G, et al. OAS1, 2, and 3: significance during active tuberculosis? J Infect Dis. 2018;217(10):1517–21.
Leisching G, et al. The association of OASL and type I Interferons in the pathogenesis and survival of intracellular replicating bacterial species. Front Cell Infect Microbiol. 2017;7:196.
Justesen J, et al. Gene structure and function of the 2′-5′-oligoadenylate synthetase family. Cell Mol Life Sci. 2000;57(11):1593–612.
Ibsen MS, et al. The 2′-5′-oligoadenylate synthetase 3 enzyme potently synthesizes the 2′-5′-oligoadenylates required for RNase L activation. J Virol. 2014;88(24):14222–31.
Li Y, et al. Activation of RNase L is dependent on OAS3 expression during infection with diverse human viruses. Proc Natl Acad Sci U S A. 2016;113(8):2241–6.
Dunn GP, et al. The immunobiology of cancer immunosurveillance and immunoediting. Immunity. 2004;21(2):137–48.
Lim JK, et al. Genetic variation in OAS1 is a risk factor for initial infection with West Nile virus in man. PLoS Pathog. 2009;5(2):e1000321.
Zhao Y, et al. Evaluate the relationship between polymorphisms of OAS1 gene and susceptibility to chronic hepatitis C with high resolution melting analysis. Clin Exp Med. 2013;13(3):171–6.
Chaaithanya IK, et al. Association of Oligoadenylate Synthetase Gene Cluster and DC-SIGN (CD209) gene polymorphisms with clinical symptoms in Chikungunya virus infection. DNA Cell Biol. 2016;35(1):44–50.
Thamizhmani R, Vijayachari P. Association of dengue virus infection susceptibility with polymorphisms of 2′-5′-oligoadenylate synthetase genes: a case-control study. Braz J Infect Dis. 2014;18(5):548–50.
Haralambieva IH, et al. Genetic polymorphisms in host antiviral genes: associations with humoral and cellular immunity to measles vaccine. Vaccine. 2011;29(48):8988–97.
Tessier MC, et al. Type 1 diabetes and the OAS gene cluster: association with splicing polymorphism or haplotype? J Med Genet. 2006;43(2):129–32.
Fedetz M, et al. OAS1 gene haplotype confers susceptibility to multiple sclerosis. Tissue Antigens. 2006;68(5):446–9.
Noguchi S, et al. Differential effects of a common splice site polymorphism on the generation of OAS1 variants in human bronchial epithelial cells. Hum Immunol. 2013;74(3):395–401.
Park YM, et al. Genome-wide detection of allelic gene expression in hepatocellular carcinoma cells using a human exome SNP chip. Gene. 2014;551(2):236–42.
Hall NB, et al. Polymorphisms in TICAM2 and IL1B are associated with TB. Genes Immun. 2015;16(2):127–33.
Sabri A, et al. Association study of genes controlling IL-12-dependent IFN-gamma immunity: STAT4 alleles increase risk of pulmonary tuberculosis in Morocco. J Infect Dis. 2014;210(4):611–8.
Chen C, et al. A rare variant at 11p13 is associated with tuberculosis susceptibility in the Han Chinese population. Sci Rep. 2016;6:24016.
He S, et al. Association of IL4, IL6, and IL10 polymorphisms with pulmonary tuberculosis in a Tibetan Chinese population. Oncotarget. 2018;9(23):16418–26.
Wang L, et al. The fifth national tuberculosis epidemiological survey in 2010. Chin J Antituber. 2012;34(8):485–508.
Schneider WM, et al. Interferon-stimulated genes: a complex web of host defenses. Annu Rev Immunol. 2014;32:513–45.
Tsutsui-Takeuchi M, et al. Roles of retinoic acid-inducible gene-I-like receptors (RLRs), toll-like receptor (TLR) 3 and 2′-5′ oligoadenylate synthetase as viral recognition receptors on human mast cells in response to viral infection. Immunol Res. 2015;61(3):240–9.
Liu X, et al. A functional variant in the OAS1 gene is associated with Sjogren’s syndrome complicated with HBV infection. Sci Rep. 2017;7(1):17571.
Boehm U, et al. Cellular responses to interferon-gamma. Annu Rev Immunol. 1997;15:749–95.
MacMicking JD, et al. Immune control of tuberculosis by IFN-gamma-inducible LRG-47. Science. 2003;302(5645):654–9.
Zhou A, et al. Interferon action and apoptosis are defective in mice devoid of 2′,5′-oligoadenylate-dependent RNase L. EMBO J. 1997;16(21):6355–63.
Corrias MV, et al. Induction of 2.5 OAS gene expression and activity is not sufficient for IFN-gamma-induced neuroblastoma cell differentiation. Int J Cancer. 1995;62(2):223–9.
Berry MP, et al. An interferon-inducible neutrophil-driven blood transcriptional signature in human tuberculosis. Nature. 2010;466(7309):973–7.
Ottenhoff TH, et al. Genome-wide expression profiling identifies type 1 interferon response pathways in active tuberculosis. PLoS One. 2012;7(9):e45839.
Bader El Din NG, et al. Impact of OAS1 exon 7 rs10774671 genetic variation on liver fibrosis progression in Egyptian HCV genotype 4 patients. Viral Immunol. 2015;28(9):509–16.
Bonnevie-Nielsen V, et al. Variation in antiviral 2′,5′-oligoadenylate synthetase (2’5’AS) enzyme activity is controlled by a single-nucleotide polymorphism at a splice-acceptor site in the OAS1 gene. Am J Hum Genet. 2005;76(4):623–33.
Gabriel SB, et al. The structure of haplotype blocks in the human genome. Science. 2002;296(5576):2225–9.
Thornton-Wells TA, et al. Genetics, statistics and human disease: analytical retooling for complexity. Trends Genet. 2004;20(12):640–7.
Cordell HJ, et al. Case/pseudocontrol analysis in genetic association studies: a unified framework for detection of genotype and haplotype associations, gene-gene and gene-environment interactions, and parent-of-origin effects. Genet Epidemiol. 2004;26(3):167–85.
Lee PH, Shatkay H. F-SNP: computationally predicted functional SNPs for disease association studies. Nucleic Acids Res. 2008;36(Database issue):D820–4.
Hamano E, et al. Polymorphisms of interferon-inducible genes OAS-1 and MxA associated with SARS in the Vietnamese population. Biochem Biophys Res Commun. 2005;329(4):1234–9.
Lee Y, et al. On the analysis of a repeated measure design in genome-wide association analysis. Int J Environ Res Public Health. 2014;11(12):12283–303.
Kim YJ, et al. Large-scale genome-wide association studies in east Asians identify new genetic loci influencing metabolic traits. Nat Genet. 2011;43(10):990–5.
The authors are grateful for all the cases and the controls who volunteered for participating in the study.
This work was supported in part by the National Natural Science Foundation of China (Beijing, China) Grants 81870015, 81072432. This funding supported the costs of the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The data analyzed during the study are not publicly available as is they are to be used for further research work on the topic. However, the data can be obtained from the authors upon request.
Ethics approval and consent to participate
The study was approved by the ethical committee of the West China Hospital of Sichuan University (the NO.131 file). The ethical approval is as follows. All participants in both groups voluntarily joined this study with written informed consents. A sample of 5 ml of peripheral blood was obtained from each participant.
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.