Association of IL-17F rs2397084 (E126G), rs11465553 (V155I) and rs763780 (H161R) variants with rheumatoid arthritis and their effects on the stability of protein.

Interleukin-17F (IL-17F), considered a pro-inflammatory cytokine, has been shown to contribute to skeletal tissue degradation and hence chronic inflammation in rheumatoid arthritis (RA). In this study we utilized bioinformatics tools to analyze the effect of three exonic SNPs (rs2397084, rs11465553, and rs763780) on the structure and function of the IL-17F gene, and evaluated their association with RA in Pakistani patients. The predicted deleterious and damaging effects of identified genetic variants were assessed through the utilization of multiple bioinformatics tools including PROVEAN, SNP&GO, SIFT, and PolyPhen2. Structural and functional effects of these variants on protein structures were evaluated through the use of additional tools such as I-Mutant, MutPred, and ConSurf. Three-dimensional (3D) models of both the wild-type and mutant proteins were constructed through the utilization of I-TASSER software, with subsequent structural comparisons between the models conducted through the use of the TM-align score. A total of 500 individuals, 250 cases and 250 controls, were genotyped through Tri-ARMS-PCR method and the resultant data was statistically analyzed using various inheritance models. Our bioinformatics analysis showed significant structural differences for wild type and mutant protein (TM-scores and RMSD values were 0.85934 and 2.34 for rs2397084 (E126G), 0.87388 and 2.49 for rs11465553 (V155I), and 0.86572 and 0.86572 for rs763780 (H161R) with decrease stability for the later. Overall, these tools enabled us to predict that these variants are crucial in causing disease phenotypes. We further tested each of these single nucleotide variants for their association with RA. Our analysis revealed a strong positive association between the genetic variant rs763780 and the risk of developing rheumatoid arthritis (RA) at both the genotypic and allelic levels. The genotypic association was statistically significant[χ2 = 111.8; P value <0.0001], as was the allelic level [OR 3.444 (2.539-4.672); P value 0.0008]. These findings suggest that the presence of this genetic variant may increase the susceptibility to RA. Similarly, we observed a significant distribution of the genetic variant rs11465553 at the genotypic level [χ2 = 25.24; P value = 0.0001]. However, this variant did not show a significant association with RA at the allelic level [OR = 1.194 (0.930-1.531); P value = 0.183]. However, the distribution of variant rs2397084 was more or less random across our sample with no significant association either at genotypic and or allelic level. Put together, our association study and in silico prediction of decreasing of IL17-F protein stabilty confirmed that two SNPs, rs11465553 and rs763780 are crucial to the suscetibility of and showed that these RA in Pakistani patients.

Exploring the Associations and Molecular Impacts of miR-146a/rs2910164 and miR-196a2/rs185070757 with Rheumatoid Arthritis in a Pakistani Population.

INTRODUCTION: MicroRNAs (miRNAs) are a new class of molecules that participate in post-transcriptional regulation of gene expression and hence have been reported to have a crucial role in the pathogenesis of rheumatoid arthritis (RA). We aimed to investigate the association of miR-146a rs2910164 (G/C) and miR-196a2 rs185070757 (T/G) with RA susceptibility in Pakistani patients and to bioinformatically predict the molecular function of these miRNAs. METHODS: A case-control study on 600 individuals was conducted, including 300 RA patients and 300 matching healthy controls. Genotyping was performed by tetra-primer amplification of refractory mutation system-polymerase chain reaction, and the association between variants and RA was statistically determined using different models. RESULTS: For the variant rs2910164 (G/C) in miR-146a, no difference in genotype distribution was observed between RA cases and controls, as determined using co-dominant (χ2 = 4.33; p = 0.114), homozygous dominant (C/C vs. G/G + C/G) (OR = 0.740 [0.531-1.032]; p = 0.091), homozygous recessive (G/G vs. C/C + G/C) (odds ratio [OR] = 01.432 [0.930-2.206]; p = 0.126), heterozygous (G/C vs. C/C + G/G) (OR = 1.084 [0.786-1.494]; p = 0.682), and additive (OR 0.778 [0.617-0.981]; p = 0.039) models. Similarly, the GT genotype in the rs185070757 (T/G) miR-196a2 variant did not differ between cases and controls with any models (p > 0.05). For the first time, we report no association of miR-146a rs2910164 (G/C) and miR-196a2 rs185070757 (T/G) with RA in a Pakistani population. A subsequent bioinformatic analysis revealed that the CC genotype of miR-146a rs2910164 might have a protective role against RA pathogenesis, with no effect observed with the miR-196a2 rs185070757. CONCLUSION: Our findings suggest that these miRNAs might have little-to-no impact on the RA pathogenesis in the Pakistani population.

Phylogenetic analysis of promoter regions of human Dolichol kinase (DOLK) and orthologous genes using bioinformatics tools.

The Dolichol kinase (DOLK) gene encodes the polytopic DOLK protein associated with the endoplasmic reticulum (ER) N-glycosylation pathway catalyzing the final step in the biosynthesis of dolichol phosphate. Dolichol phosphate is an oligosaccharide carrier required for N-glycosylation of DOLK protein, with its deficiency leading to a severe hypo glycosylation phenotype in humans which can cause congenital disorders of glycosylation and death in early infancy. The aim of the present study is to identify the phylogenetic relationship between human and ortholog species based on their conserved sequences in DOLK gene. Sequence alignment of DOLK was carried out in this study and the evolutionarily conserved regulatory sequences were identified using bioinformatics. Promoter sequence of human DOLK was compared with orthologous sequences from different organisms. Conserved non-coding sequences (CNS) and motifs in promoter regions were found by analyzing upstream promoter sequences of Homo sapiens DOLK and its orthologous genes in other organisms. Conserved sequences were predicted in the promoter regions in CNS1 and CNS2. Conserved protein sequences were also identified by alignment of the orthologous sequences. Organisms with similar gene sequences are assumed to be closely related and the ER N-glycosylation pathway is conserved in them.

Exome Sequencing Identifies a Novel FBN1 Variant in a Pakistani Family with Marfan Syndrome That Includes Left Ventricle Diastolic Dysfunction.

INTRODUCTION: Cardiomyopathies are diseases of the heart muscle and are important causes of heart failure. Dilated cardiomyopathy (DCM) is a common form of cardiomyopathy that can be acquired, syndromic or non-syndromic. The current study was conducted to explore the genetic defects in a Pakistani family with cardiac disease and features of Marfan's syndrome (MFS). METHODS: A family with left ventricle (LV) diastolic dysfunction and MFS phenotype was assessed in Pakistan. The clinical information and blood samples from the patients were collected after physical, cardiovascular, and ophthalmologic examinations. An affected individual (proband) was subjected to whole-exome sequencing (WES). The findings were further validated through Sanger sequencing in the family. RESULTS: Through WES and sanger validation, we identified a novel variant NM_000138.4; c.1402A>G in the Fibrillin-1 (FBN1) gene that segregates with LV diastolic dysfunction and MFS. Furthermore, bioinformatic evaluation suggested that the novel variant is deleterious and disease-causing. CONCLUSIONS: This study identified for the first time a novel FBN1 variant in a family with LV diastolic dysfunction and MFS in Pakistan.

Association of AFF3 Gene Polymorphism rs10865035 with Rheumatoid Arthritis: A Population-Based Case-Control Study on a Pakistani Cohort.

Rheumatoid arthritis (RA) is one of the complex diseases with the involvement of the genetic as well as environmental factors in its onset and severity. Different genome-wide association and candidate gene studies have shown the role of several genetic variants in multiple loci/genes with ethnical and geographical variations. This study was designed to detect the association of a single-nucleotide polymorphism (SNP) rs10865035 in the AFF3 gene with the genetic background of rheumatoid arthritis (RA) in the Pakistani cohort. A total of 703 individuals, including 409 RA patients and 294 healthy controls, were genotyped using TaqMan assay and Tri primer ARMS-PCR (amplification-refractory mutation system-polymerase chain reaction) methods. The association of rs10865035 with the RA was statistically determined using different models. Interestingly, besides the homozygous recessive model (G/G vs. A/G + A/A) (OR = 1.693(1.06-2.648); P = 0.025), all other models, which included the codominant (χ 2 = 5.169; P = 0.075), homozygous dominant (A/A vs. G/G + A/G) (OR = 0.867 (0.636-1.187); P = 0.41), heterozygous (A/G vs. A/A + GG) (OR = 0.491 (0.667-1.215); P = 0.49), and additive model (OR = 0.826 (0.665-1.027); P = 0.08) showed insignificant distribution of the genotypes among the cases and controls. These findings suggest that the AFF3 gene (rs10865035) has no significant role in the onset of RA in the Pakistani population.