Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Immunol.

Sec. Comparative Immunology

High-resolution single-cell RNA sequencing using canFam4 reveals novel immune subsets and checkpoint programs in healthy dogs

Provisionally accepted
Myung-Chul  KimMyung-Chul Kim1*Taeeun  GuTaeeun Gu2Hyeewon  SeoHyeewon Seo2Yewon  MoonYewon Moon2Nicholas  BorcherdingNicholas Borcherding3Ryan  KolbRyan Kolb4Youngmin  YunYoungmin Yun1Woo-Jin  SongWoo-Jin Song1Chung-Young  LeeChung-Young Lee5Hyun Je  KimHyun Je Kim2Weizhou  ZhangWeizhou Zhang4
  • 1Jeju National University, Jeju City, Republic of Korea
  • 2Seoul National University, Gwanak-gu, Republic of Korea
  • 3Washington University in St Louis School of Medicine, St. Louis, United States
  • 4University of Florida, Gainesville, United States
  • 5Kyungpook National University, Daegu, Republic of Korea

The final, formatted version of the article will be published soon.

Single-cell RNA sequencing (scRNA-seq) enables high-resolution profiling of immune heterogeneity. Although previous studies have mapped the single-cell transcriptomic atlases of peripheral leukocytes in healthy dogs, the identification and functional characterization of distinct immune subsets remain incomplete. We constructed a single-cell atlas of peripheral leukocytes from six healthy small-breed dogs using the 10x Genomics platform and the updated canFam4 genome. Analysis of 30,040 high-quality transcriptomes revealed 51 distinct immune subsets, including CD14+CD33+ monocytes, XCR1+CD1D+ dendritic cells, CEACAM1+CD24+ neutrophils, and IL32+BATF+ regulatory T cells, which were underrepresented in canFam3.1-based studies. Interferon-enriched CD14+ monocytes and CD4+ T subsets associated with myxomatous mitral valve disease were also identified. Functional enrichment analyses suggested that PDCD1 is associated with attenuated TCR signaling, whereas LAG3 was associated with malate metabolism pathways in CD4+ T cells and reduced TBX21 expression in CD8+ T cells linked to antiviral responses. CD274, which encodes PD-L1 was linked to IL-10 production in neutrophils, and CTLA4 represented an initial activation of double-negative T subsets. T cell exhaustion scores and proliferative fractions varied across cohorts, reflecting differences in environmental antigenic exposures. To our knowledge, this study represents the first comprehensive, gene-resolved single-cell analysis to reveal immunoregulatory checkpoint mechanisms underlying immune homeostasis in healthy dogs. Our dataset will serve as a valuable resource for future comparative and translational immunology research in dogs.

Keywords: canFam4, Circulating leukocytes, Reference transcriptome, ScRNA-seq, translational model

Received: 06 Aug 2025; Accepted: 24 Nov 2025.

Copyright: © 2025 Kim, Gu, Seo, Moon, Borcherding, Kolb, Yun, Song, Lee, Kim and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

* Correspondence: Myung-Chul Kim

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.