Immunogenetics

Acute myeloid leukemia (AML) is an aggressive hematologic malignancy associated with a poor prognosis and an unmet need for more effective treatment options, particularly after relapse. One potential approach to address this need, which has been used successfully in solid tumors, is adoptive transfer of T cells engineered to express specific T cell receptors (TCRs). Particularly attractive are TCRs that can target epitopes derived from oncogenes on which cancer cells depend while sparing normal cells that do not harbor the mutation. One promising oncogenic target in AML is mutated nucleophosmin 1 (mtNPM1), present in up to 30% of AML cases. Most NPM1 mutations result in a common frameshifted sequence that generates peptides presented by several class I human leukocyte antigens (HLAs), such as HLA-A*11:01. In this study, we vaccinated HLA-A*11:01-transgenic mice with the mtNPM1 peptide AVEEVSLRK and isolated six functional TCRs from mtNPM1-reactive lymphocytes. When engineered into human T cells, two of these TCRs consistently mediated robust and specific lysis of an AML cell line and primary AML cells co-expressing HLA-A*11:01 and mtNPM1. Neither TCR recognized related human peptides, HLA-A*11:01-expressing normal peripheral blood mononuclear cells or stem/progenitor cells, or cells expressing diverse allogeneic HLAs, suggesting a low risk of cross-reactivity and alloreactivity. In an AML xenograft model, both TCRs controlled AML in hematopoietic organs, although extramedullary tumors emerged during treatment, potentially reflecting the aggressive behavior of the model that may not be representative of AML in patients. Together, these data suggest that the two TCRs represent strong candidates for further preclinical development in adoptive T cell therapy for NPM1-mutated, HLA-A*11:01-positive AML. These results also support the use of HLA-transgenic mice as a platform for generating TCRs targeting oncogene-derived epitopes in hematologic malignancies.

Amoebic gill disease (AGD), caused by the parasite Neoparamoeba perurans, is a major welfare and economic burden for the Atlantic salmon aquaculture industry. However, the spatial organisation of immune functions within the gill, the site of infection, remains poorly resolved. Here, we present the first spatial transcriptomic atlas of the Atlantic salmon gill affected by AGD, generated using 10x Genomics Visium on four AGD-positive fish, confirmed by gill scoring and quantitative PCR for parasite load. As no uninfected gills were included, this study maps the anatomical and immunological organisation of AGD-affected gills rather than defining AGD-specific pathology. Unsupervised clustering of spatial gene expression resolved distinct transcriptional domains that grouped into four principal anatomical compartments of the gill, including (i) epithelial and lamellar regions defined by keratin and osmoregulatory gene expression, (ii) gill arch regions enriched for connective tissue and smooth muscle markers, (iii) an erythroid-associated vascular cluster, and (iv) a skeletal muscle cluster with signatures of glycolysis and oxidative phosphorylation. Across these domains, we identified a spatially organised and diffused immune landscape, with pattern recognition receptors, antimicrobial effectors and chemokines concentrated at the filament surface, and immunoglobulin light chain transcripts consistent with B lymphocyte activity enriched toward the gill arch. Integration of the spatial data with a published single-nucleus reference atlas using reference-based deconvolution resolved the cellular identity underlying these domains and revealed frequent co-localisation of immune and structural cell types, suggesting regulated cell to cell communication within the gill-associated lymphoid tissue. Overall, our results show that the Atlantic salmon gill is not a transcriptionally uniform surface but a regionally specialised organ, with immune functions distributed across, and specialised within, distinct anatomical compartments

The human T cell repertoire is generated through recombination of a multitude of T cell receptor gene segments, yielding a complex array of T cell clones. T cell receptor beta (TRB) V gene segment defined, T cell clonal frequency when organized in correspondence with the respective V gene segment positions on the TRB genomic locus yields a periodic, undulating curve in the spatial domain. Using the genomic distance from the TRB-D1 segment to the TRB-V1-29 segments, spectral analysis was performed utilizing Lomb-Scargle periodogram (modified Fourier analysis) to obtain spectral power curves quantifying the TRB V clonal frequencies from 12 allogeneic stem cell transplant donors (baseline) and recipients (<100 days or b/w 100-365) using a variety of analytic software. Patients either underwent HLA matched related bone marrow or HLA matched unrelated blood stem cell transplantation, using PT-CY for the former and ATG for the latter. Spectral Power curves revealed dominant spectral peaks at wavelengths ranging from 113-335 millicycles/kb in the 12 donors, with consistent frequency domain spectral patterns, supporting similar use of V segments across healthy individuals. PT-CY recipients had power spectra closely corresponding to the respective donors, however the ATG recipients had relatively dispersed spectra, with spectral centroid shifted towards higher frequencies compared to donors and a reciprocal decline in the low frequency indices. Spectral power was concentrated in the <3 kb and 3-12 kb wavelengths in both groups, consistent with the periodicity observed in the relative V gene segment usage across the repertoire. The analyses reported here demonstrate that the healthy SCT donors have a spectral signature occupying short to intermediate wavelegnths in the frequency domain, PT CY recipients recover with donor-like spectra on the average, however ATG recipients tend to shift towards higher frequencies. These findings are consistent with a normal organized distribution of TRB V segment usage in healthy individuals, with PT-CY preserving this organization more consistently than ATG. Spectral (Fourier) analysis of TRB sequencing data provides a repertoire wide summary of T cell clonal distribution.

Core fucosylation (alpha1,6-fucosylation) of N-glycans, catalyzed solely by fucosyltransferase 8 (FUT8), plays important roles in development, immunity, and antibody effector functions. Here we report the generation of Fut8 floxed mice using a CRISPR/Cas9-based Extreme Genome Editing (EGE) system, in which loxP sites flank exon 8 of the Fut8 gene and a LoxP-STOP-LoxP-RFP cassette was knocked in downstream of the promoter, rendering the targeted allele functionally null before Cre-mediated recombination. Crossing floxed mice with CMV-Cre transgenic mice yielded complete knockout mice carrying the Cre-recombined allele. A PCR-based genotyping method was established that unambiguously distinguishes the wild-type, heterozygous, and homozygous genotypes (+/+, delta/+, and delta/delta). Intercrossing of heterozygous mice revealed that homozygous knockout mice die during the early postnatal period, consistent with previous reports, whereas heterozygous mice are viable. Flow cytometric analysis of AAL binding to splenic B cells and AAL lectin blotting both demonstrated markedly reduced core fucosylation in heterozygous mice (60.0 +/- 6.8% and approximately 48% of wild-type levels, respectively). These results show that disruption of a single Fut8 allele substantially reduces core fucosylation, and the mouse model described here provides a practical tool for functional studies of core fucosylation.

During biogeographical ancestry (BGA) analysis, an unknown individual for whom ancestry informative markers (AIMs) were analyzed is often assigned to one of a few reference populations. For some individuals, such an assignment is not meaningful: the two parents of the individual may come from different populations, or the individual may come from a population not included in the reference database. In this study, we examined classification methods in combination with 104 AIMs from the VISAGE Enhanced Tool to provide one of the following three possible outcomes: (i) in-reference (IR), with an assignment to one of the reference populations; (ii) admixed (AD), with parents from two in-reference populations; and (iii) out-of-reference (OR), where the individual's ancestry is not sufficiently represented by the reference populations. We therefore compared seven classifiers in combination with class-specific acceptance thresholds for the OR decision. The best-performing method according to our selection criterion was a naive Bayes classifier over pairs of populations, which we call the recent naive Bayes (RNB) classifier. Populations held out of the reference panel were correctly rejected in $84.9\%$ of cases. We further tested the method on $800$ individuals from seven additional groups not included in the reference panel. Only $29.4\%$ of individuals from four non-admixed populations were rejected, resembling joint ancestry with reference populations. In addition, we investigated the prediction accuracy in a newly collected dataset of 28 individuals with self-reported ancestry by analyzing AIMs using massively parallel sequencing.

The allorecognition system of the hydroid Hydractinia symbiolongicarpus plays an essential role in mediating self vs. non-self recognition among cnidarian colonies, analogous to the role of the major histocompatibility complex (MHC) in bilaterian immune systems. Both complexes have contiguous stretches of similar, highly duplicated genes, suggesting a homologous relationship. Protein structural prediction methods such as AlphaFold2 also indicate structural similarity between the proteins encoded within the Hydractinia ARC and the MHC. We have determined the structure of the ARC in the related hydrozoan Podocoryna americana to determine the degree of conservation of the cnidarian allorecognition complex. Protein structural prediction methods have revealed a number of candidate allorecognition proteins that form a long contiguous complex. Phylogenies suggest homology between cnidarian allorecognition genes, as well as syntenic conservation of this gene complex across hydrozoans. Ongoing comparative genomic analyses across bilaterian species aim to clarify the relationship between the ARC and the MHC, thereby informing our understanding of the evolutionary origins of bilaterian adaptive immunity.

Hematopoiesis sustains lifelong blood production and immunity and is strongly influenced by biological sex. Males and females differ in their hematopoietic composition and susceptibility to infection and hematologic disease. Because current models cannot readily disentangle the independent contributions of sex chromosomes versus hormones, the mechanisms underlying sexual dimorphism remain poorly understood. Zebrafish have proven to be a powerful, translatable system for studying hematopoiesis, and offer a unique opportunity to isolate hormonal contributions due to their absence of sex chromosomes. However, whether sex is regulated in the zebrafish hematopoietic system remained unknown. Here, we discovered pronounced sexual dimorphism in the adult zebrafish marrow, with compositional differences peaking at reproductive age. Transcriptome analysis of marrows revealed conservation of sex-specific programs between zebrafish and human datasets. Upon immune challenge with lipopolysaccharide, all major blood lineages exhibited sex-biased changes in their composition and transcriptomes with the largest responses in erythroid and myeloid cell types. Among these sex differences was an unexpected immunoregulation of the erythroid lineage across maturation states. These findings establish biological sex as a major regulator of zebrafish hematopoiesis and immunity, paralleling humans, and lay the groundwork necessary for future mechanistic studies on sexual dimorphism in hematologic health and disease.

Two rounds of whole-genome duplication (2R-WGD) in the vertebrate lineage are thought to have driven gene diversification. Mammals have four genes encoding five nebulin superfamily proteins (LASP1, LASP2/NEBL, NEB, and NRAP) hypothesised to have arisen from invertebrate LASP via WGD. Here, we examined the evolutionary fate of the nebulin superfamily after an additional WGD, teleost-specific genome duplication (TGD). Genome and transcriptome analysis of zebrafish (Danio rerio) and medaka (Oryzias latipes), of which the complete genomes were recently determined, showed only lasp family genes retained duplicate copies. Zebrafish retained two LASP2/NEBL-derived ohnologs; one (LASP2L) had been misidentified as a nebulette gene in the Zebrafish Information Network (ZFIN) despite lacking nebulette-specific exons. Another (LASP2/NEBL) is newly identified in this study because it lies in a newly determined region of the latest genome assembly. Contrastingly, medaka retained two LASP1 ohnologs and one LASP2/NEBL. The duplicated lasp family proteins showed overlapping and specific tissue expression patterns. We confirmed all three zebrafish lasp proteins bound actin and co-localized with actin-rich structures when ectopically expressed in cultured fibroblasts. This supports that TGD-derived lasp duplicates underwent lineage-specific retention and expression divergence in teleosts whilst conserving core actin-binding properties.

Introduction: Neoepitopes form the basis of tumor-specific immune responses. Tissue biopsy, the primary source for neoepitope detection, is limited and invasive. Therefore, we aimed to identify neoepitopes by mutanome-guided immunopeptidomics from plasma of cancer patients. Methods: Mass spectrometry (MS) data analysis of HLA ligands from plasma (n = 4) was guided by patient-specific mutanomes of cell-free DNA (cfDNA) from plasma or tumor genomic DNA (tgDNA) from tissue. Matched tumor tissue and healthy donor plasma served as controls. Neoepitopes were validated with synthetic peptides, and immunogenicity was assessed using IFN-gamma ELISpot and intracellular cytokine staining. Results: Wild-type immunopeptidomes from tissue and plasma overlapped by 58%, with 91% of plasma HLA ligands rediscovered in tissue. 13 out of 15 tumor-associated HLA ligands detected in plasma were rediscovered in the matching tissue. However, no neoepitopes in plasma were identified by immunopeptidomics guided by cfDNA mutanomes, likely reflecting the limited overlap between cfDNA and tgDNA mutanomes (15%). Using the tgDNA mutanome as a complementary reference, two neoepitopes were detected in one patient's plasma, albeit at the MS detection limit. Both neoepitopes were also discovered in tissue, along with three tissue-exclusive neoepitopes. Two tissue-exclusive neoepitopes induced antigen-specific T cell responses in healthy donor PBMCs. Conclusion: In summary, plasma immunopeptidomics enables profiling of HLA ligands from wild-type proteins, including TAAs. In principle, neoepitope detection from plasma at the peptide level is feasible, but tissue remains the gold standard for variant calling and neoepitope identification. Improved detection methods may enable minimally invasive approaches in the future.

Skeletal dysplasias can occur in animals, including humans, and are attributed to both environmental and genetic causes. As our understanding of the distinguishing phenotypes of different skeletal dysplasias improves, alongside elucidation of the mechanisms and genetic networks that contribute to skeletal development, identifying and controlling genetic mutations in livestock populations has become increasingly possible, including approaches beyond selective breeding alone. Within Scotland an 'achondroplasia' syndrome that was first described in the South Country Cheviot breed of sheep (OMIA 000004-9940), has continued to occur periodically within the same geographical locations, despite culling of affected lambs over many decades. Affected lambs exhibit rounded heads, short legs, a loss of hooves, abnormal stance and short twisted tails. As the condition has been observed in flocks and as progeny from specific ram and ewe matings, and cases have been co-twinned to normal siblings, this suggests a syndrome caused by a recessive genetic condition rather than environmental factors. Through homozygosity mapping and whole genome sequence analyses, we identified a 10 bp deletion in the open reading frame of the PRICKLE1 gene that encodes prickle planar cell polarity protein 1 and is essential for normal vertebrate embryonic development. We identified additional phenotypes present in affected lambs, including a loss of phalanges, neural tube defects, vertebral and articular chondrocyte defects, hair orientation differences and polycystic kidneys, consistent with a loss of PRICKLE1 function based on comparisons to mouse mutant lines. These findings enabled the development of a simple DNA test for the identified PRICKLE1 mutation that was made commercially available to facilitate the elimination of this genetic disease from sheep flocks.

Public data sharing in immunological and clinical research is essential for transparency, reproducibility, and secondary analysis. However, many datasets in repositories such as ImmPort lack sufficient structured data to enable reuse. ImmPort's Assessments and Lab Tests templates -- which capture critical clinical information -- are often incomplete or omitted due to their complexity and the manual effort required. To address this challenge, an application with a graphical user interface requiring no programming was developed to facilitate population of these templates. The tool streamlines the submission process by reducing time and complexity for data contributors while ensuring template formatting. It integrates ImmPort files, study-specific data files, and a curated data dictionary to automate the mapping of study variables to the required fields within each template. The tool was evaluated on seven ImmPort studies with diverse designs. It successfully generated submission-ready templates that passed ImmPort's internal validation and will be incorporated into the repository. By improving the feasibility and completeness of templates, this approach enhances the quality, findability, and reusability of data within ImmPort.

Aquaporin-10 (AQP10) is an aquaglyceroporin that transports small, uncharged molecules, such as water and glycerol. In some species, such as mice and cattle, the AQP10 gene (Aqp10) is pseudogenized. In rodents, pseudogenization of Aqp10 is limited to the Myomorpha suborder, while Aqp10 remains intact in other rodent lineages. The order Cetartiodactyla encompasses diverse species, including Suina, Tylopoda, Ruminantia, and Whippomorpha. The distribution of intact Aqp10 across the entire Cetartiodactyla and analyses estimating the timing of Aqp10 pseudogenization have not been sufficiently investigated. In this study, we identified species within Cetartiodactyla that retain intact Aqp10 and analyzed the detailed patterns of pseudogenization to elucidate the specifics of pseudogenization within this order. We collected Aqp10 loci by BLAST and synteny analyses of publicly available genome databases of 51 cetartiodactyl species from 17 families. Dot plot and alignment analyses revealed that seudogenization of Aqp10 occurred in 42 Cetartiodactyla species, including the Camelidae, warthog (Suidae), Antilocapridae, okapi (Giraffidae), Cervidae, Bovidae, Hippopotamidae, and Cetacea. In contrast, nine species retained an intact Aqp10: pig and babirusa (Suidae), peccary (Tayassuidae), mouse deeres (Tragulidae), giraffes (Giraffidae), and musk deers (Moschidae). The timing of pseudogenization was estimated by comparing patterns of exon deletions, frameshifts, and nonsense mutations in Aqp10 pseudogenes. RT-PCR analysis of cattle and sheep tissues revealed that these species express Aqp3 and Aqp7 in their intestines, but not Aqp10. Collectively, these findings indicate that Aqp10 pseudogenization occurred independently in multiple cetartiodactyl lineages and provide new insights into the evolutionary history of Aqp10 loss within this order.

Coleoid cephalopods (cuttlefish, octopus, and squid) are marine mollusks with elaborate nervous systems that support a diverse repertoire of complex behaviors. These include the neural control of the color, pattern, and texture of the skin, facilitating both adaptive camouflage and innate patterning that may reflect internal state. The development of transgenic cephalopods expressing fluorescent proteins, optogenetic actuators, and reporters of neural activity would contribute a new and important technology to cephalopod biology. The generation of transgenic cephalopods, however, has remained a major challenge. Here, we report the development of stable transgenic dwarf cuttlefish (Ascarosepion bandense) expressing ubiquitous nuclear-localized mScarlet, a red fluorescent protein. We evaluated multiple strategies for transgenesis, and established cuttlefish lines using both CRISPR and the transposons Sleeping Beauty and Minos. The stable expression of transgenes enabled live imaging of cell dynamics during embryonic development. The Minos transposon emerged as the most efficient transgenesis strategy and is adaptable to promoters and transgenes of choice. These strategies now enable the generation of diverse genetic tools for mechanistic studies of cephalopod biology.

Background: Giant cell arteritis (GCA), traditionally considered T cell-driven, exhibits vascular B/plasma-cell infiltrates and altered circulating B-cell homeostasis, suggesting a pathogenic role for B cells. Relapse is common in newly diagnosed GCA (nGCA), yet predictive biomarkers remain elusive. Objective: To characterize circulating B-cell subpopulations in nGCA and assess their association with 12-month relapse. Methods: Baseline multiparametric flow cytometry was performed on peripheral blood samples from consecutive patients with nGCA and sex/age-matched healthy controls (HC). Patients received standard treatment per 2018 EULAR recommendations and were followed longitudinally. Results: 102 patients with nGCA and 102 matched HC were included. One-year outcome data were available for 55 patients, of whom 34 relapsed and 21 remained relapse-free. Baseline B-cell homeostasis differed by disease course. Non-relapsers showed reduced total CD19+ B-cell counts and contraction of most B-cell subsets, while CD27-IgD-double-negative (DN) B-cell counts remained comparable to HC. Conversely, relapsers maintained total B-cell numbers above HC despite reduced memory compartments, driven by expanded naive and DN B cells, alongside preserved plasmablast and transitional B-cell counts. At 6 months, total B-cell counts increased only in non-relapsers. DN B cells showed the greatest divergence between groups and independently predicted relapse. DN B-cell counts also correlated with circulating T peripheral helper cell numbers, supporting extrafollicular immune pathway activation. Conclusion: nGCA is characterized by outcome-dependent disturbances in B-cell homeostasis. Expansion of the DN B-cell compartment identifies a relapse-prone immunophenotype and represents a clinically relevant biomarker for early risk stratification.

The non-canonical NF-kappa B pathway, in which NFKB2 plays a central role, is critical for lymphoid organ development and immune homeostasis. Pathogenic NFKB2 variants have been strongly associated with a broad spectrum of human pathologies, notably immunodysregulation syndromes, characterized by recurrent infections, autoantibody production and endocrine abnormalities. The mechanisms underlying the resulting immune defects remain incompletely understood. In the present study, we characterized the immune phenotype of a novel rat model harboring a Nfkb2 non-sense mutation located in the c-terminal region of the gene (c.2602T>A; p.Tyr867*, Nfkb2p.Tyr867*). This variant has been identified in patients with a familial form of common variable immunodeficiency (CVID) presenting with autoimmunity manifestations. Using this model, we demonstrate that mutated NFKB2 resulted in profound alterations of the thymic microenvironment, including reduced K5+ medullary thymic epithelial cells and decreased AIRE expression, accompanied by altered thymocyte development and impaired thymic regulatory T populations. Peripheral lymphoid organization was also profoundly affected, with absence of lymph nodes, defects in the splenic marginal zone, and alterations of CD4+ and CD8+ T-cell homeostasis, including impaired Treg populations and reduced ICOS expression on CD4+ T cells. The B-cell compartment was markedly altered, with impaired B-cell maturation, characterized by accumulation of IgM+IgD- transitional cells and a marked reduction of IgD+ mature B cells. These abnormalities were associated with impaired humoral immunity and complete absence of circulating immunoglobulins. Collectively, our findings define a systemic combined immunodeficiency phenotype in the Nfkb2-mutated rat, highlighting its utility for studying the role of non-canonical NF-kB signaling in immune system organization and adaptative immunity.

Introduction: Tumor-reactive T cell receptors (TCRs) are promising candidates for the manufacturing of safe and effective adoptive TCR-T cell therapies, but their efficient identification remains a major challenge. We previously established an antigen-agnostic TCR-repertoire analysis workflow that prioritizes candidate TCRs based on intratumoral clonal expansion and high tumor-to-non-tumor (T/N) frequency ratios. Here, we investigated whether T/N-selected clonotypes show independent transcriptional and spatial features consistent with tumor association. Methods: Bulk and single-cell TCR repertoire sequencing of paired tumor and adjacent non-tumor tissues was used to identify clonotypes with high T/N ratios across 15 cancer cases, comprising six non-small cell lung cancers, five pancreatic cancers, two breast cancers, and two colorectal cancers. Candidate clonotypes were further characterized by paired single-cell TCR and gene expression profiling and, in two PDAC cases, by Xenium spatial transcriptomics to assess their cellular phenotype, spatial localization, and proximity to malignant cells. Results: High T/N ratios consistently marked T cell clonotypes with features of tumor reactivity across the studied solid tumor entities. The analysis included Cluster-A, a lead T cell cluster comprising highly similar TCRs shared among 14 HLA-A*02:01-positive patients in the present cohort. Cluster-A clonotypes were consistently enriched in tumor relative to matched non-tumor tissue, and previous functional studies demonstrated tumor recognition by Cluster-A members. In a focused analysis of a pancreatic cancer case, Cluster-A CD8 T cells were selectively detected in tumor tissue, preferentially localized near malignant epithelial cells, and displayed an antigen-experienced effector-memory phenotype with prominent expression of granzymes A and K (GZMA/GZMK). Importantly, in another PDAC case that lacked Cluster-A clonotypes, these features extended to independent clonotypes, revealing a continuous tumor-association axis in which increasing T/N ratios were associated with progressively closer localization to malignant cells and increased GZMA/GZMK expression. Thus, repertoire enrichment, spatial tumor proximity, and effector-associated transcriptional states converged across independent tumor-enriched clonotypes. Outlook: These findings support the T/N ratio as a robust primary criterion for identifying candidate tumor-reactive TCRs. Integration of single-cell and spatial transcriptomics provides supporting biological evidence of tumor association and a framework for prioritizing therapeutic TCR candidates when functional validation is limited.

The teleost whole-genome duplication (TGD) contributed to functional diversification of opsins. Some TGD paralogs, including those of parapinopsin (PP), Vertebrate Ancient (VA) opsin, and long wavelength-sensitive (LWS) opsin, show different absorption spectra and/or expression patterns. However, our knowledge of detailed evolutionary processes and mechanisms by which TGD contributed to opsin diversification is still limited. Here, we report that TGD paralogs of parietopsin (PT1 and PT2) are retained in extant species. While most species, including the zebrafish, have only PT1, several characins, including the red piranha (Pygocentrus nattereri), retain both PT1 and PT2. The Mexican tetra (Astyanax mexicanus) and catfishes have only PT2. To assess the degree of functional diversification between PT1 and PT2, we characterized spectral properties and expression patterns. Maximum absorption spectra differ slightly among PTs. Those of red piranha PT1, PT2, Mexican tetra PT2, and Japanese catfish (Silurus asotus) PT2 were located at 517 nm, 528 nm, 517 nm, and ~535 nm, respectively. Fluorescence in situ hybridization showed that (1) piranha PT1 and PT2 are expressed in the same pineal cells, and (2) they are also co-expressed with PP1. Based on these findings, functional diversification of opsins after gene duplication is discussed.

Arctic ecosystems are undergoing rapid climate warming, providing an opportunity to examine molecular evolution under environmental stress. Here, we used whole-genome resequencing to assess population structure, genetic diversity, and recent demographic history in the snowy owl (Bubo scandiacus), a highly mobile Arctic apex predator classified as Vulnerable by the IUCN. We analysed genome-wide data from individuals sampled across the species' circumpolar range using a chromosome-level reference genome. Population genomic analyses revealed an absence of detectable population structure, indicating extensive gene flow and near-panmixia across continents. Genome-wide nucleotide diversity was low ({pi}= 2.47 x 10^(-4)), while runs of homozygosity were dominated by short tracts, indicating limited standing genetic variation despite little evidence for recent close inbreeding. Our demographic inference revealed a sustained decline in effective population size over the past ~1,600 years, a pronounced reduction coinciding with the Medieval Warm Period followed by a partial recovery during the subsequent cooler interval. This direct, measurable effect of past climate warming on effective population size raises concerns for the long-term survival of the snowy owl. Our results demonstrate that high connectivity and low levels of recent inbreeding have not necessarily buffered this species against climate-driven demographic erosion of genetic diversity. This study highlights the value of using whole-genome population genomics to inform genomic monitoring and conservation decision-making in rapidly changing environments.

Introduction: Killer-cell immunoglobulin-like receptor (KIR) gene content shapes natural killer-cell responses and varies substantially among populations. We examined whether worldwide KIR diversity reflects regional structure, isolation by distance, serial founder effects, or latitudinal gradients. Data and methods: Carrier frequencies for 11 variable KIR genes were obtained from the Allele Frequency Net Database. The primary analysis included 189 locally sampled, non-diaspora, non-mixed populations from 70 countries, comprising 25,469 individuals; 161 populations had complete 11-gene profiles. Diversity indices, ordination, clustering, permutation-based variance decomposition, Mantel tests, and regression models were used to characterize geographic structure. Serial founder effects were evaluated using four distance measures from East Africa, both globally and within the Old World. Models of latitude incorporated spatial autocorrelation. Prespecified sensitivity analyses evaluated alternative population sets, gene panels, geographic classifications, weighting schemes, and multivariate outlier exclusion. Results: Geographic region explained 53.5% of variation in KIR profiles, substantially exceeding the permutation-null expectation of 6.9%, although multivariate dispersion differed among regions. Twenty-five of 28 pairwise regional comparisons remained significant after false-discovery-rate correction. Isolation by distance was modest (Mantel r=0.3003) and weakened after controlling for region (r=0.1557). KIR diversity did not decline with distance from East Africa globally but decreased strongly within the Old World (expected heterozygosity {beta}=-0.0066 and standardized gene content {beta}=-0.0789). B-content decreased nominally with absolute latitude, but latitude was collinear with region and the association survived neither spatial modelling nor sensitivity analyses. Regional differentiation remained robust across analytical specifications, whereas the Old World decline was carried by gene content once standard errors were clustered by country. Conclusion: Worldwide KIR gene-content variation shows pronounced regional differentiation, modest isolation by distance, and no independent latitudinal gradient. Diversity declines with distance from East Africa across the Old World, supporting serial founder effects, but not once the Americas and Oceania are included.

We investigated the genetic basis of a previously unclassified congenital vascular anomaly in a patient with multifocal cutaneous vascular tumors. Through genomic analysis, we identified a novel heterozygous de novo germline variant in GA-binding protein-alpha (GABPA), an ETS family transcription factor. Histopathologic evaluation demonstrated capillary-venous lesions lacking glucose transporter 1 expression, distinguishing this entity from common infantile hemangioma. Notably, the tumor vasculature exhibited a prominent alpha-SMA-positive perivascular cell layer, with nuclear localized GABPA detected in endothelial and perivascular cells. The GABPA c.509T>G (p.L170R) variant is absent from large population databases, including gnomAD, DeCAF, and RGC-MCPS databases, indicating it is extremely rare. The affected leucine residue is highly conserved across species and located within the Pointed (PNT) domain, a critical region mediating GABPA protein-protein interactions. AlphaMissense predicts this substitution to be highly damaging. Structure-based analysis using an AlphaFold-predicted approach, combined with in silico mutagenesis and molecular interaction analyses, revealed that the L170R substitution introduces new electrostatic interactions while disrupting native hydrophobic contacts within the PNT domain. Functional assessment in zebrafish demonstrated that mosaic stromal expression of GABPA-L170R caused abnormal vascular architecture compared with GABPA-wild-type controls, establishing a direct link between this variant and disrupted vascular development. Collectively, the genetic rarity, evolutionary conservation, predicted structural perturbation, and disruptive effects on vascular development establish GABPA-L170R as likely disease-causing variant. These findings define a new molecular mechanism underlying congenital vascular tumor formation and expand the role for ETS-family transcriptional regulation in human vascular anomalies.