In-Silico Computing of the Most Deleterious nsSNPs in HBA1 Gene

Background

α-Thalassemia (α-thal) is a genetic disorder caused by the substitution of single amino acid or large deletions in the HBA1 and/or HBA2 genes.

Method

Using modern bioinformatics tools as a systematic in-silico approach to predict the deleterious SNPs in the HBA1 gene and its significant pathogenic impact on the functions and structure of HBA1 protein was predicted.

Results and Discussion

A total of 389 SNPs in HBA1 were retrieved from dbSNP database, which includes: 201 non-coding synonymous (nsSNPs), 43 human active SNPs, 16 intronic SNPs, 11 mRNA 3′ UTR SNPs, 9 coding synonymous SNPs, 9 5′ UTR SNPs and other types. Structural homology-based method (PolyPhen) and sequence homology-based tool (SIFT), SNPs&Go, PROVEAN and PANTHER revealed that 2.4% of the nsSNPs are pathogenic.

Conclusions

A total of 5 nsSNPs (G60V, K17M, K17T, L92F and W15R) were predicted to be responsible for the structural and functional modifications of HBA1 protein. It is evident from the deep comprehensive in-silico analysis that, two nsSNPs such as G60Vand W15R in HBA1 are highly deleterious. These “2 pathogenic nsSNPs” can be considered for wet-lab confirmatory analysis.     

Pathway analysis of genome-wide association study on serum prostate-specific antigen levels

The wide application of prostate-specific antigen (PSA) has contributed to the early diagnosis and improved management of prostate cancer (PCa). Accumulating evidence has suggested the involvement of genetic components in regulating serum PSA levels, and several single nucleotide polymorphisms (SNPs) have been identified by genome-wide association studies (GWASs). However, the GWASs' results have the limited power to identify the causal variants and pathways. After the quality control filters, a total of 330,540 genotyped SNPs from one GWAS with 657 PCa-free Caucasian males were included for the identify candidate causal SNPs and pathways (ICSNPathway) analysis. In addition, the genotype–phenotype association analysis has been conducted with the data from HapMap database. Overall, a total of four SNPs in three genes and six pathways were identified by ICSNPathway analysis, which in total provided three hypothetical mechanisms. First, CYP26B1 rs2241057 polymorphism (nonsynonymous coding) which leads to a Leu-to-Ser amino acid shift at position 264, was implicated in the pathways including meiosis, proximal/distal pattern formation, and M phase of meiotic cell cycle. Second, CLIC5 rs3734207 and rs11752816 polymorphisms (regulatory region) to the 2 iron, 2 sulfur cluster binding pathway through regulating expression levels of CLIC5 mRNA. Third, rs4819522 polymorphism (nonsynonymous coding) leads to a Thr-to-Met transition at position 350 of TBX1 and involves in the pathways about gland and endocrine system development. In summary, our results demonstrated four candidate SNPs in three genes (CYP26B1 rs2241057, CISD1 rs2251039, rs2590370, and TBX1 rs4819522 polymorphisms), which were involved in six potential pathways to influence serum PSA levels.     

A genetic assay of three patients in the same family with Holt-Oram syndrome; a case report

Holt-Oram syndrome (HOS) is a developmental disorder inherited in an autosomal-dominant pattern. Affected organs are the heart and forelimbs with upper extremity skeletal defects and congenital heart malformation. In this study we present three cases of HOS in the same family. In one of these three individuals we detected a transition of C to T (CTG-GTT, V205V) in exon 7 of the TBX5 gene. This nucleotide change causes no amino acid change and potential pathologic effects remain unknown.Key Words: Holt-Oram syndrome, Congenital heart malformation, TBX5 gene     

Computational detection of deleterious SNPs and their effect on sequence and structural level of the <Emphasis Type="Italic">VHL</Emphasis> gene

In this work we have analyzed the genetic variation that can alter the expression and the function of the VHL gene using computational methods. Of 110 single nucleotide polymorphisms (SNPs), 33 were found to be nonsynonymous (nsSNPs) and 23 SNPs were found in untranslated regions. Of the 33 nsSNPs investigated, 36.3% were found to be deleterious by both SIFT and PolyPhen servers. An untranslated region (UTR) resource tool suggested that two SNPs in the 5' UTR region and six SNPs in the 3' UTR region might change the protein expression levels. It was found by both SIFT and PolyPhen servers that a mutation from histidine to arginine at position 115 of the native protein of the VHL gene was most deleterious. A structural analysis of this mutated protein and the native protein was performed and had a root mean square deviation (RMSD) of 2.78 A. Based on this work, we propose that the nsSNP with a SNPid of rs5030812 is an important candidate for the cause of von Hippel-Lindau syndrome via the VHL gene.     

In Silico Analysis of Functional Single Nucleotide Polymorphisms in the Human TRIM22 Gene

Tripartite motif protein 22 (TRIM22) is an evolutionarily ancient protein that plays an integral role in the host innate immune response to viruses. The antiviral TRIM22 protein has been shown to inhibit the replication of a number of viruses, including HIV-1, hepatitis B, and influenza A. TRIM22 expression has also been associated with multiple sclerosis, cancer, and autoimmune disease. In this study, multiple in silico computational methods were used to identify non-synonymous or amino acid-changing SNPs (nsSNP) that are deleterious to TRIM22 structure and/or function. A sequence homology-based approach was adopted for screening nsSNPs in TRIM22, including six different in silico prediction algorithms and evolutionary conservation data from the ConSurf web server. In total, 14 high-risk nsSNPs were identified in TRIM22, most of which are located in a protein interaction module called the B30.2 domain. Additionally, 9 of the top high-risk nsSNPs altered the putative structure of TRIM22''s B30.2 domain, particularly in the surface-exposed v2 and v3 regions. These same regions are critical for retroviral restriction by the closely-related TRIM5α protein. A number of putative structural and functional residues, including several sites that undergo post-translational modification, were also identified in TRIM22. This study is the first extensive in silico analysis of the highly polymorphic TRIM22 gene and will be a valuable resource for future targeted mechanistic and population-based studies.     

Structure analysis of deleterious nsSNPs in human PALB2 protein for functional inference

Partner and Localizer of BRCA2 or PALB2 is a typical tumor suppressor protein, that responds to DNA double stranded breaks through homologous recombination repair. Heterozygous mutations in PALB2 are known to contribute to the susceptibility of breast and ovarian cancer. However, there is no comprehensive study characterizing the structural and functional impacts of SNPs located in the PALB2 gene. Therefore, it is of interest to document a comprehensive analysis of coding and non-coding SNPs located at the PALB2 loci using in silico tools. The data for 1455 non-synonymous SNPs (nsSNPs) located in the PALB2 loci were retrieved from the dbSNP database. Comprehensive characterization of the SNPs using a combination of in silico tools such as SIFT, PROVEAN, PolyPhen, PANTHER, PhD-SNP, Pmut, MutPred 2.0 and SNAP-2, identified 28 functionally important SNPs. Among these, 16 nsSNPs were further selected for structural analysis using conservation profile and protein stability. The most deleterious nsSNPs were documented within the WD40 domain of PALB2. A general outline of the structural consequences of each variant was developed using the HOPE project data. These 16 mutant structures were further modelled using SWISS Model and three most damaging mutant models (rs78179744, rs180177123 and rs45525135) were identified. The non-coding SNPs in the 3'' UTR region of the PALB2 gene were analyzed for altered miRNA target sites. The comprehensive characterization of the coding and non-coding SNPs in the PALB2 locus has provided a list of damaging SNPs with potential disease association. Further validation through genetic association study will reveal their clinical significance.     

Association of a Polymorphism in the BIRC6 Gene with Pseudoexfoliative Glaucoma

Recently an association was observed between alleles in genes of the unfolded protein response pathway and primary open angle glaucoma (POAG). The goal of the current study is to investigate the role of these two genes, protein disulphide isomerase A member 5 (PDIA5) and baculoviral IAP repeat containing 6 (BIRC6), in different forms of glaucoma. 278 patients with POAG, 132 patients with primary angle closure glaucoma (PACG) and 135 patients with pseudoexfoliative glaucoma (PEXG) were genotyped for single nucleotide polymorphisms (SNPs) rs11720822 in PDIA5 and 471 POAG, 184 PACG and 218 PEXG patients were genotyped for rs2754511 in BIRC6. Genotyping was done by allelic discrimination PCR, and genotype and allele frequencies were calculated. Logistic regression analyses were performed using R software to determine the association of these SNPs with glaucoma. The allele and genotype frequencies of rs11720822 in PDIA5 were not associated with POAG, PACG or PEXG. The TT genotype of rs2754511 in BIRC6 was found to be protective for PEXG (p = 0.05, OR 0.42 [0.22–0.81]) in the Pakistani population, but not for POAG or PACG. This study did not confirm a previously reported association of risk alleles in PDIA5 and BIRC6 with POAG, but did demonstrate a protective role of the T allele of rs2754511 in the BIRC6 gene in PEXG. This supports a role for the unfolded protein response pathway and regulation of apoptotic cell death in the pathogenesis of PEXG.     

Connecting the Dots: Potential of Data Integration to Identify Regulatory SNPs in Late-Onset Alzheimer's Disease GWAS Findings

Late-onset Alzheimer''s disease (LOAD) is a multifactorial disorder with over twenty loci associated with disease risk. Given the number of genome-wide significant variants that fall outside of coding regions, it is possible that some of these variants alter some function of gene expression rather than tagging coding variants that alter protein structure and/or function. RegulomeDB is a database that annotates regulatory functions of genetic variants. In this study, we utilized RegulomeDB to investigate potential regulatory functions of lead single nucleotide polymorphisms (SNPs) identified in five genome-wide association studies (GWAS) of risk and age-at onset (AAO) of LOAD, as well as SNPs in LD (r²≥0.80) with the lead GWAS SNPs. Of a total 614 SNPs examined, 394 returned RegulomeDB scores of 1–6. Of those 394 variants, 34 showed strong evidence of regulatory function (RegulomeDB score <3), and only 3 of them were genome-wide significant SNPs (ZCWPW1/rs1476679, CLU/rs1532278 and ABCA7/rs3764650). This study further supports the assumption that some of the non-coding GWAS SNPs are true associations rather than tagged associations and demonstrates the application of RegulomeDB to GWAS data.     

Computational and Structural Investigation of Deleterious Functional SNPs in Breast Cancer BRCA2 Gene

In this work, we have analyzed the genetic variation that can alter the expression and the function in BRCA2 gene using computational methods. Out of the total 534 SNPs, 101 were found to be non synonymous (nsSNPs). Among the 7 SNPs in the untranslated region, 3 SNPs were found in 5′ and 4 SNPs were found in 3′ un-translated regions (UTR). Of the nsSNPs 20.7% were found to be damaging by both SIFT and PolyPhen server among the 101 nsSNPs investigated. UTR resource tool suggested that 2 SNPs in the 5′ UTR region and 4 SNPs in the 3′ UTR regions might change the protein expression levels. The mutation from asparagine to isoleucine at the position 3124 of the native protein of BRCA2 gene was most deleterious by both SIFT and PolyPhen servers. A structural analysis of this mutated protein and the native protein was made which had an RMSD value of 0.301 nm. Based on this work, we proposed that this most deleterious nsSNP with an SNPid rs28897759 is an important candidate for the cause of breast cancer by BRCA2 gene.     

Genes involved in the metabolism of poly-unsaturated fatty-acids (PUFA) and risk for Crohn's disease in children & young adults

Background and Objectives

Epidemiological evidence for the role of polyunsaturated fatty-acids (PUFA) in Crohn''s disease (CD) is unclear, although the key metabolite leucotriene B4 (LTB₄) is closely linked to the inflammatory process. We hypothesized that inherited variation in key PUFA metabolic enzymes may modify susceptibility for CD.

Methods and Principal Results

A case-control design was implemented at three pediatric gastroenterology clinics in Canada. Children ≤20 yrs diagnosed with CD and controls were recruited. 19 single nucleotide polymorphisms (SNPs) across the ALOX5 (4) CYP4F3 (5) and CYP4F2 (10) genes, were genotyped. Associations between SNPs/haplotypes and CD were examined. A total of 431 cases and 507 controls were studied. The mean (±SD) age of the cases was 12.4 (±3.3) years. Most cases were male (56.4%), had ileo-colonic disease (L3±L4, 52.7%) and inflammatory behavior (B1±p, 87%) at diagnosis. One genotyped CYP4F3 SNP (rs2683037) not in Hardy-Weinberg Equilibrium was excluded. No associations with the remaining 4 CYP4F3 SNPs with CD were evident. However haplotype analysis revealed associations with a two-marker haplotype (TG) (rs3794987 & rs1290617) (p = 0.02; permuted p = 0.08). CYP4F2 SNPs, rs3093158 (OR (recessive) = 0.56, 95% CI = 0.35–0.89; p = 0.01), rs2074902 (OR (trend) = 1.26, 95% CI = 1.00–1.60; p = 0.05), and rs2108622 (OR (recessive) = 1.6, 95% CI = 1.00–2.57; p = 0.05) were significantly associated whereas rs1272 (OR (recessive) = 0.58, 95% CI = 0.30–1.13; p = 0.10) showed suggestions for associations with CD. A haplotype comprising these 4 SNPs was significantly associated (p = 0.007, permuted p = 0.02) with CD. Associations with SNP rs3780901 in the ALOX5 gene were borderline non-significant (OR (dominant) = 1.29, 95% CI = 0.99–1.67; p = 0.056). A haplotype comprising the 4 ALOX5 SNPs (TCAA, p = 0.036) was associated with CD, but did not withstand corrections for multiple comparisons (permuted p = 0.14).

Conclusions

Inherited variation in enzymes involved in the synthesis/metabolism of LTB₄ may be associated with CD. These findings implicate PUFA metabolism as a important pathway in the CD pathogenesis.