Characterization of the canine immunoglobulin heavy chain repertoire by next generation sequencing

Application of Next-Generation Sequencing (NGS) Techniques for Selected Companion Animals

Next-Generation Sequencing (NGS) techniques have revolutionized veterinary medicine for cats and dogs, offering insights across various domains. In veterinary parasitology, NGS enables comprehensive profiling of parasite populations, aiding in understanding transmission dynamics and drug resistance mechanisms. In infectious diseases, NGS facilitates rapid pathogen identification, characterization of virulence factors, and tracking of outbreaks. Moreover, NGS sheds light on metabolic processes by elucidating gene expression patterns and metabolic pathways, essential for diagnosing metabolic disorders and designing tailored treatments. In autoimmune diseases, NGS helps identify genetic predispositions and molecular mechanisms underlying immune dysregulation. Veterinary oncology benefits from NGS through personalized tumor profiling, mutation analysis, and identification of therapeutic targets, fostering precision medicine approaches. Additionally, NGS plays a pivotal role in veterinary genetics, unraveling the genetic basis of inherited diseases and facilitating breeding programs for healthier animals. Physiological investigations leverage NGS to explore complex biological systems, unraveling gene-environment interactions and molecular pathways governing health and disease. Application of NGS in treatment planning enhances precision and efficacy by enabling personalized therapeutic strategies tailored to individual animals and their diseases, ultimately advancing veterinary care for companion animals.

Exploring the sheep (Ovis aries) immunoglobulin repertoire by next generation sequencing

Next-generation sequencing (NGS) has revolutionized the way we determine the antibody repertoires encoded by B cells in the blood or lymphoid organs and transformed our understanding of adaptive immune responses in many species. Sheep (Ovis aries) have been widely used as a host for therapeutic antibody production since the early 1980s, however, little is known about their immune repertoires or immunological processes affecting the antibody generation. The objective of this study was to employ NGS for a comprehensive analysis of immunoglobulin heavy and light chain repertoires in four healthy sheep. We obtained > 90 % complete antibody sequences and nearly 130,000, 48,000 and 218,000 unique CDR3 reads for the heavy chain (IGH), kappa chain (IGK), and lambda chain (IGL) loci, respectively. Consistent with other species, we observed biased usage of germline variable (V), diversity (D) and joining (J) genes in the heavy and kappa loci, but not in the lambda loci. Moreover, the enormous diversity of CDR3 sequences was observed through sequence clustering and convergent recombination. These data will build a foundation for future studies investigating immune repertoires in health and disease as well as contribute to further refinement of ovine-derived therapeutic antibody drugs.

Development and pharmacokinetic assessment of a fully canine anti-PD-1 monoclonal antibody for comparative translational research in dogs with spontaneous tumors

PD-1 checkpoint inhibitors have revolutionized the treatment of patients with different cancer histologies including melanoma, renal cell carcinoma, and non-small cell lung carcinoma. However, only a subset of patients show a dramatic clinical response to treatment. Despite intense biomarker discovery efforts, no single robust, prognostic correlation has emerged as a valid outcome predictor. Immune competent, pet dogs develop spontaneous tumors that share similar features to human cancers including chromosome aberrations, molecular subtypes, immune signatures, tumor heterogeneity, metastatic behavior, and chemotherapeutic response. As such, they represent a valuable parallel patient population in which to investigate predictive biomarkers of checkpoint inhibition. However, the lack of a validated, non-immunogenic, canine anti-PD-1 antibody for pre-clinical use hinders this comparative approach and prevents potential clinical benefits of PD-1 blockade being realized in the veterinary clinic. To address this, fully canine single-chain variable fragments (scFvs) that bind canine (c)PD-1 were isolated from a comprehensive canine scFv phage display library. Lead candidates were identified that bound with high affinity to cPD-1 and inhibited its interaction with canine PD-L1 (cPD-L1). The lead scFv candidate re-formatted into a fully canine IgGD reversed the inhibitory effects of cPD-1:cPD-L1 interaction on canine chimeric antigen receptor (CAR) T cell function. In vivo administration showed no toxicity and revealed favorable pharmacokinetics for a reasonable dosing schedule. These results pave the way for clinical trials with anti-cPD-1 in canine cancer patients to investigate predictive biomarkers and combination regimens to inform human clinical trials and bring a promising checkpoint inhibitor into the veterinary armamentarium.

Evolution of surrogate light chain in tetrapods and the relationship between lengths of CDR H3 and VpreB tails

In the mammalian immune system, the surrogate light chain (SLC) shapes the antibody repertoire during B cell development by serving as a checkpoint for production of functional heavy chains (HC). Structural studies indicate that tail regions of VpreB contact and cover the third complementarity-determining region of the HC (CDR H3). However, some species, particularly bovines, have CDR H3 regions that may not be compatible with this HC-SLC interaction model. With immense structural and genetic diversity in antibody repertoires across species, we evaluated the genetic origins and sequence features of surrogate light chain components. We examined tetrapod genomes for evidence of conserved gene synteny to determine the evolutionary origin of VpreB1, VpreB2, and IGLL1, as well as VpreB3 and pre-T cell receptor alpha (PTCRA) genes. We found the genes for the SLC components (VpreB1, VpreB2, and IGLL1) only in eutherian mammals. However, genes for PTCRA occurred in all amniote groups and genes for VpreB3 occurred in all tetrapod groups, and these genes were highly conserved. Additionally, we found evidence of a new VpreB gene in non-mammalian tetrapods that is similar to the VpreB2 gene of eutherian mammals, suggesting VpreB2 may have appeared earlier in tetrapod evolution and may be a precursor to traditional VpreB2 genes in higher vertebrates. Among eutherian mammals, sequence conservation between VpreB1 and VpreB2 was low for all groups except rabbits and rodents, where VpreB2 was nearly identical to VpreB1 and did not share conserved synteny with VpreB2 of other species. VpreB2 of rabbits and rodents likely represents a duplicated variant of VpreB1 and is distinct from the VpreB2 of other mammals. Thus, rabbits and rodents have two variants of VpreB1 (VpreB1-1 and VpreB1-2) but no VpreB2. Sequence analysis of VpreB tail regions indicated differences in sequence content, charge, and length; where repertoire data was available, we observed a significant relationship between VpreB2 tail length and maximum DH length. We posit that SLC components co-evolved with immunoglobulin HC to accommodate the repertoire - particularly CDR H3 length and structure, and perhaps highly unusual HC (like ultralong HC of cattle) may bypass this developmental checkpoint altogether.

Development and application of a next-generation sequencing protocol and bioinformatics pipeline for the comprehensive analysis of the canine immunoglobulin repertoire

Profiling the adaptive immune repertoire using next generation sequencing (NGS) has become common in human medicine, showing promise in characterizing clonal expansion of B cell clones through analysis of B cell receptors (BCRs) in patients with lymphoid malignancies. In contrast, most work evaluating BCR repertoires in dogs has employed traditional PCR-based approaches analyzing the IGH locus only. The objectives of this study were to: (1) describe a novel NGS protocol to evaluate canine BCRs; (2) develop a bioinformatics pipeline for processing canine BCR sequencing data; and (3) apply these methods to derive insights into BCR repertoires of healthy dogs and dogs undergoing treatment for B-cell lymphoma. RNA from peripheral blood mononuclear cells of healthy dogs (n = 25) and dogs newly diagnosed with intermediate-to-large B-cell lymphoma (n = 18) with intent to pursue chemotherapy was isolated, converted into cDNA and sequenced by NGS. The BCR repertoires were identified and quantified using a novel analysis pipeline. The IGK repertoires of the healthy dogs were far less diverse compared to IGL which, as with IGH, was highly diverse. Strong biases at key positions within the CDR3 sequence were identified within the healthy dog BCR repertoire. For a subset of the dogs with B-cell lymphoma, clonal expansion of specific IGH sequences pre-treatment and reduction post-treatment was observed. The degree of expansion and reduction correlated with the clinical outcome in this subset. Future studies employing these techniques may improve disease monitoring, provide earlier recognition of disease progression, and ultimately lead to more targeted therapeutics.

Polymorphisms in canine immunoglobulin heavy chain gene cluster: a double-edged sword for diabetes mellitus in the dog

Insulin deficiency diabetes (IDD) in dogs is an endocrine disease similar to human type 1 diabetes. There are breeds more commonly affected, such as Yorkshire Terrier and Samoyed, suggesting an underlying genetic component. However, the genetic basis for canine diabetes mellitus (DM) is not fully established. We conducted both whole-genome scans for selection signatures and GWASs to compare the genomes of 136 dogs belonging to 29 breeds previously described at low or high risk for developing DM. Candidate variants were tested in dogs with a diagnosis of IDD and controls attending the Complutense Veterinary Teaching Hospital. The only genomic region under selection (CFA8:72 700 000-74 600 000; CanFam3.1) retrieved by our analyses is included in the immunoglobulin heavy chain gene cluster, which has already been related to human human type 1 diabetes susceptibility. This region contains two non-synonymous variants, rs852072969 and rs851728071, showing significant associations with high or low risk for IDD, respectively. The first variant, rs852072969, alters a protein poorly characterised in the dog. In contrast, rs851728071 was predicted to block the synthesis of an immunoglobulin variable (V) domain in breeds at low risk for DM. Although a large and diverse V gene repertoire is thought to offer a fitness advantage, we suggest that rs851728071 prevents the formation of an auto-reactive immunoglobulin V domain probably involved in the pathophysiology of IDD and, thus, decreases the risk for the disease. These results should be interpreted with caution until the functional roles of the proposed variants have been proved in larger studies.

Topology and expressed repertoire of the Felis catus T cell receptor loci

Background

The domestic cat (Felis catus) is an important companion animal and is used as a large animal model for human disease. However, the comprehensive study of adaptive immunity in this species is hampered by the lack of data on lymphocyte antigen receptor genes and usage. The objectives of this study were to annotate the feline T cell receptor (TR) loci and to characterize the expressed repertoire in lymphoid organs of normal cats using high-throughput sequencing.

Results

The Felis catus TRG locus contains 30 genes: 12 TRGV, 12 TRGJ and 6 TRGC, the TRB locus contains 48 genes: 33 TRBV, 2 TRBD, 11 TRBJ, 2 TRBC, the TRD locus contains 19 genes: 11 TRDV, 2 TRDD, 5 TRDJ, 1 TRDC, and the TRA locus contains 127 genes: 62 TRAV, 64 TRAJ, 1 TRAC. Functional feline V genes form monophyletic clades with their orthologs, and clustering of multimember subgroups frequently occurs in V genes located at the 5′ end of TR loci. Recombination signal (RS) sequences of the heptamer and nonamer of functional V and J genes are highly conserved. Analysis of the TRG expressed repertoire showed preferential intra-cassette over inter-cassette rearrangements and dominant usage of the TRGV2–1 and TRGJ1–2 genes. The usage of TRBV genes showed minor bias but TRBJ genes of the second J-C-cluster were more commonly rearranged than TRBJ genes of the first cluster. The TRA/TRD V genes almost exclusively rearranged to J genes within their locus. The TRAV/TRAJ gene usage was relatively balanced while the TRD repertoire was dominated by TRDJ3.

Conclusions

This is the first description of all TR loci in the cat. The genomic organization of feline TR loci was similar to that of previously described jawed vertebrates (gnathostomata) and is compatible with the birth-and-death model of evolution. The large-scale characterization of feline TR genes provides comprehensive baseline data on immune repertoires in healthy cats and will facilitate the development of improved reagents for the diagnosis of lymphoproliferative diseases in cats. In addition, these data might benefit studies using cats as a large animal model for human disease.

Clonality testing in the lymph nodes from dogs with lymphadenomegaly due to Leishmania infantum infection

INTRODUCTION

In southern European countries, multicentric lymphoma and leishmaniosis are the main differential diagnoses in dogs presented with generalized lymphadenomegaly. The cytological examination is in some cases inconclusive and polymerase chain reaction (PCR) for antigen receptor rearrangement (PARR) has become a common method to confirm or rule out a lymphoproliferative neoplasia. According to the literature, leishmaniosis may lead to clonal arrangements and therefore to a false diagnosis of lymphoma, but this assumption is made from a single leishmania infected dog. Therefore, the objective of this study was to prospectively evaluate results from PARR in dogs with lymphadenomegaly due to clinical leishmaniosis at the moment of diagnosis.

MATERIALS AND METHODS

31 dogs with a diagnosis of leishmaniosis based on the LeishVet guidelines were included in the study. Samples from enlarged lymph nodes were taken for cytological examination, clonality testing and Leishmania infantum PCR.

RESULTS

All 31 dogs had medium to high positive antibody titers against Leishmania spp. and 30/31 had a positive Leishmania PCR from the lymph node. A polyclonal arrangement for B cells (immunoglobulin heavy chain gene) and T cells (T-cell receptor gamma chain gene) antigen receptors was found in 28/31 dogs. Two out of 31 dogs showed a monoclonal arrangement for Ig with high (1:2) and low (1:7) polyclonal background respectively; and one of the 31 dogs showed a monoclonal arrangement for T cell receptor with low (1:3) polyclonal background.

CONCLUSION

Infections with Leishmania infantum resulted in clonal rearrangement, and therefore in a possible false diagnosis of lymphoma, in 3 out of 31 dogs (9.7%). Although, PARR is a useful method to differentiate lymphoma from reactive lymphoid hyperplasia in dogs with leishmaniosis, mono-/biclonal results should be interpreted carefully, especially in the presence of any degree of polyclonal background, and together with other clinicopathological findings.

A review of canine B cell clonality assays and primer set optimization using large-scale repertoire data

Several molecular clonality assays have been developed to assess canine B cell proliferations. These assays were based on different sequence data, utilized different assay designs and employed different testing strategies. This has resulted in a complex body of literature and complicates evidence-based selection of primer sets. In addition, further refinement of primer sets is difficult because it is unknown how well current primer sets cover the expressed sequence repertoire. The objectives of this study were 1) to provide an overview of published IGH clonality assays that highlights key differences in assay design and testing strategy and 2) to propose a novel method for optimizing primer sets that leverages large-scale sequencing data. A review of previously published assays highlighted confounding factors that hamper a direct comparison of performance metrics between studies. These findings illustrate the need for a multi-institutional effort to harmonize veterinary clonality testing. A novel in silico analysis of primer sequences using a large dataset of expressed sequences identified shortfalls of existing primer sets and was used to guide primer optimization. Three optimized primer sets were tested and yielded qualitative sensitivity values between 80-90%. The qualitative sensitivity ranged from 1% to over 50% and was dependent on the size of the neoplastic clone and the sample DNA used. These findings illustrate that inclusion of high-throughput sequencing data for primer design can be a useful tool to guide primer design and optimization. This strategy could be applied to other antigen receptor loci or species to further improve veterinary clonality assays.