Identification of three new bovine T-cell receptor delta variable gene subgroups expressed by peripheral blood T cells

Single-cell RNA sequencing reveals the local cell landscape in mouse epididymal initial segment during aging

BACKGROUND

Morphological and functional alterations in aging reproductive organs result in decreased male fertility. The epididymis functions as the transition region for post-testicular sperm maturation. And we have previously demonstrated that the epididymal initial segment (IS), a region of the reproductive tract essential for sperm maturation and capacitation, undergoes considerable histological changes and chronic immune activation in mice during aging. However, the local aging-associated cellular and molecular changes in the aged epididymal IS are poorly understood.

RESULTS

We conducted single-cell RNA sequencing analysis on the epididymal IS of young (3-month-old) and old (21-month-old) mice. In total, 10,027 cells from the epididymal IS tissues of young and old mice were obtained and annotated. The cell composition, including the expansion of a principal cell subtype and Ms4a4b^HiMs4a6b^Hi T cells, changed with age. Aged principal cells displayed multiple functional gene expression changes associated with acrosome reaction and sperm maturation, suggesting an asynchronous process of sperm activation and maturation during epididymal transit. Meanwhile, aging-related altered pathways in immune cells, especially the "cell chemotaxis" in Cx3cr1^Hi epididymal dendritic cells (eDCs), were identified. The monocyte-specific expression of chemokine Ccl8 increased with age in eDCs. And the aged epididymal IS showed increased inflammatory cell infiltration and cytokine secretion. Furthermore, cell-cell communication analysis indicated that age increased inflammatory signaling in the epididymal IS.

CONCLUSION

Contrary to the general pattern of lower immune responses in the male proximal genital tract, we revealed an inflammaging status in mouse epididymal initial segment. These findings will allow future studies to enable the delay of male reproductive aging via immune regulation.

γδ T cells in artiodactyls: Focus on swine

Vaccination is the most effective medical strategy for disease prevention but there is a need to improve livestock vaccine efficacy. Understanding the structure of the immune system of swine, which are considered a γδ T cell "high" species, and thus, particularly how to engage their γδ T cells for immune responses, may allow for development of vaccine optimization strategies. The propensity of γδ T cells to home to specific tissues, secrete pro-inflammatory and regulatory cytokines, exhibit memory or recall responses and even function as antigen-presenting cells for αβ T cells supports the concept that they have enormous potential for priming by next generation vaccine constructs to contribute to protective immunity. γδ T cells exhibit several innate-like antigen recognition properties including the ability to recognize antigen in the absence of presentation via major histocompatibility complex (MHC) molecules enabling γδ T cells to recognize an array of peptides but also non-peptide antigens in a T cell receptor-dependent manner. γδ T cell subpopulations in ruminants and swine can be distinguished based on differential expression of the hybrid co-receptor and pattern recognition receptors (PRR) known as workshop cluster 1 (WC1). Expression of various PRR and other innate-like immune receptors diversifies the antigen recognition potential of γδ T cells. Finally, γδ T cells in livestock are potent producers of critical master regulator cytokines such as interferon (IFN)-γ and interleukin (IL)-17, whose production orchestrates downstream cytokine and chemokine production by other cells, thereby shaping the immune response as a whole. Our knowledge of the biology, receptor expression and response to infectious diseases by swine γδ T cells is reviewed here.

Special features of γδ T cells in ruminants

Ruminant γδ T cells were discovered in the mid-1980's shortly after a novel T cell receptor (TCR) gene from murine cells was described in 1984 and the murine TCRγ gene locus in 1985. It was possible to identify γδ T cell populations early in ruminants because they represent a large proportion of the peripheral blood mononuclear cells (PBMC). This null cell population, γδ T cells, was designated as such by its non-reactivity with monoclonal antibodies (mAb) against ovine and bovine CD4, CD8 and surface immunoglobulin (Ig). γδ T cells are non-conventional T cells known as innate-like cells capable of using both TCR as well as other types of receptor systems including pattern recognition receptors (PRR) and natural killer receptors (NKR). Bovine γδ T cells have been shown to respond to stimulation through toll-like receptors, NOD, and NKG2D as well as to cytokines alone, protein and non-protein antigens through their TCR, and to pathogen-infected host cells. The two main populations of γδ T cells are distinguished by the presence or absence of the hybrid co-receptor/PRR known as WC1 or T19. These two populations not only differ by their proportional representation in various tissues and organs but also by their migration into inflamed tissues. The WC1⁺ cells are found in the blood, skin and spleen while the WC1^- γδ T cells predominate in the gut, mammary gland and uterus. In ruminants, γδ T cells may produce IFNγ, IL-17, IL-10 and TGFβ, have cytotoxic activity and memory responses. The expression of particular WC1 family members controls the response to particular pathogens and correlates with differences in cytokine responses. The comparison of the WC1 gene families in cattle, sheep and goats is discussed relative to other multigenic arrays that differentiate γδ T cells by function in humans and mice.

Analysis of the Chinese Alligator TCRα/δ Loci Reveals the Evolutionary Pattern of Atypical TCRδ/TCRμ in Tetrapods

Atypical TCRδ found in sharks, amphibians, birds, and monotremes and TCRμ found in monotremes and marsupials are TCR chains that use Ig or BCR-like variable domains (VHδ/Vμ) rather than conventional TCR V domains. These unconventional TCR are consistent with a scenario in which TCR and BCR, although having diverged from each other more than 400 million years ago, continue to exchange variable gene segments in generating diversity for Ag recognition. However, the process underlying this exchange and leading to the evolution of these atypical TCR receptor genes remains elusive. In this study, we identified two TCRα/δ gene loci in the Chinese alligator (Alligator sinensis). In total, there were 144 V, 154 Jα, nine Jδ, eight Dδ, two Cα, and five Cδ gene segments in the TCRα/δ loci of the Chinese alligator, representing the most complicated TCRα/δ gene system in both genomic structure and gene content in any tetrapod examined so far. A pool of 32 VHδ genes divided into 18 subfamilies was found to be scattered over the two loci. Phylogenetic analyses revealed that these VHδ genes could be related to bird VHδ genes, VHδ/Vμ genes in platypus or opossum, or alligator VH genes. Based on these findings, a model explaining the evolutionary pattern of atypical TCRδ/TCRμ genes in tetrapods is proposed. This study sheds new light on the evolution of TCR and BCR genes, two of the most essential components of adaptive immunity.

The bovine model for elucidating the role of γδ T cells in controlling infectious diseases of importance to cattle and humans

There are several instances of co-investigation and related discoveries and achievements in bovine and human immunology; perhaps most interesting is the development of the BCG vaccine, the tuberculin skin test and the more recent interferon-gamma test that were developed first in cattle to prevent and diagnosis bovine tuberculosis and then applied to humans. There are also a number of immune-physiological traits that ruminant share with humans including the development of their immune systems in utero which increases the utility of cattle as a model for human immunology. These are reviewed here with a particular focus on the use of cattle to unravel γδ T cell biology. Based on the sheer number of γδ T cells in this γδ T cell high species, it is reasonable to expect γδ T cells to play an important role in protective immune responses. For that reason alone cattle may provide good models for elucidating at least some of the roles γδ T cells play in protective immunity in all species. This includes fundamental research on γδ T cells as well as the responses of ruminant γδ T cells to a variety of infectious disease situations including to protozoan and bacterial pathogens. The role that pattern recognition receptors (PRR) play in the activation of γδ T cells may be unique relative to αβ T cells. Here we focus on that of the γδ T cell specific family of molecules known as WC1 or T19 in ruminants, which are part of the CD163 scavenger receptor cysteine rich (SRCR) family that includes SCART1 and SCART2 expressed on murine γδ T cells. We review the evidence for WC1 being a PRR as well as an activating co-receptor and the role that γδ T cells bearing these receptors play in immunity to leptospirosis and tuberculosis. This includes the generation of memory responses to vaccines, thereby continuing the tradition of co-discovery between cattle and humans.

Spectratype analysis of the T cell receptor δ CDR3 region of bovine γδ T cells responding to leptospira

Gamma delta T cells comprise the majority of blood T cells in ruminants at birth and remain at high levels for several years with most expressing the WC1 co-receptor. A subpopulation of Bos taurus WC1(+) cells expressing a restricted set of WC1 molecules respond immediately by proliferation and interferon-γ production to leptospira following vaccination, preceding the response by CD4 T cells. Our goal is to define the γδ T cell recognition elements involved. Previously, we showed that the responding cells employed a variety of TRDV genes indicating that the CDR1 and CDR2 of TCRδ could vary and may not be principally involved in antigen specificity. Murine and human γδ T cells bind T22 and self lipids through their CDR3δ. Like mice, cattle use up to five TRDD genes in a single CDR3δ adding flexibility to length and configuration for antigen binding. Here, we used spectratyping to evaluate the CDR3δ of leptospira-responsive cells. Little or no compartmentalization of CDR3δ was found for antigen-responsive cells that incorporated TRDV1, TRDV2, or TRDV3 even though they comprise the majority of the leptospira-responding population. Compartmentalization occurred for TRDV4-containing transcripts and was maintained over time and among cattle. However, no common amino acid motif was apparent in those CDR3δ sequences, although a bias in D gene usage occurred. We hypothesize that the restricted set of WC1 co-receptors expressed by the responding cells may lend specificity to the response through their ability to bind bacteria facilitating interaction of various TCRs with bacterial components resulting in cross-linking and activation.

Genomic analysis offers insights into the evolution of the bovine TRA/TRD locus

Background

The TRA/TRD locus contains the genes for V(D)J somatic rearrangement of TRA and TRD chains expressed by αβ and γδ T cells respectively. Previous studies have demonstrated that the bovine TRA/TRD locus contains an exceptionally large number of TRAV/TRDV genes. In this study we combine genomic and transcript analysis to provide insights into the evolutionary development of the bovine TRA/TRD locus and the remarkable TRAV/TRDV gene repertoire.

Results

Annotation of the UMD3.1 assembly identified 371 TRAV/TRDV genes (distributed in 42 subgroups), 3 TRDJ, 6 TRDD, 62 TRAJ and single TRAC and TRDC genes, most of which were located within a 3.5 Mb region of chromosome 10. Most of the TRAV/TRDV subgroups have multiple members and several have undergone dramatic expansion, most notably TRDV1 (60 genes). Wide variation in the proportion of pseudogenes within individual subgroups, suggest that differential ‘birth’ and ‘death’ rates have been used to form a functional bovine TRAV/TRDV repertoire which is phylogenetically distinct from that of humans and mice. The expansion of the bovine TRAV/TRDV gene repertoire has predominantly been achieved through a complex series of homology unit (regions of DNA containing multiple gene) replications. Frequent co-localisation within homology units of genes from subgroups with low and high pseudogene proportions suggest that replication of homology units driven by evolutionary selection for the former may have led to a ‘collateral’ expansion of the latter. Transcript analysis was used to define the TRAV/TRDV subgroups available for recombination of TRA and TRD chains and demonstrated preferential usage of different subgroups by the expressed TRA and TRD repertoires, indicating that TRA and TRD selection have had distinct impacts on the evolution of the TRAV/TRDV repertoire.

Conclusion

Both TRA and TRD selection have contributed to the evolution of the bovine TRAV/TRDV repertoire. However, our data suggest that due to homology unit duplication TRD selection for TRDV1 subgroup expansion may have substantially contributed to the genomic expansion of several TRAV subgroups. Such data demonstrate how integration of genomic and transcript data can provide a more nuanced appreciation of the evolutionary dynamics that have led to the dramatically expanded bovine TRAV/TRDV repertoire.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-994) contains supplementary material, which is available to authorized users.

Characteristics of the somatic hypermutation in the Camelus dromedarius T cell receptor gamma (TRG) and delta (TRD) variable domains

In previous reports, we had shown in Camelus dromedarius that diversity in T cell receptor gamma (TRG) and delta (TRD) variable domains can be generated by somatic hypermutation (SHM). In the present paper, we further the previous finding by analyzing 85 unique spleen cDNA sequences encoding a total of 331 mutations from a single animal, and comparing the properties of the mutation profiles of dromedary TRG and TRD variable domains. The transition preference and the significant mutation frequency in the AID motifs (dgyw/wrch and wa/tw) demonstrate a strong dependence of the enzymes mediating SHM in TRG and TRD genes of dromedary similar to that of immunoglobulin genes in mammals. Overall, results reveal no asymmetry in the motifs targeting, i.e. mutations are equally distributed among g:c and a:t base pairs and replacement mutations are favored at the AID motifs, whereas neutral mutations appear to be more prone to accumulate in bases outside of the motifs. A detailed analysis of clonal lineages in TRG and TRD cDNA sequences also suggests that clonal expansion of mutated productive rearrangements may be crucial in shaping the somatic diversification in the dromedary. This is confirmed by the fact that our structural models, computed by adopting a comparative procedure, are consistent with the possibility that, irrespective of where (in the CDR-IMGT or in FR-IMGT) the diversity was generated by mutations, both clonal expansion and selection seem to be strictly related to an enhanced structural stability of the γδ subunits.

Bovine γδ T cells: cells with multiple functions and important roles in immunity

The γδ T-cell receptor (TCR)-positive lymphocytes are a major circulating lymphocyte population in cattle, especially in young calves. In contrast, human and mice have low levels of circulating γδ TCR(+) T cells (γδ T cells). The majority of the circulating γδ T cells in ruminants express the workshop cluster 1 (WC1) molecule and are of the phenotype WC1(+) CD2(-) CD4(-) CD8(-). WC1 is a 220000 molecular weight glycoprotein with homology to the scavenger receptor cysteine-rich (SRCR) family, closely related to CD163. The existence of 13 members in the bovine WC1 gene family has recently been demonstrated and although murine and human orthologues to WC1 genes exist, functional gene products have not been identified in species other than ruminants and pigs. Highly diverse TCRδ usage has been reported, with expanded variable genes in cattle compared to humans and mice. Differential γ chain usage is evident between populations of bovine γδ T cells, this may have implications for functionality. There is a growing body of evidence that WC1(+) γδ T cells are important in immune responses to mycobacteria and may have important roles in T cell regulation and antigen presentation. In this review, we will summarize recent observations in γδ T cell biology and the importance of γδ T cells in immune responses to mycobacterial infections in cattle.

Annotation and classification of the bovine T cell receptor delta genes

Background

γδ T cells differ from αβ T cells with regard to the types of antigen with which their T cell receptors interact; γδ T cell antigens are not necessarily peptides nor are they presented on MHC. Cattle are considered a "γδ T cell high" species indicating they have an increased proportion of γδ T cells in circulation relative to that in "γδ T cell low" species such as humans and mice. Prior to the onset of the studies described here, there was limited information regarding the genes that code for the T cell receptor delta chains of this γδ T cell high species.

Results

By annotating the bovine (Bos taurus) genome Btau_3.1 assembly the presence of 56 distinct T cell receptor delta (TRD) variable (V) genes were found, 52 of which belong to the TRDV1 subgroup and were co-mingled with the T cell receptor alpha variable (TRAV) genes. In addition, two genes belonging to the TRDV2 subgroup and single TRDV3 and TRDV4 genes were found. We confirmed the presence of five diversity (D) genes, three junctional (J) genes and a single constant (C) gene and describe the organization of the TRD locus. The TRDV4 gene is found downstream of the C gene and in an inverted orientation of transcription, consistent with its orthologs in humans and mice. cDNA evidence was assessed to validate expression of the variable genes and showed that one to five D genes could be incorporated into a single transcript. Finally, we grouped the bovine and ovine TRDV1 genes into sets based on their relatedness.

Conclusions

The bovine genome contains a large and diverse repertoire of TRD genes when compared to the genomes of "γδ T cell low" species. This suggests that in cattle γδ T cells play a more important role in immune function since they would be predicted to bind a greater variety of antigens.

The bovine T cell receptor alpha/delta locus contains over 400 V genes and encodes V genes without CDR2

αβ T cells and γδ T cells perform nonoverlapping immune functions. In mammalian species with a high percentage of very diverse γδ T cells, like ruminants and pigs, it is often assumed that αβ T cells are less diverse than γδ T cells. Based on the bovine genome, we have created a map of the bovine TRA/TRD locus and show that, in cattle, in addition to the anticipated >100 TRDV genes, there are also >300 TRAV or TRAV/DV genes. Among the V genes in the TRA/TRD locus, there are several genes that lack a CDR2 and are functionally rearranged and transcribed and, in some cases, have an extended CDR1. The number of bovine V genes is a multiple of the number in mice and humans and may encode T cell receptors that use a novel way of interacting with antigen.

Highly diverse TCR delta chain repertoire in bovine tissues due to the use of up to four D segments per delta chain

Tissue-specific distribution of gammadelta TCRs with limited TCR diversity is a common phenomenon in species with a low percentage of gammadelta T cells like humans and mice. We set out to investigate whether this is also the case in cattle (Bos taurus), a species with high percentages of gammadelta T cells. Using a method that was independent of variable (V) segment-specific primers, we generated 65 unique TCR delta chain sequences. We found no evidence for preferential use of certain Vdelta segments in lymph node, skin, spleen, small intestine, large intestine, and blood. The delta chain CDR3 length distribution was very wide in each tissue, which was confirmed by spectratyping. The highly variable CDR3 length was due to the use of up to four diversity (D) segments by one bovine delta chain. Human and murine delta chains contain only one or two D segments. The five functional Ddelta segments that we describe here were identified at cDNA and genomic level, and are the first ruminant D segments described. Fourteen TCR delta chain sequences used novel Vdelta1 segments, and one expressed a novel member of the Vdelta3 family. The number of known functional Vdelta segments in cattle including these new ones is 42 now, but the total number may be much higher. A high number of Vdelta segments in combination with the use of up to four out of five D segments, and the possibility of using non-template encoded (N) nucleotides on either side of these, makes the potential bovine delta chain repertoire much bigger than any known TCR chain. This situation is quite different from the situation in humans and mice, and suggests that the differences between gammadelta high and gammadelta low species in distribution, diversity, and function of gammadelta T cells may be substantial.

Comparison of gene expression by co-cultured WC1+ gammadelta and CD4+ alphabeta T cells exhibiting a recall response to bacterial antigen

Immunization of cattle with a Leptospira borgpetersenii serovar hardjo-bovis vaccine results in the development of a recall response by WC1(+) gammadelta T cells and CD4(+) alphabeta T cells characterized by proliferation and interferon-gamma production. It was hypothesized that these two T cell subpopulations had largely redundant effector functions, principally differing in their requirements for activation. To test this, gene expression in cells proliferating to antigen were compared utilizing RT-PCR and bovine microarrays. Both T cell populations had similar transcript profiles for effector molecules, including IFN-gamma, FasL and granzyme B. In contrast, transcripts for costimulatory receptors and ligands were notably different following activation, as WC1(+) T cells expressed no or lower levels of transcripts for CD28 and CD40L, while CD4(+) T cells expressed substantial levels of both. However, both cell types had high levels of CTLA-4 transcript suggesting the cells may be regulated similarly following activation but differ in their need for and ability to provide costimulation. Microarray analyses to extend the number of genes examined revealed that while both subpopulations upregulated anti-apoptotic genes as well as those involved in cell activation and protein biosynthesis, overall there were limited differences between the two antigen-activated cell populations. Those genes that did differ were involved in cell signaling, protein production and intracellular protein trafficking. These results strengthen the hypothesis that these particular activated WC1(+) and CD4(+) T cells have overlapping effector functions and therefore may differ principally with regard to how they are recruited into immune responses.

Differential TCR gene usage between WC1 − and WC1 + ruminant γδ T cell subpopulations including those responding to bacterial antigen

Ruminant gammadelta T cells are divided into subpopulations based on the presence or absence of WC1 co-receptors (scavenger-receptor-cysteine-rich family members uniquely expressed on gammadelta T cells). Evidence suggests WC1+ are inflammatory while WC1- are regulatory and that they also differ in their tissue distribution. Recently, this paradigm was refined further as cells that produce interferon-gamma and proliferate to autologous antigens, leptospira antigens, or IL-12 were largely found within the WC1+ subpopulation that bears the WC1.1 antigenic epitope but not that bearing the WC1.2 epitope. Here, the T cell receptor gene expression by these different subpopulations (WC1-, WC1.1+, and WC1.2+) was compared using flow cytometrically-purified cells and reverse transcriptase-polymerase chain reaction (RT-PCR). The WC1- gammadelta T cells had transcripts for all 11 possible combinations of the TRG subgroup V and C genes while those in both WC1+ subpopulations were restricted to TRGV3-TRGC5 and TRGV7-TRGC5. In contrast, all three subpopulations expressed transcripts from all four known bovine TRDV genes. Further analysis of the WC1+ gammadelta T cells that proliferated in leptospira antigen-stimulated cultures indicated that they do not represent a unique subpopulation within the larger WC1+ population based on their TCR gene usage. Moreover, sequencing of 65 transcripts showed that their junctional regions were diverse as TRGJ5-1, TRGJ5-2, TRDJ1, and TRDJ3 were used, and CDR3s ranged from 9 to 24 amino acids. The restricted but shared gammadelta TCR gene usage for WC1.1+, WC1.2+, and WC1(+)-antigen-responsive cells leaves open the possibility that the WC1 co-receptor is an important determining element in the activation process and subsequent response.