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.

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

  • 4,801
    views
  • 696
    downloads
  • 68
    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. 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
    2. Genetics and Genomics
    Keva Li, Nicholas Tolman ... UK Biobank Eye and Vision Consortium
    Research Article

    A glaucoma polygenic risk score (PRS) can effectively identify disease risk, but some individuals with high PRS do not develop glaucoma. Factors contributing to this resilience remain unclear. Using 4,658 glaucoma cases and 113,040 controls in a cross-sectional study of the UK Biobank, we investigated whether plasma metabolites enhanced glaucoma prediction and if a metabolomic signature of resilience in high-genetic-risk individuals existed. Logistic regression models incorporating 168 NMR-based metabolites into PRS-based glaucoma assessments were developed, with multiple comparison corrections applied. While metabolites weakly predicted glaucoma (Area Under the Curve = 0.579), they offered marginal prediction improvement in PRS-only-based models (p=0.004). We identified a metabolomic signature associated with resilience in the top glaucoma PRS decile, with elevated glycolysis-related metabolites—lactate (p=8.8E-12), pyruvate (p=1.9E-10), and citrate (p=0.02)—linked to reduced glaucoma prevalence. These metabolites combined significantly modified the PRS-glaucoma relationship (Pinteraction = 0.011). Higher total resilience metabolite levels within the highest PRS quartile corresponded to lower glaucoma prevalence (Odds Ratiohighest vs. lowest total resilience metabolite quartile=0.71, 95% Confidence Interval = 0.64–0.80). As pyruvate is a foundational metabolite linking glycolysis to tricarboxylic acid cycle metabolism and ATP generation, we pursued experimental validation for this putative resilience biomarker in a human-relevant Mus musculus glaucoma model. Dietary pyruvate mitigated elevated intraocular pressure (p=0.002) and optic nerve damage (p<0.0003) in Lmx1bV265D mice. These findings highlight the protective role of pyruvate-related metabolism against glaucoma and suggest potential avenues for therapeutic intervention.

    1. Cell Biology
    Affiong Ika Oqua, Kin Chao ... Alejandra Tomas
    Research Article

    G protein-coupled receptors (GPCRs) are integral membrane proteins which closely interact with their plasma membrane lipid microenvironment. Cholesterol is a lipid enriched at the plasma membrane with pivotal roles in the control of membrane fluidity and maintenance of membrane microarchitecture, directly impacting on GPCR stability, dynamics, and function. Cholesterol extraction from pancreatic beta cells has previously been shown to disrupt the internalisation, clustering, and cAMP responses of the glucagon-like peptide-1 receptor (GLP-1R), a class B1 GPCR with key roles in the control of blood glucose levels via the potentiation of insulin secretion in beta cells and weight reduction via the modulation of brain appetite control centres. Here, we unveil the detrimental effect of a high cholesterol diet on GLP-1R-dependent glucoregulation in vivo, and the improvement in GLP-1R function that a reduction in cholesterol synthesis using simvastatin exerts in pancreatic islets. We next identify and map sites of cholesterol high occupancy and residence time on active vs inactive GLP-1Rs using coarse-grained molecular dynamics (cgMD) simulations, followed by a screen of key residues selected from these sites and detailed analyses of the effects of mutating one of these, Val229, to alanine on GLP-1R-cholesterol interactions, plasma membrane behaviours, clustering, trafficking and signalling in INS-1 832/3 rat pancreatic beta cells and primary mouse islets, unveiling an improved insulin secretion profile for the V229A mutant receptor. This study (1) highlights the role of cholesterol in regulating GLP-1R responses in vivo; (2) provides a detailed map of GLP-1R - cholesterol binding sites in model membranes; (3) validates their functional relevance in beta cells; and (4) highlights their potential as locations for the rational design of novel allosteric modulators with the capacity to fine-tune GLP-1R responses.