Josep Lluís Gelpí

Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia

Chronic lymphocytic leukaemia (CLL), the most frequent leukaemia in adults in Western countries, is a heterogeneous disease with variable clinical presentation and evolution. Two major molecular subtypes can be distinguished, characterized respectively by a high or low number of somatic hypermutations in the variable region of immunoglobulin genes. The molecular changes leading to the pathogenesis of the disease are still poorly understood. Here we performed whole-genome sequencing of four cases of CLL and identified 46 somatic mutations that potentially affect gene function. Further analysis of these mutations in 363 patients with CLL identified four genes that are recurrently mutated: notch 1 (NOTCH1), exportin 1 (XPO1), myeloid differentiation primary response gene 88 (MYD88) and kelch-like 6 (KLHL6). Mutations in MYD88 and KLHL6 are predominant in cases of CLL with mutated immunoglobulin genes, whereas NOTCH1 and XPO1 mutations are mainly detected in patients with unmutated immunoglobulins. The patterns of somatic mutation, supported by functional and clinical analyses, strongly indicate that the recurrent NOTCH1, MYD88 and XPO1 mutations are oncogenic changes that contribute to the clinical evolution of the disease. To our knowledge, this is the first comprehensive analysis of CLL combining whole-genome sequencing with clinical characteristics and clinical outcomes. It highlights the usefulness of this approach for the identification of clinically relevant mutations in cancer.

PMUT: a web-based tool for the annotation of pathological mutations on proteins

PMUT allows the fast and accurate prediction (approximately 80% success rate in humans) of the pathological character of single point amino acidic mutations based on the use of neural networks. The program also allows the fast scanning of mutational hot spots, which are obtained by three procedures: (1) alanine scanning, (2) massive mutation and (3) genetically accessible mutations. A graphical interface for Protein Data Bank (PDB) structures, when available, and a database containing hot spot profiles for all non-redundant PDB structures are also accessible from the PMUT server.

Epigenomic analysis detects widespread gene-body DNA hypomethylation in chronic lymphocytic leukemia.

We have extensively characterized the DNA methylomes of 139 patients with chronic lymphocytic leukemia (CLL) with mutated or unmutated IGHV and of several mature B-cell subpopulations through the use of whole-genome bisulfite sequencing and high-density microarrays. The two molecular subtypes of CLL have differing DNA methylomes that seem to represent epigenetic imprints from distinct normal B-cell subpopulations. DNA hypomethylation in the gene body, targeting mostly enhancer sites, was the most frequent difference between naive and memory B cells and between the two molecular subtypes of CLL and normal B cells. Although DNA methylation and gene expression were poorly correlated, we identified gene-body CpG dinucleotides whose methylation was positively or negatively associated with expression. We have also recognized a DNA methylation signature that distinguishes new clinico-biological subtypes of CLL. We propose an epigenomic scenario in which differential methylation in the gene body may have functional and clinical implications in leukemogenesis.

Non-coding recurrent mutations in chronic lymphocytic leukaemia.

Chronic lymphocytic leukaemia (CLL) is a frequent disease in which the genetic alterations determining the clinicobiological behaviour are not fully understood. Here we describe a comprehensive evaluation of the genomic landscape of 452 CLL cases and 54 patients with monoclonal B-lymphocytosis, a precursor disorder. We extend the number of CLL driver alterations, including changes in ZNF292, ZMYM3, ARID1A and PTPN11. We also identify novel recurrent mutations in non-coding regions, including the 3′ region of NOTCH1, which cause aberrant splicing events, increase NOTCH1 activity and result in a more aggressive disease. In addition, mutations in an enhancer located on chromosome 9p13 result in reduced expression of the B-cell-specific transcription factor PAX5. The accumulative number of driver alterations (0 to ≥4) discriminated between patients with differences in clinical behaviour. This study provides an integrated portrait of the CLL genomic landscape, identifies new recurrent driver mutations of the disease, and suggests clinical interventions that may improve the management of this neoplasia.

Evidence for Transcript Networks Composed of Chimeric RNAs in Human Cells

The classic organization of a gene structure has followed the Jacob and Monod bacterial gene model proposed more than 50 years ago. Since then, empirical determinations of the complexity of the transcriptomes found in yeast to human has blurred the definition and physical boundaries of genes. Using multiple analysis approaches we have characterized individual gene boundaries mapping on human chromosomes 21 and 22. Analyses of the locations of the 5′ and 3′ transcriptional termini of 492 protein coding genes revealed that for 85% of these genes the boundaries extend beyond the current annotated termini, most often connecting with exons of transcripts from other well annotated genes. The biological and evolutionary importance of these chimeric transcripts is underscored by (1) the non-random interconnections of genes involved, (2) the greater phylogenetic depth of the genes involved in many chimeric interactions, (3) the coordination of the expression of connected genes and (4) the close in vivo and three dimensional proximity of the genomic regions being transcribed and contributing to parts of the chimeric RNAs. The non-random nature of the connection of the genes involved suggest that chimeric transcripts should not be studied in isolation, but together, as an RNA network.

A consensus view of protein dynamics

The dynamics of proteins in aqueous solution has been investigated through a massive approach based on “state of the art” molecular dynamics simulations performed for all protein metafolds using the four most popular force fields (OPLS, CHARMM, AMBER, and GROMOS). A detailed analysis of the massive database of trajectories (>1.5 terabytes of data obtained using ≈50 years of CPU) allowed us to obtain a robust-consensus picture of protein dynamics in aqueous solution.

Obligatory Amino Acid Exchange via Systems bo,+-like and y+L-like A TERTIARY ACTIVE TRANSPORT MECHANISM FOR RENAL REABSORPTION OF CYSTINE AND DIBASIC AMINO ACIDS

Mutations in the rBAT gene cause type I cystinuria, a common inherited aminoaciduria of cystine and dibasic amino acids due to their defective renal and intestinal reabsorption (Calonge, M. J., Gasparini, P., Chillarón, J., Chillón, M., Gallucci, M., Rousaud, F., Zelante, L., Testar, X., Dallapiccola, B., Di Silverio, F., Barceló, P., Estivill, X., Zorzano, A., Nunes, V., and Palacín, M. (1994) Nat. Genet. 6, 420-426; Calonge, M. J., Volipini, V., Bisceglia, L., Rousaud, F., De Sanctis, L., Beccia, E., Zelante, L., Testar, X., Zorzano, A., Estivill, X., Gasparini, P., Nunes, V., and Palacín, M. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 9667-9671). One important question that remains to be clarified is how the apparently non-concentrative system bo,+-like, associated with rBAT expression, participates in the active renal reabsorption of these amino acids. Several studies have demonstrated exchange of amino acids induced by rBAT in Xenopus oocytes. Here we offer evidence that system bo,+-like is an obligatory amino acid exchanger in oocytes and in the “renal proximal tubular” cell line OK. System bo,+-like showed a 1:1 stoichiometry of exchange, and the hetero-exchange dibasic (inward) with neutral (outward) amino acids were favored in oocytes. Obligatory exchange of amino acids via system bo,+-like fully explained the amino acid-induced current in rBAT-injected oocytes. Exchange via system bo,+-like is coupled enough to ensure a specific accumulation of substrates until the complete replacement of the internal oocyte substrates. Due to structural and functional analogies of the cell surface antigen 4F2hc to rBAT, we tested for amino acid exchange via system y+L-like. 4F2hc-injected oocytes accumulated substrates to a level higher than CAT1-injected oocytes (i.e. oocytes expressing system y+) and showed exchange of amino acids with the substrate specificity of system y+L and L-leucine-induced outward currents in the absence of extracellular sodium. In contrast to L-arginine, system y+L-like did not mediate measurable L-leucine efflux from the oocyte. We propose a role of systems bo,+-like and y+L-like in the renal reabsorption of cystine and dibasic amino acids that is based on their active tertiary transport mechanism and on the apical and basolateral localization of rBAT and 4F2hc, respectively, in the epithelial cells of the proximal tubule of the nephron.

PARMBSC1: A refined force-field for DNA simulations

We present parmbsc1, a force field for DNA atomistic simulation, which has been parameterized from high-level quantum mechanical data and tested for nearly 100 systems (representing a total simulation time of ∼140 μs) covering most of DNA structural space. Parmbsc1 provides high-quality results in diverse systems. Parameters and trajectories are available at http://mmb.irbbarcelona.org/ParmBSC1/.

Glutaryl-CoA Dehydrogenase Deficiency in Spain: Evidence of Two Groups of Patients, Genetically, and Biochemically Distinct

Glutaryl-CoA dehydrogenase (GCDH) deficiency causes glutaric aciduria type I (GA I), an inborn error of metabolism that is characterized clinically by dystonia and dyskinesia and pathologically by neural degeneration of the caudate and putamen. Studies of metabolite excretion allowed us to categorize 43 GA I Spanish patients into two groups: group 1 (26 patients), those presenting with high excretion of both glutarate and 3-hydroxyglutarate, and group 2 (17 patients) , those who might not be detected by routine urine organic acid analysis because glutarate might be normal and 3-hydroxyglutarate only slightly higher than controls. Single-strand conformation polymorphism (SSCP) screening and sequence analysis of the 11 exons and the corresponding intron boundaries of the GCDH gene allowed us to identify 13 novel and 10 previously described mutations. The most frequent mutations in group 1 were A293T and R402W with an allele frequency of 30% and 28%, respectively. These two mutations were also found in group 2, but always in heterozygosity, in particular in combination with mutations V400M or R227P. Interestingly, mutations V400M and R227P were only found in group 2, and at least one of these mutations was found in 11 of 15 unrelated alleles, accounting together for 53% of the mutant alleles in group 2. Therefore, it seems clear that two genetically and biochemically distinct groups of patients exist. The severity of the clinical phenotype seems to be closely linked to the development of encephalopathic crises rather than to residual enzyme activity or genotype. Comparison of GCDH protein with other acyl-CoA dehydrogenases (whose x-ray crystal structure has been determined) reveals that most of the mutations identified in GCDH protein seem to affect folding and tetramerization, as has been described for a number of mutations affecting mitochondrial β-oxidation acyl-CoA dehydrogenases.

MoDEL (Molecular Dynamics Extended Library): A Database of Atomistic Molecular Dynamics Trajectories

More than 1700 trajectories of proteins representative of monomeric soluble structures in the protein data bank (PDB) have been obtained by means of state-of-the-art atomistic molecular dynamics simulations in near-physiological conditions. The trajectories and analyses are stored in a large data warehouse, which can be queried for dynamic information on proteins, including interactions. Here, we describe the project and the structure and contents of our database, and provide examples of how it can be used to describe the global flexibility properties of proteins. Basic analyses and trajectories stripped of solvent molecules at a reduced resolution level are available from our web server.