A novel missense compound heterozygous variant in TLR1 gene is associated with susceptibility to rheumatoid arthritis - structural perspective and functional annotations

Introduction Besides human leukocyte antigen (HLA-DRB1) locus, more than 100 loci across the genome have been identified and linked with the onset, expression and/or progression of rheumatoid arthritis (RA). However, there are still grey areas in our understanding of the key genetic contributors of the disease, particularly in familial cases.Methods In the present study, we have performed the whole exome sequencing (WES) of RA patients from two consanguineous families of Pakistan in a quest to identify novel, high-impact, RA-susceptibility genetic variants.Results Through stepwise filtering, around 17,000 variants (common in the affected members) were recognized, out of which 2651 were predicted to be deleterious. Of these, 196 had direct relevance to RA. When selected for homozygous recessive mode of inheritance, two novel pathogenic variants (c.1324T>C, p.Cys(442)?Arg(442); c.2036T>C, p.Ile(679)?Thr(679)) in the TLR1 gene displayed the role of compound heterozygosity in modulating the phenotypic expression and penetrance of RA. The structural and functional consequences of the TLR1 missense single nucleotide mutations (Cys(442)?Arg(442); Ile(679)?Thr(679)) were evaluated through molecular dynamic simulation (MDS) studies. Analysis showed domain's rigidification, conferring stability to mutant TLR1-TIR/TIRAP-TIR complex with concomitant increase in molecular interactions with pro-inflammatory cytokines, compared to the wild-type conformation. Gene co-expression network analysis highlighted interlinked partnering genes along with interleukin-6 production of TLR1 (corrected p-value 2.98e-4) and acetylcholine receptor activity of CHRNG (corrected p-value 6.12e-2) as highly enriched associated functions.Conclusion The results, validated through case-control study subjects, suggested that the variants identified through WES and confirmed through Sanger sequencing and MDS are the novel disease variants and are likely to confer RA-susceptibility, independently and/or in a family-specific context.