Systematic characterization of immunoglobulin loci and deep sequencing of the expressed repertoire in the Atlantic cod (Gadus morhua).

BACKGROUND: The Atlantic cod is a prolific species in the Atlantic, despite its inconsistent specific antibody response. It presents a peculiar case within vertebrate immunology due to its distinct immune system, characterized by the absence of MHCII antigen presentation pathway, required for T cell-dependent antibody responses. Thorough characterisation of immunoglobulin loci and analysis of the antibody repertoire is necessary to further our understanding of the Atlantic cod's immune response on a molecular level. RESULTS: A comprehensive search of the cod genome (gadmor3.0) identified the complete set of IgH genes organized into three sequential translocons on chromosome 2, while IgL genes were located on chromosomes 2 and 5. The Atlantic cod displayed a moderate germline V gene diversity, comprising four V gene families for both IgH and IgL, each with distinct chromosomal locations and organizational structures. 5'RACE sequencing revealed a diverse range of heavy chain CDR3 sequences and relatively limited CDR3 diversity in light chains. The analysis highlighted a differential impact of V-gene germline CDR3 length on receptor CDR3 length between heavy and light chains, underlining different recombination processes. CONCLUSIONS: This study reveals that the Atlantic cod, despite its inconsistent antibody response, maintains a level of immunoglobulin diversity comparable to other fish species. The findings suggest that the extensive recent duplications of kappa light chain genes do not result in increased repertoire diversity. This research provides a comprehensive view of the Atlantic cod's immunoglobulin gene organization and repertoire, necessary for future studies of antibody responses at the molecular level.

High-throughput sequencing of insect specimens with sub-optimal DNA preservation using a practical, plate-based Illumina-compatible Tn5 transposase library preparation method.

Entomological sampling and storage conditions often prioritise efficiency, practicality and conservation of morphological characteristics, and may therefore be suboptimal for DNA preservation. This practice can impact downstream molecular applications, such as the generation of high-throughput genomic libraries, which often requires substantial DNA input amounts. Here, we use a practical Tn5 transposase tagmentation-based library preparation method optimised for 96-well plates and low yield DNA extracts from insect legs that were stored under sub-optimal conditions for DNA preservation. The samples were kept in field vehicles for extended periods of time, before long-term storage in ethanol in the freezer, or dry at room temperature. By reducing DNA input to 6ng, more samples with sub-optimal DNA yields could be processed. We matched this low DNA input with a 6-fold dilution of a commercially available tagmentation enzyme, significantly reducing library preparation costs. Costs and workload were further suppressed by direct post-amplification pooling of individual libraries. We generated medium coverage (>3-fold) genomes for 88 out of 90 specimens, with an average of approximately 10-fold coverage. While samples stored in ethanol yielded significantly less DNA compared to those which were stored dry, these samples had superior sequencing statistics, with longer sequencing reads and higher rates of endogenous DNA. Furthermore, we find that the efficiency of tagmentation-based library preparation can be improved by a thorough post-amplification bead clean-up which selects against both short and large DNA fragments. By opening opportunities for the use of sub-optimally preserved, low yield DNA extracts, we broaden the scope of whole genome studies of insect specimens. We therefore expect these results and this protocol to be valuable for a range of applications in the field of entomology.

Gluten-Free Diet Induces Rapid Changes in Phenotype and Survival Properties of Gluten-Specific T Cells in Celiac Disease.

BACKGROUND & AIMS: The treatment of celiac disease (CeD) with gluten-free diet (GFD) normalizes gut inflammation and disease-specific antibodies. CeD patients have HLA-restricted, gluten-specific T cells persisting in the blood and gut even after decades of GFD, which are reactivated and disease driving upon gluten exposure. Our aim was to examine the transition of activated gluten-specific T cells into a pool of persisting memory T cells concurrent with normalization of clinically relevant biomarkers during the first year of treatment. METHODS: We followed 17 CeD patients during their initial GFD year, leading to disease remission. We assessed activation and frequency of gluten-specific CD4+ blood and gut T cells with HLA-DQ2.5:gluten tetramers and flow cytometry, disease-specific serology, histology, and symptom scores. We assessed gluten-specific blood T cells within the first 3 weeks of GFD in 6 patients and serology in an additional 9 patients. RESULTS: Gluten-specific CD4+ T cells peaked in blood at day 14 while up-regulating Bcl-2 and down-regulating Ki-67 and then decreased in frequency within 10 weeks of GFD. CD38, ICOS, HLA-DR, and Ki-67 decreased in gluten-specific cells within 3 days. PD-1, CD39, and OX40 expression persisted even after 12 months. IgA-transglutaminase 2 decreased significantly within 4 weeks. CONCLUSIONS: GFD induces rapid changes in the phenotype and number of gluten-specific CD4+ blood T cells, including a peak of nonproliferating, nonapoptotic cells at day 14. Subsequent alterations in T-cell phenotype associate with the quiescent but chronic nature of treated CeD. The rapid changes affecting gluten-specific T cells and disease-specific antibodies offer opportunities for clinical trials aiming at developing nondietary treatments for patients with newly diagnosed CeD.

Identification of a pharyngeal mucosal lymphoid organ in zebrafish and other teleosts: Tonsils in fish?

The constant exposure of the fish branchial cavity to aquatic pathogens causes local mucosal immune responses to be extremely important for their survival. Here, we used a marker for T lymphocytes/natural killer (NK) cells (ZAP70) and advanced imaging techniques to investigate the lymphoid architecture of the zebrafish branchial cavity. We identified a sub-pharyngeal lymphoid organ, which we tentatively named "Nemausean lymphoid organ" (NELO). NELO is enriched in T/NK cells, plasma/B cells, and antigen-presenting cells embedded in a network of reticulated epithelial cells. The presence of activated T cells and lymphocyte proliferation, but not V(D)J recombination or hematopoiesis, suggests that NELO is a secondary lymphoid organ. In response to infection, NELO displays structural changes including the formation of T/NK cell clusters. NELO and gill lymphoid tissues form a cohesive unit within a large mucosal lymphoid network. Collectively, we reveal an unreported mucosal lymphoid organ reminiscent of mammalian tonsils that evolved in multiple teleost fish families.

Lymphocyte subsets in Atlantic cod (Gadus morhua) interrogated by single-cell sequencing.

Atlantic Cod (Gadus morhua) has lost the major histocompatibility complex class II presentation pathway. We recently identified CD8-positive T cells, B cells, and plasma cells in cod, but further characterisation of lymphocyte subsets is needed to elucidate immune adaptations triggered by the absence of CD4-positive T lymphocytes. Here, we use single-cell RNA sequencing to examine the lymphocyte heterogeneity in Atlantic cod spleen. We describe five T cell subsets and eight B cell subsets and propose a B cell trajectory of differentiation. Notably, we identify a subpopulation of T cells that are CD8-negative. Most of the CD8-negative T lymphocytes highly express the homologue of monocyte chemotactic protein 1b, and another subset of CD8-negative T lymphocytes express the homologue of the scavenger receptor m130. Uncovering the multiple lymphocyte cell sub-clusters reveals the different immune states present within the B and T cell populations, building a foundation for further work.

Stereotyped B-cell responses are linked to IgG constant region polymorphisms in multiple sclerosis.

Clonally related B cells infiltrate the brain, meninges, and cerebrospinal fluid of MS patients, but the mechanisms driving the B-cell response and shaping the immunoglobulin repertoires remain unclear. Here, we used single-cell full-length RNA-seq and BCR reconstruction to simultaneously assess the phenotypes, isotypes, constant region polymorphisms, and the paired heavy- and light-chain repertoires in intrathecal B cells. We detected extensive clonal connections between the memory B cell and antibody-secreting cell (ASC) compartments and observed clonally related cells of different isotypes including IgM/IgG1, IgG1/IgA1, IgG1/IgG2, and IgM/IgA1. There was a strong dominance of the G1m1 allotype constant region polymorphisms in ASCs, but not in memory B cells. Tightly linked to the G1m1 allotype, we found a preferential pairing of the immunoglobulin heavy-chain variable (IGHV)4 gene family with the κ variable (IGKV)1 gene family. The IGHV4-39 gene was most used and showed the highest frequency of pairing with IGKV1-5 and IGKV1(D)-33. These results link IgG constant region polymorphisms to stereotyped B-cell responses in MS and indicate that the intrathecal B-cell response in these patients could be directed against structurally similar epitopes.

TCRpower: quantifying the detection power of T-cell receptor sequencing with a novel computational pipeline calibrated by spike-in sequences.

T-cell receptor (TCR) sequencing has enabled the development of innovative diagnostic tests for cancers, autoimmune diseases and other applications. However, the rarity of many T-cell clonotypes presents a detection challenge, which may lead to misdiagnosis if diagnostically relevant TCRs remain undetected. To address this issue, we developed TCRpower, a novel computational pipeline for quantifying the statistical detection power of TCR sequencing methods. TCRpower calculates the probability of detecting a TCR sequence as a function of several key parameters: in-vivo TCR frequency, T-cell sample count, read sequencing depth and read cutoff. To calibrate TCRpower, we selected unique TCRs of 45 T-cell clones (TCCs) as spike-in TCRs. We sequenced the spike-in TCRs from TCCs, together with TCRs from peripheral blood, using a 5' RACE protocol. The 45 spike-in TCRs covered a wide range of sample frequencies, ranging from 5 per 100 to 1 per 1 million. The resulting spike-in TCR read counts and ground truth frequencies allowed us to calibrate TCRpower. In our TCR sequencing data, we observed a consistent linear relationship between sample and sequencing read frequencies. We were also able to reliably detect spike-in TCRs with frequencies as low as one per million. By implementing an optimized read cutoff, we eliminated most of the falsely detected sequences in our data (TCR α-chain 99.0% and TCR ß-chain 92.4%), thereby improving diagnostic specificity. TCRpower is publicly available and can be used to optimize future TCR sequencing experiments, and thereby enable reliable detection of disease-relevant TCRs for diagnostic applications.

Innate lymphoid cell characterization in the rat and their correlation to gut commensal microbes.

Innate lymphoid cells (ILCs) are important for tissue immune homeostasis, and are thoroughly characterized in mice and humans. Here, we have performed in-depth characterization of rat ILCs. Rat ILCs were identified based on differential expression of transcription factors and lack of lineage markers. ILC3s represented the major ILC population of the small intestine, while ILC2s were infrequent but most prominent in liver and adipose tissue. Two major subsets of group 1 ILCs were defined. Lineage- T-bet+ Eomes+ cells were identified as conventional NK cells, while lineage- T-bet+ Eomes- cells were identified as the probable rat counterpart of ILC1s based on their selective expression of the ILC marker CD200R. Rat ILC1s were particularly abundant in liver and intestinal tissues, and were functionally similar to NK cells. Single-cell transcriptomics of spleen and liver cells confirmed the main division of NK cells and ILC1-like cells, and demonstrated Granzyme A as an additional ILC1 marker. We further report differential distributions of NK cells and ILCs along the small and large intestines, and the association of certain bacterial taxa to frequencies of ILCs. In conclusion, we provide a framework for future studies of ILCs in diverse rat experimental models, and novel data on the potential interplay between commensals and intestinal ILCs.

Differential expression profile of gluten-specific T cells identified by single-cell RNA-seq.

Gluten-specific CD4+ T cells drive the pathogenesis of celiac disease and circulating gluten-specific T cells can be identified by staining with HLA-DQ:gluten tetramers. In this first single-cell RNA-seq study of tetramer-sorted T cells from untreated celiac disease patients blood, we found that gluten-specific T cells showed distinct transcriptomic profiles consistent with activated effector memory T cells that shared features with Th1 and follicular helper T cells. Compared to non-specific cells, gluten-specific T cells showed differential expression of several genes involved in T-cell receptor signaling, translational processes, apoptosis, fatty acid transport, and redox potentials. Many of the gluten-specific T cells studied shared T-cell receptor with each other, indicating that circulating gluten-specific T cells belong to a limited number of clones. Moreover, the transcriptional profiles of cells that shared the same clonal origin were transcriptionally more similar compared with between clonally unrelated gluten-specific cells.

Frequency of Gluten-Reactive T Cells in Active Celiac Lesions Estimated by Direct Cell Cloning.

Chronic inflammation of the small intestine in celiac disease is driven by activation of CD4+ T cells that recognize gluten peptides presented by disease-associated HLA-DQ molecules. We have performed direct cell cloning of duodenal biopsies from five untreated and one refractory celiac disease patients, and three non-celiac disease control subjects in order to assess, in an unbiased fashion, the frequency of gluten-reactive T cells in the disease-affected tissue as well as the antigen fine specificity of the responding T cells. From the biopsies of active disease lesions of five patients, 19 T-cell clones were found to be gluten-reactive out of total 1,379 clones tested. This gave an average of 1.4% (range 0.7% - 1.9%) of gluten-reactive T cells in lamina propria of active celiac lesions. Interestingly, also the patient with refractory celiac disease had gluten-reactive T cell clones in the lamina propria (5/273; 1.8%). In comparison, we found no gluten-reactive T cells in any of the total 984 T-cell clones screened from biopsies from three disease control donors. Around two thirds of the gluten-reactive clones were reactive to a panel of peptides representing known gluten T-cell epitopes, of which two thirds were reactive to the immunodominant DQ2.5-glia-α1/DQ2.5-glia-α2 and DQ2.5-glia-ω1/DQ2.5-glia-ω2 epitopes. This study shows that gluten-reactive T cells in the inflamed duodenal tissue are prevalent in the active disease lesion, and that many of these T cells are reactive to T-cell epitopes that are not yet characterized. Knowledge of the prevalence and epitope specificity of gluten-specific T cells is a prerequisite for therapeutic efforts that target disease-specific T cells in celiac disease.

Comprehensive Analysis of CDR3 Sequences in Gluten-Specific T-Cell Receptors Reveals a Dominant R-Motif and Several New Minor Motifs.

Gluten-specific CD4+ T cells are drivers of celiac disease (CeD). Previous studies of gluten-specific T-cell receptor (TCR) repertoires have found public TCRs shared across multiple individuals, biased usage of particular V-genes and conserved CDR3 motifs. The CDR3 motifs within the gluten-specific TCR repertoire, however, have not been systematically investigated. In the current study, we analyzed the largest TCR database of gluten-specific CD4+ T cells studied so far consisting of TCRs of 3122 clonotypes from 63 CeD patients. We established a TCR database from CD4+ T cells isolated with a mix of HLA-DQ2.5:gluten tetramers representing four immunodominant gluten epitopes. In an unbiased fashion we searched by hierarchical clustering for common CDR3 motifs among 2764 clonotypes. We identified multiple CDR3α, CDR3ß, and paired CDR3α:CDR3ß motif candidates. Among these, a previously known conserved CDR3ß R-motif used by TRAV26-1/TRBV7-2 TCRs specific for the DQ2.5-glia-α2 epitope was the most prominent motif. Furthermore, we identified the epitope specificity of altogether 16 new CDR3α:CDR3ß motifs by comparing with TCR sequences of 231 T-cell clones with known specificity and TCR sequences of cells sorted with single HLA-DQ2.5:gluten tetramers. We identified 325 public TCRα and TCRß sequences of which 145, 102 and 78 belonged to TCRα, TCRß and paired TCRαß sequences, respectively. While the number of public sequences was depended on the number of clonotypes in each patient, we found that the proportion of public clonotypes from the gluten-specific TCR repertoire of given CeD patients appeared to be stable (median 37%). Taken together, we here demonstrate that the TCR repertoire of CD4+ T cells specific to immunodominant gluten epitopes in CeD is diverse, yet there is clearly biased V-gene usage, presence of public TCRs and existence of conserved motifs of which R-motif is the most prominent.

Potential impact of celiac disease genetic risk factors on T cell receptor signaling in gluten-specific CD4+ T cells.

Celiac disease is an auto-immune disease in which an immune response to dietary gluten leads to inflammation and subsequent atrophy of small intestinal villi, causing severe bowel discomfort and malabsorption of nutrients. The major instigating factor for the immune response in celiac disease is the activation of gluten-specific CD4+ T cells expressing T cell receptors that recognize gluten peptides presented in the context of HLA-DQ2 and DQ8. Here we provide an in-depth characterization of 28 gluten-specific T cell clones. We assess their transcriptional and epigenetic response to T cell receptor stimulation and link this to genetic factors associated with celiac disease. Gluten-specific T cells have a distinct transcriptional profile that mostly resembles that of Th1 cells but also express cytokines characteristic of other types of T-helper cells. This transcriptional response appears not to be regulated by changes in chromatin state, but rather by early upregulation of transcription factors and non-coding RNAs that likely orchestrate the subsequent activation of genes that play a role in immune pathways. Finally, integration of chromatin and transcription factor binding profiles suggest that genes activated by T cell receptor stimulation of gluten­specific T cells may be impacted by genetic variation at several genetic loci associated with celiac disease.

Immunobiology and conflicting roles of the human CD161 receptor in T cells.

Human C-type lectin-like CD161 is a type-II transmembrane protein expressed on the surface of various lymphocytes across innate and adaptive immune systems. CD161+ T cells displayed enhanced ability to produce cytokines and were shown to be enriched in the gut. Independently of function, CD161 was used as marker of innate-like T cells and marker of IL-17-producing cells. The function of CD161 is still not fully understood. In T cells, CD161 was proposed to act as co-signalling receptor that influence T-cell receptor-dependent responses. However, conflicting studies were published demonstrating lack of agreement over the role of CD161 during T-cell activation. In this review, we outline phenotypical and functional consequences of CD161 expression in T cells. We provide critical discussion over the most pressing issues including in depth evaluation of the literature concerning CD161 putative co-signalling properties.

Characterization of T-cell receptor transgenic mice recognizing immunodominant HLA-DQ2.5-restricted gluten epitopes.

We created a TCR transgenic mouse with CD4+ T cells recognizing the immunodominant DQ2.5-glia-ω2 gluten epitope. We show that these cells respond to deamidated gluten feed in vivo and compare them to previously published α2- and γ1-specific mice. These mice may help enlighten key aspects of celiac disease pathogenesis.

C-type lectin-like CD161 is not a co-signalling receptor in gluten-reactive CD4 + T cells.

C-type lectin-like CD161, a class II transmembrane protein, is a surface receptor expressed by NK cells and T cells. In coeliac disease, CD161 was expressed more frequently on gluten-reactive CD4 + T cells compared to other memory CD4 + T cells isolated from the same tissue compartment. CD161 is a putative co-signalling molecule that was proposed to act as co-stimulatory receptor in the context of signalling through TCR, but contradicting results were published. In order to understand the role of CD161 in gluten-reactive CD4 + T cells, we combined T cell stimulation assays or T cell proliferation assays with ligation of CD161 and intracellular cytokine staining. We found that CD161 ligation provided neither co-stimulatory nor co-inhibitory signals to modulate proliferation and IFN-γ or IL-21 production by gluten-reactive CD4 + T cell clones. Thus, we suggest that CD161 does not function as a co-signalling receptor in the context of gluten-reactive CD4 + T cells.

T cell receptor repertoire as a potential diagnostic marker for celiac disease.

An individual's T cell repertoire is skewed towards some specificities as a result of past antigen exposure and subsequent clonal expansion. Identifying T cell receptor signatures associated with a disease is challenging due to the overall complexity of antigens and polymorphic HLA allotypes. In celiac disease, the antigen epitopes are well characterised and the specific HLA-DQ2-restricted T-cell repertoire associated with the disease has been explored in depth. By investigating T cell receptor repertoires of unsorted lamina propria T cells from 15 individuals, we provide the first proof-of-concept study showing that it could be possible to infer disease state by matching against a priori known disease-associated T cell receptor sequences.

Longevity, clonal relationship, and transcriptional program of celiac disease-specific plasma cells.

Disease-specific plasma cells (PCs) reactive with transglutaminase 2 (TG2) or deamidated gluten peptides (DGPs) are abundant in celiac disease (CeD) gut lesions. Their contribution toward CeD pathogenesis is unclear. We assessed expression of markers associated with PC longevity in 15 untreated and 26 treated CeD patients in addition to 13 non-CeD controls and performed RNA sequencing with clonal inference and transcriptomic analysis of 3,251 single PCs. We observed antigen-dependent V-gene selection and stereotypic antibodies. Generation of recombinant DGP-specific antibodies revealed a key role of a heavy chain residue that displays polymorphism, suggesting that immunoglobulin gene polymorphisms may influence CeD-specific antibody responses. We identified transcriptional differences between CeD-specific and non-disease-specific PCs and between short-lived and long-lived PCs. The short-lived CD19+CD45+ phenotype dominated in untreated and short-term-treated CeD, in particular among disease-specific PCs but also in the general PC population. Thus, the disease lesion of untreated CeD is characterized by massive accumulation of short-lived PCs that are not only directed against disease-specific antigens.

Single-Cell Transcriptome Profiling of Immune Cell Repertoire of the Atlantic Cod Which Naturally Lacks the Major Histocompatibility Class II System.

The Atlantic cod's unusual immune system, entirely lacking the Major Histocompatibility class II pathway, has prompted intriguing questions about what mechanisms are used to combat bacterial infections and how immunological memory is generated. By single-cell RNA sequencing we here report an in-depth characterisation of cell types found in immune tissues, the spleen and peripheral blood leukocytes of Atlantic cod. Unbiased transcriptional clustering revealed eleven distinct immune cell signatures. Resolution at the single cell level enabled characterisation of the major cell subsets including the cytotoxic T cells, B cells, erythrocytes, thrombocytes, neutrophils, and macrophages. Additionally, to our knowledge we are the first to uncover cell subsets in Atlantic cod which may represent dendritic cells, natural killer-like cells, and a population of cytotoxic cells expressing GATA-3, a master transcription factor of T helper 2 cells. We further identify putative gene markers for each cluster and describe the relative proportions of each cell type in the spleen and peripheral blood leukocytes. Of the major haematopoietic cell populations, the lymphocytes make up 55 and 68% of the spleen and peripheral blood leukocytes respectively, while the myeloid cells make up 45 and 32%. By single-cell analysis, this study provides the most detailed molecular and cellular characterisation of the immune system of the Atlantic cod so far.

A molecular basis for the T cell response in HLA-DQ2.2 mediated celiac disease.

The highly homologous human leukocyte antigen (HLA)-DQ2 molecules, HLA-DQ2.5 and HLA-DQ2.2, are implicated in the pathogenesis of celiac disease (CeD) by presenting gluten peptides to CD4+ T cells. However, while HLA-DQ2.5 is strongly associated with disease, HLA-DQ2.2 is not, and the molecular basis underpinning this differential disease association is unresolved. We here provide structural evidence for how the single polymorphic residue (HLA-DQ2.5-Tyr22α and HLA-DQ2.2-Phe22α) accounts for HLA-DQ2.2 additionally requiring gluten epitopes possessing a serine at the P3 position of the peptide. In marked contrast to the biased T cell receptor (TCR) usage associated with HLA-DQ2.5-mediated CeD, we demonstrate with extensive single-cell sequencing that a diverse TCR repertoire enables recognition of the immunodominant HLA-DQ2.2-glut-L1 epitope. The crystal structure of two CeD patient-derived TCR in complex with HLA-DQ2.2 and DQ2.2-glut-L1 (PFSEQEQPV) revealed a docking strategy, and associated interatomic contacts, which was notably distinct from the structures of the TCR:HLA-DQ2.5:gliadin epitope complexes. Accordingly, while the molecular surfaces of the antigen-binding clefts of HLA-DQ2.5 and HLA-DQ2.2 are very similar, differences in the nature of the peptides presented translates to differences in responding T cell repertoires and the nature of engagement of the respective antigen-presenting molecules, which ultimately is associated with differing disease penetrance.