1. Chromosomes and Gene Expression
  2. Developmental Biology

Allele-specific gene expression can underlie altered transcript abundance in zebrafish mutants

  1. Richard J White
  2. Eirinn Mackay
  3. Stephen W Wilson
  4. Elisabeth M Busch-Nentwich  Is a corresponding author
  1. University of Cambridge, United Kingdom
  2. University College London, United Kingdom
https://doi.org/10.7554/eLife.72825
Share this article
Cite this article
  1. Richard J White
  2. Eirinn Mackay
  3. Stephen W Wilson
  4. Elisabeth M Busch-Nentwich
(2022)
Allele-specific gene expression can underlie altered transcript abundance in zebrafish mutants
eLife 11:e72825.
https://doi.org/10.7554/eLife.72825
  2. Download BibTeX
  3. Download .RIS

Abstract

In model organisms, RNA sequencing is frequently used to assess the effect of genetic mutations on cellular and developmental processes. Typically, animals heterozygous for a mutation are crossed to produce offspring with different genotypes. Resultant embryos are grouped by genotype to compare homozygous mutant embryos to heterozygous and wild-type siblings. Genes that are differentially expressed between the groups are assumed to reveal insights into the pathways affected by the mutation. Here we show that in zebrafish, differentially expressed genes are often overrepresented on the same chromosome as the mutation due to different levels of expression of alleles from different genetic backgrounds. Using an incross of haplotype-resolved wild-type fish, we found evidence of widespread allele-specific expression, which appears as differential expression when comparing embryos homozygous for a region of the genome to their siblings. When analysing mutant transcriptomes, this means that the differential expression of genes on the same chromosome as a mutation of interest may not be caused by that mutation. Typically, the genomic location of a differentially expressed gene is not considered when interpreting its importance with respect to the phenotype. This could lead to pathways being erroneously implicated or overlooked due to the noise of spurious differentially expressed genes on the same chromosome as the mutation. These observations have implications for the interpretation of RNA-seq experiments involving outbred animals and non-inbred model organisms.

Data availability

Sequencing data have been deposited in ENA under the accessions shown in the Materials and Methods. Differentially expressed gene lists for all the experiments are available at doi.org/10.6084/m9.figshare.15082239.

The following data sets were generated
The following previously published data sets were used

Article and author information

Author details

  1. Richard J White

    Department of Medicine, University of Cambridge, Cambridge, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.

  2. Eirinn Mackay

    Department of Cell and Developmental Biology, University College London, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.

  3. Stephen W Wilson

    Department of Cell and Developmental Biology, University College London, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8557-5940

  4. Elisabeth M Busch-Nentwich

    Department of Medicine, University of Cambridge, Cambridge, United Kingdom
    For correspondence
    e.busch-nentwich@qmul.ac.uk
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6450-744X

Funding

Medical Research Council (MR/L003775/1)

  • Stephen W Wilson

Medical Research Council (MR/T020164/1)

  • Stephen W Wilson

Wellcome Trust (095722/Z/11/Z)

  • Stephen W Wilson

Wellcome Trust (206194)

  • Richard J White
  • Elisabeth M Busch-Nentwich

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2022, White et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 2,980
    views
  • 220
    downloads
  • 10
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Richard J White
  2. Eirinn Mackay
  3. Stephen W Wilson
  4. Elisabeth M Busch-Nentwich
(2022)
Allele-specific gene expression can underlie altered transcript abundance in zebrafish mutants
eLife 11:e72825.
https://doi.org/10.7554/eLife.72825

Share this article

https://doi.org/10.7554/eLife.72825

Categories and tags

Research organism

Further reading

    1. Chromosomes and Gene Expression
    2. Genetics and Genomics

    Regulatory genome annotation of 33 insect species

    Hasiba Asma, Ellen Tieke ... Marc S Halfon
    Tools and Resources

    Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules—e.g., enhancers and silencers—that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

    1. Chromosomes and Gene Expression

    Transcriptional inhibition after irradiation occurs preferentially at highly expressed genes in a manner dependent on cell cycle progression

    Zulong Chen, Xin Wang ... Jessica K Tyler
    Research Article

    In response to DNA double-strand damage, ongoing transcription is inhibited to facilitate accurate DNA repair while transcriptional recovery occurs after DNA repair is complete. However, the mechanisms at play and the identity of the transcripts being regulated in this manner are unclear. In contrast to the situation following UV damage, we found that transcriptional recovery after ionizing radiation (IR) occurs in a manner independent of the HIRA histone chaperone. Sequencing of the nascent transcripts identified a programmed transcriptional response, where certain transcripts and pathways are rapidly downregulated after IR, while other transcripts and pathways are upregulated. Specifically, most of the loss of nascent transcripts occurring after IR is due to inhibition of transcriptional initiation of the highly transcribed histone genes and the rDNA. To identify factors responsible for transcriptional inhibition after IR in an unbiased manner, we performed a whole genome gRNA library CRISPR/Cas9 screen. Many of the top hits on our screen were factors required for protein neddylation. However, at short times after inhibition of neddylation, transcriptional inhibition still occurred after IR, even though neddylation was effectively inhibited. Persistent inhibition of neddylation blocked transcriptional inhibition after IR, and it also leads to cell cycle arrest. Indeed, we uncovered that many inhibitors and conditions that lead to cell cycle arrest in G1 or G2 phase also prevent transcriptional inhibition after IR. As such, it appears that transcriptional inhibition after IR occurs preferentially at highly expressed genes in cycling cells.

    1. Chromosomes and Gene Expression

    An anciently diverged family of RNA binding proteins maintain correct splicing of a class of ultra-long exons through cryptic splice site repression

    Chileleko Siachisumo, Sara Luzzi ... David J Elliott
    Research Advance

    Previously, we showed that the germ cell-specific nuclear protein RBMXL2 represses cryptic splicing patterns during meiosis and is required for male fertility (Ehrmann et al., 2019). Here, we show that in somatic cells the similar yet ubiquitously expressed RBMX protein has similar functions. RBMX regulates a distinct class of exons that exceed the median human exon size. RBMX protein-RNA interactions are enriched within ultra-long exons, particularly within genes involved in genome stability, and repress the selection of cryptic splice sites that would compromise gene function. The RBMX gene is silenced during male meiosis due to sex chromosome inactivation. To test whether RBMXL2 might replace the function of RBMX during meiosis we induced expression of RBMXL2 and the more distantly related RBMY protein in somatic cells, finding each could rescue aberrant patterns of RNA processing caused by RBMX depletion. The C-terminal disordered domain of RBMXL2 is sufficient to rescue proper splicing control after RBMX depletion. Our data indicate that RBMX and RBMXL2 have parallel roles in somatic tissues and the germline that must have been conserved for at least 200 million years of mammalian evolution. We propose RBMX family proteins are particularly important for the splicing inclusion of some ultra-long exons with increased intrinsic susceptibility to cryptic splice site selection.