2′-5′-Oligoadenylate synthetase 1 polymorphisms are associated with tuberculosis: a case-control study

Background

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).

Methods

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.

Results

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.

Conclusions

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.

     

Potential association between <Emphasis Type="Italic">ITPKC</Emphasis> genetic variations and Hirschsprung disease

Hirschsprung disease (HSCR) is a congenital and complex disorder characterized by intestinal obstruction due to the absence of enteric neurons along variable lengths of the hindgut. Our recent genome-wide association study (GWAS) has revealed regional associations with HSCR at several loci of inositol-trisphosphate 3-kinase C (ITPKC). For fine mapping, we additionally selected and genotyped a total of 12 single nucleotide polymorphisms (SNPs) of ITPKC in 187 HSCR patients and 283 unaffected controls, and performed a further combined imputation analysis based on genotype data from this second stage of fine mapping and our previous GWAS stage, totaling 902 subjects (187 HSCR cases and 715 controls). As a result, several SNPs (minimum P?=?0.004) and a haplotype (P?=?0.02) were found to be significantly associated with HSCR. In further in silico analyses to ascertain the potential functions of the significant variants, the change from the common allele to the rare allele of the highly conserved nonsynonymous rs76785336 showed a difference in mRNA folding structure. In the case of intronic SNPs, rs2607420 with a high consensus value was predicted to be a new splice site. Although this study has limitations (such as lack of functional evaluations, small number of cases, and further need of replication in other cohorts), our findings suggest that genetic variants of ITPKC may have a potential association with HSCR susceptibility and/or developmental diseases related to enteric nervous system development.     

Association study of genetic markers of schizophrenia and its cognitive endophenotypes

A replicative analysis of associations of 15 SNPs located in the regions of 11 genes (TCF4, VRK2, NOTCH4, ZNF804A, AGBL1, RELN, ZFP64P1, KCNB2, CSMD1, CPVL, NRIP1) and three intergenic regions (SLCO6A1/LINCOO491, LOC105376248/LOC105376249, SPA17/NRGN) with schizophrenia was conducted in the Russian population of the Siberian region. These SNPs were previously identified in genome-wide association studies (GWAS) of schizophrenia and cognitive abnormalities. The present study confirmed associations of KCNB2 rs2247572, CSMD1 rs2616984, and intergenic rs12807809 located in SPA17/NRGN with schizophrenia. It was established that the frequency of the CSMD1 rs2616984 G/G genotype was higher in patients compared to the control group (OR = 1.73; CI: 1.14–2.62; р = 0.0337). The frequencies of the KCNB2 rs2247572 TT genotype (OR = 0.41; CI: 0.20–0.87; р = 0.0485) and intergenic rs12807809 CT genotype located in SPA17/NRGN (OR = 0.70; CI: 0.53–0.94; р = 0.0464) were significantly decreased in patients compared to the control group.     

Association of three single nucleotide polymorphisms at the SIX1-SIX6 locus with primary open angle glaucoma in the Chinese population

This study investigated the genetic association of three single nucleotide polymorphisms(SNPs; rs10483727, rs33912345, and rs146737847) at the SIX1-SIX6 locus with primary open angle glaucoma(POAG) in the Chinese population. A total of 866 subjects with POAG(685 high-tension glaucoma(HTG) and 181 normal-tension glaucoma(NTG)) and 266 control individuals were included. Significant genetic association was identified for rs10483727 in HTG(P=0.02; odds ratio(OR)=1.31), NTG(P=7.41×10~(-6); OR=2.71), and POAG(i.e., HTG and NTG combined; P=0.001; OR=1.44). rs33912345 was also significantly associated with HTG(P=0.008; OR=1.36), NTG(P=2.72×10~(-6); OR=2.27), and POAG(P=3.84×10~(-4); OR=1.49). The rare SIX6 mutation, rs146737847, was not found in the subjects enrolled in this study. Stratification by patient age identified that both rs10483727 and rs33912345 were significantly associated with NTG in patients aged above 40 years(P=2.08×10~(-5); OR=2.28), whereas in patients aged between 20–40 years, rs33912345 was significantly associated with NTG(P=0.017; OR=2.06). In HTG, the genetic associations for both rs10483727 and rs33912345 were significant in patients aged between 20–40 years(P=0.006; OR=1.56) but not in those aged above 40 years(P=0.118, OR=1.21 and P=0.042, OR=1.29, respectively). This study replicated the association of POAG with two SNPs at the SIX1-SIX6 locus and demonstrated that SNPs, rs10483727 and rs33912345, are significantly associated with POAG, especially with NTG in patients aged above 40 years.     

Association analysis of ILVBL gene polymorphisms with aspirin-exacerbated respiratory disease in asthma

Background

We previously reported that the ILVBL gene on chromosome 19p13.1 was associated with the risk for aspirin-exacerbated respiratory disease (AERD) and the percent decline of forced expired volume in one second (FEV1) after an oral aspirin challenge test. In this study, we confirmed the association between polymorphisms and haplotypes of the ILVBL gene and the risk for AERD and its phenotype.

Methods

We recruited 141 AERD and 995 aspirin-tolerant asthmatic (ATA) subjects. All study subjects underwent an oral aspirin challenge (OAC). Nine single nucleotide polymorphisms (SNPs) with minor allele frequencies above 0.05, which were present in the region from 2 kb upstream to 0.5 kb downstream of ILVBL in Asian populations, were selected and genotyped.

Results

In an allelic association analysis, seven of nine SNPs were significantly associated with the risk for AERD after correction for multiple comparisons. In a codominant model, the five SNPs making up block2 (rs2240299, rs7507755, rs1468198, rs2074261, and rs13301) showed significant associations with the risk for AERD (corrected P?=?0.001–0.004, OR?=?0.59–0.64). Rs1468198 was also significantly associated with the percent decline in FEV1 in OAC tests after correction for multiple comparisons in the codominant model (corrected P?=?0.033), but the other four SNPs in hapblock2 were not.

Conclusion

To the best of our knowledge, this is the first report of an association between SNPs on ILVBL and AERD. SNPs on ILVBL could be promising genetic markers of this condition.

     

Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached “suggestive” genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.     

Effects of <Emphasis Type="Italic">HLA</Emphasis>-<Emphasis Type="Italic">DRB1</Emphasis>/<Emphasis Type="Italic">DQB1</Emphasis> Genetic Variants on Neuroimaging in Healthy,Mild Cognitive Impairment,and Alzheimer’s Disease Cohorts

Alzheimer’s disease (AD) is the most common form of dementia and exhibits a considerable level of heritability. Previous association studies gave evidence for the associations of HLA-DRB1/DQB1 alleles with AD. However, how and when the gene variants in HLA-DRB1/DQB1 function in AD pathogenesis has yet to be determined. Here, we firstly investigated the association of gene variants in HLA-DRB1/DQB1 alleles and AD related brain structure on magnetic resonance imaging (MRI) in a large sample from the Alzheimer’s Disease Neuroimaging Initiative (ADNI). We selected hippocampus, subregion, parahippocampus, posterior cingulate, precuneus, middle temporal, entorhinal cortex, and amygdala as regions of interest (ROIs). Twelve SNPs in HLA-DRB1/DQB1 were identified in the dataset following quality control measures. In the total group hybrid population analysis, our study (rs35445101, rs1130399, and rs28746809) were associated with the smaller baseline volume of the left posterior cingulate and rs2854275 was associated with the larger baseline volume of the left posterior cingulate. Furthermore, we detected the above four associations in mild cognitive impairment (MCI) sub-group analysis, and two risk loci (rs35445101 and rs1130399) were also the smaller baseline volume of the left posterior cingulate in (NC) sub-group analysis. Our study suggested that HLA-DRB1/DQB1 gene variants appeared to modulate the alteration of the left posterior cingulate volume, hence modulating the susceptibility of AD.     

Field Synopsis and Re-analysis of Systematic Meta-analyses of Genetic Association Studies in Multiple Sclerosis: a Bayesian Approach

To provide an up-to-date summary of multiple sclerosis-susceptible gene variants and assess the noteworthiness in hopes of finding true associations, we investigated the results of 44 meta-analyses on gene variants and multiple sclerosis published through December 2016. Out of 70 statistically significant genotype associations, roughly a fifth (21%) of the comparisons showed noteworthy false-positive rate probability (FPRP) at a statistical power to detect an OR of 1.5 and at a prior probability of 10^?6 assumed for a random single nucleotide polymorphism. These associations (IRF8/rs17445836, STAT3/rs744166, HLA/rs4959093, HLA/rs2647046, HLA/rs7382297, HLA/rs17421624, HLA/rs2517646, HLA/rs9261491, HLA/rs2857439, HLA/rs16896944, HLA/rs3132671, HLA/rs2857435, HLA/rs9261471, HLA/rs2523393, HLA-DRB1/rs3135388, RGS1/rs2760524, PTGER4/rs9292777) also showed a noteworthy Bayesian false discovery probability (BFDP) and one additional association (CD24 rs8734/rs52812045) was also noteworthy via BFDP computation. Herein, we have identified several noteworthy biomarkers of multiple sclerosis susceptibility. We hope these data are used to study multiple sclerosis genetics and inform future screening programs.     

MicroRNA Biogenesis Pathway Gene Polymorphisms Are Associated with Breast Cancer Risk

MicroRNAs (miRNAs) play an important role as epigenetic regulators in cancer initiation and progression. One of the mechanisms of miRNA dysregulation is altered functioning of proteins involved in miRNA processing machinery. It has been suggested that single nucleotide polymorphisms (SNPs) within miRNA gene regions, miRNA target genes, and miRNA machinery genes may affect the miRNAs regulation. We selected 25 SNPs in the key genes of miRNA biosynthesis, including DROSHA/RNASEN, DGCR8, DICER1, XPO5, RAN, PIWIL1/HIWI, AGO1/EIF2C1, AGO2, GEMIN4, GEMIN3/DDX20, and DDX5, and investigated the association between these SNPs and the risk of breast cancer. The total number of breast cancer cases and cancer-free controls enrolled in the investigation were 778 (417 breast cancer patients and 361 healthy women). We found that rs11060845 and rs10773771 in the PIWIL1 gene, rs3809142/RAN, rs10719/DROSHA, rs1640299/DGCR8, rs563002/DDX20, rs595055/AGO1, and rs2740348/GEMIN4 were associated with breast cancer risk in Russians.     

Genetic Association of <Emphasis Type="Italic">HLA</Emphasis> Gene Variants with MRI Brain Structure in Alzheimer’s Disease

There is accumulating evidence that the human leukocyte antigen (HLA) gene variants are associated with Alzheimer’s disease (AD). However, how they affect AD occurrence is still unknown. In this study, we firstly investigated the association of gene variants in HLA gene variants and brain structures on MRI in a large sample from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) to explore the effects of HLA on AD pathogenesis. We selected hippocampus, hippocampus CA1 subregion, parahippocampus, posterior cingulate, precuneus, middle temporal, entorhinal cortex, and amygdala as regions of interest (ROIs). According to the previous association studies of HLA variants and AD, 12 SNPs in HLA were identified in the dataset following quality control measures. In total group analysis, our results showed that TNF-α SNPs at rs2534672 and rs2395488 were significantly positively associated with the volume of the left middle temporal lobe (rs2534672: P?=?0.00035, Pc?=?0.004; rs2395488: P?=?0.0038, Pc?=?0.023) at baseline. In the longitudinal study, HFE rs1800562 was remarkably correlated with the lower atrophy rate of right middle temporal lobe (P?=?0.0003, Pc?=?0.003) and RAGE rs2070600 was associated with the atrophy rate of right hippocampus substructure-CA1 over 2 years (P?=?0.003, Pc?=?0.035). Furthermore, we detected the above four associations in mild cognitive impairment (MCI) subgroup analysis, as well as the association of rs2534672 with the baseline volume of the left middle temporal lobe in normal cognition (NC) subgroup analysis. Our study provided preliminary evidences that HLA gene variants might participate in the structural alteration of AD associated brain regions, hence modulating the susceptibility of AD.     

Haplotype-specific <Emphasis Type="Italic">MAPT</Emphasis> exon 3 expression regulated by common intronic polymorphisms associated with Parkinsonian disorders

Background

Genome wide association studies have identified microtubule associated protein tau (MAPT) H1 haplotype single nucleotide polymorphisms (SNPs) as leading common risk variants for Parkinson’s disease, progressive supranuclear palsy and corticobasal degeneration. The MAPT risk variants fall within a large 1.8 Mb region of high linkage disequilibrium, making it difficult to discern the functionally important risk variants. Here, we leverage the strong haplotype-specific expression of MAPT exon 3 to investigate the functionality of SNPs that fall within this H1 haplotype region of linkage disequilibrium.

Methods

In this study, we dissect the molecular mechanisms by which haplotype-specific SNPs confer allele-specific effects on the alternative splicing of MAPT exon 3. Firstly, we use haplotype-hybrid whole-locus genomic MAPT vectors studies to identify functional SNPs. Next, we characterise the RNA-protein interactions at two loci by mass spectrometry. Lastly, we knockdown candidate splice factors to determine their effect on MAPT exon 3 using a novel allele-specific qPCR assay.

Results

Using whole-locus genomic DNA expression vectors to express MAPT haplotype variants, we demonstrate that rs17651213 regulates exon 3 inclusion in a haplotype-specific manner. We further investigated the functionality of this region using RNA-electrophoretic mobility shift assays to show differential RNA-protein complex formation at the H1 and H2 sequence variants of SNP rs17651213 and rs1800547 and subsequently identified candidate trans-acting splicing factors interacting with these functional SNPs sequences by RNA-protein pull-down experiment and mass spectrometry. Finally, gene knockdown of candidate splice factors identified by mass spectrometry demonstrate a role for hnRNP F and hnRNP Q in the haplotype-specific regulation of exon 3 inclusion.

Conclusions

We identified common splice factors hnRNP F and hnRNP Q regulating the haplotype-specific splicing of MAPT exon 3 through intronic variants rs1800547 and rs17651213. This work demonstrates an integrated approach to characterise the functionality of risk variants in large regions of linkage disequilibrium.

     

Association of <Emphasis Type="Italic">ACYP2</Emphasis> and <Emphasis Type="Italic">TSPYL6</Emphasis> Genetic Polymorphisms with Risk of Ischemic Stroke in Han Chinese Population

The development of ischemic stroke is associated with advanced age. Telomere length, as a marker of biological aging, has been reported to influence the risk of several age-related diseases, including ischemic stroke. Recent studies have identified the genetic variant within ACYP2 and TSPYL6 associated with shorter telomere length. The objective of this study is to investigate the putative association of ischemic stroke with common polymorphisms in ACYP2 and TSPYL6 genes in a Chinese Han population. We found that the risk alleles of six single nucleotide polymorphisms (SNPs), including rs11125529, rs12615793, rs843711, rs11896604, and rs843706 within both ACYP2 and TSPYL6, and rs17045754 in ACYP2 gene, were related with increased risk of ischemic stroke according to both allelic and genotype association analyses. The significant correlations between ACYP2 and TSPYL6 SNPs and ischemic stroke risk were also observed in dominant, recessive, and additive models, respectively. Two blocks in high linkage disequilibrium were identified in this study, and two haplotypes were associated with higher ischemic stroke susceptibility. In conclusion, the genetic polymorphisms of ACYP2 and TSPYL6 are associated with increased risk of developing ischemic stroke. Further studies with larger sample sizes are required to validate our findings.     

<Emphasis Type="Italic">PPARGC1B</Emphasis> gene is associated with Kashin-Beck disease in Han Chinese

Kashin-Beck disease (KBD) is a chronic osteochondropathy. The genetic basis of KBD remains elusive now. To investigate the relationship between PPARGC1B gene polymorphism and KBD, we conducted a two-stage association study using 2743 unrelated Han Chinese subjects. In the first stage, three SNPs rs1078324, rs4705372, and rs11743128 of PPARGC1B gene were genotyped in 559 KBD patients and 467 health controls using Sequenom MassARRAY platform. In the second stage, the association analysis results of PPARGC1B with KBD were replicated using an independent sample of 1717 subjects. SNP association analysis was conducted by PLINK software. Genotype imputation was conducted by IMPUTE 2.0 against the reference panel of the 1000 genome project. Bonferroni multiple testing correction was performed. We observed a significant association signal at rs4705372 (P?=?0.0160) and a suggestive association signal at rs11743128 (P?=?0.0290). Further replication study confirmed the association signals of rs4705372 (P?=?0.0026) and rs11743128 (P?=?0.0387) in the independent validation sample. Our study results suggest that PPARGC1B is a novel susceptibility gene of KBD.     

The analysis of association between type 2 diabetes and polymorphic markers in the <Emphasis Type="Italic">CDKAL1</Emphasis> gene and in the <Emphasis Type="Italic">HHEX/IDE</Emphasis> locus

The increase in diabetes was noted at the turn of the 21st century. Patients with type 2 diabetes (T2DM) make up the majority of patients. Diabetes is a multifactorial disease. It arises from adverse effects of environmental factors on the body of genetically susceptible peoples. According to modern concepts, T2DM is a polygenic disease. Each of the involved genes contributes to the risk of developing of this disease. In our study, the association between polymorphic genetic markers rs7756992, rs9465871, rs7754840, and rs10946398 in the CDKAL1 gene and rs1111875 in the HHEX/IDE locus and T2DM in the Russian population were studied. Four hundred forty patients with type 2 diabetes and 264 healthy individuals without any signs of the disease were examined. The comparative analysis of distribution of genotypes and allele frequencies points to an association between polymorphic genetic markers rs7756992, rs9465871, and rs10946398 in the CDKAL1 gene and this disease. For the other polymorphic genetic markers (rs7754840 in the CDKAL1 gene and rs1111875 in the HHEX/IDE locus), no statistically significant associations are found. On the basis of these data, we can conclude that the CDKAL1 gene is associated with development of T2DM. For the HHEX/IDE locus, such an association is absent.     

Association of SNPs in dopamine and serotonin pathway genes and their interacting genes with temperament traits in Charolais cows

Cattle temperament is a complex trait, and molecular studies aimed at defining this trait are scarce. We used an interaction networks approach to identify new genes (interacting genes) and to estimate their effects and those of 19 dopamine- and serotonin-related genes on the temperament traits of Charolais cattle. The genes proopiomelanocortin (POMC), neuropeptide Y (NPY), solute carrier family 18, member 2 (SLC18A2) and FBJ murine osteosarcoma viral oncogene homologue (FOSFBJ) were identified as new candidates. Their potential to be associated with temperament was estimated according to their reported biological activities, which included interactions with neural activity, receptor function, targeting or synthesis of neurotransmitters and association with behaviour. Pen score (PS) and exit velocity (EV) measures were determined from 412 Charolais cows to calculate their temperament score (TS). Based on the TS, calm (n?=?55; TS, 1.09?±?0.33) and temperamental (n?=?58; TS, 2.27?±?0.639) cows were selected and genotyped using a 248 single-nucleotide variation (SNV) panel. Of the 248 variations in the panel, only 151 were confirmed to be polymorphic (single-nucleotide polymorphisms; SNPs) in the tested population. Single-marker association analyses between genotypes and temperament measures (EV, PS and/or TS) indicated significant associations of six SNPs from four candidate genes. The markers rs109576799 and rs43696138, located in the DRD3 and HTR2A genes, respectively, were significantly associated with both EV and TS traits. Four markers, rs110365063 and rs137756569 from the POMC gene and rs110365063 and rs135155082 located in SLC18A2 and DRD2, respectively, were associated with PS. The variant rs110365063 located in bovine SLC18A2 causes a change in the amino acid sequence from Ala to Thr. Further studies are needed to confirm the association of genetic profile with cattle temperament; however, our study represents important progress in understanding the regulation of cattle temperament by different genes with divergent functions.     

Association of the apolipoprotein M gene polymorphisms and serum lipid levels

The association of rs707921 and rs707922 SNPs in the apolipoprotein M (APOM) gene and serum lipid levels is still controversial. This study aimed to detect the association of the APOM rs707921 and rs707922 SNPs and several environmental factors with serum lipid profiles. Genotyping of rs707921 and rs707922 was performed in 703 of Mulao’s and 707 of Han’s participants. The serum levels of TG in Mulao, and TG and HDL-C in Han were different between the A and C allele carriers of rs707921 (P < 0.05–0.01); while the serum levels of TG in both Mulao and Han were different between the T and G allele carriers of rs707922 (P < 0.05–0.01). According to the gender-subgroup analysis, the levels of TC in Mulao females, TG and ApoB in Han males, and HDL-C in Han females were associated with the genotypes of rs707921 (P < 0.05 for each); whereas the levels of TG in Mulao males, and TG and ApoB in Han males were correlated with the genotypes of rs707922 (P < 0.05 for each). Serum lipid parameters were also associated with several environmental factors (P < 0.05–0.001). The APOM gene rs707921 and rs707922 SNPs are associated with some serum lipid parameters in the two ethnic groups, but the trends of association suggest that the two SNPs might have racial/ethnic- and/or gender- specificity.     

Association of <Emphasis Type="Italic">Toll like receptor 2</Emphasis> and <Emphasis Type="Italic">9</Emphasis> gene variants with pulmonary tuberculosis: exploration in a northern Indian population

Tuberculosis (TB) is a disease of global importance. There is an increasing recognition of the role of Toll like receptors, important pattern recognition receptors of host immune system, in determining the susceptibility or resistance to TB in various populations. In an attempt to examine the importance of Toll like receptors in immune response to Mycobacterium tuberculosis infection, we explored two variants each of TLR2 and TLR9 in a population residing in Uttar Pradesh, India. Genotyping was performed to detect -196 to -174 del polymorphism and G2258A SNP (Arg753Gln, rs5743708) in TLR2 gene and -T1237C (rs5743836) and G2848A (rs352140) SNP in TLR9 gene in patients with pulmonary TB and healthy controls. The A allele of G2848A SNP in TLR9 gene was found with a marginally higher frequency among TB patients as compared to healthy controls, suggesting that A allele at position 2848 of TLR9 gene may be associated with susceptibility to TB in North Indian population [p?=?0.05, Mantel–Haenszel OR?=?1.34, 95% CI (1.0–1.82)].     

The association of AKNA gene polymorphisms with knee osteoarthritis suggests the relevance of this immune response regulator in the disease genetic susceptibility

Recent studies have identified AKNA as a potential susceptibility gene for several inflammatory diseases. Here, we aimed to assess the potential association of AKNA polymorphisms with knee osteoarthritis (KOA) susceptibility in a Mexican population, following STREGA recommendations. From a DNA bank of 181 KOA patients and 140 healthy controls, two AKNA SNPs were genotyped using TaqMan probes. The association between KOA susceptibility and AKNA polymorphisms genotypes was evaluated by multivariated logistic regression analysis. Information regarding patients’ inflammatory biomarkers levels was obtained and their association with AKNA polymorphisms genotypes was assessed by lineal regression. We found a positive association with the recessive inheritance model of both AKNA polymorphisms (A/A genotype for both) and KOA susceptibility adjusting by age, body mass index (BMI), gender and place of birth (OR?=?2.48, 95% CI 1.09–5.65 for rs10817595 polymorphism; and OR?=?4.96; 95% CI 2.421–10.2 for rs3748176 polymorphism). Additionally these associations were also seen after stratifying patients by KOA severity and age. Furthermore the total leukocyte count was positively associated with rs10817595 AKNA polymorphism (β?=?1.39; 95% CI 0.44–2.34) adjusting by age, BMI, gender, place of birth and disease severity. We suggest that regulatory and coding polymorphisms of the inflammatory modulator gene AKNA can influence the development of KOA. Further structural and functional studies might reveal the role of AKNA in OA and other rheumatic diseases.