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₄₄₂→Arg₄₄₂; c.2036T>C, p.Ile₆₇₉→Thr₆₇₉) 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₄₄₂→Arg₄₄₂; Ile₆₇₉→Thr₆₇₉) 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. Key Points • Exploration of population based/ethno-specific big data is imperative to identify novel causal variants of RA. • Two new deleterious missense mutations in mutational hotspot exon 4 of TLR1 gene have been identified in Pakistani RA patients. • MD simulation data provides evidence for domain’s rigidification, conferring stability to mutant TLR1-TIR/TIRAP-TIR complex, with concomitant increase in production of pro-inflammatory cytokines, thus adding to the onset/erosive outcome of RA.