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,969
    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. Immunology and Inflammation

    Dynamic chromatin architecture identifies new autoimmune-associated enhancers for IL2 and novel genes regulating CD4+ T cell activation

    Matthew C Pahl, Prabhat Sharma ... Andrew D Wells
    Research Article

    Genome-wide association studies (GWAS) have identified hundreds of genetic signals associated with autoimmune disease. The majority of these signals are located in non-coding regions and likely impact cis-regulatory elements (cRE). Because cRE function is dynamic across cell types and states, profiling the epigenetic status of cRE across physiological processes is necessary to characterize the molecular mechanisms by which autoimmune variants contribute to disease risk. We localized risk variants from 15 autoimmune GWAS to cRE active during TCR-CD28 co-stimulation of naïve human CD4+ T cells. To characterize how dynamic changes in gene expression correlate with cRE activity, we measured transcript levels, chromatin accessibility, and promoter–cRE contacts across three phases of naive CD4+ T cell activation using RNA-seq, ATAC-seq, and HiC. We identified ~1200 protein-coding genes physically connected to accessible disease-associated variants at 423 GWAS signals, at least one-third of which are dynamically regulated by activation. From these maps, we functionally validated a novel stretch of evolutionarily conserved intergenic enhancers whose activity is required for activation-induced IL2 gene expression in human and mouse, and is influenced by autoimmune-associated genetic variation. The set of genes implicated by this approach are enriched for genes controlling CD4+ T cell function and genes involved in human inborn errors of immunity, and we pharmacologically validated eight implicated genes as novel regulators of T cell activation. These studies directly show how autoimmune variants and the genes they regulate influence processes involved in CD4+ T cell proliferation and activation.

    1. Chromosomes and Gene Expression
    2. Developmental Biology

    OVO positively regulates essential maternal pathways by binding near the transcriptional start sites in the Drosophila female germline

    Leif Benner, Savannah Muron ... Brian Oliver
    Research Article

    Differentiation of female germline stem cells into a mature oocyte includes the expression of RNAs and proteins that drive early embryonic development in Drosophila. We have little insight into what activates the expression of these maternal factors. One candidate is the zinc-finger protein OVO. OVO is required for female germline viability and has been shown to positively regulate its own expression, as well as a downstream target, ovarian tumor, by binding to the transcriptional start site (TSS). To find additional OVO targets in the female germline and further elucidate OVO’s role in oocyte development, we performed ChIP-seq to determine genome-wide OVO occupancy, as well as RNA-seq comparing hypomorphic and wild type rescue ovo alleles. OVO preferentially binds in close proximity to target TSSs genome-wide, is associated with open chromatin, transcriptionally active histone marks, and OVO-dependent expression. Motif enrichment analysis on OVO ChIP peaks identified a 5’-TAACNGT-3’ OVO DNA binding motif spatially enriched near TSSs. However, the OVO DNA binding motif does not exhibit precise motif spacing relative to the TSS characteristic of RNA polymerase II complex binding core promoter elements. Integrated genomics analysis showed that 525 genes that are bound and increase in expression downstream of OVO are known to be essential maternally expressed genes. These include genes involved in anterior/posterior/germ plasm specification (bcd, exu, swa, osk, nos, aub, pgc, gcl), egg activation (png, plu, gnu, wisp, C(3)g, mtrm), translational regulation (cup, orb, bru1, me31B), and vitelline membrane formation (fs(1)N, fs(1)M3, clos). This suggests that OVO is a master transcriptional regulator of oocyte development and is responsible for the expression of structural components of the egg as well as maternally provided RNAs that are required for early embryonic development.

    1. Chromosomes and Gene Expression
    2. Genetics and Genomics

    BRCA1/BRC-1 and SMC-5/6 regulate DNA repair pathway engagement during Caenorhabditis elegans meiosis

    Erik Toraason, Alina Salagean ... Diana E Libuda
    Research Article Updated

    The preservation of genome integrity during sperm and egg development is vital for reproductive success. During meiosis, the tumor suppressor BRCA1/BRC-1 and structural maintenance of chromosomes 5/6 (SMC-5/6) complex genetically interact to promote high fidelity DNA double strand break (DSB) repair, but the specific DSB repair outcomes these proteins regulate remain unknown. Using genetic and cytological methods to monitor resolution of DSBs with different repair partners in Caenorhabditis elegans, we demonstrate that both BRC-1 and SMC-5 repress intersister crossover recombination events. Sequencing analysis of conversion tracts from homolog-independent DSB repair events further indicates that BRC-1 regulates intersister/intrachromatid noncrossover conversion tract length. Moreover, we find that BRC-1 specifically inhibits error prone repair of DSBs induced at mid-pachytene. Finally, we reveal functional interactions of BRC-1 and SMC-5/6 in regulating repair pathway engagement: BRC-1 is required for localization of recombinase proteins to DSBs in smc-5 mutants and enhances DSB repair defects in smc-5 mutants by repressing theta-mediated end joining (TMEJ). These results are consistent with a model in which some functions of BRC-1 act upstream of SMC-5/6 to promote recombination and inhibit error-prone DSB repair, while SMC-5/6 acts downstream of BRC-1 to regulate the formation or resolution of recombination intermediates. Taken together, our study illuminates the coordinated interplay of BRC-1 and SMC-5/6 to regulate DSB repair outcomes in the germline.