@article {611, title = {Association of a single nucleotide polymorphism in the ubxn6 gene with long-term non-progression phenotype in HIV-positive individuals.}, journal = {Clin Microbiol Infect}, volume = {26}, year = {2020}, month = {2020 Jan}, pages = {107-114}, abstract = {

OBJECTIVES: The long-term non-progressors (LTNPs) are a heterogeneous group of HIV-positive individuals characterized by their ability to maintain high CD4 T-cell counts and partially control viral replication for years in the absence of antiretroviral therapy. The present study aims to identify host single nucleotide polymorphisms (SNPs) associated with non-progression in a cohort of 352 individuals.

METHODS: DNA microarrays and exome sequencing were used for genotyping about 240~000 functional polymorphisms throughout more than 20~000 human genes. The allele frequencies of 85 LTNPs were compared with a control population. SNPs associated with LTNPs were confirmed in a population of typical progressors. Functional analyses in the affected gene were carried out through knockdown experiments in HeLa-P4, macrophages and dendritic cells.

RESULTS: Several SNPs located within the major histocompatibility complex region previously related to LTNPs were confirmed in this new cohort. The SNP rs1127888 (UBXN6) surpassed the statistical significance of these markers after Bonferroni correction (q~=~2.11~{\texttimes}~10). An uncommon allelic frequency of rs1127888 among LTNPs was confirmed by comparison with typical progressors and other publicly available populations. UBXN6 knockdown experiments caused an increase in CAV1 expression and its accumulation in the plasma membrane. In~vitro infection of different cell types with HIV-1 replication-competent recombinant viruses caused a reduction of the viral replication capacity compared with their corresponding wild-type cells expressing UBXN6.

CONCLUSIONS: A higher prevalence of Ala31Thr in UBXN6 was found among LTNPs within its N-terminal region, which is crucial for UBXN6/VCP protein complex formation. UBXN6 knockdown affected CAV1 turnover and HIV-1 replication capacity.

}, keywords = {Adaptor Proteins, Vesicular Transport, Autophagy-Related Proteins, Caveolin 1, Cohort Studies, Dendritic Cells, Disease Progression, Gene Frequency, Gene Knockdown Techniques, Genetic Association Studies, HeLa Cells, HIV Infections, HIV Long-Term Survivors, HIV-1, Humans, Macrophages, Oligonucleotide Array Sequence Analysis, Phenotype, Polymorphism, Single Nucleotide, whole exome sequencing}, issn = {1469-0691}, doi = {10.1016/j.cmi.2019.05.015}, author = {D{\'\i}ez-Fuertes, F and De La Torre-Tarazona, H E and Calonge, E and Pernas, M and Bermejo, M and Garc{\'\i}a-P{\'e}rez, J and {\'A}lvarez, A and Capa, L and Garc{\'\i}a-Garc{\'\i}a, F and Saumoy, M and Riera, M and Boland-Auge, A and L{\'o}pez-Gal{\'\i}ndez, C and Lathrop, M and Dopazo, J and Sakuntabhai, A and Alcam{\'\i}, J} } @article {422, title = {Differential metabolic activity and discovery of therapeutic targets using summarized metabolic pathway models.}, journal = {NPJ Syst Biol Appl}, volume = {5}, year = {2019}, month = {2019}, pages = {7}, abstract = {

In spite of the increasing availability of genomic and transcriptomic data, there is still a gap between the detection of perturbations in gene expression and the understanding of their contribution to the molecular mechanisms that ultimately account for the phenotype studied. Alterations in the metabolism are behind the initiation and progression of many diseases, including cancer. The wealth of available knowledge on metabolic processes can therefore be used to derive mechanistic models that link gene expression perturbations to changes in metabolic activity that provide relevant clues on molecular mechanisms of disease and drug modes of action (MoA). In particular, pathway modules, which recapitulate the main aspects of metabolism, are especially suitable for this type of modeling. We present Metabolizer, a web-based application that offers an intuitive, easy-to-use interactive interface to analyze differences in pathway metabolic module activities that can also be used for class prediction and in silico prediction of knock-out (KO) effects. Moreover, Metabolizer can automatically predict the optimal KO intervention for restoring a diseased phenotype. We provide different types of validations of some of the predictions made by Metabolizer. Metabolizer is a web tool that allows understanding molecular mechanisms of disease or the MoA of drugs within the context of the metabolism by using gene expression measurements. In addition, this tool automatically suggests potential therapeutic targets for individualized therapeutic interventions.

}, keywords = {Computational Biology, Computer Simulation, Drug discovery, Gene Regulatory Networks, Humans, Internet, Metabolic Networks and Pathways, Models, Biological, Neoplasms, Phenotype, Software, Transcriptome}, issn = {2056-7189}, doi = {10.1038/s41540-019-0087-2}, author = {Cubuk, Cankut and Hidalgo, Marta R and Amadoz, Alicia and Rian, Kinza and Salavert, Francisco and Pujana, Miguel A and Mateo, Francesca and Herranz, Carmen and Carbonell-Caballero, Jos{\'e} and Dopazo, Joaquin} } @article {610, title = {Exploring the druggable space around the Fanconi anemia pathway using machine learning and mechanistic models.}, journal = {BMC Bioinformatics}, volume = {20}, year = {2019}, month = {2019 Jul 02}, pages = {370}, abstract = {

BACKGROUND: In spite of the abundance of genomic data, predictive models that describe phenotypes as a function of gene expression or mutations are difficult to obtain because they are affected by the curse of dimensionality, given the disbalance between samples and candidate genes. And this is especially dramatic in scenarios in which the availability of samples is difficult, such as the case of rare diseases.

RESULTS: The application of multi-output regression machine learning methodologies to predict the potential effect of external proteins over the signaling circuits that trigger Fanconi anemia related cell functionalities, inferred with a mechanistic model, allowed us to detect over 20 potential therapeutic targets.

CONCLUSIONS: The use of artificial intelligence methods for the prediction of potentially causal relationships between proteins of interest and cell activities related with disease-related phenotypes opens promising avenues for the systematic search of new targets in rare diseases.

}, keywords = {Databases, Factual, Fanconi Anemia, Genomics, Humans, Machine Learning, Phenotype, Proteins, Signal Transduction}, issn = {1471-2105}, doi = {10.1186/s12859-019-2969-0}, author = {Esteban-Medina, Marina and Pe{\~n}a-Chilet, Maria and Loucera, Carlos and Dopazo, Joaquin} } @article {405, title = {Gene Expression Integration into Pathway Modules Reveals a Pan-Cancer Metabolic Landscape.}, journal = {Cancer Res}, volume = {78}, year = {2018}, month = {2018 11 01}, pages = {6059-6072}, abstract = {

Metabolic reprogramming plays an important role in cancer development and progression and is a well-established hallmark of cancer. Despite its inherent complexity, cellular metabolism can be decomposed into functional modules that represent fundamental metabolic processes. Here, we performed a pan-cancer study involving 9,428 samples from 25 cancer types to reveal metabolic modules whose individual or coordinated activity predict cancer type and outcome, in turn highlighting novel therapeutic opportunities. Integration of gene expression levels into metabolic modules suggests that the activity of specific modules differs between cancers and the corresponding tissues of origin. Some modules may cooperate, as indicated by the positive correlation of their activity across a range of tumors. The activity of many metabolic modules was significantly associated with prognosis at a stronger magnitude than any of their constituent genes. Thus, modules may be classified as tumor suppressors and oncomodules according to their potential impact on cancer progression. Using this modeling framework, we also propose novel potential therapeutic targets that constitute alternative ways of treating cancer by inhibiting their reprogrammed metabolism. Collectively, this study provides an extensive resource of predicted cancer metabolic profiles and dependencies. Combining gene expression with metabolic modules identifies molecular mechanisms of cancer undetected on an individual gene level and allows discovery of new potential therapeutic targets. .

}, keywords = {Cell Line, Tumor, Cluster Analysis, Disease Progression, Gene Expression Profiling, Gene Expression Regulation, Neoplastic, Gene Regulatory Networks, Humans, Kaplan-Meier Estimate, Metabolome, mutation, Neoplasms, Oncogenes, Phenotype, Prognosis, RNA, Small Interfering, Sequence Analysis, RNA, Transcriptome, Treatment Outcome}, issn = {1538-7445}, doi = {10.1158/0008-5472.CAN-17-2705}, author = {Cubuk, Cankut and Hidalgo, Marta R and Amadoz, Alicia and Pujana, Miguel A and Mateo, Francesca and Herranz, Carmen and Carbonell-Caballero, Jos{\'e} and Dopazo, Joaquin} } @article {384, title = {Genomic expression differences between cutaneous cells from red hair color individuals and black hair color individuals based on bioinformatic analysis.}, journal = {Oncotarget}, volume = {8}, year = {2017}, month = {2017 Feb 14}, pages = {11589-11599}, abstract = {

The MC1R gene plays a crucial role in pigmentation synthesis. Loss-of-function MC1R variants, which impair protein function, are associated with red hair color (RHC) phenotype and increased skin cancer risk. Cultured cutaneous cells bearing loss-of-function MC1R variants show a distinct gene expression profile compared to wild-type MC1R cultured cutaneous cells. We analysed the gene signature associated with RHC co-cultured melanocytes and keratinocytes by Protein-Protein interaction (PPI) network analysis to identify genes related with non-functional MC1R variants. From two detected networks, we selected 23 nodes as hub genes based on topological parameters. Differential expression of hub genes was then evaluated in healthy skin biopsies from RHC and black hair color (BHC) individuals. We also compared gene expression in melanoma tumors from individuals with RHC versus BHC. Gene expression in normal skin from RHC cutaneous cells showed dysregulation in 8 out of 23 hub genes (CLN3, ATG10, WIPI2, SNX2, GABARAPL2, YWHA, PCNA and GBAS). Hub genes did not differ between melanoma tumors in RHC versus BHC individuals. The study suggests that healthy skin cells from RHC individuals present a constitutive genomic deregulation associated with the red hair phenotype and identify novel genes involved in melanocyte biology.

}, keywords = {Adult, Coculture Techniques, Computational Biology, gene expression, Genetic Predisposition to Disease, Genomics, Hair Color, Humans, Keratinocytes, Melanocytes, Middle Aged, Phenotype, Receptor, Melanocortin, Type 1}, issn = {1949-2553}, doi = {10.18632/oncotarget.14140}, url = {http://www.impactjournals.com/oncotarget/index.php?journal=oncotarget\&page=article\&op=view\&path\%5B\%5D=14140\&path\%5B\%5D=45094}, author = {Puig-Butille, Joan Anton and Gimenez-Xavier, Pol and Visconti, Alessia and Nsengimana, J{\'e}r{\'e}mie and Garcia-Garcia, Francisco and Tell-Marti, Gemma and Escamez, Maria Jos{\'e} and Newton-Bishop, Julia and Bataille, Veronique and Del Rio, Marcela and Dopazo, Joaquin and Falchi, Mario and Puig, Susana} } @article {433, title = {Mutations in TRAPPC11 are associated with a congenital disorder of glycosylation.}, journal = {Hum Mutat}, volume = {38}, year = {2017}, month = {2017 02}, pages = {148-151}, abstract = {

Congenital disorders of glycosylation (CDG) are a heterogeneous and rapidly growing group of diseases caused by abnormal glycosylation of proteins and/or lipids. Mutations in genes involved in the homeostasis of the endoplasmic reticulum (ER), the Golgi apparatus (GA), and the vesicular trafficking from the ER to the ER-Golgi intermediate compartment (ERGIC) have been found to be associated with CDG. Here, we report a patient with defects in both N- and O-glycosylation combined with a delayed vesicular transport in the GA due to mutations in TRAPPC11, a subunit of the TRAPPIII complex. TRAPPIII is implicated in the anterograde transport from the ER to the ERGIC as well as in the vesicle export from the GA. This report expands the spectrum of genetic alterations associated with CDG, providing new insights for the diagnosis and the understanding of the physiopathological mechanisms underlying glycosylation disorders.

}, keywords = {Abnormalities, Multiple, Alleles, Amino Acid Substitution, Brain, Congenital Disorders of Glycosylation, Genotype, Humans, Magnetic Resonance Imaging, Male, mutation, Phenotype, Vesicular Transport Proteins, Whole Genome Sequencing}, issn = {1098-1004}, doi = {10.1002/humu.23145}, author = {Matalonga, Leslie and Bravo, Miren and Serra-Peinado, Carla and Garc{\'\i}a-Pelegr{\'\i}, Elisabeth and Ugarteburu, Olatz and Vidal, Silvia and Llambrich, Maria and Quintana, Ester and Fuster-Jorge, Pedro and Gonzalez-Bravo, Maria Nieves and Beltran, Sergi and Dopazo, Joaquin and Garcia-Garcia, Francisco and Foulquier, Fran{\c c}ois and Matthijs, Gert and Mills, Philippa and Ribes, Antonia and Egea, Gustavo and Briones, Paz and Tort, Frederic and Gir{\'o}s, Marisa} } @article {452, title = {Improving the management of Inherited Retinal Dystrophies by targeted sequencing of a population-specific gene panel.}, journal = {Sci Rep}, volume = {6}, year = {2016}, month = {2016 Apr 01}, pages = {23910}, abstract = {

Next-generation sequencing (NGS) has overcome important limitations to the molecular diagnosis of Inherited Retinal Dystrophies (IRD) such as the high clinical and genetic heterogeneity and the overlapping phenotypes. The purpose of this study was the identification of the genetic defect in 32 Spanish families with different forms of IRD. With that aim, we implemented a custom NGS panel comprising 64 IRD-associated genes in our population, and three disease-associated intronic regions. A total of 37 pathogenic mutations (14 novels) were found in 73\% of IRD patients ranging from 50\% for autosomal dominant cases, 75\% for syndromic cases, 83\% for autosomal recessive cases, and 100\% for X-linked cases. Additionally, unexpected phenotype-genotype correlations were found in 6 probands, which led to the refinement of their clinical diagnoses. Furthermore, intra- and interfamilial phenotypic variability was observed in two cases. Moreover, two cases unsuccessfully analysed by exome sequencing were resolved by applying this panel. Our results demonstrate that this hypothesis-free approach based on frequently mutated, population-specific loci is highly cost-efficient for the routine diagnosis of this heterogeneous condition and allows the unbiased analysis of a miscellaneous cohort. The molecular information found here has aid clinical diagnosis and has improved genetic counselling and patient management.

}, keywords = {Alleles, Computer Simulation, DNA Copy Number Variations, DNA Mutational Analysis, Eye Proteins, Gene Library, Genetic Association Studies, Genetic Heterogeneity, Genetic Therapy, High-Throughput Nucleotide Sequencing, Humans, mutation, Phenotype, Retinal Dystrophies}, issn = {2045-2322}, doi = {10.1038/srep23910}, author = {Bravo-Gil, Nereida and M{\'e}ndez-Vidal, Cristina and Romero-P{\'e}rez, Laura and Gonz{\'a}lez-del Pozo, Mar{\'\i}a and Rodr{\'\i}guez-de la R{\'u}a, Enrique and Dopazo, Joaquin and Borrego, Salud and Anti{\v n}olo, Guillermo} } @article {449, title = {Mutations in the MORC2 gene cause axonal Charcot-Marie-Tooth disease.}, journal = {Brain}, volume = {139}, year = {2016}, month = {2016 Jan}, pages = {62-72}, abstract = {

Charcot-Marie-Tooth disease (CMT) is a complex disorder with wide genetic heterogeneity. Here we present a new axonal Charcot-Marie-Tooth disease form, associated with the gene microrchidia family CW-type zinc finger 2 (MORC2). Whole-exome sequencing in a family with autosomal dominant segregation identified the novel MORC2 p.R190W change in four patients. Further mutational screening in our axonal Charcot-Marie-Tooth disease clinical series detected two additional sporadic cases, one patient who also carried the same MORC2 p.R190W mutation and another patient that harboured a MORC2 p.S25L mutation. Genetic and in silico studies strongly supported the pathogenicity of these sequence variants. The phenotype was variable and included patients with congenital or infantile onset, as well as others whose symptoms started in the second decade. The patients with early onset developed a spinal muscular atrophy-like picture, whereas in the later onset cases, the initial symptoms were cramps, distal weakness and sensory impairment. Weakness and atrophy progressed in a random and asymmetric fashion and involved limb girdle muscles, leading to a severe incapacity in adulthood. Sensory loss was always prominent and proportional to disease severity. Electrophysiological studies were consistent with an asymmetric axonal motor and sensory neuropathy, while fasciculations and myokymia were recorded rather frequently by needle electromyography. Sural nerve biopsy revealed pronounced multifocal depletion of myelinated fibres with some regenerative clusters and occasional small onion bulbs. Morc2 is expressed in both axons and Schwann cells of mouse peripheral nerve. Different roles in biological processes have been described for MORC2. As the silencing of Charcot-Marie-Tooth disease genes have been associated with DNA damage response, it is tempting to speculate that a deregulation of this pathway may be linked to the axonal degeneration observed in MORC2 neuropathy, thus adding a new pathogenic mechanism to the long list of causes of Charcot-Marie-Tooth disease.

}, keywords = {Adult, Aged, Animals, Axons, Charcot-Marie-Tooth Disease, Female, gene expression, Humans, Infant, Male, Mice, Middle Aged, mutation, Pedigree, Phenotype, Sciatic Nerve, Sural Nerve, Transcription Factors, Young Adult}, issn = {1460-2156}, doi = {10.1093/brain/awv311}, author = {Sevilla, Teresa and Lupo, Vincenzo and Mart{\'\i}nez-Rubio, Dolores and Sancho, Paula and Sivera, Rafael and Chumillas, Mar{\'\i}a J and Garc{\'\i}a-Romero, Mar and Pascual-Pascual, Samuel I and Muelas, Nuria and Dopazo, Joaquin and V{\'\i}lchez, Juan J and Palau, Francesc and Espin{\'o}s, Carmen} } @article {450, title = {Screening of CD96 and ASXL1 in 11 patients with Opitz C or Bohring-Opitz syndromes.}, journal = {Am J Med Genet A}, volume = {170A}, year = {2016}, month = {2016 Jan}, pages = {24-31}, abstract = {

Opitz C trigonocephaly (or Opitz C syndrome, OTCS) and Bohring-Opitz syndrome (BOS or C-like syndrome) are two rare genetic disorders with phenotypic overlap. The genetic causes of these diseases are not understood. However, two genes have been associated with OTCS or BOS with dominantly inherited de novo mutations. Whereas CD96 has been related to OTCS (one case) and to BOS (one case), ASXL1 has been related to BOS only (several cases). In this study we analyze CD96 and ASXL1 in a group of 11 affected individuals, including 2 sibs, 10 of them were diagnosed with OTCS, and one had a BOS phenotype. Exome sequences were available on six patients with OTCS and three parent pairs. Thus, we could analyze the CD96 and ASXL1 sequences in these patients bioinformatically. Sanger sequencing of all exons of CD96 and ASXL1 was carried out in the remaining patients. Detailed scrutiny of the sequences and assessment of variants allowed us to exclude putative pathogenic and private mutations in all but one of the patients. In this patient (with BOS) we identified a de novo mutation in ASXL1 (c.2100dupT). By nature and location within the gene, this mutation resembles those previously described in other BOS patients and we conclude that it may be responsible for the condition. Our results indicate that in 10 of 11, the disease (OTCS or BOS) cannot be explained by small changes in CD96 or ASXL1. However, the cohort is too small to make generalizations about the genetic etiology of these diseases.

}, keywords = {Adolescent, Antigens, CD, Child, Child, Preschool, Craniosynostoses, Exome, Female, High-Throughput Nucleotide Sequencing, Humans, Infant, Intellectual Disability, Male, mutation, Pedigree, Phenotype, Prognosis, Repressor Proteins}, issn = {1552-4833}, doi = {10.1002/ajmg.a.37418}, author = {Urreizti, Roser and Roca-Ayats, Neus and Trepat, Judith and Garcia-Garcia, Francisco and Alem{\'a}n, Alejandro and Orteschi, Daniela and Marangi, Giuseppe and Neri, Giovanni and Opitz, John M and Dopazo, Joaquin and Cormand, Bru and Vilageliu, Llu{\"\i}sa and Balcells, Susana and Grinberg, Daniel} } @article {456, title = {The EGR2 gene is involved in axonal Charcot-Marie-Tooth disease.}, journal = {Eur J Neurol}, volume = {22}, year = {2015}, month = {2015 Dec}, pages = {1548-55}, abstract = {

BACKGROUND AND PURPOSE: A three-generation family affected by axonal Charcot-Marie-Tooth disease (CMT) was investigated with the aim of discovering genetic defects and to further characterize the phenotype.

METHODS: The clinical, nerve conduction studies and muscle magnetic resonance images of the patients were reviewed. A whole exome sequencing was performed and the changes were investigated by genetic studies, in silico analysis and luciferase reporter assays.

RESULTS: A novel c.1226G>A change (p.R409Q) in the EGR2 gene was identified. Patients presented with a typical, late-onset axonal CMT phenotype with variable severity that was confirmed in the ancillary tests. The in silico studies showed that the residue R409 is an evolutionary conserved amino acid. The p.R409Q mutation, which is predicted as probably damaging, would alter the conformation of the protein slightly and would cause a decrease of gene expression.

CONCLUSIONS: This is the first report of an EGR2 mutation presenting as an axonal CMT phenotype with variable severity. This study broadens the phenotype of the EGR2-related neuropathies and suggests that the genetic testing of patients suffering from axonal CMT should include the EGR2 gene.

}, keywords = {Adult, Aged, Aged, 80 and over, Axons, Charcot-Marie-Tooth Disease, Early Growth Response Protein 2, Exome, Female, Humans, Male, Middle Aged, mutation, Pedigree, Phenotype, Severity of Illness Index, Young Adult}, issn = {1468-1331}, doi = {10.1111/ene.12782}, author = {Sevilla, T and Sivera, R and Mart{\'\i}nez-Rubio, D and Lupo, V and Chumillas, M J and Calpena, E and Dopazo, J and V{\'\i}lchez, J J and Palau, F and Espin{\'o}s, C} } @article {490, title = {Novel RP1 mutations and a recurrent BBS1 variant explain the co-existence of two distinct retinal phenotypes in the same pedigree.}, journal = {BMC Genet}, volume = {15}, year = {2014}, month = {2014 Dec 14}, pages = {143}, abstract = {

BACKGROUND: Molecular diagnosis of Inherited Retinal Dystrophies (IRD) has long been challenging due to the extensive clinical and genetic heterogeneity present in this group of disorders. Here, we describe the clinical application of an integrated next-generation sequencing approach to determine the underlying genetic defects in a Spanish family with a provisional clinical diagnosis of autosomal recessive Retinitis Pigmentosa (arRP).

RESULTS: Exome sequencing of the index patient resulted in the identification of the homozygous BBS1 p.M390R mutation. Sanger sequencing of additional members of the family showed lack of co-segregation of the p.M390R variant in some individuals. Clinical reanalysis indicated co-ocurrence of two different phenotypes in the same family: Bardet-Biedl syndrome in the individual harboring the BBS1 mutation and non-syndromic arRP in extended family members. To identify possible causative mutations underlying arRP, we conducted disease-targeted gene sequencing using a panel of 26 IRD genes. The in-house custom panel was validated using 18 DNA samples known to harbor mutations in relevant genes. All variants were redetected, indicating a high mutation detection rate. This approach allowed the identification of two novel heterozygous null mutations in RP1 (c.4582_4585delATCA; p.I1528Vfs*10 and c.5962dupA; p.I1988Nfs*3) which co-segregated with the disease in arRP patients. Additionally, a mutational screening in 96 patients of our cohort with genetically unresolved IRD revealed the presence of the c.5962dupA mutation in one unrelated family.

CONCLUSIONS: The combination of molecular findings for RP1 and BBS1 genes through exome and gene panel sequencing enabled us to explain the co-existence of two different retinal phenotypes in a family. The identification of two novel variants in RP1 suggests that the use of panels containing the prevalent genes of a particular population, together with an optimized data analysis pipeline, is an efficient and cost-effective approach that can be reliably implemented into the routine diagnostic process of diverse inherited retinal disorders. Moreover, the identification of these novel variants in two unrelated families supports the relatively high prevalence of RP1 mutations in Spanish population and the role of private mutations for commonly mutated genes, while extending the mutational spectrum of RP1.

}, keywords = {Bardet-Biedl Syndrome, Base Sequence, Case-Control Studies, DNA Mutational Analysis, Eye Proteins, Genes, Recessive, Genetic Association Studies, Humans, Microsatellite Repeats, Microtubule-Associated Proteins, Mutation, Missense, Pedigree, Phenotype, Retina, Retinitis pigmentosa}, issn = {1471-2156}, doi = {10.1186/s12863-014-0143-2}, author = {M{\'e}ndez-Vidal, Cristina and Bravo-Gil, Nereida and Gonz{\'a}lez-del Pozo, Mar{\'\i}a and Vela-Boza, Alicia and Dopazo, Joaquin and Borrego, Salud and Anti{\v n}olo, Guillermo} } @article {494, title = {The role of the interactome in the maintenance of deleterious variability in human populations.}, journal = {Mol Syst Biol}, volume = {10}, year = {2014}, month = {2014 Sep 26}, pages = {752}, abstract = {

Recent genomic projects have revealed the existence of an unexpectedly large amount of deleterious variability in the human genome. Several hypotheses have been proposed to explain such an apparently high mutational load. However, the mechanisms by which deleterious mutations in some genes cause a pathological effect but are apparently innocuous in other genes remain largely unknown. This study searched for deleterious variants in the 1,000 genomes populations, as well as in a newly sequenced population of 252 healthy Spanish individuals. In addition, variants causative of monogenic diseases and somatic variants from 41 chronic lymphocytic leukaemia patients were analysed. The deleterious variants found were analysed in the context of the interactome to understand the role of network topology in the maintenance of the observed mutational load. Our results suggest that one of the mechanisms whereby the effect of these deleterious variants on the phenotype is suppressed could be related to the configuration of the protein interaction network. Most of the deleterious variants observed in healthy individuals are concentrated in peripheral regions of the interactome, in combinations that preserve their connectivity, and have a marginal effect on interactome integrity. On the contrary, likely pathogenic cancer somatic deleterious variants tend to occur in internal regions of the interactome, often with associated structural consequences. Finally, variants causative of monogenic diseases seem to occupy an intermediate position. Our observations suggest that the real pathological potential of a variant might be more a systems property rather than an intrinsic property of individual proteins.

}, keywords = {Alleles, Exome, Gene Library, Genetic Variation, Genetics, Population, Genome, Human, Genomics, Humans, Models, Genetic, mutation, Phenotype, Protein Conformation, Protein Interaction Maps, Sequence Analysis, DNA, Whites}, issn = {1744-4292}, doi = {10.15252/msb.20145222}, author = {Garc{\'\i}a-Alonso, Luz and Jim{\'e}nez-Almaz{\'a}n, Jorge and Carbonell-Caballero, Jos{\'e} and Vela-Boza, Alicia and Santoyo-L{\'o}pez, Javier and Anti{\v n}olo, Guillermo and Dopazo, Joaquin} } @article {501, title = {Role of CPI-17 in restoring skin homoeostasis in cutaneous field of cancerization: effects of topical application of a film-forming medical device containing photolyase and UV filters.}, journal = {Exp Dermatol}, volume = {22}, year = {2013}, month = {2013 Jul}, pages = {494-6}, abstract = {

Cutaneous field of cancerization (CFC) is caused in part by the carcinogenic effect of the cyclobutane pyrimidine dimers CPD and 6-4 photoproducts (6-4PPs). Photoreactivation is carried out by photolyases which specifically recognize and repair both photoproducts. The study evaluates the molecular effects of topical application of a film-forming medical device containing photolyase and UV filters on the precancerous field in AK from seven patients. Skin improvement after treatment was confirmed in all patients by histopathological and molecular assessment. A gene set analysis showed that skin recovery was associated with biological processes involved in tissue homoeostasis and cell maintenance. The CFC response was associated with over-expression of the CPI-17 gene, and a dependence on the initial expression level was observed (P~=~0.001). Low CPI-17 levels were directly associated with pro-inflammatory genes such as TNF (P~=~0.012) and IL-1B (P~=~0.07). Our results suggest a role for CPI-17 in restoring skin homoeostasis in CFC lesions.

}, keywords = {Administration, Topical, Adult, Aged, Aged, 80 and over, Biopsy, Deoxyribodipyrimidine Photo-Lyase, Female, Gene Expression Profiling, Gene Expression Regulation, Enzymologic, Gene Expression Regulation, Neoplastic, Homeostasis, Humans, Inflammation, Intracellular Signaling Peptides and Proteins, Liposomes, Male, Middle Aged, Muscle Proteins, Phenotype, Phosphoprotein Phosphatases, Reactive Oxygen Species, Skin, Skin Neoplasms, Ultraviolet Rays}, issn = {1600-0625}, doi = {10.1111/exd.12177}, author = {Puig-Butille, Joan Anton and Malvehy, Josep and Potrony, Miriam and Trullas, Carles and Garcia-Garcia, Francisco and Dopazo, Joaquin and Puig, Susana} } @article {519, title = {SNPeffect 4.0: on-line prediction of molecular and structural effects of protein-coding variants.}, journal = {Nucleic Acids Res}, volume = {40}, year = {2012}, month = {2012 Jan}, pages = {D935-9}, abstract = {

Single nucleotide variants (SNVs) are, together with copy number variation, the primary source of variation in the human genome and are associated with phenotypic variation such as altered response to drug treatment and susceptibility to disease. Linking structural effects of non-synonymous SNVs to functional outcomes is a major issue in structural bioinformatics. The SNPeffect database (http://snpeffect.switchlab.org) uses sequence- and structure-based bioinformatics tools to predict the effect of protein-coding SNVs on the structural phenotype of proteins. It integrates aggregation prediction (TANGO), amyloid prediction (WALTZ), chaperone-binding prediction (LIMBO) and protein stability analysis (FoldX) for structural phenotyping. Additionally, SNPeffect holds information on affected catalytic sites and a number of post-translational modifications. The database contains all known human protein variants from UniProt, but users can now also submit custom protein variants for a SNPeffect analysis, including automated structure modeling. The new meta-analysis application allows plotting correlations between phenotypic features for a user-selected set of variants.

}, keywords = {Databases, Protein, Humans, Internet, Meta-Analysis as Topic, Phenotype, Polymorphism, Single Nucleotide, Protein Conformation, Proteins}, issn = {1362-4962}, doi = {10.1093/nar/gkr996}, author = {De Baets, Greet and Van Durme, Joost and Reumers, Joke and Maurer-Stroh, Sebastian and Vanhee, Peter and Dopazo, Joaquin and Schymkowitz, Joost and Rousseau, Frederic} } @article {542, title = {Functional analysis of multiple genomic signatures demonstrates that classification algorithms choose phenotype-related genes.}, journal = {Pharmacogenomics J}, volume = {10}, year = {2010}, month = {2010 Aug}, pages = {310-23}, abstract = {

Gene expression signatures of toxicity and clinical response benefit both safety assessment and clinical practice; however, difficulties in connecting signature genes with the predicted end points have limited their application. The Microarray Quality Control Consortium II (MAQCII) project generated 262 signatures for ten clinical and three toxicological end points from six gene expression data sets, an unprecedented collection of diverse signatures that has permitted a wide-ranging analysis on the nature of such predictive models. A comprehensive analysis of the genes of these signatures and their nonredundant unions using ontology enrichment, biological network building and interactome connectivity analyses demonstrated the link between gene signatures and the biological basis of their predictive power. Different signatures for a given end point were more similar at the level of biological properties and transcriptional control than at the gene level. Signatures tended to be enriched in function and pathway in an end point and model-specific manner, and showed a topological bias for incoming interactions. Importantly, the level of biological similarity between different signatures for a given end point correlated positively with the accuracy of the signature predictions. These findings will aid the understanding, and application of predictive genomic signatures, and support their broader application in predictive medicine.

}, keywords = {Algorithms, Databases, Genetic, Endpoint Determination, Gene Expression Profiling, Genomics, Humans, Neural Networks, Computer, Oligonucleotide Array Sequence Analysis, Phenotype, Predictive Value of Tests, Proteins, Quality Control}, issn = {1473-1150}, doi = {10.1038/tpj.2010.35}, author = {Shi, W and Bessarabova, M and Dosymbekov, D and Dezso, Z and Nikolskaya, T and Dudoladova, M and Serebryiskaya, T and Bugrim, A and Guryanov, A and Brennan, R J and Shah, R and Dopazo, J and Chen, M and Deng, Y and Shi, T and Jurman, G and Furlanello, C and Thomas, R S and Corton, J C and Tong, W and Shi, L and Nikolsky, Y} } @article {578, title = {Exploring the antimicrobial action of a carbon monoxide-releasing compound through whole-genome transcription profiling of Escherichia coli.}, journal = {Microbiology (Reading)}, volume = {155}, year = {2009}, month = {2009 Mar}, pages = {813-824}, abstract = {

We recently reported that carbon monoxide (CO) has bactericidal activity. To understand its mode of action we analysed the gene expression changes occurring when Escherichia coli, grown aerobically and anaerobically, is treated with the CO-releasing molecule CORM-2 (tricarbonyldichlororuthenium(II) dimer). Microarray analysis shows that the E. coli CORM-2 response is multifaceted, with a high number of differentially regulated genes spread through several functional categories, namely genes involved in inorganic ion transport and metabolism, regulators, and genes implicated in post-translational modification, such as chaperones. CORM-2 has a higher impact in E. coli cells grown anaerobically, as judged by the repression of genes belonging to eight functional classes which are not seen in the response of aerobically CORM-2-treated cells. The biological relevance of the variations caused by CORM-2 was substantiated by studying the CORM-2 sensitivity of selected E. coli mutants. The results show that the deletion of redox-sensing regulators SoxS and OxyR increased the sensitivity to CORM-2 and suggest that while SoxS plays an important role in protection against CORM-2 under both growth conditions, OxyR seems to participate only in the aerobic CORM-2 response. Under anaerobic conditions, we found that the heat-shock proteins IbpA and IbpB contribute to CORM-2 defence since the deletion of these genes increases the sensitivity of the strain. The induction of several met genes and the hypersensitivity to CORM-2 of the DeltametR, DeltametI and DeltametN mutant strains suggest that CO has effects on the methionine metabolism of E. coli. CORM-2 also affects the transcription of several E. coli biofilm-related genes and increases biofilm formation in E. coli. In particular, the absence of tqsA or bhsA increases the resistance of E. coli to CORM-2, and deletion of tsqA leads to a strain that has lost its capacity to form biofilm upon treatment with CORM-2. In spite of the relatively stable nature of the CO molecule, our results show that CO is able to trigger a significant alteration in the transcriptome of E. coli which necessarily has effects in several key metabolic pathways.

}, keywords = {Biofilms, Carbon Monoxide, Escherichia coli, Escherichia coli Proteins, Gene Expression Profiling, Gene Expression Regulation, Bacterial, Genes, Bacterial, Genes, Regulator, Genetic Complementation Test, Methionine, Microbial Viability, mutation, Oligonucleotide Array Sequence Analysis, Organometallic Compounds, Phenotype, RNA, Bacterial}, issn = {1350-0872}, doi = {10.1099/mic.0.023911-0}, author = {Nobre, L{\'\i}gia S and Al-Shahrour, F{\'a}tima and Dopazo, Joaquin and Saraiva, L{\'\i}gia M} } @article {590, title = {Direct functional assessment of the composite phenotype through multivariate projection strategies.}, journal = {Genomics}, volume = {92}, year = {2008}, month = {2008 Dec}, pages = {373-83}, abstract = {

We present a novel approach for the analysis of transcriptomics data that integrates functional annotation of gene sets with expression values in a multivariate fashion, and directly assesses the relation of functional features to a multivariate space of response phenotypical variables. Multivariate projection methods are used to obtain new correlated variables for a set of genes that share a given function. These new functional variables are then related to the response variables of interest. The analysis of the principal directions of the multivariate regression allows for the identification of gene function features correlated with the phenotype. Two different transcriptomics studies are used to illustrate the statistical and interpretative aspects of the methodology. We demonstrate the superiority of the proposed method over equivalent approaches.

}, keywords = {Breast Neoplasms, Computational Biology, Databases, Genetic, Female, Gene Expression Profiling, Humans, Mathematical Computing, Multivariate Analysis, Phenotype}, issn = {1089-8646}, doi = {10.1016/j.ygeno.2008.05.015}, author = {Conesa, Ana and Bro, Rasmus and Garcia-Garcia, Francisco and Prats, Jos{\'e} Manuel and G{\"o}tz, Stefan and Kjeldahl, Karin and Montaner, David and Dopazo, Joaquin} }