The effects of death and post-mortem cold ischemia on human tissue transcriptomes

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.Peer ReviewedPostprint (published version

Gene expression variation and constraint across organs and species

Mice are the premier model organisms to study human biology and disease, but there is still debate about the extent to which molecular mechanisms are conserved between human and mouse. With the advance of next-generation sequencing technologies, comparative transcriptomics can be carried out at unprecedented resolution. In this thesis we confirm findings that human and mouse transcriptomes are globally conserved and we identify and characterize the properties of a core set of genes with constrained expression between the two species. Additionally, we show that clustering of gene expression across humans, mice and other vertebrates across multiple tissues depends on which genes and samples are included. Finally, we analyze gene expression of primary cells in human to understand how functional units of organs contribute to the expression of an organ as a whole. Indeed, profiling entire organs constitutes one of the main limitations of current comparative studies.Los ratones son los principales organismos modelos para estudiar la biología y las enfermedades humanas, pero aún está en debate el nivel de conservación molecular entre humanos y ratones. Con el progreso de las tecnologías de secuenciación masiva, la transcriptómica comparativa ha llegado a una resolución sin precedentes. En esta tesis confirmamos que los transcriptomas de humano y ratón están globalmente conservados y identificamos y caracterizamos las propiedades de un conjunto de genes con expresión parecida entre las dos especies. Además, demostramos que diferentes tejidos de humanos, ratones y otros vertebrados se agrupan en base a su expresión génica segón los genes y las muestras incluidas en el análisis. Finalmente, analizamos la expresión génica de líneas celulares primarias humanas para investigar cómo las unidades funcionales de los órganos afectan la expresión de todo un órgano entero. De hecho, los estudios comparativos actuales tienen como limitación que se basan en datos de órganos enteros

ggsashimi: Sashimi plot revised for browser- and annotation-independent splicing visualization

We present ggsashimi, a command-line tool for the visualization of splicing events across multiple samples. Given a specified genomic region, ggsashimi creates sashimi plots for individual RNA-seq experiments as well as aggregated plots for groups of experiments, a feature unique to this software. Compared to the existing versions of programs generating sashimi plots, it uses popular bioinformatics file formats, it is annotation-independent, and allows the visualization of splicing events even for large genomic regions by scaling down the genomic segments between splice sites. ggsashimi is freely available at https://github.com/guigolab/ggsashimi. It is implemented in python, and internally generates R code for plotting

ggsashimi: Sashimi plot revised for browser- and annotation-independent splicing visualization

We present ggsashimi, a command-line tool for the visualization of splicing events across multiple samples. Given a specified genomic region, ggsashimi creates sashimi plots for individual RNA-seq experiments as well as aggregated plots for groups of experiments, a feature unique to this software. Compared to the existing versions of programs generating sashimi plots, it uses popular bioinformatics file formats, it is annotation-independent, and allows the visualization of splicing events even for large genomic regions by scaling down the genomic segments between splice sites. ggsashimi is freely available at https://github.com/guigolab/ggsashimi. It is implemented in python, and internally generates R code for plotting

Transcriptional diversity during lineage commitment of human blood progenitors

Blood cells derive from hematopoietic stem cells through stepwise fating events. To characterize gene expression programs driving lineage choice, we sequenced RNA from eight primary human hematopoietic progenitor populations representing the major myeloid commitment stages and the main lymphoid stage. We identified extensive cell type-specific expression changes: 6711 genes and 10,724 transcripts, enriched in non-protein-coding elements at early stages of differentiation. In addition, we found 7881 novel splice junctions and 2301 differentially used alternative splicing events, enriched in genes involved in regulatory processes. We demonstrated experimentally cell-specific isoform usage, identifying nuclear factor I/B (NFIB) as a regulator of megakaryocyte maturation-the platelet precursor. Our data highlight the complexity of fating events in closely related progenitor populations, the understanding of which is essential for the advancement of transplantation and regenerative medicine.The work described in this manuscript was primarily supported by the European Commission Seventh Framework Program through the BLUEPRINT grant with code HEALTH-F5-2011-282510 (DH, FB, GC, JHAM, KD, LC, MF, SC, SF and SPG). Research in the Ouwehand laboratory is further supported by program grants from the National Institute for Health Research (NIHR, http://www.nihr.ac.uk; to AA, MK, PP, SBGJ, SN, and WHO); and the British Heart Foundation under numbers RP-PG-0310-1002 and RG/09/12/28096 (http://www.bhf.org.uk; to AR and WJA). KF and MK were supported by Marie Curie funding from the NETSIM FP7 program funded by the European Commission. The Cambridge BioResource (http://www.cambridgebioresource.org.uk), the Cell Phenotyping Hub, and the Cambridge Translational GenOmics laboratory (http://www.catgo.org.uk) are supported by an NIHR grant to the Cambridge NIHR Biomedical Research Centre (BRC). Research in the Soranzo laboratory (LV, NS and SW) is further supported by the Wellcome Trust (Grant Codes WT098051 and WT091310) and the EU FP7 EPIGENESYS initiative (Grant Code 257082). Research in the Cvejic laboratory (AC and CL) is funded by the Cancer Research UK under grant number C45041/A14953. SJS is funded by NIHR. MEF is supported by a British Heart Foundation Clinical Research Training Fellowship, number FS/12/27/29405. EBM is supported by a Wellcome Trust grant, number 084183/Z/07/Z. FAC, CL and SW are supported by MRC Clinical Training Fellowships and TB by a British Society of Haematology/NHS Blood and Transplant grant. RJR is a Principal Research Fellow of the Wellcome Trust, grant No. 082961/Z/07/Z. Research in the Flicek laboratory is also supported by the Wellcome Trust (grant number 095908) and EMBL. Research in the Bertone laboratory is supported by EMBL. KF and CvG are supported by FWO-Vlaanderen through grant G.0B17.13

Gene-specific patterns of expression variation across organs and species

Background: A comparison of transcriptional profiles derived from different tissues in a given species or among different species assumes that commonalities reflect evolutionarily conserved programs and that differences reflect species or tissue responses to environmental conditions or developmental program staging. Apparently conflicting results have been published regarding whether organ-specific transcriptional patterns dominate over species-specific patterns, or vice versa, making it unclear to what extent the biology of a given organism can be extrapolated to another. These studies have in common that they treat the transcriptomes monolithically, implicitly ignoring that each gene is likely to have a specific pattern of transcriptional variation across organs and species. Results: We use linear models to quantify this pattern. We find a continuum in the spectrum of expression variation: the expression of some genes varies considerably across species and little across organs, and simply reflects evolutionary distance. At the other extreme are genes whose expression varies considerably across organs and little across species; these genes are much more likely to be associated with diseases than are genes whose expression varies predominantly across species. Conclusions: Whether transcriptomes, when considered globally, cluster preferentially according to one component or the other may not be a property of the transcriptomes, but rather a consequence of the dominant behavior of a subset of genes. Therefore, the values of the components of the variance of expression for each gene could become a useful resource when planning, interpreting, and extrapolating experimental data from mouse to humans.This project was supported by awards U54HG007004 and U41HG007234 from the National Human Genome Research Institute of the National Institutes of Health, as well as from the Spanish Ministry of Economy and Competitiveness, Centro de Excelencia Severo Ochoa 2013–2017, SEV-2012-0208, and Programa de Ayudas FPI del Ministerio de Economia y Competitividad, BES-2012-055848. We would also like to acknowledge support from the European Research Council (ERC) under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement 294653

bsAS, an antisense long non-coding RNA, essential for correct wing development through regulation of blistered/DSRF isoform usage

Natural Antisense Transcripts (NATs) are long non-coding RNAs (lncRNAs) that overlap coding genes in the opposite strand. NATs roles have been related to gene regulation through different mechanisms, including post-transcriptional RNA processing. With the aim to identify NATs with potential regulatory function during fly development, we generated RNA-Seq data in Drosophila developing tissues and found bsAS, one of the most highly expressed lncRNAs in the fly wing. bsAS is antisense to bs/DSRF, a gene involved in wing development and neural processes. bsAS plays a crucial role in the tissue specific regulation of the expression of the bs/DSRF isoforms. This regulation is essential for the correct determination of cell fate during Drosophila development, as bsAS knockouts show highly aberrant phenotypes. Regulation of bs isoform usage by bsAS is mediated by specific physical interactions between the promoters of these two genes, which suggests a regulatory mechanism involving the collision of RNA polymerases transcribing in opposite directions. Evolutionary analysis suggests that bsAS NAT emerged simultaneously to the long-short isoform structure of bs, preceding the emergence of wings in insects.This work was supported by the European Community under the FP7 program (ERC-2011-AdG-294653-RNA-MAPS to R.G.), by the Spanish Ministry of Economy and Competitiveness (MEC) (BIO2011-26205 to R.G), by the Centro de Excelencia Severo Ochoa, from the CERCA Programme (Generalitat de Catalunya), and from the Spanish Ministry of Science and Innovation to the EMBL partnership. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

The histone code reader PHD finger protein 7 controls sex-linked disparities in gene expression and malignancy in Drosophila

The notable male predominance across many human cancer types remains unexplained. Here, we show that Drosophila l(3)mbt brain tumors are more invasive and develop as malignant neoplasms more often in males than in females. By quantitative proteomics, we have identified a signature of proteins that are differentially expressed between male and female tumor samples. Prominent among them is the conserved chromatin reader PHD finger protein 7 (Phf7). We show that Phf7 depletion reduces sex-dependent differences in gene expression and suppresses the enhanced malignant traits of male tumors. Our results identify potential regulators of sex-linked tumor dimorphism and show that these genes may serve as targets to suppress sex-linked malignant traits.Part of this study was funded by ERC AdG 2011 294603 and BFU2015-66304-P from the MINECO, Spain

The histone code reader PHD finger protein 7 controls sex-linked disparities in gene expression and malignancy in Drosophila

The notable male predominance across many human cancer types remains unexplained. Here, we show that Drosophila l(3)mbt brain tumors are more invasive and develop as malignant neoplasms more often in males than in females. By quantitative proteomics, we have identified a signature of proteins that are differentially expressed between male and female tumor samples. Prominent among them is the conserved chromatin reader PHD finger protein 7 (Phf7). We show that Phf7 depletion reduces sex-dependent differences in gene expression and suppresses the enhanced malignant traits of male tumors. Our results identify potential regulators of sex-linked tumor dimorphism and show that these genes may serve as targets to suppress sex-linked malignant traits.Part of this study was funded by ERC AdG 2011 294603 and BFU2015-66304-P from the MINECO, Spain

Enhanced transcriptome maps from multiple mouse tissues reveal evolutionary constraint in gene expression

Mice have been a long-standing model for human biology and disease. Here we characterize, by RNA sequencing, the transcriptional profiles of a large and heterogeneous collection of mouse tissues, augmenting the mouse transcriptome with thousands of novel transcript candidates. Comparison with transcriptome profiles in human cell lines reveals substantial conservation of transcriptional programmes, and uncovers a distinct class of genes with levels of expression that have been constrained early in vertebrate evolution. This core set of genes captures a substantial fraction of the transcriptional output of mammalian cells, and participates in basic functional and structural housekeeping processes common to all cell types. Perturbation of these constrained genes is associated with significant phenotypes including embryonic lethality and cancer. Evolutionary constraint in gene expression levels is not reflected in the conservation of the genomic sequences, but is associated with conserved epigenetic marking, as well as with characteristic post-transcriptional regulatory programme, in which sub-cellular localization and alternative splicing play comparatively large roles.This work was supported by the National Human Genome Research Institute (NHGRI) grants number 1U54HG007004-1, U41HG007000, U01HG004695, U54HG004555 and U41HG007234, by the Spanish Plan Nacional grants number BIO2011-26205 and BFU2011-28575, the ERC grant number 294653, LaCaixa and the EU-FP7 quantomics project