@article {768, title = {Editorial: Critical assessment of massive data analysis (CAMDA) annual conference 2021.}, journal = {Front Genet}, volume = {14}, year = {2023}, month = {2023}, pages = {1154398}, issn = {1664-8021}, doi = {10.3389/fgene.2023.1154398}, author = {{\L}abaj, Pawe{\l} P and Dopazo, Joaquin and Xiao, Wenzhong and Kreil, David P} } @article {801, title = {Evaluation of a combined detection of SARS-CoV-2 and its variants using real-time allele-specific PCR strategy: an advantage for clinical practice.}, journal = {Epidemiol Infect}, volume = {151}, year = {2023}, month = {2023 Nov 24}, pages = {e201}, abstract = {

This study aimed to assess the ability of a real-time reverse transcription polymerase chain reaction (RT-PCR) with multiple targets to detect SARS-CoV-2 and its variants in a single test. Nasopharyngeal specimens were collected from patients in Granada, Spain, between January 2021 and December 2022. Five allele-specific RT-PCR kits were used sequentially, with each kit designed to detect a predominant variant at the time. When the Alpha variant was dominant, the kit included the HV69/70 deletion, E and N genes. When Delta replaced Alpha, the kit incorporated the L452R mutation in addition to E and N genes. When Omicron became dominant, L452R was replaced with the N679K mutation. Before incorporating each variant kit, a comparative analysis was carried out with SARS-CoV-2 whole genome sequencing (WGS). The results demonstrated that RT-PCR with multiple targets can provide rapid and effective detection of SARS-CoV-2 and its variants in a single test. A very high degree of agreement (96.2\%) was obtained between the comparison of RT-PCR and WGS. Allele-specific RT-PCR assays make it easier to implement epidemiological surveillance systems for effective public health decision making.

}, keywords = {Alleles, COVID-19, COVID-19 Testing, Humans, Real-Time Polymerase Chain Reaction, SARS-CoV-2, Sensitivity and Specificity}, issn = {1469-4409}, doi = {10.1017/S095026882300184X}, author = {Chaves-Blanco, Luc{\'\i}a and de Salazar, Adolfo and Fuentes, Ana and Vi{\~n}uela, Laura and Perez-Florido, Javier and Dopazo, Joaquin and Garc{\'\i}a, Federico} } @article {800, title = {Evidence of the association between increased use of direct oral anticoagulants and a reduction in the rate of atrial fibrillation-related stroke and major bleeding at the population level (2012-2019).}, journal = {Med Clin (Barc)}, year = {2023}, month = {2023 Nov 20}, abstract = {

BACKGROUND: The introduction of direct-acting oral anticoagulants (DOACs) has shown to decrease atrial fibrillation (AF)-related stroke and bleeding rates in clinical studies, but there is no certain evidence about their effects at the population level. Our aim was to assess changes in AF-related stroke and major bleeding rates between 2012 and 2019 in Andalusia (Spain), and the association between DOACs use and events rates at the population level.

METHODS: All patients with an AF diagnosis from 2012 to 2019 were identified using the Andalusian Health Population Base, that provides clinical information on all Andalusian people. Annual ischemic and hemorrhagic stroke, major bleeding rates, and used antithrombotic treatments were determined. Marginal hazard ratios (HR) were calculated for each treatment.

RESULTS: A total of 95,085 patients with an AF diagnosis were identified. Mean age was 76.1{\textpm}10.2 years (49.7\% women). An increase in the use of DOACs was observed throughout the study period in both males and females (p<0.001). The annual rate of ischemic stroke decreased by one third, while that of hemorrhagic stroke and major bleeding decreased 2-3-fold from 2012 to 2019. Marginal HR was lower than 0.50 for DOACs compared to VKA for all ischemic or hemorrhagic events.

CONCLUSIONS: In this contemporary population-based study using clinical and administrative databases in Andalusia, a significant reduction in the incidence of AF-related ischemic and hemorrhagic stroke and major bleeding was observed between 2012 and 2019. The increased use of DOACs seems to be associated with this reduction.

}, issn = {1578-8989}, doi = {10.1016/j.medcli.2023.10.008}, author = {Loucera, Carlos and Carmona, Rosario and Bostelmann, Gerrit and Mu{\~n}oyerro-Mu{\~n}iz, Dolores and Villegas, Rom{\'a}n and Gonzalez-Manzanares, Rafael and Dopazo, Joaquin and Anguita, Manuel} } @article {763, title = {Endoglin and MMP14 Contribute to Ewing Sarcoma Spreading by Modulation of Cell-Matrix Interactions.}, journal = {Int J Mol Sci}, volume = {23}, year = {2022}, month = {2022 Aug 04}, abstract = {

Endoglin (ENG) is a mesenchymal stem cell (MSC) marker typically expressed by active endothelium. This transmembrane glycoprotein is shed by matrix metalloproteinase 14 (MMP14). Our previous work demonstrated potent preclinical activity of first-in-class anti-ENG antibody-drug conjugates as a nascent strategy to eradicate Ewing sarcoma (ES), a devastating rare bone/soft tissue cancer with a putative MSC origin. We also defined a correlation between ENG and MMP14 expression in ES. Herein, we show that ENG expression is significantly associated with a dismal prognosis in a large cohort of ES patients. Moreover, both ENG/MMP14 are frequently expressed in primary ES tumors and metastasis. To deepen in their functional relevance in ES, we conducted transcriptomic and proteomic profiling of in vitro ES models that unveiled a key role of ENG and MMP14 in cell mechano-transduction. Migration and adhesion assays confirmed that loss of ENG disrupts actin filament assembly and filopodia formation, with a concomitant effect on cell spreading. Furthermore, we observed that ENG regulates cell-matrix interaction through activation of focal adhesion signaling and protein kinase C expression. In turn, loss of MMP14 contributed to a more adhesive phenotype of ES cells by modulating the transcriptional extracellular matrix dynamics. Overall, these results suggest that ENG and MMP14 exert a significant role in mediating correct spreading machinery of ES cells, impacting the aggressiveness of the disease.

}, keywords = {Bone Neoplasms, Endoglin, Humans, Matrix Metalloproteinase 14, Proteomics, Receptors, Growth Factor, Sarcoma, Ewing, Signal Transduction}, issn = {1422-0067}, doi = {10.3390/ijms23158657}, author = {Puerto-Camacho, Pilar and Diaz-Martin, Juan and Olmedo-Pelayo, Joaqu{\'\i}n and Bolado-Carrancio, Alfonso and Salguero-Aranda, Carmen and Jord{\'a}n-P{\'e}rez, Carmen and Esteban-Medina, Marina and Alamo-Alvarez, Inmaculada and Delgado-Bellido, Daniel and Lobo-Selma, Laura and Dopazo, Joaquin and Sastre, Ana and Alonso, Javier and Gr{\"u}newald, Thomas G P and Bernabeu, Carmelo and Byron, Adam and Brunton, Valerie G and Amaral, Ana Teresa and de Alava, Enrique} } @article {729, title = {The ELIXIR Human Copy Number Variations Community: building bioinformatics infrastructure for research.}, journal = {F1000Res}, volume = {9}, year = {2020}, month = {2020}, chapter = {1229}, abstract = {

Copy number variations (CNVs) are major causative contributors both in the genesis of genetic diseases and human neoplasias. While "High-Throughput" sequencing technologies are increasingly becoming the primary choice for genomic screening analysis, their ability to efficiently detect CNVs is still heterogeneous and remains to be developed. The aim of this white paper is to provide a guiding framework for the future contributions of ELIXIR{\textquoteright}s recently established with implications beyond human disease diagnostics and population genomics. This white paper is the direct result of a strategy meeting that took place in September 2018 in Hinxton (UK) and involved representatives of 11 ELIXIR Nodes. The meeting led to the definition of priority objectives and tasks, to address a wide range of CNV-related challenges ranging from detection and interpretation to sharing and training. Here, we provide suggestions on how to align these tasks within the ELIXIR Platforms strategy, and on how to frame the activities of this new ELIXIR Community in the international context.

}, keywords = {Computational Biology, DNA Copy Number Variations, High-Throughput Nucleotide Sequencing, Humans}, issn = {2046-1402}, doi = {10.12688/f1000research.24887.1}, author = {Salgado, David and Armean, Irina M and Baudis, Michael and Beltran, Sergi and Capella-Gut{\'\i}errez, Salvador and Carvalho-Silva, Denise and Dominguez Del Angel, Victoria and Dopazo, Joaquin and Furlong, Laura I and Gao, Bo and Garcia, Leyla and Gerloff, Dietlind and Gut, Ivo and Gyenesei, Attila and Habermann, Nina and Hancock, John M and Hanauer, Marc and Hovig, Eivind and Johansson, Lennart F and Keane, Thomas and Korbel, Jan and Lauer, Katharina B and Laurie, Steve and Lesko{\v s}ek, Brane and Lloyd, David and Marqu{\'e}s-Bonet, Tom{\'a}s and Mei, Hailiang and Monostory, Katalin and Pi{\~n}ero, Janet and Poterlowicz, Krzysztof and Rath, Ana and Samarakoon, Pubudu and Sanz, Ferran and Saunders, Gary and Sie, Daoud and Swertz, Morris A and Tsukanov, Kirill and Valencia, Alfonso and Vidak, Marko and Yenyxe Gonz{\'a}lez, Cristina and Ylstra, Bauke and B{\'e}roud, Christophe} } @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 {397, title = {The effects of death and post-mortem cold ischemia on human tissue transcriptomes.}, journal = {Nat Commun}, volume = {9}, year = {2018}, month = {2018 Feb 13}, pages = {490}, abstract = {

Post-mortem tissues samples are a key resource for investigating patterns of gene expression. However, the processes triggered by death and the post-mortem interval (PMI) can significantly alter physiologically normal RNA levels. We investigate the impact of PMI on gene expression using data from multiple tissues of post-mortem donors obtained from the GTEx project. We find that many genes change expression over relatively short PMIs in a tissue-specific manner, but this potentially confounding effect in a biological analysis can be minimized by taking into account appropriate covariates. By comparing ante- and post-mortem blood samples, we identify the cascade of transcriptional events triggered by death of the organism. These events do not appear to simply reflect stochastic variation resulting from mRNA degradation, but active and ongoing regulation of transcription. Finally, we develop a model to predict the time since death from the analysis of the transcriptome of a few readily accessible tissues.

}, keywords = {Blood, Cold Ischemia, Death, Female, gene expression, Humans, Models, Biological, Postmortem Changes, RNA, Messenger, Stochastic Processes, Transcriptome}, issn = {2041-1723}, doi = {10.1038/s41467-017-02772-x}, author = {Ferreira, Pedro G and Mu{\~n}oz-Aguirre, Manuel and Reverter, Ferran and S{\'a} Godinho, Caio P and Sousa, Abel and Amadoz, Alicia and Sodaei, Reza and Hidalgo, Marta R and Pervouchine, Dmitri and Carbonell-Caballero, Jos{\'e} and Nurtdinov, Ramil and Breschi, Alessandra and Amador, Raziel and Oliveira, Patr{\'\i}cia and Cubuk, Cankut and Curado, Jo{\~a}o and Aguet, Fran{\c c}ois and Oliveira, Carla and Dopazo, Joaquin and Sammeth, Michael and Ardlie, Kristin G and Guig{\'o}, Roderic} } @article {559, title = {Extension of human lncRNA transcripts by RACE coupled with long-read high-throughput sequencing (RACE-Seq)}, journal = {Nature Communications}, volume = {7}, year = {2016}, month = {Jan-11-2016}, doi = {10.1038/ncomms12339}, url = {http://www.nature.com/articles/ncomms12339http://www.nature.com/articles/ncomms12339.pdfhttp://www.nature.com/articles/ncomms12339.pdfhttp://www.nature.com/articles/ncomms12339}, author = {Lagarde, Julien and Uszczynska-Ratajczak, Barbara and Santoyo-L{\'o}pez, Javier and Gonzalez, Jose Manuel and Tapanari, Electra and Mudge, Jonathan M. and Steward, Charles A. and Wilming, Laurens and Tanzer, Andrea and Howald, C{\'e}dric and Chrast, Jacqueline and Vela-Boza, Alicia and Rueda, Antonio and Lopez-Domingo, Francisco J. and Dopazo, Joaquin and Reymond, Alexandre and Guig{\'o}, Roderic and Harrow, Jennifer} } @article {1211, title = {Extension of human lncRNA transcripts by RACE coupled with long-read high-throughput sequencing (RACE-Seq).}, journal = {Nature communications}, volume = {7}, year = {2016}, month = {2016}, pages = {12339}, abstract = {Long non-coding RNAs (lncRNAs) constitute a large, yet mostly uncharacterized fraction of the mammalian transcriptome. Such characterization requires a comprehensive, high-quality annotation of their gene structure and boundaries, which is currently lacking. Here we describe RACE-Seq, an experimental workflow designed to address this based on RACE (rapid amplification of cDNA ends) and long-read RNA sequencing. We apply RACE-Seq to 398 human lncRNA genes in seven tissues, leading to the discovery of 2,556 on-target, novel transcripts. About 60\% of the targeted loci are extended in either 5{\textquoteright} or 3{\textquoteright}, often reaching genomic hallmarks of gene boundaries. Analysis of the novel transcripts suggests that lncRNAs are as long, have as many exons and undergo as much alternative splicing as protein-coding genes, contrary to current assumptions. Overall, we show that RACE-Seq is an effective tool to annotate an organism{\textquoteright}s deep transcriptome, and compares favourably to other targeted sequencing techniques.}, issn = {2041-1723}, doi = {10.1038/ncomms12339}, url = {http://www.nature.com/articles/ncomms12339}, author = {Lagarde, Julien and Uszczynska-Ratajczak, Barbara and Santoyo-L{\'o}pez, Javier and Gonzalez, Jose Manuel and Tapanari, Electra and Mudge, Jonathan M and Steward, Charles A and Wilming, Laurens and Tanzer, Andrea and Howald, C{\'e}dric and Chrast, Jacqueline and Vela-Boza, Alicia and Antonio Rueda and L{\'o}pez-Domingo, Francisco J and Dopazo, Joaquin and Reymond, Alexandre and Guig{\'o}, Roderic and Harrow, Jennifer} } @article {471, title = {Exome sequencing reveals a high genetic heterogeneity on familial Hirschsprung disease}, journal = {Scientific Reports}, volume = {5}, year = {2015}, month = {Jan-12-2015}, doi = {10.1038/srep16473}, url = {http://www.nature.com/articles/srep16473http://www.nature.com/articles/srep16473.pdfhttp://www.nature.com/articles/srep16473.pdfhttp://www.nature.com/articles/srep16473}, author = {Luz{\'o}n-Toro, Berta and Gui, Hongsheng and Ruiz-Ferrer, Macarena and Sze-Man Tang, Clara and Fern{\'a}ndez, Raquel M. and Sham, Pak-Chung and Torroglosa, Ana and Kwong-Hang Tam, Paul and Espino-Pais{\'a}n, Laura and Cherny, Stacey S. and Bleda, Marta and Enguix-Riego, Mar{\'\i}a Del Valle and Dopazo, Joaquin and Anti{\v n}olo, Guillermo and Garcia-Barcel{\'o}, Maria-Merc{\`e} and Borrego, Salud} } @article {489, title = {Exome sequencing reveals novel and recurrent mutations with clinical significance in inherited retinal dystrophies.}, journal = {PLoS One}, volume = {9}, year = {2014}, month = {2014}, pages = {e116176}, abstract = {

This study aimed to identify the underlying molecular genetic cause in four Spanish families clinically diagnosed of Retinitis Pigmentosa (RP), comprising one autosomal dominant RP (adRP), two autosomal recessive RP (arRP) and one with two possible modes of inheritance: arRP or X-Linked RP (XLRP). We performed whole exome sequencing (WES) using NimbleGen SeqCap EZ Exome V3 sample preparation kit and SOLID 5500xl platform. All variants passing filter criteria were validated by Sanger sequencing to confirm familial segregation and the absence in local control population. This strategy allowed the detection of: (i) one novel heterozygous splice-site deletion in RHO, c.937-2_944del, (ii) one rare homozygous mutation in C2orf71, c.1795T>C; p.Cys599Arg, not previously associated with the disease, (iii) two heterozygous null mutations in ABCA4, c.2041C>T; p.R681* and c.6088C>T; p.R2030*, and (iv) one mutation, c.2405-2406delAG; p.Glu802Glyfs*31 in the ORF15 of RPGR. The molecular findings for RHO and C2orf71 confirmed the initial diagnosis of adRP and arRP, respectively, while patients with the two ABCA4 mutations, both previously associated with Stargardt disease, presented symptoms of RP with early macular involvement. Finally, the X-Linked inheritance was confirmed for the family with the RPGR mutation. This latter finding allowed the inclusion of carrier sisters in our preimplantational genetic diagnosis program.

}, keywords = {Adolescent, Adult, Amino Acid Sequence, Base Sequence, Child, Chromosome Segregation, DNA Mutational Analysis, Exome, Family, Female, Humans, Inheritance Patterns, Male, Middle Aged, Molecular Sequence Data, mutation, Pedigree, Retinal Dystrophies, Rhodopsin}, issn = {1932-6203}, doi = {10.1371/journal.pone.0116176}, author = {Gonz{\'a}lez-del Pozo, Mar{\'\i}a and M{\'e}ndez-Vidal, Cristina and Bravo-Gil, Nereida and Vela-Boza, Alicia and Dopazo, Joaquin and Borrego, Salud and Anti{\v n}olo, Guillermo} } @article {566, title = {Exome sequencing identifies a new mutation in SERAC1 in a patient with 3-methylglutaconic aciduria.}, journal = {Mol Genet Metab}, volume = {110}, year = {2013}, month = {2013 Sep-Oct}, pages = {73-7}, abstract = {

3-Methylglutaconic aciduria (3-MGA-uria) is a heterogeneous group of syndromes characterized by an increased excretion of 3-methylglutaconic and 3-methylglutaric acids. Five types of 3-MGA-uria (I to V) with different clinical presentations have been described. Causative mutations in TAZ, OPA3, DNAJC19, ATP12, ATP5E, and TMEM70 have been identified. After excluding the known genetic causes of 3-MGA-uria we used exome sequencing to investigate a patient with Leigh syndrome and 3-MGA-uria. We identified a homozygous variant in SERAC1 (c.202C>T; p.Arg68*), that generates a premature stop codon at position 68 of SERAC1 protein. Western blot analysis in patient{\textquoteright}s fibroblasts showed a complete absence of SERAC1 that was consistent with the prediction of a truncated protein and supports the pathogenic role of the mutation. During the course of this project a parallel study identified mutations in SERAC1 as the genetic cause of the disease in 15 patients with MEGDEL syndrome, which was compatible with the clinical and biochemical phenotypes of the patient described here. In addition, our patient developed microcephaly and optic atrophy, two features not previously reported in MEGDEL syndrome. We highlight the usefulness of exome sequencing to reveal the genetic bases of human rare diseases even if only one affected individual is available.

}, keywords = {Adolescent, Adult, Carboxylic Ester Hydrolases, Child, Exome, Female, High-Throughput Nucleotide Sequencing, Humans, Infant, Male, Metabolism, Inborn Errors, mutation}, issn = {1096-7206}, doi = {10.1016/j.ymgme.2013.04.021}, author = {Tort, Frederic and Garc{\'\i}a-Silva, Mar{\'\i}a Teresa and Ferrer-Cort{\`e}s, X{\`e}nia and Navarro-Sastre, Aleix and Garcia-Villoria, Judith and Coll, Maria Josep and Vidal, Enrique and Jim{\'e}nez-Almaz{\'a}n, Jorge and Dopazo, Joaquin and Briones, Paz and Elpeleg, Orly and Ribes, Antonia} } @article {514, title = {Extensive translatome remodeling during ER stress response in mammalian cells.}, journal = {PLoS One}, volume = {7}, year = {2012}, month = {2012}, pages = {e35915}, abstract = {

In this work we have described the translatome of two mammalian cell lines, NIH3T3 and Jurkat, by scoring the relative polysome association of \~{}10,000 mRNA under normal and ER stress conditions. We have found that translation efficiencies of mRNA correlated poorly with transcript abundance, although a general tendency was observed so that the highest translation efficiencies were found in abundant mRNA. Despite the differences found between mouse (NIH3T3) and human (Jurkat) cells, both cell types share a common translatome composed by \~{}800-900 mRNA that encode proteins involved in basic cellular functions. Upon stress, an extensive remodeling in translatomes was observed so that translation of \~{}50\% of mRNA was inhibited in both cell types, this effect being more dramatic for those mRNA that accounted for most of the cell translation. Interestingly, we found two subsets comprising 1000-1500 mRNA whose translation resisted or was induced by stress. Translation arrest resistant class includes many mRNA encoding aminoacyl tRNA synthetases, ATPases and enzymes involved in DNA replication and stress response such as BiP. This class of mRNA is characterized by high translation rates in both control and stress conditions. Translation inducible class includes mRNA whose translation was relieved after stress, showing a high enrichment in early response transcription factors of bZIP and zinc finger C2H2 classes. Unlike yeast, a general coordination between changes in translation and transcription upon stress (potentiation) was not observed in mammalian cells. Among the different features of mRNA analyzed, we found a relevant association of translation efficiency with the presence of upstream ATG in the 5{\textquoteright}UTR and with the length of coding sequence of mRNA, and a looser association with other parameters such as the length and the G+C content of 5{\textquoteright}UTR. A model for translatome remodeling during the acute phase of stress response in mammalian cells is proposed.

}, keywords = {Animals, Endoplasmic Reticulum Stress, Humans, Jurkat Cells, Mice, NIH 3T3 Cells, Oligonucleotide Array Sequence Analysis, Protein Biosynthesis, RNA, Messenger, Transcription, Genetic}, issn = {1932-6203}, doi = {10.1371/journal.pone.0035915}, author = {Ventoso, Iv{\'a}n and Kochetov, Alex and Montaner, David and Dopazo, Joaquin and Santoyo, Javier} } @article {531, title = {Early peroxisome proliferator-activated receptor gamma regulated genes involved in expansion of pancreatic beta cell mass.}, journal = {BMC Med Genomics}, volume = {4}, year = {2011}, month = {2011 Dec 30}, pages = {86}, abstract = {

BACKGROUND: The progression towards type 2 diabetes depends on the allostatic response of pancreatic beta cells to synthesise and secrete enough insulin to compensate for insulin resistance. The endocrine pancreas is a plastic tissue able to expand or regress in response to the requirements imposed by physiological and pathophysiological states associated to insulin resistance such as pregnancy, obesity or ageing, but the mechanisms mediating beta cell mass expansion in these scenarios are not well defined. We have recently shown that ob/ob mice with genetic ablation of PPARγ2, a mouse model known as the POKO mouse failed to expand its beta cell mass. This phenotype contrasted with the appropriate expansion of the beta cell mass observed in their obese littermate ob/ob mice. Thus, comparison of these models islets particularly at early ages could provide some new insights on early PPARγ dependent transcriptional responses involved in the process of beta cell mass expansion

RESULTS: Here we have investigated PPARγ dependent transcriptional responses occurring during the early stages of beta cell adaptation to insulin resistance in wild type, ob/ob, PPARγ2 KO and POKO mice. We have identified genes known to regulate both the rate of proliferation and the survival signals of beta cells. Moreover we have also identified new pathways induced in ob/ob islets that remained unchanged in POKO islets, suggesting an important role for PPARγ in maintenance/activation of mechanisms essential for the continued function of the beta cell.

CONCLUSIONS: Our data suggest that the expansion of beta cell mass observed in ob/ob islets is associated with the activation of an immune response that fails to occur in POKO islets. We have also indentified other PPARγ dependent differentially regulated pathways including cholesterol biosynthesis, apoptosis through TGF-β signaling and decreased oxidative phosphorylation.

}, keywords = {Animals, Cell Proliferation, Cell Survival, Cholesterol, Down-Regulation, Female, Gene Expression Regulation, Gene Knockout Techniques, Insulin Resistance, Insulin-Secreting Cells, Mice, obesity, Oxidation-Reduction, Phosphorylation, PPAR gamma, Signal Transduction, Transcription, Genetic, Transforming Growth Factor beta}, issn = {1755-8794}, doi = {10.1186/1755-8794-4-86}, author = {Vivas, Yurena and Martinez-Garcia, Cristina and Izquierdo, Adriana and Garcia-Garcia, Francisco and Callejas, Sergio and Velasco, Ismael and Campbell, Mark and Ros, Manuel and Dopazo, Ana and Dopazo, Joaquin and Vidal-Puig, Antonio and Medina-Gomez, Gema} } @article {532, title = {Early transcriptional defense responses in Arabidopsis cell suspension culture under high-light conditions.}, journal = {Plant Physiol}, volume = {156}, year = {2011}, month = {2011 Jul}, pages = {1439-56}, abstract = {

The early transcriptional defense responses and reactive oxygen species (ROS) production in Arabidopsis (Arabidopsis thaliana) cell suspension culture (ACSC), containing functional chloroplasts, were examined at high light (HL). The transcriptional analysis revealed that most of the ROS markers identified among the 449 transcripts with significant differential expression were transcripts specifically up-regulated by singlet oxygen ((1)O(2)). On the contrary, minimal correlation was established with transcripts specifically up-regulated by superoxide radical or hydrogen peroxide. The transcriptional analysis was supported by fluorescence microscopy experiments. The incubation of ACSC with the (1)O(2) sensor green reagent and 2{\textquoteright},7{\textquoteright}-dichlorofluorescein diacetate showed that the 30-min-HL-treated cultures emitted fluorescence that corresponded with the production of (1)O(2) but not of hydrogen peroxide. Furthermore, the in vivo photodamage of the D1 protein of photosystem II indicated that the photogeneration of (1)O(2) took place within the photosystem II reaction center. Functional enrichment analyses identified transcripts that are key components of the ROS signaling transduction pathway in plants as well as others encoding transcription factors that regulate both ROS scavenging and water deficit stress. A meta-analysis examining the transcriptional profiles of mutants and hormone treatments in Arabidopsis showed a high correlation between ACSC at HL and the fluorescent mutant family of Arabidopsis, a producer of (1)O(2) in plastids. Intriguingly, a high correlation was also observed with ABA deficient1 and more axillary growth4, two mutants with defects in the biosynthesis pathways of two key (apo)carotenoid-derived plant hormones (i.e. abscisic acid and strigolactones, respectively). ACSC has proven to be a valuable system for studying early transcriptional responses to HL stress.

}, keywords = {Arabidopsis, Blotting, Western, Cell Culture Techniques, Cells, Cultured, Chloroplasts, Cluster Analysis, Gene Expression Profiling, Gene Expression Regulation, Plant, Hydrogen Peroxide, Light, mutation, Oligonucleotide Array Sequence Analysis, Photosystem II Protein Complex, Plant Growth Regulators, Reproducibility of Results, Reverse Transcriptase Polymerase Chain Reaction, RNA, Messenger, Signal Transduction, Stress, Physiological, Transcription, Genetic}, issn = {1532-2548}, doi = {10.1104/pp.111.177766}, author = {Gonz{\'a}lez-P{\'e}rez, Sergio and Guti{\'e}rrez, Jorge and Garcia-Garcia, Francisco and Osuna, Daniel and Dopazo, Joaquin and Lorenzo, Oscar and Revuelta, Jos{\'e} L and Arellano, Juan B} } @article {533, title = {Evidence for short-time divergence and long-time conservation of tissue-specific expression after gene duplication.}, journal = {Brief Bioinform}, volume = {12}, year = {2011}, month = {2011 Sep}, pages = {442-8}, abstract = {

Gene duplication is one of the main mechanisms by which genomes can acquire novel functions. It has been proposed that the retention of gene duplicates can be associated to processes of tissue expression divergence. These models predict that acquisition of divergent expression patterns should be acquired shortly after the duplication, and that larger divergence in tissue expression would be expected for paralogs, as compared to orthologs of a similar age. Many studies have shown that gene duplicates tend to have divergent expression patterns and that gene family expansions are associated with high levels of tissue specificity. However, the timeframe in which these processes occur have rarely been investigated in detail, particularly in vertebrates, and most analyses do not include direct comparisons of orthologs as a baseline for the expected levels of tissue specificity in absence of duplications. To assess the specific contribution of duplications to expression divergence, we combine here phylogenetic analyses and expression data from human and mouse. In particular, we study differences in spatial expression among human-mouse paralogs, specifically duplicated after the radiation of mammals, and compare them to pairs of orthologs in the same species. Our results show that gene duplication leads to increased levels of tissue specificity and that this tends to occur promptly after the duplication event.

}, keywords = {Animals, Conserved Sequence, Evolution, Molecular, Gene Duplication, gene expression, Genome, Humans, Mice, Organ Specificity}, issn = {1477-4054}, doi = {10.1093/bib/bbr022}, author = {Huerta-Cepas, Jaime and Dopazo, Joaquin and Huynen, Martijn A and Gabald{\'o}n, Toni} } @article {21731483, title = {An evolutionary trade-off between protein turnover rate and protein aggregation favors a higher aggregation propensity in fast degrading proteins.}, journal = {PLoS computational biology}, volume = {7}, year = {2011}, month = {2011 Jun}, pages = {e1002090}, abstract = {

We previously showed the existence of selective pressure against protein aggregation by the enrichment of aggregation-opposing {\textquoteright}gatekeeper{\textquoteright} residues at strategic places along the sequence of proteins. Here we analyzed the relationship between protein lifetime and protein aggregation by combining experimentally determined turnover rates, expression data, structural data and chaperone interaction data on a set of more than 500 proteins. We find that selective pressure on protein sequences against aggregation is not homogeneous but that short-living proteins on average have a higher aggregation propensity and fewer chaperone interactions than long-living proteins. We also find that short-living proteins are more often associated to deposition diseases. These findings suggest that the efficient degradation of high-turnover proteins is sufficient to preclude aggregation, but also that factors that inhibit proteasomal activity, such as physiological ageing, will primarily affect the aggregation of short-living proteins.

}, doi = {10.1371/journal.pcbi.1002090}, author = {De Baets, Greet and Reumers, Joke and Delgado Blanco, Javier and Dopazo, Joaquin and Schymkowitz, Joost and Rousseau, Frederic} } @article {548, title = {ETE: a python Environment for Tree Exploration.}, journal = {BMC Bioinformatics}, volume = {11}, year = {2010}, month = {2010 Jan 13}, pages = {24}, abstract = {

BACKGROUND: Many bioinformatics analyses, ranging from gene clustering to phylogenetics, produce hierarchical trees as their main result. These are used to represent the relationships among different biological entities, thus facilitating their analysis and interpretation. A number of standalone programs are available that focus on tree visualization or that perform specific analyses on them. However, such applications are rarely suitable for large-scale surveys, in which a higher level of automation is required. Currently, many genome-wide analyses rely on tree-like data representation and hence there is a growing need for scalable tools to handle tree structures at large scale.

RESULTS: Here we present the Environment for Tree Exploration (ETE), a python programming toolkit that assists in the automated manipulation, analysis and visualization of hierarchical trees. ETE libraries provide a broad set of tree handling options as well as specific methods to analyze phylogenetic and clustering trees. Among other features, ETE allows for the independent analysis of tree partitions, has support for the extended newick format, provides an integrated node annotation system and permits to link trees to external data such as multiple sequence alignments or numerical arrays. In addition, ETE implements a number of built-in analytical tools, including phylogeny-based orthology prediction and cluster validation techniques. Finally, ETE{\textquoteright}s programmable tree drawing engine can be used to automate the graphical rendering of trees with customized node-specific visualizations.

CONCLUSIONS: ETE provides a complete set of methods to manipulate tree data structures that extends current functionality in other bioinformatic toolkits of a more general purpose. ETE is free software and can be downloaded from http://ete.cgenomics.org.

}, keywords = {Computational Biology, Databases, Genetic, Phylogeny, Software}, issn = {1471-2105}, doi = {10.1186/1471-2105-11-24}, author = {Huerta-Cepas, Jaime and Dopazo, Joaquin and Gabald{\'o}n, Toni} } @article {549, title = {Exploring the link between germline and somatic genetic alterations in breast carcinogenesis.}, journal = {PLoS One}, volume = {5}, year = {2010}, month = {2010 Nov 22}, pages = {e14078}, abstract = {

Recent genome-wide association studies (GWASs) have identified candidate genes contributing to cancer risk through low-penetrance mutations. Many of these genes were unexpected and, intriguingly, included well-known players in carcinogenesis at the somatic level. To assess the hypothesis of a germline-somatic link in carcinogenesis, we evaluated the distribution of somatic gene labels within the ordered results of a breast cancer risk GWAS. This analysis suggested frequent influence on risk of genetic variation in loci encoding for "driver kinases" (i.e., kinases encoded by genes that showed higher somatic mutation rates than expected by chance and, therefore, whose deregulation may contribute to cancer development and/or progression). Assessment of these predictions using a population-based case-control study in Poland replicated the association for rs3732568 in EPHB1 (odds ratio (OR) = 0.79; 95\% confidence interval (CI): 0.63-0.98; P(trend) = 0.031). Analyses by early age at diagnosis and by estrogen receptor α (ERα) tumor status indicated potential associations for rs6852678 in CDKL2 (OR = 0.32, 95\% CI: 0.10-1.00; P(recessive) = 0.044) and rs10878640 in DYRK2 (OR = 2.39, 95\% CI: 1.32-4.30; P(dominant) = 0.003), and for rs12765929, rs9836340, rs4707795 in BMPR1A, EPHA3 and EPHA7, respectively (ERα tumor status P(interaction)<0.05). The identification of three novel candidates as EPH receptor genes might indicate a link between perturbed compartmentalization of early neoplastic lesions and breast cancer risk and progression. Together, these data may lay the foundations for replication in additional populations and could potentially increase our knowledge of the underlying molecular mechanisms of breast carcinogenesis.

}, keywords = {Adult, Bone Morphogenetic Protein Receptors, Type I, Breast, Breast Neoplasms, Calcium-Calmodulin-Dependent Protein Kinases, Case-Control Studies, Cyclin-Dependent Kinases, Disease Progression, Estrogen Receptor alpha, Female, Gene Frequency, Genetic Predisposition to Disease, Genome-Wide Association Study, Genotype, Germ-Line Mutation, Humans, Odds Ratio, Poland, Polymorphism, Single Nucleotide, Protein Serine-Threonine Kinases, Protein-Tyrosine Kinases, Receptor Protein-Tyrosine Kinases, Receptor, EphA3, Receptor, EphA7, Receptor, EphB1, Risk Factors}, issn = {1932-6203}, doi = {10.1371/journal.pone.0014078}, author = {Bonifaci, N{\'u}ria and G{\'o}rski, Bohdan and Masoj{\'c}, Bartlomiej and Woko{\l}orczyk, Dominika and Jakubowska, Anna and D{\k e}bniak, Tadeusz and Berenguer, Antoni and Serra Musach, Jordi and Brunet, Joan and Dopazo, Joaquin and Narod, Steven A and Lubi{\'n}ski, Jan and L{\'a}zaro, Conxi and Cybulski, Cezary and Pujana, Miguel Angel} } @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 {591, title = {Expression and microarrays.}, journal = {Methods Mol Biol}, volume = {453}, year = {2008}, month = {2008}, pages = {245-55}, abstract = {

High throughput methodologies have increased by several orders of magnitude the amount of experimental microarray data available. Nevertheless, translating these data into useful biological knowledge remains a challenge. There is a risk of perceiving these methodologies as mere factories that produce never-ending quantities of data if a proper biological interpretation is not provided. Methods of interpreting these data are continuously evolving. Typically, a simple two-step approach has been used, in which genes of interest are first selected based on thresholds for the experimental values, and then enrichment in biologically relevant terms in the annotations of these genes is analyzed in a second step. For various reasons, such methods are quite poor in terms of performance and new procedures inspired by systems biology that directly address sets of functionally related genes are currently under development.

}, keywords = {Animals, Computational Biology, gene expression, Gene Expression Profiling, Humans, Oligonucleotide Array Sequence Analysis}, issn = {1064-3745}, doi = {10.1007/978-1-60327-429-6_12}, author = {Dopazo, Joaquin and Al-Shahrour, F{\'a}tima} } @article {604, title = {Evidence for systems-level molecular mechanisms of tumorigenesis.}, journal = {BMC Genomics}, volume = {8}, year = {2007}, month = {2007 Jun 20}, pages = {185}, abstract = {

BACKGROUND: Cancer arises from the consecutive acquisition of genetic alterations. Increasing evidence suggests that as a consequence of these alterations, molecular interactions are reprogrammed in the context of highly connected and regulated cellular networks. Coordinated reprogramming would allow the cell to acquire the capabilities for malignant growth.

RESULTS: Here, we determine the coordinated function of cancer gene products (i.e., proteins encoded by differentially expressed genes in tumors relative to healthy tissue counterparts, hereafter referred to as "CGPs") defined as their topological properties and organization in the interactome network. We show that CGPs are central to information exchange and propagation and that they are specifically organized to promote tumorigenesis. Centrality is identified by both local (degree) and global (betweenness and closeness) measures, and systematically appears in down-regulated CGPs. Up-regulated CGPs do not consistently exhibit centrality, but both types of cancer products determine the overall integrity of the network structure. In addition to centrality, down-regulated CGPs show topological association that correlates with common biological processes and pathways involved in tumorigenesis.

CONCLUSION: Given the current limited coverage of the human interactome, this study proposes that tumorigenesis takes place in a specific and organized way at the molecular systems-level and suggests a model that comprises the precise down-regulation of groups of topologically-associated proteins involved in particular functions, orchestrated with the up-regulation of specific proteins.

}, keywords = {Cell Transformation, Neoplastic, Gene Expression Profiling, Gene Expression Regulation, Neoplastic, Humans, Male, Models, Biological, Models, Genetic, Models, Statistical, Neoplasm Proteins, Neoplasms, Prostatic Neoplasms, Protein Interaction Mapping, RNA, Messenger, Signal Transduction, Systems biology}, issn = {1471-2164}, doi = {10.1186/1471-2164-8-185}, author = {Hern{\'a}ndez, Pilar and Huerta-Cepas, Jaime and Montaner, David and Al-Shahrour, F{\'a}tima and Valls, Joan and G{\'o}mez, Laia and Capell{\`a}, Gabriel and Dopazo, Joaquin and Pujana, Miguel Angel} }