Immunogenetics

Turbot (Scophthalmus maximus) is an important aquaculture species, but the genomic organization and expressed diversity of its antibody repertoire remain incompletely characterized. In this study, we annotated the immunoglobulin heavy chain (IGH) locus using the haplotype resolved fScoMax1.1 genome assembly, and we used this as a reference to profile the expressed turbot IgM, IgD and IgT repertoires in skin and spleen. The primary IGH locus was located on chromosome 19, spanned approximately 72 kb, and contained 25 IGHV genes, including 24 functional genes and one pseudogene, together with three IGHD, seven IGHJ and three IGHC genes corresponding to IgT, IgM and IgD. Comparison with the alternate fScoMax1.1 haplotype and a second turbot genome assembly showed conserved IGHD, IGHJ and IGHC content, whereas IGHV gene number differed among assemblies. High throughput 5RACE repertoire sequencing revealed isotype and tissue associated differences in expressed IGH diversity. IgM represented the dominant productive repertoire in both skin and spleen and showed the highest clonotypic diversity, particularly in spleen. IgD displayed an intermediate profile, whereas IgT was more enriched in skin and exhibited the strongest clonal restriction. IGHV subgroup usage was dominated by IGHV3 in IgM and IgD, whereas IgT showed a distinct profile characterized by preferential use of IGHV4, especially in skin. Gene level analysis further showed broad IGHV-IGHJ pairing in IgM and IgD, with preferential use IGHJ3 segment, while IgT sequences paired exclusively with IGHJT. Clonotype sharing between skin and spleen was isotype dependent, being strongest for IgT, intermediate for IgM, and negligible for IgD, suggesting that clonal expansion did not necessarily predict inter tissue trafficking. Together, these results provide a curated genomic and expressed repertoire framework for turbot IGH genes and reveal isotype specific organization of antibody diversity, with IgT displaying a particular repertoire pattern.

TAPBP is a key chaperone of the peptide-loading complex that facilitates peptide loading onto major histocompatibility complex class I (MHC I) molecules. This study characterized TAPBP alleles in Korean Native Chickens (KNCs), identified novel variants, and evaluated haplotypic associations with BF2. Thirty-six samples representing six KNC lines were genotyped using LEI0258 and the MHC-B SNP panel, and individuals homozygous at both markers were classified into 16 groups. The same samples were subjected to Sanger sequencing of TAPBP exons 3-8. Sequences were assembled and aligned against MHC-B reference haplotypes and the Red Junglefowl reference. Additional comparisons with "tapasin allele" datasets enabled the identification of novel variants. Six novel nucleotide variants were detected across exons 3-6, including one nonsynonymous substitution in exon 4 (D251H). This residue corresponds to position Q265 in human TAPBP and lies adjacent to residues involved in MHC I interaction, suggesting potential functional relevance. Furthermore, TAPBP exhibited high haplotype diversity (Hd = 0.93) and moderate nucleotide diversity ({pi} = 0.00892), with exon 5 showing the highest diversity ({pi} = 0.01). B9 was the most frequent haplotype at the nucleotide level, whereas B6/B24 predominated at the amino acid level. Comparison with BF2 data revealed haplotype-dependent pairing patterns: BF2-B9 consistently matched TAPBP-B9, whereas BF2-B6 was associated with distinct TAPBP nucleotide variants, indicating allelic diversification within a shared haplotypic background. Homozygosity at LEI0258 and the SNP panel corresponded with TAPBP homozygosity, supporting marker-based prediction. These findings highlight potential BF2-TAPBP associations and provide a foundation for understanding variation in MHC I peptide loading.

NLRP3 is a cytosolic pattern recognition receptor that controls the formation of an inflammasome, a multimolecular complex that cleaves the pro-inflammatory cytokines interleukin-1{beta} and interleukin-18 into their bioactive forms. NLRP3 has been widely assumed to be conserved across vertebrates, suggesting it plays an indispensable role within the vertebrate innate immune system. Here, we used gene synteny, phylogeny, and structural analysis to examine the evolutionary origins of NLRP3 in greater detail. Our analysis revealed that modern NLRP3, defined by gene synteny and structural features, is unique to the eutherian lineage. Non-eutherian NLRP genes do not share synteny with the eutherian NLRP3 locus and lack conservation of key features including disc forming residues, cage interfaces, and membrane binding regions. NLRP3s characteristic regulatory architecture therefore appears to have evolved after eutherians split from marsupials 160 million years ago. These findings have important implications for understanding the beneficial roles of NLRP3 signalling and suggest that enhanced regulatory control of NLRP3 activity arose in response to distinct aspects of eutherian physiology.

Naturally occurring variation in the bovine mannose receptor C-type 1 gene (MRC1) may shape macrophage responses to Mycobacterium (M.) bovis, a key driver of bovine tuberculosis (bTB). We identified four coding region SNPs in MRC1 across Bos taurus (Holstein Friesian, Brown Swiss) and Bos indicus (Boran, Sahiwal) cattle breeds, including a non-synonymous variant, rs380943118 (c.2963G>A; Ser988Asn) in C-type lectin-like domain (CTLD) 6, most prevalent in Sahiwal cattle. Structural modelling suggested that the S988N substitution, which is spatially separated from the monosaccharide binding site of CTLD4, might indirectly affect glycan binding, perhaps through a conformational change in the receptor. Monocyte-derived macrophages upregulated MR expression during differentiation, with heterozygous (G/A) animals showing higher MR expression and increased uptake of GFP-M. bovis BCG, although differences were not statistically significant. Anti-CD206 blockade did not inhibit BCG internalization, either indicating that this specific antibody did not bind to a CTLD involved in ligand binding or that MR is not the sole entry receptor. These results highlight naturally occurring MRC1 polymorphisms that may influence MR structure and macrophage function, providing a foundation for future studies to assess their role in bTB susceptibility.

Black-footed ferrets (Mustela nigripes) are highly susceptible to sylvatic plague caused by Yersinia pestis, but the genetic basis of this vulnerability remains poorly understood. Here, comparative immunogenomic analyses across Carnivora species identified a conserved class of immunoglobulin lambda variable (IGLV) genes with unusually long antigen-binding sites (CDRL1) that are common among Caniformia species but absent in Feliformia species. First discovered in the domestic ferret (Mustela putorius furo), these genes encode tyrosine-rich and anionic motifs resembling the chemokine receptor CCR5 and contain experimentally validated sulfotyrosines previously associated with pathogen-interacting interfaces. Evolutionary analyses revealed distinct selective pressures across Caniformia lineages and showed strong purifying selection acting on long-CDRL1 IGLV genes in mustelids and bears. Antibody repertoire sequencing demonstrated that these genes are actively utilized in expressed repertoires and that their usage correlates with evolutionary conservation. Functional analyses of monoclonal antibodies derived from the long-CDRL1 IGLV gene identified an antibody that significantly reduced intracellular Y. pestis survival in macrophages and revealed a positive correlation between anti-plague activity and sulfotyrosine signal. Notably, all analyzed black-footed ferrets carried a frameshifting deletion in the long-CDRL1 IGLV gene resulting in loss of its expression in antibody repertoires. Together, these findings uncover a germline-encoded immunoglobulin feature conserved across dog-like carnivores and suggest a potential link between antibody germline variation and immune responses to plague.

We generated the first reference-level nuclear genome assembly of the keystone Caribbean long-spined black sea urchin species, Diadema antillarum (Philippi, 1845). Using whole-genome sequencing data from PacBio HiFi, Oxford Nanopore, and Illumina platforms, we employed multiple assembly strategies to generate a high-quality, near-complete genome. The final assembly spans 1.73 Gbp, consists of 2,964 scaffolds, and has an N50 of 1.56 Mbp. BUSCO analysis (metazoa_odb10) indicates 98.4% completeness. The genome displays a heterozygosity rate of 2.52% and contains 42.85% repetitive elements, of which 29.96% are unclassified. Coverage analysis reveals that while most of the genome was assembled at 11x depth, certain regions exhibit up to 530x coverage. Notably, regions exceeding 33x coverage account for 30.53% of the repetitive content, suggesting localized expansion of repeats. Duplication analysis of the assembled contigs shows that approximately 66% of contigs have duplicated, which supports segmental genome duplication in the past, and is further evidenced by the moderate level of heterozygosity of the assembly. While these characteristics contribute to the complexity of the genome, they do not diminish the quality of our assembly. Despite this complexity, our assembly maintains high completeness and contiguity. Our assembly provides a valuable resource for future genetic studies and serves as a critical framework for conservation, monitoring, and restoration of D. antillarum populations across the Caribbean.

Hemoglobinopathies are inherited disorders of hemoglobin, most notably sickle cell anemia and thalassemia. These conditions result from mutations in the globin genes, leading either to structural abnormalities in the globin chains or to reduced synthesis of normal globin chains. Hemoglobinopathies is a worldwide health problem according to the World Health Organization; it affects mostly the indigenous Tharu groups in Nepal. Both the global and local rates of illness and death associated with these diseases are on the rise. The objective of this study was to assess the presence of hemoglobinopathies and common mutations of the beta-globin gene within the Tharu population in western Nepal. A cross-sectional study of 1,400 Tharu individuals was conducted among individuals obtained through hospitals within the Banke district, Bardiya district, and Kailali district in western Nepal. A thorough hematological analysis was done with the use of a Sysmex XN-350 analyzer. Hemoglobin variants were detected via high-performance liquid chromatography (HPLC). The molecular characterization of the seven most common mutations of {beta}-thalassemia was performed on a subset of 20 confirmed cases by using a real-time PCR kit.The total number of cases diagnosed with hemoglobinopathies was 14.43% (n=202 out of 1,400). Sickle cell trait (HbAS) was reported as the most prevalent type of Hemoglobinopathies (8.50% of population), followed by {beta}-thalassemia trait (4.00%). In addition to these disorders were sickle cell disease (HbSS), HbE trait, and compound heterozygous states. Hematological parameters differed significantly across types of hemoglobinopathies, and the patterns of microcytic, hypochromic, and hemolytic anemia were also distinct. Commonly documented symptoms included fatigue and joint pain (42.5% and 23.1%, respectively). Molecular characterization of {beta}-thalassemia cases demonstrated that most individuals were compound heterozygotes with IVS1-6 (T>C) as the most prevalent variant. The research identified that the Tharu population in western Nepal has a significant burden of hemoglobinopathies (especially sickle cell trait and {beta}-thalassemia), highlighting the requirement for appropriate screening programs, genetic counseling and public health strategies to help manage and prevent these conditions within this particular region.

Defining HLA restriction of T cell epitopes is essential for understanding immune responses in infectious disease, autoimmunity, and vaccine design. Current bioinformatic programs, including the IEDB RATE tool, enable inference of single-HLA restrictions from immune response data of HLA-typed donors. However, T cell epitopes are frequently presented by multiple HLA alleles, a phenomenon termed promiscuous restriction, limiting the utility of single-allele approaches. To address this limitation, we developed COMBO-RATE, an extension of RATE that systematically evaluates combinations of HLA alleles to identify multi-allelic restriction patterns. Analysis of three independent datasets spanning distinct antigen systems and different epitope discovery strategies revealed that promiscuous restriction is a near-universal feature of immunodominant epitopes. Focusing on the 43 immunodominant CD4 T cell epitopes identified in a B. pertussis genome-wide screen, COMBO-RATE outperformed conventional RATE, identifying restrictions for 35 of 43 epitopes, compared to 24 by RATE alone, and uncovered 64 additional allele restrictions, including 29 unique alleles. Experimental validation using single-HLA transfected cell lines and antigen presentation assays confirmed COMBO-RATE-inferred restrictions, demonstrating that a single epitope can be independently presented by distinct HLA alleles. Overall, COMBO-RATE provides a robust and scalable framework for defining complete HLA restriction profiles from existing population response data, with important implications for the design of vaccines requiring broad HLA coverage across genetically diverse populations. This pipeline is available as both a Python package and a user-friendly web application.

Complex multi-tissue regeneration capacity varies across vertebrates. Mammals are amongst the least regenerative species, while salamanders can regenerate complex tissues such as limbs and tails throughout life. Previous studies have shown that innate and adaptive immune cells are present during salamander limb regeneration. While innate immune cells have been shown to promote limb regeneration, it is unknown whether adaptive immunity is responsive to amputation or plays a role in appendage regeneration. Here we show that during limb regeneration in axolotls, the immune response is characterized by a coordinated immunoregulatory signature including the downregulation of antigen presentation, cytokine secretion, and T cell activation. To test the role of adaptive immune cells in regeneration, we generated Recombination activating gene 1 deficient (Rag1-/-) newts. Rag1-/- newts lack antigen receptor recombination and show a marked reduction of adaptive immune cells. We find that Rag1-/- newts do not reject allografts, confirming their functional immunodeficiency. Finally, we demonstrate that both larval and adult newts regenerate appendages in the absence of adaptive immunity. Our work demonstrates that the adaptive arm of the immune system is not required for appendage regeneration and establishes an important model for novel experimental approaches in comparative immunology and regenerative biology.

Despite extensive pretransplant serological screening and HLA matching, 10-15% of kidney allografts experience acute rejection within the first year. Currently, risk stratification for transplantation relies primarily on antibody reactivity to HLA molecules, with no assessment of the T cell compartment before or after transplantation. In our previously established longitudinal cohort of 54 patients, T cell receptor {beta} (TCR{beta}) sequencing was performed on paired kidney biopsy and peripheral blood samples. Here, we further analyzed the data to construct a comprehensive set of sequence-similarity networks and quantify over 30 network metrics. After adjusting for repertoire size, graft status was the strongest signal for the underlying differences in network metrics. Individuals who rejected the kidney graft generally exhibited more fragmented and less connected networks at baseline, with fewer interconnect T cell clones and more isolated sequences. Notably, pre-transplant peripheral blood mononuclear cell (PBMC) network topology alone predicted non-stable outcomes with an area under the curve (AUC) of 0.81, sensitivity of 76%, and specificity of 76%. The performance of this prediction model was independent of HLA mismatch, while changes in network topology at three months post-transplantation further improved prediction to an AUC of 0.88 (permutation p = 0.009). Collectively, TCR sequencing and network analysis represent a potential novel, non-invasive approach for pre-transplant risk stratification and immune monitoring, capturing functional immunological risk that may not be accessible through HLA genotyping or serology.

Currently, the genetic architecture of Middle Eastern populations is underrepresented in global genomic databases. This gap increases the rate of Variants of Uncertain Significance (VUSs) and clinical misinterpretations of genomic data especially in Middle Eastern populations. Whole exome sequencing was conducted on 90 healthy individuals from Jordan and the data were analysed using Principal Component Analysis (PCA) and multi-computational filtering. PCA revealed a double ancestry (EUR-AFR) admixture rather than a triple admixture (EUR-AFR-AMR). More than 3,500 populations-specific variants (PSVs) were identified, of which 72% were singletons. Additionally, 19 variants were significantly enriched compared to the maximum allele frequencies in public global databases (Fisher's exact test with Benjamini-Hochberg false discovery rate correction, p-value < 0.05). Consequently, the results suggest the reclassification of variants of Uncertain Significance (VUS) which reside in the ECE2 gene to likely benign and the variants of Conflicting Classification of Pathogenicity in the genes IL1RN and THPO to benign based on the significant allele frequency (AF=0.0389, p-value < 0.05). Furthermore, a pathogenic ClinVar variant was identified in a healthy individual, warranting careful interpretation. The findings underscore the importance of identifying PSVs in order to minimize or even prevent clinical misdiagnosis and highlight the unique genetic signature in Jordan. The study serves as a foundational resource for precision medicine in the region.

The mucosal chemoattractant C-X-C motif chemokine ligand 17 (CXCL17) was recently identified as a ligand for the orphan G protein-coupled receptor 25 (GPR25). Although CXCL17 orthologs have been identified in fishes, amphibians, and mammals, their presence in reptiles and birds remains unclear. In this study, we employed bioinformatic searches based on gene synteny and sequence features to identify CXCL17 orthologs in public databases. We identified functional CXCL17 orthologs in 46 reptilian species, including lizards, snakes, turtles, and alligators. In contrast, we found only non-functional gene relics in 22 bird species, suggesting the avian lineage lost functional CXCL17 during evolution. A recombinant reptilian CXCL17 from the loggerhead turtle (Caretta caretta), termed Cc-CXCL17, directly bound to and efficiently activated its corresponding receptor, Cc-GPR25, in a C-terminal fragment-dependent manner. Activation of Cc-GPR25 by Cc-CXCL17 also induced chemotactic movement of transfected human embryonic kidney (HEK) 293T cells. In cross-species activity assays, CXCL17s from human and tropical clawed frog could activate Cc-GPR25 albeit with lower potency, but fish orthologs lacked this activity; all tested CXCL17s had no detectable activity towards chicken GPR25, but Cc-CXCL17 had low activity towards mallard GPR25. Our findings demonstrate the presence of functional CXCL17 orthologs in extant reptiles and provide evidence for their evolutionary loss in birds, offering new insights into the phylogenetic distribution of the newly identified CXCL17-GPR25 signaling system.

During meiosis, homologous chromosomes pair to form synaptonemal complexes (SCs) and exchange genetic material through a process known as meiotic recombination. First, programmed DNA double-strand breaks form, followed by the assembly of recombination foci on SCs. These foci mark the sites of recombination intermediates and future crossovers. Distributions of recombination foci along SCs have been studied in many eukaryotes, revealing the interplay between recombination patterns and genome evolution. However, in fish, data on recombination patterns are scarce, and, for the majority of groups, completely absent. Here, we measure the positions of MLH1 foci in 3,504 SCs from 219 male meiotic cells of an African annual killifish Nothobranchius virgatus, a representative of a genus with remarkable karyotype and genome diversity, and present a detailed statistical analysis of its recombination patterns. We found that, in contrast to the several other fish species characterised to date, recombination in N. virgatus occurs across almost entire chromosome arms, excluding (peri)centromeres and telomeres. In the longest SCs, we observed a proximal and a distal peak of the recombination focus frequency and explained the peaks by chromosome pairing dynamics. We also revealed the typical positions of focus pairs, demonstrated interference between foci, with the minimal interfocus distance of 4 m, and described regions of the total recombination suppression near centromeres and telomeres. In sum, our study provides a detailed analysis of recombination patterns in a killifish with a fully acrocentric karyotype and contributes to cytogenomic and statistical methodology for future exploration of meiotic recombination patterns.

Junctional adhesion molecule-like (JAML) is an adhesion molecule known to promote T cell activation and T cell-mediated tumor rejection. In the current study, we show that JAML expression is enriched on mouse intratumoral NK cells compared with splenic NK cells. JAML+ NK cells were associated with tissue residency and co-expressed the immune checkpoints PD-1 and LAG3. JAML expression could be induced on splenic NK cells by IL-2 and further enhanced by IL-21. JAML levels were inversely correlated with inhibitory signaling, as NK cells expressing self-recognizing Ly49 receptors had reduced JAML expression, suggesting regulation of JAML expression by MHC class I molecules. Interaction with the JAML ligand CXADR also reduced JAML surface expression, indicating that tumor-mediated membrane stripping may represent a mechanism of immunoediting. Although JAML RNA transcripts were detectable in human NK cells, JAML protein was found only intracellularly. Together, these findings identify the JAML-CXADR interaction as a potential regulatory pathway in NK cell-mediated killing of tumors.

The phylum Mollusca constitutes one of the most taxonomically and morphologically diverse animal clades; however, the genomic exploration of this group has been hampered by fragmented and taxonomically incomplete transcriptomic resources. To address this fundamental limitation, we present MolluscaGenes, a centralized database that unifies transcriptomes from 299 molluscan species spanning all eight recognized classes, encompassing a broad array of tissues and developmental stages. MolluscaGenes provides searchable databases via BLAST and DIAMOND alongside a suite of 196 molluscan-optimized Hidden Markov Models (HMMs) for sensitive protein family identification. To demonstrate the utility of this resource, we performed a comprehensive phylum-wide characterization of the nicotinic acetylcholine receptor (nAChR) superfamily, recovering 3,586 sequences from over 190 species and resolving 15 distinct phylogenetic clades. This analysis revealed substantial lineage-specific expansions across multiple molluscan classes, the identification of novel clades with substitutions in canonical ligand-binding residues, and the evolutionary placement of chemotactile receptors (CRs) and CR-like sequences as predominantly cephalopod clades within the broader nAChR phylogeny. MolluscaGenes constitutes a foundational resource that will accelerate the elucidation of the unique biology and evolutionary history of Mollusca.

Early complications following allogeneic hematopoietic stem cell transplantation (allo-HSCT), including graft-versus-host disease (GVHD) and viral reactivation, remain major causes of post-transplant morbidity, but whether immune perturbations precede these events remains unclear. We performed longitudinal TCR{beta} repertoire profiling in 108 allo-HSCT recipients and their corresponding donors at baseline and three early post-transplant time points to characterize immune reconstitution dynamics. Reduced baseline TCR diversity was most strongly associated with subsequent Epstein-Barr virus (EBV) reactivation, whereas cytomegalovirus (CMV) reactivation was more closely linked to post-transplant repertoire remodeling characterized by clonal expansion and reduced donor-recipient repertoire similarity. Sequence-based predictive modeling demonstrated meaningful discrimination, with fusion models achieving area under the curve (AUC) values of 0.745 for CMV, 0.819 for EBV, and 0.834 for GVHD. Temporal analyses further revealed complication-specific predictive windows. These findings indicate that major post-transplant complications are preceded by detectable immune perturbations and support the potential utility of TCR repertoire monitoring for early risk stratification after transplantation.

The Toll/interleukin-1 receptor (TIR) domain is important for immune signaling across bacteria, plants, and animals. In human innate immunity, TIR domains are known to function as adaptors mediating protein-protein interactions, yet studies in bacteria and plants revealed that TIR domains often act as enzymes that produce immune signaling molecules. Here, we show that TIR domains from evolutionarily diverse animals have conserved active sites, implying that they can function as enzymes. In vitro experiments with animal TIRs show that the TIR domain of several Toll-like receptors (TLRs), including that of human TLR4, can produce cyclic ADP-ribose (cADPR), revealing an enzymatic activity previously unknown for TLR TIRs. We show that production of cADPR is a conserved feature of TIR domains across the animal tree of life, implying a role for this molecule in animal TIR signaling. Finally, we report a TIR domain from green algae that synthesizes 3'cADPR, suggesting conservation of 3'cADPR signaling between bacteria and eukaryotes. Our results reveal that the catalytic activity of TIR domains is widespread in animals and conserved across the tree of life.

Neoantigens are cancer-specific antigens arising from genomic alterations. Single Amino Acid Variants (SAAVs) represent a primary class of these neoantigens. To evaluate the therapeutic potential of Neurofibromin 1 (NF1)-derived SAAVs - given that NF1 is frequently mutated in malignant brain tumors - we prioritized the 40 NF1 SAAVs determined to be HLA-A*02:01 binders using computational prediction coupled with experimental validation. To validate these predicted neoepitopes, we employed a two-tiered experimental approach in HLA-A*02:01 homozygous U87-MG cells. We first synthesized minigene constructs encoding the predicted neoepitopes, introduced them via lentiviral transfection and confirmed their expression by mass spectrometry (MS). Subsequently, we performed endogenous validation using pan-HLA immunoprecipitation mass spectrometry (IP-MS), confirming 4 (10 neoepitopes) of the 40 candidate SAAVs. We observed a discrepancy between in silico predictions and the observed sequences. Our endogenous peptidomics further revealed conserved peptide motifs and demonstrated that peptide selection for HLA presentation is transient. While our study substantiates the therapeutic feasibility of T-cell immunotherapies targeting NF1 mutations, these results underscore a limitation in current computational prediction. Our study highlights the necessity of experimental validation to refine neoantigen prioritization strategies.

The thymic medulla provides the microenvironment for negative selection, late thymocyte maturation, and thymocyte egress, and is generally characterized by widespread distribution of medullary thymic epithelial cells (mTECs). In contrast, rat thymic medulla contains medullary epithelium-free areas (mEFAs), but the cellular composition and functional significance of these regions remain unclear. Here, we combined spatial transcriptomics and scRNA-seq, using robust cell-type decomposition (RCTD) to characterize mEFAs in Lewis-strain rat thymus. These analyses revealed that more mature-phenotypes of CD4SP, CD8SP, and regulatory T-cell-lineage thymocytes were preferentially localized in mEFAs, whereas immature SP subsets were enriched in medullary epithelium-containing areas. Newly found rat thymic mesenchymal cell-3 and -4 (TMC3 and TMC4) subsets were also enriched in mEFAs. These subsets were broadly similar to mouse medullary fibroblasts but displayed distinct predicted interactions with SP thymocytes, including costimulatory molecule- receptor, chemokine-receptor, and ECM-integrin axes. In addition, the venous endothelial cells (vECs) expressing portal endothelial cell markers were accumulated in mEFAs. The S1P transporter gene Spns2 was preferentially expressed in both TMC4 and vEC subsets, suggesting increased local concentration in mEFAs. These findings indicate that rat mEFAs are specialized medullary niches linking stromal organization, thymocyte maturation, and thymic egress.

Behcets disease (BD) is a systemic inflammatory disease. It is considered as an autoinflammatory disease triggered by innate immunity rather than adaptive immunity. Human leukocyte antigen-B51 (HLA-B51) is the strongest genetic factor associated with BD. This study investigated how HLA class 1 molecules interact with innate immune cells and induce cytokine secretion. For this purpose, 293T cells transfected with a plasmid encoding HLA-B51 were cultured with natural killer (NK) cells obtained from healthy human donors. Within 24 h, the concentrations of interleukin-4 (IL-4), IL-8, and interferon-{gamma} (IFN-{gamma}) in the medium increased, indicating that NK cells secreted cytokines without undergoing cellular expansion for cytolysis. NK cells stimulated by nonself HLA-B51 produced IFN-{gamma} levels comparable to those produced by NK cells stimulated by self HLA-B51. NK cells carrying HLA-B51 were accurately recognized by overexpressing HLA-B51 on 293T cells. Moreover, ample intracellular IFN-{gamma} levels were detected in NK cells after stimulation with phorbol 12-myristate-13-acetate (PMA) plus ionomycin. KLRK1 (CD314)-positive cells mainly primarily accounted for IFN-{gamma}-producing cells, whereas KLRK1-negative cells did not. In contrast, both NCR1 (CD335)-positive and -negative cells contributed to IFN-{gamma} production. We next investigated whether HLA-B51 on the surface of 293T cells stimulates KLRK1 as a ligand causing IFN-{gamma} secretion. In masking experiments using anti-KLRK1 antibodies, NK cells with high levels of cell surface KLRK1 decreased the production of IFN-{gamma}. Conversely, human NK cell line KHYG1 cells also produced IFN-{gamma} in culture with 293T cells, but did not increase IFN-{gamma} through HLA-B51 stimulation. The mRNA expression of the signal adaptor protein HCST (DAP10) in KHYG1 cells was lower than that in NK cells, whereas the relative expression of IL-2RA in KHYG1 cells was higher than that in NK cells. These findings suggest that HLA-B51 can interact with KLRK1 on the NK cells inducing IFN-{gamma} secretion, whereas IL-2 signals outweigh HLA-51 stimulation in KHYG1 cells.