Investigating the role of single point mutations in the human proteome: a computational study

In the post genomic era with the massive production of biological data the understanding of factors affecting protein stability is one of the most important and challenging tasks for highlighting the role of mutations in relation to human maladies. The problem is at the basis of what is referred to as molecular medicine with the underlying idea that pathologies can be detailed at a molecular level. To this purpose scientific efforts focus on characterising mutations that hamper protein functions and by these affect biological processes at the basis of cell physiology. New techniques have been developed with the aim of detailing single nucleotide polymorphisms (SNPs) at large in all the human chromosomes and by this information in specific databases are exponentially increasing. Eventually mutations that can be found at the DNA level, when occurring in transcribed regions may then lead to mutated proteins and this can be a serious medical problem, largely affecting the phenotype. Bioinformatics tools are urgently needed to cope with the flood of genomic data stored in database and in order to analyse the role of SNPs at the protein level. In principle several experimental and theoretical observations are suggesting that protein stability in the solvent-protein space is responsible of the correct protein functioning. Then mutations that are found disease related during DNA analysis are often assumed to perturb protein stability as well. However so far no extensive analysis at the proteome level has investigated whether this is the case. Also computationally methods have been developed to infer whether a mutation is disease related and independently whether it affects protein stability. Therefore whether the perturbation of protein stability is related to what it is routinely referred to as a disease is still a big question mark. In this work we have tried for the first time to explore the relation among mutations at the protein level and their relevance to diseases with a large-scale computational study of the data from different databases. To this aim in the first part of the thesis for each mutation type we have derived two probabilistic indices (for 141 out of 150 possible SNPs): the perturbing index (Pp), which indicates the probability that a given mutation effects protein stability considering all the “in vitro” thermodynamic data available and the disease index (Pd), which indicates the probability of a mutation to be disease related, given all the mutations that have been clinically associated so far. We find with a robust statistics that the two indexes correlate with the exception of all the mutations that are somatic cancer related. By this each mutation of the 150 can be coded by two values that allow a direct comparison with data base information. Furthermore we also implement computational methods that starting from the protein structure is suited to predict the effect of a mutation on protein stability and find that overpasses a set of other predictors performing the same task. The predictor is based on support vector machines and takes as input protein tertiary structures. We show that the predicted data well correlate with the data from the databases. All our efforts therefore add to the SNP annotation process and more importantly found the relationship among protein stability perturbation and the human variome leading to the diseasome.

CYP17A1 polymorphisms and clinical outcome of patients with metastatic castration-resistant prostate cancer treated with abiraterone

Background. Abiraterone acetate is a potent inhibitor of cytochrome P450 17 α-hydrolase (CYP17A1) that causes a reduction in the synthesis of testosterone in the adrenal glands, testes and tumor microenvironment. Blocking androgen production, abiraterone has been shown to prolong progression-free survival (PFS) and overall survival (OS) in patients with metastatic castration-resistant prostate cancer (CRPC) previously submitted to chemotherapy. The aim of our study was to verify the role of single nucleotide polymorphisms (SNPs) in predicting clinical outcome in CRPC patients treated with abiraterone after chemotherapy. Methods. We analyzed 48 CRPC consecutive patients treated with abiraterone after at least one chemotherapeutic regimen with docetaxel. DNA was extracted from peripheral blood and genotyped for four polymorphisms in the CYP17A1 gene (rs743572, rs10883783, rs17115100, rs284849). PFS and OS survival curves were used to identify statistical associations between haplotypes and clinical outcome. Results. Forty-eight Caucasian patients with metastatic CRPC treated with abiraterone were genotyped for polymorphisms in the CYP17A1 gene. All samples were evaluable for both sequencing and TaqMan Genotyping assay. The CRPC patients treated with abiraterone had a median PFS and OS of 7.6 months (95% CI: 4.3-10.5) and 17.6 months (95% CI: 10.5-19.0), respectively Statistical analyses highlighted a difference approaching statistical significance (log-rank test p = 0.0534) between rs10883783 and PFS. Other polymorphisms were not associated with a benefit from treatment with abiraterone. Conclusions. In our case series of 48 treated patients, rs10883783 only was identified as a possible predictive marker, results showing a trend toward statistical significance. Further analysis of this polymorphism is needed in larger series of patients to confirm our findings.

Identification of SNPs in the cystic fibrosis interactome influencing pulmonary progression in cystic fibrosis

There is growing evidence that the great phenotypic variability in patients with cystic fibrosis (CF) not only depends on the genotype, but apart from a combination of environmental and stochastic factors predominantly also on modifier gene effects. It has been proposed that genes interacting with CF transmembrane conductance regulator (CFTR) and epithelial sodium channel (ENaC) are potential modifiers. Therefore, we assessed the impact of single-nucleotide polymorphisms (SNPs) of several of these interacters on CF disease outcome. SNPs that potentially alter gene function were genotyped in 95 well-characterized p.Phe508del homozygous CF patients. Linear mixed-effect model analysis was used to assess the relationship between sequence variants and the repeated measurements of lung function parameters. In total, we genotyped 72 SNPs in 10 genes. Twenty-five SNPs were used for statistical analysis, where we found strong associations for one SNP in PPP2R4 with the lung clearance index (P ≤ 0.01), the specific effective airway resistance (P ≤ 0.005) and the forced expiratory volume in 1 s (P ≤ 0.005). In addition, we identified one SNP in SNAP23 to be significantly associated with three lung function parameters as well as one SNP in PPP2R1A and three in KRT19 to show a significant influence on one lung function parameter each. Our findings indicate that direct interacters with CFTR, such as SNAP23, PPP2R4 and PPP2R1A, may modify the residual function of p.Phe508del-CFTR while variants in KRT19 may modulate the amount of p.Phe508del-CFTR at the apical membrane and consequently modify CF disease.

Influence of host interleukin-10 polymorphisms on development of traveler's diarrhea due to heat-labile enterotoxin-producing Escherichia coli in travelers from the United States who are visiting Mexico.

Up to 60% of U.S. visitors to Mexico develop traveler's diarrhea (TD), mostly due to enterotoxigenic Escherichia coli (ETEC) strains that produce heat-labile (LT) and/or heat-stable (ST) enterotoxins. Distinct single-nucleotide polymorphisms (SNPs) within the interleukin-10 (IL-10) promoter have been associated with high, intermediate, or low production of IL-10. We conducted a prospective study to investigate the association of SNPs in the IL-10 promoter and the occurrence of TD in ETEC LT-exposed travelers. Sera from U.S. travelers to Mexico collected on arrival and departure were studied for ETEC LT seroconversion by using cholera toxin as the antigen. Pyrosequencing was performed to genotype IL-10 SNPs. Stools from subjects who developed diarrhea were also studied for other enteropathogens. One hundred twenty-one of 569 (21.3%) travelers seroconverted to ETEC LT, and among them 75 (62%) developed diarrhea. Symptomatic seroconversion was more commonly seen in subjects who carried a genotype producing high levels of IL-10; it was seen in 83% of subjects with the GG genotype versus 54% of subjects with the AA genotype at IL-10 gene position -1082 (P, 0.02), in 71% of those with the CC genotype versus 33% of those with the TT genotype at position -819 (P, 0.005), and in 71% of those with the CC genotype versus 38% of those with the AA genotype at position -592 (P, 0.02). Travelers with the GCC haplotype were more likely to have symptomatic seroconversion than those with the ATA haplotype (71% versus 38%; P, 0.002). Travelers genetically predisposed to produce high levels of IL-10 were more likely to experience symptomatic ETEC TD.

Prognostic significance of xCT polymorphisms and expression in patients with advanced pancreatic cancer treated with chemotherapy

The plasma membrane xc- cystine/glutamate transporter mediates cellular uptake of cystine in exchange for intracellular glutamate and is highly expressed by pancreatic cancer cells. The xCT gene, encoding the cystine-specific xCT protein subunit of xc-, is important in regulating intracellular glutathione (GSH) levels, critical for cancer cell protection against oxidative stress, tumor growth and resistance to chemotherapeutic agents including platinum. We examined 4 single nucleotide polymorphisms (SNPs) of the xCT gene in 269 advanced pancreatic cancer patients who received first line gemcitabine with or without cisplatin or oxaliplatin. Genotyping was performed using Taqman real-time PCR assays. A statistically significant correlation was noted between the 3' untranslated region (UTR) xCT SNP rs7674870 and overall survival (OS): Median survival time (MST) was 10.9 and 13.6 months, respectively, for the TT and TC/CC genotypes (p = 0.027). Stratified analysis showed the genotype effect was significant in patients receiving gemcitabine in combination with platinum therapy (n = 145): MST was 10.5 versus 14.1 months for the TT and TC/CC genotypes, respectively (p = 0.013). The 3' UTR xCT SNP rs7674870 may correlate with OS in pancreatic cancer patients receiving gemcitabine and platinum combination therapy. Paraffin-embedded core and surgical biopsy tumor specimens from 98 patients with metastatic pancreatic adenocarcinoma were analyzed by immunohistochemistry using an xCT specific antibody. xCT protein IHC expression scores were analyzed in relation to overall survival in 86 patients and genotype in 12 patients and no statistically significant association was found between the level of xCT IHC expression score and overall survival (p = 0.514). When xCT expression was analyzed in terms of treatment response, no statistically significant associations could be determined (p = 0.908). These data suggest that polymorphic variants of xCT may have predictive value, and that the xc- transporter may represent an important target for therapy in pancreatic cancer.

IL1B and DEFB1 Polymorphisms Increase Susceptibility to Invasive Mold Infection After Solid Organ Transplantation.

BACKGROUND  Single nucleotide polymorphisms (SNPs) in immune genes have been associated with susceptibility to invasive mold infection (IMI) among hematopoietic stem cell (HSCT) but not solid organ transplant (SOT) recipients. METHODS  24 SNPs from systematically selected genes were genotyped among 1101 SOT recipients (715 kidneys, 190 liver, 102 lungs, 79 hearts, 15 other) from the Swiss Transplant Cohort Study. Association between SNPs and the endpoint were assessed by log-rank test and Cox regression models. Cytokine production upon Aspergillus stimulation was measured by ELISA in PBMCs from healthy volunteers and correlated with relevant genotypes. RESULTS  Mold colonization (N=45) and proven/probable IMI (N=26) were associated with polymorphisms in interleukin-1 beta (IL1B, rs16944; log-rank test, recessive mode, colonization P=0.001 and IMI P=0.00005), interleukin-1 receptor antagonist (IL1RN, rs419598; P=0.01 and P=0.02) and β-defensin-1 (DEFB1, rs1800972; P=0.001 and P=0.0002, respectively). The associations with IL1B and DEFB1 remained significant in a multivariate regression model (IL1B rs16944 P=0.002; DEFB1 rs1800972 P=0.01). Presence of two copies of the rare allele of rs16944 or rs419598 was associated with reduced Aspergillus-induced IL-1β and TNFα secretion by PBMCs. CONCLUSIONS  Functional polymorphisms in IL1B and DEFB1 influence susceptibility to mold infection in SOT recipients. This observation may contribute to individual risk stratification.

Polymorphisms within the canine MLPH gene are associated with dilute coat color in dogs

BACKGROUND Pinschers and other dogs with coat color dilution show a characteristic pigmentation phenotype. The fur colors are a lighter shade, e.g. silvery grey (blue) instead of black and a sandy color (Isabella fawn) instead of red or brown. In some dogs the coat color dilution is sometimes accompanied by hair loss and recurrent skin inflammation, the so called color dilution alopecia (CDA) or black hair follicular dysplasia (BHFD). In humans and mice a comparable pigmentation phenotype without any documented hair loss is caused by mutations within the melanophilin gene (MLPH). RESULTS We sequenced the canine MLPH gene and performed a mutation analysis of the MLPH exons in 6 Doberman Pinschers and 5 German Pinschers. A total of 48 sequence variations was identified within and between the breeds. Three families of dogs showed co-segregation for at least one polymorphism in an MLPH exon and the dilute phenotype. No single polymorphism was identified in the coding sequences or at splice sites that is likely to be causative for the dilute phenotype of all dogs examined. In 18 German Pinschers a mutation in exon 7 (R199H) was consistently associated with the dilute phenotype. However, as this mutation was present in homozygous state in four dogs of other breeds with wildtype pigmentation, it seems unlikely that this mutation is truly causative for coat color dilution. In Doberman Pinschers as well as in Large Munsterlanders with BHFD, a set of single nucleotide polymorphisms (SNPs) around exon 2 was identified that show a highly significant association to the dilute phenotype. CONCLUSION This study provides evidence that coat color dilution is caused by one or more mutations within or near the MLPH gene in several dog breeds. The data on polymorphisms that are strongly associated with the dilute phenotype will allow the genetic testing of Pinschers to facilitate the breeding of dogs with defined coat colors and to select against Large Munsterlanders carrying BHFD.

Genetic polymorphisms associated with the inflammatory response in bacterial meningitis.

BACKGROUND Bacterial meningitis (BM) is an infectious disease that results in high mortality and morbidity. Despite efficacious antibiotic therapy, neurological sequelae are often observed in patients after disease. Currently, the main challenge in BM treatment is to develop adjuvant therapies that reduce the occurrence of sequelae. In recent papers published by our group, we described the associations between the single nucleotide polymorphisms (SNPs) AADAT +401C > T, APEX1 Asn148Glu, OGG1 Ser326Cys and PARP1 Val762Ala and BM. In this study, we analyzed the associations between the SNPs TNF -308G > A, TNF -857C > T, IL-8 -251A > T and BM and investigated gene-gene interactions, including the SNPs that we published previously. METHODS The study was conducted with 54 BM patients and 110 healthy volunteers (as the control group). The genotypes were investigated via primer-introduced restriction analysis-polymerase chain reaction (PIRA-PCR) or polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) analysis. Allelic and genotypic frequencies were also associated with cytokine and chemokine levels, as measured with the x-MAP method, and cell counts. We analyzed gene-gene interactions among SNPs using the generalized multifactor dimensionality reduction (GMDR) method. RESULTS We did not find significant association between the SNPs TNF -857C > T and IL-8 -251A > T and the disease. However, a higher frequency of the variant allele TNF -308A was observed in the control group, associated with changes in cytokine levels compared to individuals with wild type genotypes, suggesting a possible protective role. In addition, combined inter-gene interaction analysis indicated a significant association between certain genotypes and BM, mainly involving the alleles APEX1 148Glu, IL8 -251 T and AADAT +401 T. These genotypic combinations were shown to affect cyto/chemokine levels and cell counts in CSF samples from BM patients. CONCLUSIONS In conclusion, this study revealed a significant association between genetic variability and altered inflammatory responses, involving important pathways that are activated during BM. This knowledge may be useful for a better understanding of BM pathogenesis and the development of new therapeutic approaches.

CCDC26, CDKN2BAS, RTEL1 and TERT Polymorphisms in pediatric brain tumor susceptibility.

The role of genetic polymorphisms in pediatric brain tumor (PBT) etiology is poorly understood. We hypothesized that single nucleotide polymorphisms (SNPs) identified in genome-wide association studies (GWAS) on adult glioma would also be associated with PBT risk. The study is based on the Cefalo study, a population-based multicenter case-control study. Saliva DNA from 245 cases and 489 controls, aged 7-19 years at diagnosis/reference date, was extracted and genotyped for 29 SNPs reported by GWAS to be significantly associated with risk of adult glioma. Data were analyzed using unconditional logistic regression. Stratified analyses were performed for two histological subtypes: astrocytoma alone and the other tumor types combined. The results indicated that four SNPs, CDKN2BAS rs4977756 (p = 0.036), rs1412829 (p = 0.037), rs2157719 (p = 0.018) and rs1063192 (p = 0.021), were associated with an increased susceptibility to PBTs, whereas the TERT rs2736100 was associated with a decreased risk (p = 0.018). Moreover, the stratified analyses showed a decreased risk of astrocytoma associated with RTEL1 rs6089953, rs6010620 and rs2297440 (p trend = 0.022, p trend = 0.042, p trend = 0.029, respectively) as well as an increased risk of this subtype associated with RTEL1 rs4809324 (p trend = 0.033). In addition, SNPs rs10464870 and rs891835 in CCDC26 were associated with an increased risk of non-astrocytoma tumor subtypes (p trend = 0.009, p trend = 0.007, respectively). Our findings indicate that SNPs in CDKN2BAS, TERT, RTEL1 and CCDC26 may be associated with the risk of PBTs. Therefore, we suggest that pediatric and adult brain tumors might share common genetic risk factors and similar etiological pathways.

Association of p21/p27 polymorphisms with risk of HPV16-associated oral squamous cell carcinoma

The increasing incidence of oral squamous cell carcinoma (OSCC) among young adults has been associated with sexually transmitted infection of human papillomavirus (HPV), particularly HPV16. Given the roles of p21 (WAF1/Cip1/CDKN1A) and p27 (Kip1/CDKNIB) in cell-cycle regulation and of HPV16 E6 and E7 oncoproteins in p53 degradation and pRb inactivation, the effect of HPV16 L1 seropositivity and three putatively functional single-nucleotide polymorphisms (SNPs) of p21 (p21 C70T and p21 C98A) and p27 (p27 T109G), individually and in combination, on the risk of OSCC was evaluated in a hospital-based case-control study of 327 cases and 401 cancer-free controls who were frequency-matched on age, gender and smoking status. Individuals with HPV16 L1 seropositivity had an overall 3-fold increased risk of having OSCC than those with HPV16 seronegativity. The increased risk of HPV16-associated OSCC was particularly found among younger people (aged ≤ 50 years), males, never smokers, never drinkers and oropharynx cancer patients. None of three p21 and p27 polymorphisms alone was significantly associated with risk of OSCC. Individuals with variant genotypes for both p21 polymorphisms were more likely to have OSCC than individuals with wild-type genotypes or variant genotypes for either one of the p21 polymorphisms (adjusted OR, 1.4; 95% CI, 0.9-2.1). There was a borderline significant or significant interaction between the p21 C70T, combined p21 and combined p21/p27 genotypes and HPV16 L1 seropositivity on risk of OSCC. The three studied p21 and p27 polymorphisms, individually or in combination, did not appear to have an effect on HPV16-related clinical outcomes (overall and disease-free survival and tumor recurrence). Despite the fact that the exact biological mechanism remains to be explored, these findings suggest possible involvement of p21variants, particularly the p21 C70T variant genotypes (CT/TT), in the etiology of HPV16-associated OPSCC. Further large and functional studies are required to validate the findings.^

Role of obesity, weight gain, physical activity, and polymorphisms in the mTOR pathway and the risk of renal cell carcinoma

The interplay between obesity, physical activity, weight gain and genetic variants in mTOR pathway have not been studied in renal cell carcinoma (RCC). We examined the associations between obesity, weight gain, physical activity and RCC risk. We also analyzed whether genetic variants in the mTOR pathway could modify the association. Incident renal cell carcinoma cases and healthy controls were recruited from the University of Texas MD Anderson Cancer Center in Houston, Texas. Cases and controls were frequency-matched by age (±5 years), ethnicity, sex, and county of residence. Epidemiologic data were collected via in-person interview. A total of 577 cases and 593 healthy controls (all white) were included. One hundred ninety-two (192) SNPs from 22 genes were available and their genotyping data were extracted from previous genome-wide association studies. Logistic regression and regression spline were performed to obtain odds ratios. Obesity at age 20, 40, and 3 years prior to diagnosis/recruitment, and moderate and large weight gain from age 20 to 40 were each significantly associated with increased RCC risk. Low physical activity was associated with a 4.08-fold (95% CI: 2.92-5.70) increased risk. Five single nucleotide polymorphisms (SNPs) were significantly associated with RCC risk and their cumulative effect increased the risk by up to 72% (95% CI: 1.20-2.46). Strata specific effects for weight change and genotyping cumulative groups were observed. However, no interaction was suggested by our study. In conclusion, energy balance related risk factors and genetic variants in the mTOR pathway may jointly influence susceptibility to RCC. ^

Polymorphisms in the trace amine receptor 4 (TRAR4) gene on chromosome 6q23.2 are associated with susceptibility to schizophrenia

Several linkage studies across multiple population groups provide convergent support for a susceptibility locus for schizophrenia - and, more recently, for bipolar disorder - on chromosome 6q13-q26. We genotyped 192 European-ancestry and African American (AA) pedigrees with schizophrenia from samples that previously showed linkage evidence to 6q13-q26, focusing on the MOXD1-STX7-TRARs gene cluster at 6q23.2, which contains a number of prime candidate genes for schizophrenia. Thirty-one screening single-nucleotide polymorphisms (SNPs) were selected, providing a minimum coverage of at least 1 SNP/20 kb. The association observed with rs4305745 (P = .0014) within the TRAR4 (trace amine receptor 4) gene remained significant after correction for multiple testing. Evidence for association was proportionally stronger in the smaller AA sample. We performed database searches and sequenced genomic DNA in a 30-proband subsample to obtain a high-density map of 23 SNPs spanning 21.6 kb of this gene. Single-SNP analyses and also haplotype analyses revealed that rs4305745 and/or two other polymorphisms in perfect linkage disequilibrium (LD) with rs4305745 appear to be the most likely variants underlying the association of the TRAR4 region with schizophrenia. Comparative genomic analyses further revealed that rs4305745 and/or the associated polymorphisms in complete LD with rs4305745 could potentially affect gene expression. Moreover, RT-PCR studies of various human tissues, including brain, confirm that TRAR4 is preferentially expressed in those brain regions that have been implicated in the pathophysiology of schizophrenia. These data provide strong preliminary evidence that TRAR4 is a candidate gene for schizophrenia; replication is currently being attempted in additional clinical samples.

Association between polymorphisms in the progesterone receptor gene and endometriosis

The progesterone receptor (PR) is a candidate gene for the development of endometriosis, a complex disease with strong hormonal features, common in women of reproductive age. We typed the 306 base pair Alu insertion (AluIns) polymorphism in intron G of PR in 101 individuals, estimated linkage disequilibrium (LD) between five single-nucleotide polymorphisms (SNPs) across the PR locus in 980 Australian triads (endometriosis case and two parents) and used transmission disequilibrium testing (TDT) for association with endometriosis. The five SNPs showed strong pairwise LD, and the AluIns was highly correlated with proximal SNPs rs1042839 ({Delta}2 = 0.877, D9 = 1.00, P < 0.0001) and rs500760 ({Delta}2 = 0.438, D9 = 0.942, P < 0.0001). TDT showed weak evidence of allelic association between endometriosis and rs500760 (P = 0.027) but not in the expected direction. We identified a common susceptibility haplotype GGGCA across the five SNPs (P = 0.0167) in the whole sample, but likelihood ratio testing of haplotype transmission and non-transmission of the AluIns and flanking SNPs showed no significant pattern. Further, analysis of our results pooled with those from two previous studies suggested that neither the T2 allele of the AluIns nor the T1/T2 genotype was associated with endometriosis.

Polymorphisms in the vitamin D receptor and their associations with risk of schizophrenia and selected anthropometric measures

The association between vitamin D levels and skeletal growth has long been recognized. However, exposure to low levels of vitamin D during early life is also known to alter brain development, and is a candidate risk factor for schizophrenia. This study examines the association between four polymorphisms in the vitamin D receptor (VDR) and 1) risk of schizophrenia, and 2) three anthropometric variables (height, head size, and head shape). Four single-nucleotide polymorphisms (SNPs; rs10735810/FokI, rsl544410/BsmI, rs7975232/ApaI, and rs731236/TaqI) in the VDR gene were genotyped in 179 individuals with schizophrenia and 189 healthy controls. No significant associations were detected between any of the four VDR SNPs and risk of schizophrenia. Patients were slightly but significantly shorter compared to controls. Of the four SNPs, only rs10735810/FokI was associated with any of the anthropometric measures: the M4 isoform of this SNP was significantly associated with larger head size (P = 0.002). In light of the evidence demonstrating a role for vitamin D during brain development, the association between polymorphisms in VDR and brain development warrants closer scrutiny.

A Test of Association between Qualitative Trait and a Set of SNPs

2000 Mathematics Subject Classification: 62P10, 92D10, 92D30, 62F03