A single-cell profile of the porcine immune repertoire: concurrent profiling of T cell receptor and B cell receptor V(D)J sequences in health

Abstract

Adaptive immune diversity is shaped by the complex organization of T cell receptor (TCR) and B cell receptor (BCR) repertoires across tissues. However, a systematic cross-tissue immune repertoire profiling in healthy pigs is still lacking. In this study, we used 3 14-day-old Duroc × Landrace × Yorkshire (DLY) crossbred piglets, all of which were male. After sample purification, a total of 71,731 cells were retained. We established single-cell profiling integrating paired BCR and TCR repertoires across 11 immune-associated tissues using single-cell RNA sequencing combined with V(D)J sequencing based on DLY crossbred piglets. Diverse populations of immune cells were identified, with significant tissue-specific heterogeneity being revealed in lymphatic, mucosal, and circulatory organs. Immune repertoire analysis demonstrated biased V(D)J gene usage, conserved complementarity-determining region 3 (CDR3) length distributions, preferential receptor pairing patterns, and extensive clonotype diversity. Cross-tissue analyses further revealed coordinated clonotype sharing and distinct clonal expansion patterns among immune cell subsets. In addition, using the Pig-specific V(D)J primers we designed, the sequencing results have achieved a modest improvement in the detection of immunoglobulin alleles compared to previous studies. It is expected to provide a more reliable technical solution for related research. Collectively, this work provides a comprehensive reference profiling of the healthy DLY crossbred piglet adaptive immune system and offers valuable resources for future studies of porcine immunology, infectious diseases, and vaccine responses.

Keywords

ScRNA-seq / V(D)J sequencing / Immune repertoire / Pig

References

[1] Baumgartner CK, Malherbe LP. Antigen-driven T-cell repertoire selection during adaptive immune responses. Immunol Cell Biol, 2011, 89: 54-59.

[2] Berruezo-Llacuna M, Kabrani E, Di Virgilio M. The B cell dilemma: diversity or fidelity?. DNA Repair, 2025, 153. 103888

[3] Bram S, et al. Parallel single B cell transcriptomics to elucidate pig B cell repertoire. Sci Rep, 2024, 14. 15997

[4] Briney B, et al. Commonality Despite Exceptional Diversity in the Baseline Human Antibody Repertoire. Nature, 2019, 566: 393-397.

[5] Buggert M, et al. The identity of human tissue-emigrant CD8+T cells. Cell, 2020, 183: 1946-1961.e15.

[6] Burton OT, et al. The tissue-resident regulatory T cell pool is shaped by transient multi-tissue migration and a conserved residency program. Immunity, 2024, 57715861602.e10.

[7] Butler JE, et al. Porcine IgG: structure, genetics, and evolution. Immunogenetics, 2009, 61: 209-230.

[8] Butler JE, Wertz N, Sinkora M. Antibody Repertoire Development in Swine. Annual Review of Animal Biosciences, 2017, 5: 255-279.

[9] Butler A, et al. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol, 2018, 36: 411-420.

[10] Chao A. Nonparametric-Estimation of the Number of Classes in A Population. Scand J Stat, 1984, 11: 265-270

[11] Chao A, Lee SM. Estimating the Number of Classes Via Sample Coverage. J Am Stat Assoc, 1992, 87417210217.

[12] Chen N, Trible BRR. Amplification and selection of PRRSV-activated VDJ repertoires in pigs secreting distinct neutralizing antiboidies. Vet Immunol Immunopathol, 2017, 189: 53-57.

[13] Chi H, Pepper M, Thomas PG. Principles and therapeutic applications of adaptive immunity. Cell, 2024, 187920522078.

[14] Dunn-Walters D, et al. Immunoglobulin gene analysis as a tool for investigating human immune responses. Immunol Rev, 2018, 2841132147.

[15] Franzen O, Gan LM, Bjorkegren JLM. PanglaoDB: a web server for exploration of mouse and human single-cell RNA sequencing data. Database, 2019, 20191baz046

[16] Graf R, et al. BCR-dependent lineage plasticity in mature B cells. Science, 2019, 363: 748-753.

[17] Gu W, et al. Single-cell antigen receptor sequencing in pigs with influenza. Communications Biology, 2025, 811108.

[18] Hafemeister C, Satija R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biol, 2019, 20. 296

[19] Hammer SE, et al. Development of a RACE-based RNA-Seq approach to characterize the T-cell receptor repertoire of porcine Γδ T cells. Dev Comp Immunol, 2019, 105: 103575-103575.

[20] Herrera-Uribe J, et al. Reference transcriptomes of porcine peripheral immune cells created through bulk and single-cell RNA sequencing. Front Genet, 2021, 12. 689406

[21] Hill MO. Diversity and evenness: a unifying notation and its consequences. Ecology, 1973, 54: 427-432.

[22] Hobeika E, Maity PC, Jumaa H. Control of B cell responsiveness by isotype and structural elements of the antigen receptor. Trends Immunol, 2016, 375310320.

[23] Hong Y, et al. In-depth analysis of human neonatal and adult IgM antibody repertoires. Front Immunol, 2018, 9. 128

[24] Jiang S, et al. Generic Diagramming Platform (GDP): a Comprehensive Database of High-Quality Biomedical Graphics. Nucleic Acids Res, 2025, 53: D1670D1676.

[25] Kekarainen T, et al. Immune responses and vaccine-induced immunity against Porcine circovirus type 2. Vet Immunol Immunopathol, 2010, 1363–4185193.

[26] Kitaura T, et al. A new high-throughput sequencing method for determining diversity and similarity of T cell receptor (TCR) α and β repertoires and identifying potential new invariant TCR α chains. BMC Immunol., 2016, 17138.

[27] Knox RV. Worldwide perspective for swine production and reproduction for the next 20 years. Theriogenology, 2025, 234: 24-33.

[28] Le Page L, Baldwin CL, Telfer JC. γδ T cells in artiodactyls: Focus on swine. Dev Comp Immunol, 2022, 128. 104334

[29] Lepore A, et al. Parallel T-cell cloning and deep sequencing of human MAIT cells reveal stable oligoclonal TCRβ repertoire. Nat Commun, 2014, 5: 3866.

[30] Li J, et al. Single-cell transcriptomic analysis reveals transcriptional and cell subpopulation differences between human and pig immune cells. Genes Genomics, 2024, 463303322.

[31] Mandric I, et al. Profiling immunoglobulin repertoires across multiple human tissues using RNA sequencing. Nat Commun, 2020, 1113126

[32] Mark M, et al. A hierarchy of selection pressures determines the organization of the T cell receptor repertoire. Front Immunol, 2022, 13. 939394

[33] Mika J, et al. A comprehensive evaluation of diversity measures for TCR repertoire profiling. BMC Biol, 2025, 231133

[34] Muñoz-Rojas AR, Mathis D. Tissue regulatory T cells: regulatory chameleons. Nat Rev Immunol, 2021, 219597611.

[35] Pai JA, Satpathy AT. High-throughput and single-cell T cell receptor sequencing technologies. Nat Methods, 2021, 188881892.

[36] Partners, C.N.M.A.Database resources of the National Genomics Data Center, China National Center for Bioinformation in 2026. Nucleic Acids Res, 2026, 54: D28D47.

[37] Patil V, et al. Characterization of a novel functional porcine CD3(+)CD4(low)CD8α(+)CD8β(+) T-helper/memory lymphocyte subset in the respiratory tract lymphoid tissues of swine influenza A virus vaccinated pigs. Vet Immunol Immunopathol, 2024, 274. 110785

[38] Pelissier A, et al. Convergent evolution and B-cell recirculation in germinal centers in a human lymph node. Life Sci Alliance, 2023, 611e202301959

[39] Qiu Y, et al. Differential analysis of immunoglobulin gene expression pattern in chickens of distinct breeds and developmental periods. J Anim Sci, 2024, 102. skae111

[40] Roth DB. V(D)J recombination: mechanism, errors, and fidelity. Microbiol Spectr, 2014, 26MDNA3-0041-2014

[41] Schad SE, et al. Tumor-induced double positive T cells display distinct lineage commitment mechanisms and functions. J Exp Med, 2022, 2196e20212169

[42] Schwartz JC, Murtaugh MP. Viral Infection Skews Immunoglobulin Light Chain Repertoire Diversity. BioRxiv, 2018.

[43] Scott MC, et al. Deep profiling deconstructs features associated with memory CD8+ T cell tissue residence. Immunity, 2024, 58: 162-181.e10.

[44] Seo K, Choi JK. Comprehensive Analysis of TCR and BCR Repertoires: Insights into Methodologies, Challenges, and Applications. Genomics & Informatics, 2025, 2316.

[45] Silva SL, et al. Human naïve regulatory T-cells feature high steady-state turnover and are maintained by IL-7. Oncotarget, 2016, 7111216312175.

[46] Silver J, et al. Stochasticity enables BCR-independent germinal center initiation and antibody affinity maturation. J Exp Med, 2017, 215: 77-90.

[47] Stepanova K, et al. The order of immunoglobulin light chain Κ and Λ usage in primary and secondary lymphoid tissues of germ-free and conventional piglets. Dev Comp Immunol, 2022, 131. 104392

[48] Stuart T, et al. Comprehensive integration of single cell data. Cell, 2019, 177: 1888-1902.e21.

[49] Sureshchandra S, et al. Tissue Determinants of the Human T Cell Receptor Repertoire. BioRxiv, 2024.

[50] Thai C, et al. CANTAO: guiding clustering and annotation in single-cell RNA sequencing using average overlap. Mol Syst Biol, 2026, 223461475.

[51] VanderWaal K, Deen J. Global trends in infectious diseases of swine. Proc Natl Acad Sci U S A, 2018, 115451149511500.

[52] Weisberg SP, Ural BB, Farber DL. Tissue-specific immunity for a changing world. Cell, 2021, 184615171529.

[53] Wiarda J, et al. A single-cell immune atlas of primary and secondary lymphoid organs in pigs. Front Immunol, 2026, 17: 1704257.

[54] Wu M, et al. Immune repertoire: revealing the "real-time" adaptive immune response in autoimmune diseases. Autoimmunity, 2021, 5426175.

[55] Wu X, et al. Single-cell Sequencing to Multi-Omics: Technologies and Applications. Biomarker Research, 2024, 12: 110.

[56] Xia Y, et al. Development of a CRISPR-Cas12a based assay for the detection of swine enteric coronaviruses in pig herds in China. Adv Biotechnol, 2024, 2: 7.

[57] Yan S, et al. A Huntingtin Knockin Pig Model Recapitulates Features of Selective Neurodegeneration in Huntington's Disease. Cell, 2018, 17349891002.e13.

[58] Yan Z, et al. Impact of African swine fever epidemic on the cost intensity of pork production in China. Agriculture, 2023, 132497

[59] Yang Q, et al. Screpertoire 2: enhanced and efficient toolkit for single-cell immune profiling. PLoS Comput Biol, 2025, 21. e1012760

[60] Zhang X, et al. CellMarker: a manually curated resource of cell markers in human and mouse. Nucleic Acids Res, 2018, 47: D721D728.

[61] Zhang M, et al. Revisiting the pig IGHC gene locus in different breeds uncovers nine distinct IGHG genes. J Immunol, 2020, 205821372145.

[62] Zhang Z, et al. The Economic Impact of Porcine Reproductive and Respiratory Syndrome Outbreak in Four Chinese Farms: Based on Cost and Revenue Analysis. Frontiers in Veterinary Science, 2022, 9: 1024720.

[63] Zhang S, et al. The GSA family in 2025: a broadened sharing platform for multi-omics and multimodal data. Genomics Proteomics Bioinformatics, 2025, 234qzaf072

[64] Zhu T, et al. Single-cell omics in plant biology: mechanistic insights and applications for crop improvement. Adv Biotechnol, 2025, 3. 20

[65] Zou Q, et al. A single-cell transcriptome atlas of pig skin characterizes anatomical positional heterogeneity. Elife, 2023, 12. e86504

Funding

Key Technologies Research and Development Program(2025YFD1800901)

Rights & permissions

The Author(s)