An improved zebrafish transcriptome annotation for sensitive and comprehensive detection of cell type-specific genes

  1. Nathan D Lawson  Is a corresponding author
  2. Rui Li
  3. Masahiro Shin
  4. Ann Grosse
  5. Onur Yukselen
  6. Oliver A Stone
  7. Alper Kucukural
  8. Lihua Zhu
  1. University of Massachusetts Medical School, United States
  2. University of Oxford, United Kingdom

Abstract

The zebrafish is ideal for studying embryogenesis and is increasingly applied to model human disease. In these contexts, RNA-sequencing (RNA-seq) provides mechanistic insights by identifying transcriptome changes between experimental conditions. Application of RNA-seq relies on accurate transcript annotation for a genome of interest. Here, we find discrepancies in analysis from RNA-seq datasets quantified using Ensembl and RefSeq zebrafish annotations. These issues were due, in part, to variably annotated 3' untranslated regions and thousands of gene models missing from each annotation. Since these discrepancies could compromise downstream analyses and biological reproducibility, we built a more comprehensive zebrafish transcriptome annotation that addresses these deficiencies. Our annotation improves detection of cell type-specific genes in both bulk and single cell RNA-seq datasets, where it also improves resolution of cell clustering. Thus, we demonstrate that our new transcriptome annotation can outperform existing annotations, providing an important resource for zebrafish researchers.

Data availability

All data generated in this study are available in accompanying source data files. Transcriptome annotation files described in this study are available for download at zf-transcriptome.umassmed.edu. Raw and processed RNA-seq data generated in this study are available at GEO (GSE152759).

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

Article and author information

Author details

  1. Nathan D Lawson

    Department of Molecular, Cell, and Cancer Biology, University of Massachusetts Medical School, Worcester, United States
    For correspondence
    nathan.lawson@umassmed.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7788-9619
  2. Rui Li

    Department of Molecular, Cell, and Cancer Biology, University of Massachusetts Medical School, Worcester, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Masahiro Shin

    Department of Molecular, Cell, and Cancer Biology, University of Massachusetts Medical School, Worcester, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Ann Grosse

    Department of Molecular, Cell, and Cancer Biology, University of Massachusetts Medical School, Worcester, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Onur Yukselen

    Department of Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, United States
    Competing interests
    The authors declare that no competing interests exist.
  6. Oliver A Stone

    Department of Physiology, Anatomy and Genetics, University of Oxford, Oxford, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  7. Alper Kucukural

    Department of Bioinformatic and Integrative Biology, University of Massachusetts Medical School, Worcester, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9983-394X
  8. Lihua Zhu

    Department of Molecular, Cell and Cancer Biology, University of Massachusetts Medical School, Worcester, United States
    Competing interests
    The authors declare that no competing interests exist.

Funding

National Heart, Lung, and Blood Institute (R35HL140017)

  • Nathan D Lawson

National Human Genome Research Institute (U01HG007910)

  • Onur Yukselen
  • Alper Kucukural

National Center for Advancing Translational Sciences (UL1TR001453)

  • Onur Yukselen
  • Alper Kucukural

National Institute of Neurological Disorders and Stroke (R21NS105654)

  • Nathan D Lawson

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

Ethics

Animal experimentation: Zebrafish studies were performed in accordance with protocols #A2613 and #A2632 approved by the University of Massachusetts institutional animal care and use committee (IACUC).

Copyright

© 2020, Lawson 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

  • 10,566
    views
  • 891
    downloads
  • 115
    citations

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

Download links

Share this article

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

Further reading

    1. Cell Biology
    2. Developmental Biology
    Pavan K Nayak, Arul Subramanian, Thomas F Schilling
    Research Article

    Mechanical forces play a critical role in tendon development and function, influencing cell behavior through mechanotransduction signaling pathways and subsequent extracellular matrix (ECM) remodeling. Here we investigate the molecular mechanisms by which tenocytes in developing zebrafish embryos respond to muscle contraction forces during the onset of swimming and cranial muscle activity. Using genome-wide bulk RNA sequencing of FAC-sorted tenocytes we identify novel tenocyte markers and genes involved in tendon mechanotransduction. Embryonic tendons show dramatic changes in expression of matrix remodeling associated 5b (mxra5b), matrilin1 (matn1), and the transcription factor kruppel-like factor 2a (klf2a), as muscles start to contract. Using embryos paralyzed either by loss of muscle contractility or neuromuscular stimulation we confirm that muscle contractile forces influence the spatial and temporal expression patterns of all three genes. Quantification of these gene expression changes across tenocytes at multiple tendon entheses and myotendinous junctions reveals that their responses depend on force intensity, duration and tissue stiffness. These force-dependent feedback mechanisms in tendons, particularly in the ECM, have important implications for improved treatments of tendon injuries and atrophy.

    1. Developmental Biology
    2. Genetics and Genomics
    Menglei Yang, Hafiz Muhammad Jafar Hussain ... Baolu Shi
    Research Article

    Asthenoteratozoospermia, a prevalent cause of male infertility, lacks a well-defined etiology. DNAH12 is a special dynein featured by the absence of a microtubule-binding domain, however, its functions in spermatogenesis remain largely unknown. Through comprehensive genetic analyses involving whole-exome sequencing and subsequent Sanger sequencing on infertile patients and fertile controls from six distinct families, we unveiled six biallelic mutations in DNAH12 that co-segregate recessively with male infertility in the studied families. Transmission electron microscopy (TEM) revealed pronounced axonemal abnormalities, including inner dynein arms (IDAs) impairment and central pair (CP) loss in sperm flagella of the patients. Mouse models (Dnah12-/- and Dnah12mut/mut) were generated and recapitulated the reproductive defects in the patients. Noteworthy, DNAH12 deficiency did not show effects on cilium organization and function. Mechanistically, DNAH12 was confirmed to interact with two other IDA components DNALI1 and DNAH1, while disruption of DNAH12 leads to failed recruitment of DNALI1 and DNAH1 to IDAs and compromised sperm development. Furthermore, DNAH12 also interacts with radial spoke head proteins RSPH1, RSPH9, and DNAJB13 to regulate CP stability. Moreover, the infertility of Dnah12-/- mice could be overcome by intracytoplasmic sperm injection (ICSI) treatment. Collectively, DNAH12 plays a crucial role in the proper organization of axoneme in sperm flagella, but not cilia, by recruiting DNAH1 and DNALI1 in both humans and mice. These findings expand our comprehension of dynein component assembly in flagella and cilia and provide a valuable marker for genetic counseling and diagnosis of asthenoteratozoospermia in clinical practice.