Immunogenetics

T cell repertoire sequencing has unveiled a vast, complex array of T cells responsible for the human immune responses. Traditional analytic methodology fails to fully characterize and quantify the diversity of T cell receptors constituting the T cell repertoire. T cell receptor clonal frequency measured in terms of T cell receptor beta (TRB) V gene segment usage when arrayed in correspondence with the respective V gene segment positions on the TRB loci yields a periodic, undulating curve in the spatial domain of the TRB genomic locus. Using the genomic distance from the TRB-D1 segment to the TRB-V1-29 segments, Fourier analysis was performed utilizing Lomb-Scargle periodogram to obtain Spectral Power curves quantifying the TRB V clonal frequencies from 6 allogeneic stem cell transplant donors (baseline) and recipients (>/=100 days) using a variety of analytic software. Spectral Power curves revealed dominant spectral peaks at wavelengths ranging from 4-9 kb (113-252 millicycles/kb) in the six donors, with consistent frequency domain spectral patterns. This is consistent with similar use of V segments across healthy individuals. Recipients on the other hand demonstrated more dispersed spectra, with a spectral centroid shifted towards higher frequencies compared to donors (260 vs. 247 millicycles/kb). Consistent with this observation, the Low Frequency Index was lower in the recipients (0.18 vs 0.20). Power was concentrated in the <3 kb and 3-12 kb wavelengths in both groups. The analyses reported here demonstrate that the healthy SCT donors have a remarkably similar spectral signature occupying short to intermediate wavelegnths in the frequency domain, whereas recipients tend to shift towards higher frequencies. These findings are consistent with a normal organized distribution of TRB V segment usage in healthy individuals (by analogy other loci), and a more diffuse and disorderly usage in recipients, consistent with the notion of T cell responses constituting a dynamical system which evolves as a function of time. Fourier analysis of TRB (and potentially TRA) sequencing data provides a repertoire wide summary of T cell clonal distribution.

Advances in immunophenotyping of tumour-infiltrating lymphocytes (TILs) have improved our understanding of prognostic biomarkers and immune targets in human melanoma. However, whether the tumour-immune landscape in canine oral malignant melanoma (COMM) is concordant with that of human melanoma has not been properly defined. To address this gap, we developed a 16-colour spectral flow cytometry panel to characterise TIL phenotypic and functional profiles in COMM. Validation using mitogen-stimulated peripheral blood mononuclear cells from healthy dogs (n = 5) demonstrated robust identification of major T cell lineages, including regulatory T cells (Tregs) and memory subsets, and reliable evaluation of their activation and exhaustion status. COMM patients (n = 8) displayed a distinct protumour microenvironment, characterised by an increased proportion of Tregs, enrichment of tumour-specific exhausted-like T cells co-expressing programmed cell death protein 1 (PD-1) and tumour necrosis factor receptor 2 (TNFR2), and a reduction in cytotoxic CD8 and natural killer T (NKT) cell populations compared with tissue-resident (n = 4) and circulating (n = 8) lymphocytes. Analysis of additional solid tumours, including a mast cell tumour, nerve sheath tumour and adrenal cortical carcinoma, further supported the capability of the panel to identify similar patterns of immune dysregulation across diverse canine tumour landscapes. Collectively, this work describes the first detailed evaluation of canine TIL immunophenotyping using spectral flow cytometry and provides insights into the immunosuppressive mechanisms shaping the tumour microenvironment in COMM. These findings not only increase our understanding of canine tumour immunology but also identify potential immune targets and support ongoing comparative immuno-oncology efforts.

Introduction: How much of the variable (V) and joining (J) gene identity of a T-cell receptor is recoverable from its third complementarity-determining region (CDR3) amino-acid sequence alone? Immune repertoire studies often report the CDR3 with V and J annotation that is missing, low-confidence, or inconsistent, so what the CDR3 alone can and cannot fix is both a basic question about the receptor and a practical one for reading those repertoires. Methods: For each of 118,096 pooled human rearrangements (37,687 and 80,409 {beta}) we computed the posterior distribution over candidate genes under a generative model of V(D)J recombination and under its post-selection counterpart, and measured recoverability by conditional entropy, the candidate-list size needed to contain the annotated gene, the fraction of sequences admitting a high-confidence single-gene call, and the structure of gene-by-gene confusion. Results: The J gene was nearly determined by the CDR3 in both chains. The V gene was only partially recoverable, and behaved as a group rather than a gene: junctional trimming and non-templated insertion, together with the loss of synonymous codon information in translation, leave sets of mutually confusable V genes whose grouping departs sharply from germline family nomenclature (adjusted Rand index 0.05 for and 0.21 for {beta}). Selection sharpened the V posterior modestly (usage-controlled entropy shift -0.06 nats for and -0.28 for {beta}) and redistributed which V gene was most probable, a locus-scale rewrite in {beta} against a mild reweight in . Both the recoverability measurements and the confusion grouping reproduced in two held-out tumor cohorts. Discussion: V identity is an emergent, system-level property of the repertoire, set jointly by recombination and selection and invisible in any single rearrangement, so it should be reported as a calibrated group rather than a single gene. We also release the pipeline with a computational tool which can output a set of candidate genes with confidence values given a CDR3 sequence.

Production of aberrant glycans by cancer cells constitutes a key immunosuppressive strategy to avoid destruction by immune cells. Although sialic acid-containing glycans are known to dampen the activation of lymphocytes, including Natural Killer (NK) cells, the role of fucose-containing glycans remains poorly characterized. In this work, we explored the role of Lewis X (LeX) in cancer cell-NK cell interactions. We induced ectopic expression of FUT9, an 1-3/4-fucosyltransferase, in two colorectal cancer cell lines and showed this enzyme only synthesized LeX structures but not sialyl-LeX. FUT9 introduction was not associated with altered MHC class I surface expression, nor with CD2 (which has been proposed as a receptor for LeX) binding to cancer cells. By inhibiting fucosylation we could demonstrate that CD2 binding was furthermore independent of surface fucosylated glycans in three independent cell lines. Lastly, FUT9/LeX had a limited role in cancer cell destruction and expression of activation markers by NK cells. Overall, our study suggests that, unlike sialylated glycans, 1-3/4-fucosylated glycans have limited impact on cancer cell evasion of NK cell-mediated destruction.

Background and objectivesMyasthenia gravis associated with antibodies against muscle-specific kinase (MuSK-MG) is a well-characterized IgG4-autoimmune disease, however, the mechanisms driving IgG4 predominance remain poorly understood. This study investigated whether promoter DNA methylation of cytokine genes involved in IgG4 class switching is associated with this immune response. MethodsPeripheral blood mononuclear cells were isolated from MuSK-MG patients (n=36), acetylcholine receptor myasthenia gravis (AChR-MG) patients as disease controls (n=7), and sex-matched healthy controls (n=12). Promoter DNA methylation of IL4, IL10, and IL13 was assessed by methylation-sensitive high-resolution melting and relative cytokine mRNA expression by qPCR. Associations with clinical variables, and antibody levels were subsequently evaluated. ResultsMuSK-MG patients showed lower median IL13 promoter methylation compared with healthy controls (p = 0.004). Median IL4 promoter methylation was also reduced in MuSK-MG compared with healthy controls (p < 0.001) and AChR-MG disease controls (p < 0.001), whereas no differences were observed for IL10 promoter methylation. Relative mRNA expression of IL4 (p = 0.0005), IL10 (p = 0.0462), and IL13 (p = 0.0002) was increased in MuSK-MG compared with AChR-MG. Compared with healthy controls, only IL4 expression remained significantly increased (p < 0.0001). Promoter methylation was inversely correlated with relative mRNA expression for IL4 (p < 0.0001), while IL13 showed a similar but non-significant trend (p = 0.054), no association was observed for IL10. Multivariable analysis demonstrated that treatment at sampling was independently associated with lower IL10 and IL13 promoter methylation, whereas no associations were observed with age, sex, disease phase, or disease duration. Promoter methylation did not correlate with total serum IgG4 or anti-MuSK IgG4 levels. DiscussionMuSK-MG is associated with selective hypomethylation of IL4 and IL13 promoters accompanied by increased cytokine gene expression, while IL10 promoter methylation remains unchanged. The association between treatment and IL10 and IL13 promoter methylation suggests that immunosuppressive therapy may influence epigenetic regulation in MuSK-MG. Together, these findings support a role for epigenetic dysregulation of Th2-associated cytokines in the immunological environment associated with IgG4 subclass switch. To our knowledge, this is the first study investigating IL4, IL10, and IL13 promoter DNA methylation in MuSK-MG.

EU has reclassified the sika deer (Cervus nippon) as an undesirable invasive species based on reports that hybridization with the indigenous red deer (C. elaphus) may produce fertile offspring. Since sika-derived DNA previosuly introduced into the red deer population (introgression) cannot be removed, the crucial question is whether new (F1) hybridisation occur. To address this, we analysed the chromosomes in 56 sika and 22 red deer. All red deer had a chromosome number 2n=68. In contrast, the chromosome number in sika ranged from 64 to 67, due to the variable presence of two sika-specific Robertsonian translocations (ROB1,ROB2). In the free-ranging sika population in Jutland, >90% of the sika deer were homozygote for at least one of these ROBs, excluding that they could be F1-hybrids. Moreover, ROB2 was in Hardy-Weinberg equilibrium, further supporting the absence of gene flow between the two species. In contrast, ROB1 was in Hardy-Weinberg disequilibrium, suggesting negative fitness of heterozygotes, including potential F1-hybrids. In Jaegersborg Deer Park, the eight examined sika deer had the same genotype (absence of ROB1, homozygosity of ROB2), supporting that it is a founder population which may have been isolated for [~]100 years. Again, none of these can be F1-hybrids due to the homozygosity of ROB2. We conclude that F1-hybridisation between sika and red deer either does not occur or occur very rarely in Denmark. The study establish the Danish sika-populations as unique models for adressing important biological questions: What underlies the absence of hybridisation? Why are ROBs frequent in sika deer but not in the closely related red deer? How fast do new species/subspecies develop in isolated founder populations? Which factors determine, that some ROBs have little heterozygous effects, whereas others are selected against, with implications for the role of ROBs as genetic barriers promoting speciation, and for fertility problems in some human ROB carriers.

Glycosylation of proteins is arguably the most diverse post-translational modification that is altered in almost all cancer types, which has been demonstrated to play a crucial role in creating an immunosuppressive microenvironment that promotes immune tolerance and evasion. However, the biosynthesis of N-linked glycan is mediated by a series of enzymatic reactions catalysed by glycosyltransferases and glycosidases in a template-independent manner, hindering our understanding of specific structure-function relationships and roles of specific glycans on specific proteins. Here we use HLA class I-negative cell line K562, a known reference target for NK-mediated cytolysis, to establish a model investigating how cell surface glycan dynamics influence its susceptibility to cytolysis mediated by NK-92 cells. Treatment of K562 cells with kifunensine, swainsonine, 2F-peracetyl-fucose, or 3Fax-peracetyl Neu5Ac, inhibitors of N-linked glycan processing, resulted in drastic alterations in cell surface carbohydrate phenotype, as could be shown by flow cytometric analysis of the lectiNbinding properties of the cells. Despite these clear changes in carbohydrate phenotype, only K562 cells treated with either kifunensine or 3Fax-peracetyl Neu5Ac exhibited higher susceptibility to the cytolysis medidated by NK-92 cells accompanied with an increased CD107a expression by NK-92 cells. Although K562 cells overexpressing gene MGAT3 exhibited a decreased NK-susceptibility, we further found that this decrease was not exclusively determined by the overexpression of gene MGAT3 product bisecting {beta}1,4-GlcNAc, because the treatment of 3Fax-peracetyl Neu5Ac reversed the resistance of K562 cell against NK-92 cell in despite of expressing higher levels of bisecting {beta}1,4-GlcNAc. Expressing HLA-G on cell surface as extravillous trophoblast did not change the NK-susceptibility of K562 cells, despite evidence that HLA-G molecules expressed by K562 cells can bind to inhibitory receptor ILT2 expressed on NK-92 cell surface. These findings suggest that the level of terminal sialylation, outweighing other components in N-linked glycan, determines the NK-susceptibility of K562 cell, offering a new strategy to weaken the resistance of cancer cells so that the immune system can maximise the elimination.

NKG2D (encoded by KLRK1 in humans and Klrk1 in mice) is an activating receptor expressed by cytotoxic lymphocytes. In humans, NKG2D signaling is regulated post-transcriptionally: activated T cells retain intron 4 of KLRK1 to generate NKG2D-TR, a truncated dominant-negative isoform that limits receptor signaling. Whether mice, the principal preclinical model for NKG2D-directed therapies, possess an analogous regulatory mechanism remains unknown. Across four murine RNA-seq datasets comprising 50 samples and spanning T-cell differentiation, graft-versus-host disease, and acute and chronic LCMV infection, we examined the retained-intron isoform Klrk1-203. The transcript retains the canonical start codon, while its predicted stop codon lies within the terminal exon downstream of the final exon-exon junction, suggesting that it may escape nonsense-mediated decay. If translated, Klrk1-203 is predicted to encode a truncated product that retains the cytoplasmic and transmembrane domains but lacks most of the ligand-binding ectodomain. This predicted architecture resembles human NKG2D-TR, although the murine product contains a short C-terminal sequence encoded by the retained intron. Klrk1-203 was below the detection limit in unchallenged naive and early-effector T cells but was induced in differentiated effector and memory populations, reaching approximately one-fifth of total Klrk1 transcripts in one effector-memory sample. Read-level analyses independently demonstrated increased intron 4 retention with differentiation; however, short-read sequencing cannot fully distinguish Klrk1-203 from the co-retained Klrk1-204 transcript, making isoform-specific abundance dependent on model-based quantification. An independent coding-potential algorithm classified Klrk1-203 as non-coding, providing an important counterpoint to the structural predictions. Together, these findings identify Klrk1-203 as a candidate NMD-resistant, differentiation-associated regulator of murine NKG2D and a potential counterpart of human NKG2D-TR that warrants experimental validation.

Bone marrow chimeras are widely used to study immune development and function. As the field moves from cesium-137 (137Cs)-based irradiators to X-ray irradiators for safety reasons, it is essential to determine if there are differences in immune system reconstitution after irradiating mice with one of these two radiation sources. Here, we performed a comprehensive immunological comparison of mice lethally irradiated with 137Cs or one of two different X-ray platforms and reconstituted with congenic bone marrow. Mice received 12 Gy total body radiation in two 6 Gy sessions followed by intravenous transfer of donor hematopoietic stem cells and were analyzed eight weeks post-transplant. We assessed mouse survival, donor chimerism, immune cell subset distribution, and activation states across multiple organs (bone marrow, spleen, lymph nodes, liver, and lung). All groups exhibited comparable survival and high levels of donor chimerism, with expected organ-specific reconstitution patterns. Immune lineage distributions, CD4/CD8 ratios, and activation states did not differ by irradiation type. Host-derived radioresistant cells were also similar across all irradiation groups and were predominantly composed of T cells skewed toward an activated phenotype. Overall, our data show that X-ray irradiation with proper filters and energy levels (225 KVp and 320 KVp) can yield equivalent immunological outcomes, including immune reconstitution and activation states, as compared to the same radiation dose from 137Cs-based irradiation in bone marrow chimera models. These results support the continued adoption of X-ray irradiation systems in place of 137Cs for generating bone marrow chimeras to be used across a wide range of immunologic studies.

Japanese macaque (Macaca fuscata) is used in biomedical and neurobiology research, yet transcriptomic resources for the brain are limited. We present a hybrid RNA sequencing dataset and a prefrontal cortex transcriptome assembly from two healthy 6-year-old animals. Short-read Illumina ({approx}70 million paired-end reads per sample) and long-read Oxford Nanopore direct RNA sequencing ({approx}2.5 million reads per sample) were combined. Reads were quality controlled, aligned to the macFus_1.0 reference genome, and assembled with StringTie2. Transcripts were annotated using Trinotate and eggNOG-mapper, and open reading frames were predicted with TransDecoder. The released data package includes raw reads (NCBI SRA BioProject PRJNA1295993), transcript sequences and structural annotation files, predicted coding sequences and proteins, functional annotation tables, and transcript abundance estimates (TPM). Technical validation includes read-level QC and protein-level comparisons to expressed gene sets from human, rhesus macaque and chimpanzee prefrontal cortex. These resources enable reuse for transcript-level expression studies, isoform characterization and comparative primate neurogenomics.

CD8 T cells play a central role in immune responses to infection and cancer. However, the diversity of T cell receptor (TCR) specificities makes it challenging to study the mechanisms that regulate T cell activation, differentiation, and effector function. Beyond TCR transgenic mouse models, various complex genome-editing approaches have been employed to overcome this challenge. However, these strategies are often technically demanding, time-intensive, and difficult to adapt. Investigators who are interested in testing de novo TCRs under their chosen experimental conditions would benefit from a standardized and accessible method. Here, we describe a protocol that combines ribonucleoprotein (RNP)-based CRISPR-Cas9 editing with retroviral transduction to enable efficient genetic manipulation of murine CD8 T cells. We show that T cells engineered via this protocol can be generated at sufficient scale for downstream in vitro assays and in vivo adoptive transfer experiments. We expect this method will be useful for investigators who require a standardized and accessible way to study how TCR specificity impacts CD8 T cell responses.

Rosa gallica L., the French rose, is a perennial, tetraploid, heterozygous species that naturally propagates by seed and sucker. It occurs in the wild, primarily in Europe, and also exists as cultivated varieties. R. gallica cultivars were extensively bred and cultivated in France at the beginning of the 19th century. Although several hypotheses have been proposed regarding the species expansion based on historical records, none have been assessed using molecular data. Indeed, its genetic diversity has so far been investigated only at local or regional scales, hindering the identification of the evolutionary factors shaping its present-day distribution. Using 29 sequenced microsatellites, we genotyped a comprehensive sample of 1618 individuals, including wild R. gallica from 219 sites across the species range, rose cultivars, and specimens from other Rosa species. We then detected clonal lineages and characterized the range-wide genetic diversity and structure, aiming to disentangle the roles of natural and human factors in shaping the distribution of R. gallica, with particular focus on France. French diversity appears particularly structured compared to the rest of the range, suggesting multiple origins within France. Populations in South Alps, Central Eastern Europe, and Eastern France appear to have recolonized naturally from a single southern glacial refugium. In contrast, populations in the western part of France likely resulted from more recent natural or human-mediated dispersal. Finally, clonal lineages containing both wild and cultivated individuals were predominantly found in France, highlighting the role of human-mediated dispersal in 28 of the 98 French sites studied. These findings show that the present-day natural range of R. gallica was shaped primarily by post-glacial recolonization, but also reveal a contribution of human activities to its recent dispersal, particularly in France, where cultivated varieties were intensively bred and exchanged.

Antibodies exhibit species-specific sequence and structural features that influence their antigen-recognition properties. Although several studies have investigated porcine antibodies, their repertoire and structural characteristics remain less well characterized than those of several other mammalian species. In this study, we analyzed public porcine heavy-chain repertoire sequencing data together with available antibody structural data to identify characteristic features of porcine antibodies. We found several residues enriched in porcine antibody framework regions, particularly at the base of heavy-chain complementarity-determining region 3 (CDR-H3). In particular, Arg101 and Glu123 were closely positioned in available structures and may influence CDR-H3 conformation at its base, whereas Pro120 may help constrain local backbone conformation. We also observed non-canonical cysteine usage in both framework region 1 and CDR-H3, which may contribute to structural diversity in the porcine repertoire. Finally, we evaluated the humanization potential of a porcine antibody using a human antibody language model and found that human-likeness increased after model-guided substitutions, although the resulting sequences did not exceed the T20 score threshold. Overall, these results indicate that porcine antibodies possess distinct sequence and structural features that may influence CDR-H3 properties and should be considered in future antibody analysis and engineering.

Background: Viral myocarditis presents a significant burden of disease, particularly among children and young adults. However, clinical guidelines and treatment strategies for pediatric patients are derived from those for adult patients due to a lack of pediatric data. Current animal models of viral myocarditis use adult mice, so conclusions from these models cannot necessarily be extrapolated to the pediatric population. We sought to develop a juvenile mouse model of myocarditis to examine differences between these two distinct clinical populations. Methods: Male and female BALB/c 3-4-week-old 'juvenile' and 8-week-old 'adult' mice were infected intraperitoneally with 103 PFU of heart-passaged coxsackievirus B3. Sera was used to evaluate testosterone and estradiol levels. Cardiac histological evaluations included overall inflammation, fibrosis, and specific cell-type infiltration. RNA was extracted from cardiac tissue and evaluated for changes in gene expression of cell-type markers, complement components, and NLRP3 inflammasome components. Results: Juvenile mice exhibited more severe inflammation than adult mice but no sex differences in overall inflammation. Juvenile mice demonstrated increased infiltration of CD11b+ cells, F4/80+ cells, and CD3+ T-cells vs. adults. Inflammasome genes NLRP3 and caspase-1 were significantly increased in juvenile compared with adult myocarditis. Conclusions: This paper is the first to describe a juvenile mouse model of coxsackievirus B3 myocarditis and provides a direct comparison to a translational adult mouse model. Juvenile mice had greater cardiac inflammation than adults. This model replicates clinical populations and provides a valuable tool to study age as a factor in the pathogenesis of myocarditis.

Extravillous trophoblasts (EVTs) are a trophoblast subpopulation critical for feto-maternal tolerance during early pregnancy, primarily using HLA-G to exert immunomodulatory effect, and possessing N-glycomic profiles distinct from other trophoblast subpopulations. However, whether the N-glycosylation confers distinct immunological properties to EVTs remains poorly understood. To investigate this, we employed JEG-3, a human choriocarcinoma cell line having the capacity to produce pregnancy-related hormones and expressing both HLA-C and HLA-G resembling placental EVTs, as an in vitro EVT model, alongside cell line JAR which exhibits villous trophoblast phenotypes distinct from JEG-3. Both cell lines were treated with kifunensine or swainsonine, inhibitors of -mannosidases, to remodel their N-glycosylation patterns. This led to significant remodelling of their N-glycomic profiles, with JEG-3 cells showing an increased level of polylactosamine chains and decreased levels of -2,6-sialylation and core -1,6-fucosylation. Western blot analysis showed that inhibiting -mannosidases altered only the composition of N-glycans on cell-surface HLA-G, without affecting the overall abundance of cell-surface HLA-G. In kifunensine-treated JEG-3 cells that predominantly express oligomannose type N-glycans, an intracellular accumulation of unfolded HLA-G fragments, increased hCG secretion, and down-regulations of EVT markers GATA3 and KRT7 were observed compared to untreated control, while swainsonine treatment did not impact N-glycan expression. Cytotoxicity assays using NK-92 as effector cells showed that the de-sialylation of JEG-3 by neuraminidase treatment led to increased NK-92 mediated killing. JEG-3 cell sustained its EVT immunological properties through generating bioequivalent N-glycans, exemplified by NK-92 cells pre-conditioned with used culture media of kifunensine-treated JEG-3 cells displaying reduced cytotoxicity toward NK-sensitive lymphoblast cell line K562, an effect not observed with swainsonine-treated JEG-3 cells. This model suggests that EVTs immunological properties are dependent on specific N-glycomic profiles that are maintained by unique N-glycosylation homeostasis, and overall improves our understanding of how EVTs maintain their immunomodulatory effect at the maternal-fetal interface.

Understanding the genetic diversity and population structure of breeding materials is essential for developing stress-resilient cultivars. In tropical maize, where drought and low soil nitrogen (low N) severely limit productivity, continuous development of tolerant varieties remains a priority. This study assessed the genetic diversity and population structure of 250 doubled haploid lines (DHLs) derived from five drought- and low N-tolerant tropical populations. Genotyping was performed using mid-density DArTseq markers, yielding 3,305 high-quality SNPs for analysis. Results revealed a moderate level of diversity among the DHLs, with an average genetic distance of 0.39, a polymorphism information content (PIC) of 0.33, and a minor allele frequency (MAF) of 0.29. These values reflect substantial allelic variation, important for identifying complementary parental combinations in hybrid development. Discriminant analysis of principal components (DAPC) grouped the DHLs into five distinct clusters, largely corresponding to their source populations, although some admixture was observed. This indicates that while the genetic backgrounds of the source populations were mostly retained, recombination introduced useful variation. Overall, the clear population structure and high diversity observed among these DHLs provide a strong genetic foundation for future maize improvement. These lines represent valuable resources for heterotic group formation, hybrid development, and recurrent selection schemes aimed at enhancing drought and low nitrogen tolerance in tropical maize.

Background: The gradual decline in the functionality of the immune system during aging is commonly referred to as immunosenescence. This study aims to investigate changes in the B-cell receptor immunoglobulin heavy chain (BCR IGH) gene repertoire during natural aging and evaluate the dynamics of emergent low-level BCR IGH clonality in the elderly. We conduct a longitudinal study nested within the Doetinchem Cohort study, comprising 98 participants aged between 31 and 59 years old at study entry who had repeated blood samples drawn at 5 year intervals over a 30 year period (n=548 samples). We sequenced the IGH gene repertoire and evaluated the impact of aging on IGH gene repertoire clonality and diversity using linear mixed effects modeling. Results: Participants older than 60 years exhibited increased BCR IGH clonality and reduced IGH gene repertoire diversity. In a multivariable model, IGH gene repertoire diversity was significantly decreased for individuals with a dominant clonotype ratio greater than 10 (Beta=-0.57, p < 0.001). Additionally, a trend toward reduced IGH gene repertoire diversity was observed in participants aged 60-70 years (Beta =-0.19, p = 0.1) and those aged 70 years or older (Beta =-0.20, p = 0.13). IGH gene repertoire diversity was determined primarily by the naive and transitional B-cell pool, while BCR IGH clonality was influenced by switched memory and age-associated B-cell counts. Conclusions: In summary, our study indicates that IGH gene repertoire diversity decreases significantly after age 60, which coincides with an increased incidence of low-level BCR IGH clonality. This clonality may be driven largely by switched memory and age-associated B-cells. By providing deep insights into age-related dynamic changes in the IGH gene repertoire, these findings lay the groundwork for the molecular assessment and monitoring of incident clonality by clinicians and researchers alike.

We describe here an efficient method for gene editing in mouse T cells, based on well-established, high-efficiency retroviral transduction protocols. Our platform allows analysis of mutant phenotypes in primary murine T cells in vitro and in vivo. This approach uses a single retroviral vector to simultaneously knockout an endogenous gene and ectopically express its mutant version. This knockout/re-expression vector can be used as the only plasmid to transduce Cas9-expressing T cells, or used together with a Cas9 retroviral vector to transduce T cells from any mouse strain. We validated the system for analysis of murine T cells by targeting key molecules in proximal T cell signaling, i.e. CD3{gamma} and Zap70. We obtain high knockout and re-expression efficiencies in both Cas9-expressing and non-Cas9 T cells. Knockout efficiencies can be further improved by gRNA multiplexing. Endogenous proteins compete with their ectopically expressed mutants or tagged versions for cellular location, protein interactions and cellular functions. Here, we quantified the incorporation of CD3{gamma}-GFP into surface T cell receptor (TCR) complexes. Our data shows that the knockout and re-expression platform improves integration of CD3{gamma}-GFP into the TCR. Therefore, eliminating competition between endogenous and ectopic proteins benefits analyses of protein assemblies and signaling pathways in primary T cells. Furthermore, we validated advantages of our system for mutant analysis using wild-type and mutant Zap70s. Zap70 mutants deficient in TCR binding or kinase activity show their phenotypes only in the absence of endogenous protein, further validating our knockout/re-expression approach. Most importantly, this system can be used to generate gene-edited primary T cells for in vivo studies, such as the quantification of anti-tumor responses. Our knockout and re-expression platform provides a useful gene editing tool for primary T cells in fundamental research and immunotherapy development.

Understanding the evolutionary history of species within a geographic context is key to historical biogeography, as it reveals how geological and climatic changes shaped biodiversity. This is especially important in ecologically significant regions like the Amazon and Orinoco basins. Together, they host the worlds greatest freshwater fish diversity ([~]3,500 species), sharing a common but not yet fully understood evolutionary history. The gilded catfish (Brachyplatystoma rousseauxii), an ancient species widely distributed as a metapopulation in Neotropics, is an important model for studying past connectivity, divergence, and historical processes shaping fish diversity between these basins. This study analyzed the genetic structure, connectivity routes, and demographic history of B. rousseauxii using nuclear (microsatellite and ddRADseq) and mitochondrial DNA. Population structure analyses and coalescent models indicate that B. rousseauxii populations from the Orinoco and Amazon basins are genetically distinct, with no evidence of current gene flow. However, our results support the occurrence of a possible secondary contact event after the divergence, with the Boa Vista population retaining the genetic signal of this process. The ancestral population split occurred at the Rupununi Portal around 2.54 Ma (ddRAD) or 1.31 Ma (mtDNA). Then, the species colonized the Branco and Orinoco Rivers [~]1.90 Ma (ddRAD) or 0.6 Ma (mtDNA), rapidly expanding in the Orinoco (>1.3 or >0.29 Ma), while colonization of the Amazon from the Branco River was more recent ([&le;]1.0 or [&le;]0.15 Ma). Population expansion signal was detected in the Orinoco ([~]0.20 Ma), whereas the Amazon remained stable. Our findings suggest that the rise of the Vaupes Arch in the Late Miocene does not explain the observed genetic divergence. Likewise, the Casiquiare Canal and Japura-Guaviare headwaters are not connectivity routes between basins. Instead, the Rupununi Portal, including the recent capture of the Branco River by the Negro River, was the last point of connection and played a key role in shaping B. rousseauxiis distribution. These findings provide insights into Neotropical fish biogeography and the historical configuration of the Orinoco and Amazon basins.

Interactions between Human Leukocyte Antigen (HLA) molecules and their cognate immunoreceptors are essential for regulating innate and adaptive immune cell functions. Leukocyte Immunoglobulin-like Receptors (LILRs) are key regulators of HLA-mediated immune responses, owing to their broad expression across immune cell populations and their ability to modulate both immune activation and tolerance. Among these, LILRB1-HLA interactions are increasingly recognised as important in transplantation, chronic infection and cancer therapies. Unlike other HLA-binding receptors, which recognise epitopes specific to HLA subsets, LILRB1 primarily engages the relatively conserved 3 and {beta}2-microglobulin components of HLA molecules, supporting its role as a broad regulator of pan-HLA class I-mediated functions. Nonetheless, there have been conflicting findings regarding the breadth of LILRB1-HLA-I interactions. While direct affinity studies on a limited subset of HLA-I molecules have revealed no significant differences in LILRB1 binding, broader analyses using single-antigen bead arrays suggest underlying variability. Here, we show through a broad binding assay that, while LILRB1 is a broad HLA-I-binding receptor, it exhibits differential preferences across HLA-I allotypes. Molecular dynamics analyses of the HLA-I-LILRB1 interface suggest that HLA-3 domain dynamism underlies these binding differences. We further determined the crystal structure of LILRB1 and used it to highlight intrinsic structural flexibility within its domains. Finally, these structural insights were leveraged to refine our understanding of the binding modalities of therapeutic monoclonal antibodies currently described. Together, our findings establish structural and mechanistic bases for differential HLA-I recognition by LILRB1 and provide insights into immunotherapeutic targeting of LILRB1.