Single-cell RNA-seq analysis reveals penaeid shrimp hemocyte subpopulations and cell differentiation process

  1. Keiichiro Koiwai  Is a corresponding author
  2. Takashi Koyama
  3. Soichiro Tsuda
  4. Atsushi Toyoda
  5. Kiyoshi Kikuchi
  6. Hiroaki Suzuki
  7. Ryuji Kawano
  1. Tokyo University of Marine Science and Technology, Japan
  2. Nagasaki University, Japan
  3. bitBiome Inc., Japan
  4. National Institute of Genetics, Japan
  5. University of Tokyo, Japan
  6. Chuo University, Japan
  7. Tokyo University of Agriculture and Technology, Japan

Abstract

Crustacean aquaculture is expected to be a major source of fishery commodities in the near future. Hemocytes are key players of the immune system in shrimps; however, their classification, maturation, and differentiation are still under debate. To date, only discrete and inconsistent information on the classification of shrimp hemocytes has been reported, showing that the morphological characteristics are not sufficient to resolve their actual roles. Our present study using single-cell RNA sequencing, revealed six types of hemocytes of Marsupenaeus japonicus based on their transcriptional profiles. We identified markers of each subpopulation and predicted the differentiation pathways involved in their maturation. We also predicted cell growth factors that might play crucial roles in hemocyte differentiation. Different immune roles among these subpopulations were suggested from the analysis of differentially expressed immune-related genes. These results provide a unified classification of shrimp hemocytes, which improves the understanding of its immune system.

Data availability

Sequencing data have been deposited in DDBJ under accession codes DRA010948, DRA010949, DRA010950, DRA010951, and DRA010952.Digital expression data of Drop-seq from three shrimp were archived in DDBJ under accession code E-GEAD-403.Data code can be accessed at https://github.com/KeiichiroKOIWAI/Drop-seq_on_shrimp.All data generated or analyzed during this study are included in the manuscript and supporting files.All tables are provided as Supplementary files.Source data files are provided to support the sinaplots in Figure 2, Figure 2-figure supplement 1.Source data files are provided to support thepercentage of cell state in Figure 5.Source data file is provided to support the bar graph in Figure 9F.

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

Article and author information

Author details

  1. Keiichiro Koiwai

    Laboratory of genome Science, Tokyo University of Marine Science and Technology, Tokyo, Japan
    For correspondence
    koiwai@kaiyodai.ac.jp
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8890-8229
  2. Takashi Koyama

    Graduate School of Fisheries and Environmental Sciences, Nagasaki University, Nagasaki, Japan
    Competing interests
    No competing interests declared.
  3. Soichiro Tsuda

    Research and Development Division, bitBiome Inc., Tokyo, Japan
    Competing interests
    Soichiro Tsuda, S.T. is an employee of bitBiome inc., and does not own any stock of the company..
  4. Atsushi Toyoda

    Department of Genomics and Evolutionary Biology, National Institute of Genetics, Mishima, Japan
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0728-7548
  5. Kiyoshi Kikuchi

    Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, University of Tokyo, Hamamatsu, Japan
    Competing interests
    No competing interests declared.
  6. Hiroaki Suzuki

    Department of Precision Mechanics, Chuo University, Tokyo, Japan
    Competing interests
    No competing interests declared.
  7. Ryuji Kawano

    Department of Biotechnology and Life Science, Tokyo University of Agriculture and Technology, Tokyo, Japan
    Competing interests
    No competing interests declared.

Funding

Japan Society for the Promotion of Science (20K15603)

  • Keiichiro Koiwai

Japan Society for the Promotion of Science (19J00539)

  • Keiichiro Koiwai

Japan Society for the Promotion of Science (17H06425)

  • Kiyoshi Kikuchi

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

Reviewing Editor

  1. Jiwon Shim, Hanyang University, Republic of Korea

Version history

  1. Received: January 28, 2021
  2. Accepted: June 15, 2021
  3. Accepted Manuscript published: June 16, 2021 (version 1)
  4. Version of Record published: July 8, 2021 (version 2)

Copyright

© 2021, Koiwai 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

  • 3,880
    Page views
  • 569
    Downloads
  • 46
    Citations

Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, PubMed Central.

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. Keiichiro Koiwai
  2. Takashi Koyama
  3. Soichiro Tsuda
  4. Atsushi Toyoda
  5. Kiyoshi Kikuchi
  6. Hiroaki Suzuki
  7. Ryuji Kawano
(2021)
Single-cell RNA-seq analysis reveals penaeid shrimp hemocyte subpopulations and cell differentiation process
eLife 10:e66954.
https://doi.org/10.7554/eLife.66954

Share this article

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

Further reading

    1. Cell Biology
    Kazuki Hanaoka, Kensuke Nishikawa ... Kouichi Funato
    Research Article

    Membrane contact sites (MCSs) are junctures that perform important roles including coordinating lipid metabolism. Previous studies have indicated that vacuolar fission/fusion processes are coupled with modifications in the membrane lipid composition. However, it has been still unclear whether MCS-mediated lipid metabolism controls the vacuolar morphology. Here, we report that deletion of tricalbins (Tcb1, Tcb2, and Tcb3), tethering proteins at endoplasmic reticulum (ER)–plasma membrane (PM) and ER–Golgi contact sites, alters fusion/fission dynamics and causes vacuolar fragmentation in the yeast Saccharomyces cerevisiae. In addition, we show that the sphingolipid precursor phytosphingosine (PHS) accumulates in tricalbin-deleted cells, triggering the vacuolar division. Detachment of the nucleus–vacuole junction (NVJ), an important contact site between the vacuole and the perinuclear ER, restored vacuolar morphology in both cells subjected to high exogenous PHS and Tcb3-deleted cells, supporting that PHS transport across the NVJ induces vacuole division. Thus, our results suggest that vacuolar morphology is maintained by MCSs through the metabolism of sphingolipids.

    1. Cell Biology
    2. Chromosomes and Gene Expression
    Monica Salinas-Pena, Elena Rebollo, Albert Jordan
    Research Article

    Histone H1 participates in chromatin condensation and regulates nuclear processes. Human somatic cells may contain up to seven histone H1 variants, although their functional heterogeneity is not fully understood. Here, we have profiled the differential nuclear distribution of the somatic H1 repertoire in human cells through imaging techniques including super-resolution microscopy. H1 variants exhibit characteristic distribution patterns in both interphase and mitosis. H1.2, H1.3, and H1.5 are universally enriched at the nuclear periphery in all cell lines analyzed and co-localize with compacted DNA. H1.0 shows a less pronounced peripheral localization, with apparent variability among different cell lines. On the other hand, H1.4 and H1X are distributed throughout the nucleus, being H1X universally enriched in high-GC regions and abundant in the nucleoli. Interestingly, H1.4 and H1.0 show a more peripheral distribution in cell lines lacking H1.3 and H1.5. The differential distribution patterns of H1 suggest specific functionalities in organizing lamina-associated domains or nucleolar activity, which is further supported by a distinct response of H1X or phosphorylated H1.4 to the inhibition of ribosomal DNA transcription. Moreover, H1 variants depletion affects chromatin structure in a variant-specific manner. Concretely, H1.2 knock-down, either alone or combined, triggers a global chromatin decompaction. Overall, imaging has allowed us to distinguish H1 variants distribution beyond the segregation in two groups denoted by previous ChIP-Seq determinations. Our results support H1 variants heterogeneity and suggest that variant-specific functionality can be shared between different cell types.