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.

Detailed analyses of mutations affecting structural formation of neuromuscular synapses in Drosophila melanogaster

The establishment of appropriate synapses between neurons and their target cells is an essential requirement for the formation of functional neuronal circuits. However, there is very little insight into the mechanisms underlying de novo formation of synapses and synaptic terminals. To identify novel genes involved in signalling or structural aspects of these processes I capitalised on possibilities provided by the model organism Drosophila. Thus, I contributed to a screen of a collection of third chromosomal mutations (Salzberg et al., 1997, Genetics 147, 1723ff.) selecting those mutant strains displaying structural defects of Drosophila neuromuscular junctions (NMJ). Carrying out genetic mapping experiments, I could assign 7 genes to interesting candidate mutations. All 7 mutations selected in this process cause size alterations of the embryonic NMJ, and one shows additional disturbances in the distribution of synaptic markers. 4 of these turned out to be transcription factors, not falling into the remit of this project. Only for one of these, the neuronal transcription factor Castor, I could show that its overgrown mutant NMJ phenotype is due to an increase in the number of motorneurons. The remaining genes encode a potential nitrophenylphosphatase, the translation initiation factor eIF4AIII, and a novel protein Waharan. Unfortunately, the nitophenylphosphatase gene was identified too late to carry out functional studies in the context of this project, but potential roles are discussed. eIF4AIII promotes NMJ size tempting to speculate that local translation at the NMJ is affected. I found that the synaptic scaffolding molecule Discs large (Dlg; orthologue of PSD95) is upregulated at eIF4AIII mutant NMJs. Targeted upregulation of Dlg can not mimic the eIF4AIII mutant phenotype, but dlg mutations suppress it. Therefore, Dlg function is required but not sufficient in this context. My findings are discussed in detail, pointing out future directions. The main focus of this work is the completely novel gene waharan (wah), an orthologue of the human gene KIAA1267 encoding a big brain protein of likewise unknown structure and function. My studies show that mutations or RNAi knock-down of wah cause NMJ overgrowth and reveal additional crucial roles in the patterning of wing imginal discs. RNAi studies suggest Wah to be required pre- and postsynaptically at NMJs and, consistently, wah is transcribed in the nervous system and muscles. Anti-Wah antisera were produced but could no longer be tested here, but preliminary studies with newly generated HA-targeted constructs suggest that Wah localises at NMJs and in neuronal nuclei. In silico analyses predict Wah to be structurally related to the Rad23-family of proteins, likely to target ubiquitinated proteins to the proteasome for degradation (Chen et al., 2002, Mol Cell Biol 22, 4902ff.) . In agreement with this prediction, poly-ubiquitinated proteins were found to accumulate in the absence of wah function, and wah-like mutant phenotypes were induced in NMJs and wing discs by knocking down proteasome function. My analysis further revealed that poly-ubiquitinated proteins are reduced in nuclei of wah mutant neurons and muscles, suggesting that Wah may play additional roles in ubiquitin-mediated nuclear import. Taken together, this study has uncovered a number of interesting candidate genes required for the de novo formation of Drosophila NMJs. 3 of these genes fell into the focus of this project. As discussed in detail, discovery of these genes and insights gained into their function have high potential to be translatable into vertebrate systems.

High sensitivity analysis of BRAF mutations in neoplastic and non-neoplastic thyroid lesions

The clonal distribution of BRAFV600E in papillary thyroid carcinoma (PTC) has been recently debated. No information is currently available about precursor lesions of PTCs. My first aim was to establish whether the BRAFV600E mutation occurs as a subclonal event in PTCs. My second aim was to screen BRAF mutations in histologically benign tissue of cases with BRAFV600E or BRAFwt PTCs in order to identify putative precursor lesions of PTCs. Highly sensitive semi-quantitative methods were used: Allele Specific LNA quantitative PCR (ASLNAqPCR) and 454 Next-Generation Sequencing (NGS). For the first aim 155 consecutive formalin-fixed and paraffin-embedded (FFPE) specimens of PTCs were analyzed. The percentage of mutated cells obtained was normalized to the estimated number of neoplastic cells. Three groups of tumors were identified: a first had a percentage of BRAF mutated neoplastic cells > 80%; a second group showed a number of BRAF mutated neoplastic cells < 30%; a third group had a distribution of BRAFV600E between 30-80%. The large presence of BRAFV600E mutated neoplastic cell sub-populations suggests that BRAFV600E may be acquired early during tumorigenesis: therefore, BRAFV600E can be heterogeneously distributed in PTC. For the second aim, two groups were studied: one consisted of 20 cases with BRAFV600E mutated PTC, the other of 9 BRAFwt PTCs. Seventy-five and 23 histologically benign FFPE thyroid specimens were analyzed from the BRAFV600E mutated and BRAFwt PTC groups, respectively. The screening of BRAF mutations identified BRAFV600E in “atypical” cell foci from both groups of patients. “Unusual” BRAF substitutions were observed in histologically benign thyroid associated with BRAFV600E PTCs. These mutations were very uncommon in the group with BRAFwt PTCs and in BRAFV600E PTCs. Therefore, lesions carrying BRAF mutations may represent “abortive” attempts at cancer development: only BRAFV600E boosts neoplastic transformation to PTC. BRAFV600E mutated “atypical foci” may represent precursor lesions of BRAFV600E mutated PTCs.

Next-Generation Sequencing-Based Mutations Scanning Strategy of the BCR-ABL Kinase Domain in Patients with PhiladelPhia-Chromosome Positive Leukemias Treated with Tyrosine Kinase Inhibitors

In chronic myeloid leukemia and Philadelphia-positive acute lymphoblastic leukemia patients resistant to tyrosine kinase inhibitors (TKIs), BCR-ABL kinase domain mutation status is an essential component of the therapeutic decision algorithm. The recent development of Ultra-Deep Sequencing approach (UDS) has opened the way to a more accurate characterization of the mutant clones surviving TKIs conjugating assay sensitivity and throughput. We decided to set-up and validated an UDS-based for BCR-ABL KD mutation screening in order to i) resolve qualitatively and quantitatively the complexity and the clonal structure of mutated populations surviving TKIs, ii) study the dynamic of expansion of mutated clones in relation to TKIs therapy, iii) assess whether UDS may allow more sensitive detection of emerging clones, harboring critical 2GTKIs-resistant mutations predicting for an impending relapse, earlier than SS. UDS was performed on a Roche GS Junior instrument, according to an amplicon sequencing design and protocol set up and validated in the framework of the IRON-II (Interlaboratory Robustness of Next-Generation Sequencing) International consortium.Samples from CML and Ph+ ALL patients who had developed resistance to one or multiple TKIs and collected at regular time-points during treatment were selected for this study. Our results indicate the technical feasibility, accuracy and robustness of our UDS-based BCR-ABL KD mutation screening approach. UDS was found to provide a more accurate picture of BCR-ABL KD mutation status, both in terms of presence/absence of mutations and in terms of clonal complexity and showed that BCR-ABL KD mutations detected by SS are only the “tip of iceberg”. In addition UDS may reliably pick 2GTKIs-resistant mutations earlier than SS in a significantly greater proportion of patients.The enhanced sensitivity as well as the possibility to identify low level mutations point the UDS-based approach as an ideal alternative to conventional sequencing for BCR-ABL KD mutation screening in TKIs-resistant Ph+ leukemia patients

Analysis of Copy Number Variants identifies new candidate genes for Autism Spectrum Disorder and Intellectual Disability

Autism spectrum disorder (ASD) and Intellectual Disability (ID) are complex neuropsychiatric disorders characterized by extensive clinical and genetic heterogeneity and with overlapping risk factors. The aim of my project was to further investigate the role of Copy Numbers Variants (CNVs), identified through genome-wide studies performed by the Autism Geome Project (AGP) and the CHERISH consortium in large cohorts of ASD and ID cases, respectively. Specifically, I focused on four rare genic CNVs, selected on the basis of their impact on interesting ASD/ID candidate genes: a) a compound heterozygous deletion involving CTNNA3, predicted to cause the lack of functional protein; b) a 15q13.3 duplication containing CHRNA7; c) a 2q31.1 microdeletion encompassing KLHL23, SSB and METTL5; d) Lastly, I investigated the putative imprinting regulation of the CADPS2 gene, disrupted by a maternal deletion in two siblings with ASD and ID. This study provides further evidence for the role of CTNNA3, CHRNA7, KLHL23 and CADPS2 as ASD and/or ID susceptibility genes, and highlights that rare genetic variation contributes to disease risk in different ways: some rare mutations, such as those impacting CTNNA3, act in a recessive mode of inheritance, while other CNVs, such as those occurring in the 15q13.3 region, are implicated in multiple developmental and/or neurological disorders possibly interacting with other susceptibility variants elsewhere in the genome. On the other hand, the discovery of a tissue-specific monoallelic expression for the CADPS2 gene, implicates the involvement of epigenetic regulatory mechanisms as risk factors conferring susceptibility to ASD/ID.

Growth hormone (GH)-releasing hormone increases the expression of the dominant-negative GH isoform in cases of isolated GH deficiency due to GH splice-site mutations

An autosomal dominant form of isolated GH deficiency (IGHD II) can result from heterozygous splice site mutations that weaken recognition of exon 3 leading to aberrant splicing of GH-1 transcripts and production of a dominant-negative 17.5-kDa GH isoform. Previous studies suggested that the extent of missplicing varies with different mutations and the level of GH expression and/or secretion. To study this, wt-hGH and/or different hGH-splice site mutants (GH-IVS+2, GH-IVS+6, GH-ISE+28) were transfected in rat pituitary cells expressing human GHRH receptor (GC-GHRHR). Upon GHRH stimulation, GC-GHRHR cells coexpressing wt-hGH and each of the mutants displayed reduced hGH secretion and intracellular GH content when compared with cells expressing only wt-hGH, confirming the dominant-negative effect of 17.5-kDa isoform on the secretion of 22-kDa GH. Furthermore, increased amount of 17.5-kDa isoform produced after GHRH stimulation in cells expressing GH-splice site mutants reduced production of endogenous rat GH, which was not observed after GHRH-induced increase in wt-hGH. In conclusion, our results support the hypothesis that after GHRH stimulation, the severity of IGHD II depends on the position of splice site mutation leading to the production of increasing amounts of 17.5-kDa protein, which reduces the storage and secretion of wt-GH in the most severely affected cases. Due to the absence of GH and IGF-I-negative feedback in IGHD II, a chronic up-regulation of GHRH would lead to an increased stimulatory drive to somatotrophs to produce more 17.5-kDa GH from the severest mutant alleles, thereby accelerating autodestruction of somatotrophs in a vicious cycle.

Mutations affecting disulphide bonds contribute to a fairly common prevalence of F13B gene defects: results of a genetic study in 14 families with factor XIII B deficiency

Severe factor XIII (FXIII) deficiency is a rare autosomal recessive coagulation disorder affecting one in two million individuals. The aim of the present study was to screen for and analyse F13B gene defects in the German population. A total of 150 patients presenting with suspected FXIII deficiency and one patient with severe (homozygous) FXIII deficiency were screened for mutations in F13A and F13B genes. Twenty-five individuals presented with detectable heterozygous mutations, 12 of them in the F13A gene and 13 of them in the F13B gene. We report on the genotype-phenotype correlations of the individuals showing defects in the F13B gene. Direct sequencing revealed 12 unique mutations including seven missense mutations (Cys5Arg, Ile81Asn, Leu116Phe, Val217Ile, Cys316Phe, Val401Glu, Pro428Ser), two splice site mutations (IVS2-1G>C, IVS3-1G>C), two insertions (c.1155_1158dupACTT, c.1959insT) and one in-frame deletion (c.471-473delATT). Two of the missense mutations (Cys5Arg, Cys316Phe) eliminated disulphide bonds (Cys5-Cys56, Cys316-Cys358). Another three missense mutations, (Leu116Phe, Val401Glu, Pro428Ser) were located proximal to other cysteine disulphide bonds, therefore indicating that the region in and around these disulphide bonds is prone to functionally relevant mutations in the FXIII-B subunit. The present study reports on a fairly common prevalence of F13B gene defects in the German population. The regions in and around the cysteine disulphide bonds in the FXIII-B protein may be regions prone to frequent mutations.

Identification of eight novel coagulation factor XIII subunit A mutations: implied consequences for structure and function

Severe hereditary coagulation factor XIII deficiency is a rare homozygous bleeding disorder affecting one person in every two million individuals. In contrast, heterozygous factor XIII deficiency is more common, but usually not associated with severe hemorrhage such as intracranial bleeding or hemarthrosis. In most cases, the disease is caused by F13A gene mutations. Causative mutations associated with the F13B gene are rarer.

Brain metabolite composition in relation to cognitive function and dystrophin mutations in boys with Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is a hereditary X-linked recessive disorder affecting the synthesis of dystrophin, a protein essential for structural stability in muscle. Dystrophin also occurs in the central nervous system, particularly in the neocortex, hippocampus and cerebellum. Quantitative metabolic analysis by localized (1) H MRS was performed in the cerebellum (12 patients and 15 controls) and a temporo-parietal location (eight patients and 15 controls) in patients with DMD and healthy controls to investigate possible metabolic differences. In addition, the site of individual mutations on the dystrophin gene was analyzed and neuropsychological cognitive functions were examined. Cognitive deficits in the patient group were found in line with earlier investigations, mainly concerning verbal short-term memory, visuo-spatial long-term memory and verbal fluency, but also the full-scale IQ. Causal mutations were identified in all patients with DMD. Quantitative MRS showed consistent choline deficits, in both cerebellar white matter and temporo-parietal cortex, as well as small, but significant, metabolic abnormalities for glutamate and total N-acetyl compounds in the temporo-parietal region. Compartment water analysis did not reveal any abnormalities. In healthy subjects, choline levels were age related in the cerebellum. The choline deficit contrasts with earlier findings in DMD, where a surplus of choline was postulated for the cerebellum. In patients, total N-acetyl compounds in the temporo-parietal region were related to verbal IQ and verbal short-term memory. However, choline, the putative main metabolic abnormality, was not found to be associated with cognitive deficits. Furthermore, in contrast with the cognitive performance, the metabolic brain composition did not depend significantly on whether or not gene mutations concerned the expression of the dystrophin isoform Dp140, leading to the conclusion that the effect of the missing Dp140 isoform on cognitive performance is not mediated through the observed metabolite composition, or is caused by local effects beyond the resolution accessible to MRS investigations.

Neuropsychological impairments and the impact of dystrophin mutations on general cognitive functioning of patients with Duchenne muscular dystrophy

Mutations in the dystrophin gene have long been recognised as a cause of mental retardation. However, for reasons that are unclear, some boys with dystrophin mutations do not show general cognitive deficits. To investigate the relationship between dystrophin mutations and cognition, the general intellectual abilities of a group of 25 boys with genetically confirmed Duchenne muscular dystrophy were evaluated. Furthermore, a subgroup underwent additional detailed neuropsychological assessment. The results showed a mean full scale intelligence quotient (IQ) of 88 (standard deviation 24). Patients performed very poorly on various neuropsychological tests, including arithmetics, digit span tests and verbal fluency. No simple relationship between dystrophin mutations and cognitive functioning could be detected. However, our analysis revealed that patients who lack the dystrophin isoform Dp140 have significantly greater cognitive problems.

Desmosomal gene analysis in arrhythmogenic right ventricular dysplasia/cardiomyopathy: spectrum of mutations and clinical impact in practice

Five desmosomal genes have been recently implicated in arrhythmogenic right ventricular dysplasia/cardiomyopathy (ARVD/C) but the clinical impact of genetics remains poorly understood. We wanted to address the potential impact of genotyping.

Uromodulin is expressed in renal primary cilia and UMOD mutations result in decreased ciliary uromodulin expression

Uromodulin (UMOD) mutations are responsible for three autosomal dominant tubulo-interstitial nephropathies including medullary cystic kidney disease type 2 (MCKD2), familial juvenile hyperuricemic nephropathy and glomerulocystic kidney disease. Symptoms include renal salt wasting, hyperuricemia, gout, hypertension and end-stage renal disease. MCKD is part of the 'nephronophthisis-MCKD complex', a group of cystic kidney diseases. Both disorders have an indistinguishable histology and renal cysts are observed in either. For most genes mutated in cystic kidney disease, their proteins are expressed in the primary cilia/basal body complex. We identified seven novel UMOD mutations and were interested if UMOD protein was expressed in the primary renal cilia of human renal biopsies and if mutant UMOD would show a different expression pattern compared with that seen in control individuals. We demonstrate that UMOD is expressed in the primary cilia of renal tubules, using immunofluorescent studies in human kidney biopsy samples. The number of UMOD-positive primary cilia in UMOD patients is significantly decreased when compared with control samples. Additional immunofluorescence studies confirm ciliary expression of UMOD in cell culture. Ciliary expression of UMOD is also confirmed by electron microscopy. UMOD localization at the mitotic spindle poles and colocalization with other ciliary proteins such as nephrocystin-1 and kinesin family member 3A is demonstrated. Our data add UMOD to the group of proteins expressed in primary cilia, where mutations of the gene lead to cystic kidney disease.

Reduction in hepatic drug metabolizing CYP3A4 activities caused by P450 oxidoreductase mutations identified in patients with disordered steroid metabolism

Cytochrome P450 3A4 (CYP3A4), the major P450 present in human liver metabolizes approximately half the drugs in clinical use and requires electrons supplied from NADPH through NADPH-P450 reductase (POR, CPR). Mutations in human POR cause a rare form of congenital adrenal hyperplasia from diminished activities of steroid metabolizing P450s. In this study we examined the effect of mutations in POR on CYP3A4 activity. We used purified preparations of wild type and mutant human POR and in vitro reconstitution with purified CYP3A4 to perform kinetic studies. We are reporting that mutations in POR identified in patients with disordered steroidogenesis/Antley-Bixler syndrome (ABS) may reduce CYP3A4 activity, potentially affecting drug metabolism in individuals carrying mutant POR alleles. POR mutants Y181D, A457H, Y459H, V492E and R616X had more than 99% loss of CYP3A4 activity, while POR mutations A287P, C569Y and V608F lost 60-85% activity. Loss of CYP3A4 activity may result in increased risk of drug toxicities and adverse drug reactions in patients with POR mutations.

Altered heme catabolism by heme oxygenase-1 caused by mutations in human NADPH cytochrome P450 reductase

Human heme oxygenase-1 (HO-1) carries out heme catabolism supported by electrons supplied from the NADPH through NADPH P450 reductase (POR, CPR). Previously we have shown that mutations in human POR cause a rare form of congenital adrenal hyperplasia. In this study, we have evaluated the effects of mutations in POR on HO-1 activity. We used purified preparations of wild type and mutant human POR and in vitro reconstitution with purified HO-1 to measure heme degradation in a coupled assay using biliverdin reductase. Here we show that mutations in POR found in patients may reduce HO-1 activity, potentially influencing heme catabolism in individuals carrying mutant POR alleles. POR mutants Y181D, A457H, Y459H, V492E and R616X had total loss of HO-1 activity, while POR mutations A287P, C569Y and V608F lost 50-70% activity. The POR variants P228L, R316W and G413S, A503V and G504R identified as polymorphs had close to WT activity. Loss of HO-1 activity may result in increased oxidative neurotoxicity, anemia, growth retardation and iron deposition. Further examination of patients affected with POR deficiency will be required to assess the metabolic effects of reduced HO-1 activity in affected individuals.

Risk profiles and penetrance estimations in multiple endocrine neoplasia type 2A caused by germline RET mutations located in exon 10

Multiple endocrine neoplasia type 2 is characterized by germline mutations in RET. For exon 10, comprehensive molecular and corresponding phenotypic data are scarce. The International RET Exon 10 Consortium, comprising 27 centers from 15 countries, analyzed patients with RET exon 10 mutations for clinical-risk profiles. Presentation, age-dependent penetrance, and stage at presentation of medullary thyroid carcinoma (MTC), pheochromocytoma, and hyperparathyroidism were studied. A total of 340 subjects from 103 families, age 4-86, were registered. There were 21 distinct single nucleotide germline mutations located in codons 609 (45 subjects), 611 (50), 618 (94), and 620 (151). MTC was present in 263 registrants, pheochromocytoma in 54, and hyperparathyroidism in 8 subjects. Of the patients with MTC, 53% were detected when asymptomatic, and among those with pheochromocytoma, 54%. Penetrance for MTC was 4% by age 10, 25% by 25, and 80% by 50. Codon-associated penetrance by age 50 ranged from 60% (codon 611) to 86% (620). More advanced stage and increasing risk of metastases correlated with mutation in codon position (609?620) near the juxtamembrane domain. Our data provide rigorous bases for timing of premorbid diagnosis and personalized treatment/prophylactic procedure decisions depending on specific RET exon 10 codons affected.