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

Local chicken breeds exhibit abundant TCR-V segments but similar repertoire diversity

The thymus-derived lymphocytes of jawed vertebrates have four T-cell receptor (TCR) chains that play a significant role in immunity. As chickens have commercial value, their immune systems require a great deal of attention. Local chicken breeds are an essential part of poultry genetic resources in China. Here, we used high-throughput sequencing to analyze the TCRα and TCRβ repertoires and their relative expression levels in the native chicken breeds Baier Buff, Longyou Partridge, Xiaoshan, and Xianju. We found that TCR Vα and TCR Vβ were expressed and included 17, 19, 17, and six segments of the Vα2, Vα3, Vβ1, and Vβ2 subgroups, respectively. V-J pairing was biased; Jα11 was utilized by nearly all Vα segments and was the most commonly used. Breed-specific V segments and V-J pairings were detected as well. The results of the principal coordinate analysis (PCoA) as well as the V-J pairing and CDR3 diversity analyses suggested that the four local chicken breeds did not significantly differ in terms of TCR diversity. Hence, they expressed not significant differentiation, and they are rich genetic resources for the development and utilization of immune-related poultry breeding.

Integrated mRNA-miRNA transcriptome profiling of blood immune responses potentially related to pulmonary fibrosis in forest musk deer

Background

Forest musk deer (FMD, Moschus Berezovskii) is a critically endangered species world-widely, the death of which can be caused by pulmonary disease in the farm. Pulmonary fibrosis (PF) was a huge threat to the health and survival of captive FMD. MicroRNAs (miRNAs) and messenger RNAs (mRNAs) have been involved in the regulation of immune genes and disease development. However, the regulatory profiles of mRNAs and miRNAs involved in immune regulation of FMD are unclear.

Methods

In this study, mRNA-seq and miRNA-seq in blood were performed to constructed coexpression regulatory networks between PF and healthy groups of FMD. The hub immune- and apoptosis-related genes in the PF blood of FMD were explored through Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Further, protein-protein interaction (PPI) network of immune-associated and apoptosis-associated key signaling pathways were constructed based on mRNA-miRNA in the PF blood of the FMD. Immune hub DEGs and immune hub DEmiRNAs were selected for experimental verification using RT-qPCR.

Results

A total of 2744 differentially expressed genes (DEGs) and 356 differentially expressed miRNAs (DEmiRNAs) were identified in the PF blood group compared to the healthy blood group. Among them, 42 DEmiRNAs were negatively correlated with 20 immune DEGs from a total of 57 correlations. The DEGs were significantly associated with pathways related to CD molecules, immune disease, immune system, cytokine receptors, T cell receptor signaling pathway, Th1 and Th2 cell differentiation, cytokine-cytokine receptor interaction, intestinal immune network for IgA production, and NOD-like receptor signaling pathway. There were 240 immune-related DEGs, in which 186 immune-related DEGs were up-regulated and 54 immune-related DEGs were down-regulated. In the protein-protein interaction (PPI) analysis of immune-related signaling pathway, TYK2, TLR2, TLR4, IL18, CSF1, CXCL13, LCK, ITGB2, PIK3CB, HCK, CD40, CD86, CCL3, CCR7, IL2RA, TLR3, and IL4R were identified as the hub immune genes. The mRNA-miRNA coregulation analysis showed that let-7d, miR-324-3p, miR-760, miR-185, miR-149, miR-149-5p, and miR-1842-5p are key miRNAs that target DEGs involved in immune disease, immune system and immunoregulation.

Conclusion

The development and occurrence of PF were significantly influenced by the immune-related and apoptosis-related genes present in PF blood. mRNAs and miRNAs associated with the development and occurrence of PF in the FMD.

Subgenomic T cell receptor alpha and delta (TRA/TRD) loci in common carp

In jawed vertebrates, the T cell receptor alpha (TRA) and delta (TRD) genes, which encode the TRα and TRδ chains, respectively, are located as a nested structure on a single chromosome. To date, no animal has been reported to harbor multiple TRA/TRD loci on different chromosomes. Therefore, herein, we describe the first full annotation of the TRA/TRD genomic regions of common carp, an allo-tetraploid fish species that experiences cyprinid-specific whole-genome duplication (WGD) in evolution. Fine genomic maps of TRA/TRD genomic regions 1 and 2, on LG30 and LG22, respectively, were constructed using the annotations of complete sets of TRA and TRD genes, including TRA/TRD variable (V), TRA junction (J), and constant (C), TRD diversity (D), and the J and C genes. The structure and synteny of the TRA/TRD genomic regions were highly conserved in zebrafish, indicating that these regions are on individual chromosomes. Furthermore, analysis of the variable regions of the TRA and TRD genes in a monoclonal T cell line revealed that both subgenomic regions 1 and 2 were indeed rearranged. Although carp TRAV and TRDV genes were phylogenetically divided into different lineages, they were mixed and organized into the TRA/TRD V gene clusters on the genome, similar to that in other vertebrates. Notably, 285 potential TRA/TRD V genes were detected in the TRA/TRD genomic regions, which is the most abundant number of genes in vertebrates and approximately two-fold that in zebrafish. The recombination signal sequences (RSSs) at the end of each V gene differed between TRAV and TRDV, suggesting that RSS variations might separate each V gene into a TRα or TRδ chain. This study is the first to describe subgenomic TRA/TRD loci in animals. Our findings provide fundamental insights to elucidate the impact of WGD on the evolution of immune repertoire.

New insights into the germline genes and CDR3 repertoire of the TCRβ chain in Chiroptera

Introduction

Bats are recognized as natural reservoirs for many viruses, and their unique immune system enables them to coexist with these viruses without frequently exhibiting disease symptoms. However, the current understanding of the bat adaptive immune system is limited due to the lack of a database or tool capable of processing T-cell receptor (TCR) sequences for bats.

Methods

We performed germline gene annotation in three bat species using homologous genes and RSSs (Recombinational Signal Sequences) scanning method. Then we used the conserved C gene to construct the TCRβ chain receptor library of the Intermediate Horseshoe Bat. Bats' TCRβ data will be analyzed using MiXCR and constructed reference library.

Results

Regarding the annotation results, we found that the Pale Spear-nosed Bat has 37 members in the TRBV12 family, which is more than the total number of TRBV genes in the Greater Horseshoe Bat. The average number of unique TCRβ chain receptor sequences in each Intermediate Horseshoe Bat sample reached 24,904.

Discussion

The distinct variations in the distribution of TRBV genes among the three types of bats could have a direct impact on the diversity of the TCR repertoire, as evidenced by the presence of conserved amino acids that indicate the T-cell recognition of antigens in bats is MHC-restricted. The bats’ TCRβ repertoire is formed through the rearrangement of the V-D-J-C genes, with D-J/V-D deletions and insertions resulting in high diversity.

Blood transcriptome analysis revealed the immune changes and immunological adaptation of wildness training giant pandas

The giant panda (Ailuropoda melanoleuca) is a global flagship species for biodiversity conservation. As the time for captive giant pandas to be released into the wild matures, wildness training is provided to allow adaptation to their natural environment. It is assumed that changes in the immune system would be integral in this adaptation from captive to wild, where many more pathogens would be encountered in their natural habitats. Therefore, this study aims to determine the expression changes of immune-related genes and their potential as immunoassay markers for adaptation monitoring in wildness training giant pandas, and then to understand the adaptation strategy of wildness training giant pandas to the wild environment, thereby improving the success rate of panda reintroduction. We obtained 300 differentially expressed genes (DEGs) by RNA-seq, with 239 up-regulated and 61 down-regulated DEGs in wildness training giant pandas compared to captive pandas. Functional enrichment analysis indicated that up-regulated DEGs were enriched in several immune-related terms and pathways. There were 21 immune-related DEGs, in which most of them were up-regulated in wildness training giant pandas, including several critical innate and cellular immune genes. IL1R2 was the most significantly up-regulated gene and is a signature of homeostasis within the immune system. In the protein-protein interaction (PPI) analysis, CXCL8, CXCL10, and CCL5 were identified as the hub immune genes. Our results suggested that wildness training giant pandas have stronger innate and cellular immunity than captive giant pandas, and we proposed that a gene set of CXCL8, CXCL10, CCL5, CD3D, NFKBIA, TBX21, IL12RB2, and IL1R2 may serve as potential immunoassay markers to monitor and assess the immune status of wildness training giant pandas. Our study offers the first insight into immune alterations of wildness training giant pandas, paving the way for monitoring and evaluating the immune status of giant pandas when reintroducing them into the wild.

TR Locus Annotation and Characteristics of Rhinolophus ferrumequinum

T-cell antigen receptors (TRs) in vertebrates can be divided into αβ or γδ, encoded by TRA/D, TRG, or TRB loci. TRs play a central role in mammal cellular immunity, which occurs by rearrangement of V, D, J, and C genes in the loci. The bat is the only mammal with flying ability and is considered the main host of zoonotic viruses, an important public health concern. However, at present, little is known about the composition of bat TR genes. Based on the whole genome sequence of the greater horseshoe bat (Rhinolophus ferrumequinum) and referring to the TR/IG annotation rules formulated by the international ImMunoGeneTics information system (IMGT), we present a complete annotation of TRA/D, TRG, and TRB loci of R. ferrumequinum. A total of 128 V segments, three D segments, 85 J segments, and 6 C segments were annotated and compared with other known mammalian data. The characteristics of the TR locus and germline genes of R. ferrumequinum are analyzed.

Identification and Phenotype of MAIT Cells in Cattle and Their Response to Bacterial Infections

Mucosal-associated invariant T (MAIT) cells are a population of innate-like T cells that utilize a semi-invariant T cell receptor (TCR) α chain and are restricted by the highly conserved antigen presenting molecule MR1. MR1 presents microbial riboflavin biosynthesis derived metabolites produced by bacteria and fungi. Consistent with their ability to sense ligands derived from bacterial sources, MAIT cells have been associated with the immune response to a variety of bacterial infections, such as Mycobacterium spp., Salmonella spp. and Escherichia coli. To date, MAIT cells have been studied in humans, non-human primates and mice. However, they have only been putatively identified in cattle by PCR based methods; no phenotypic or functional analyses have been performed. Here, we identified a MAIT cell population in cattle utilizing MR1 tetramers and high-throughput TCR sequencing. Phenotypic analysis of cattle MAIT cells revealed features highly analogous to those of MAIT cells in humans and mice, including expression of an orthologous TRAV1-TRAJ33 TCR α chain, an effector memory phenotype irrespective of tissue localization, and expression of the transcription factors PLZF and EOMES. We determined the frequency of MAIT cells in peripheral blood and multiple tissues, finding that cattle MAIT cells are enriched in mucosal tissues as well as in the mesenteric lymph node. Cattle MAIT cells were responsive to stimulation by 5-OP-RU and riboflavin biosynthesis competent bacteria in vitro. Furthermore, MAIT cells in milk increased in frequency in cows with mastitis. Following challenge with virulent Mycobacterium bovis, a causative agent of bovine tuberculosis and a zoonosis, peripheral blood MAIT cells expressed higher levels of perforin. Thus, MAIT cells are implicated in the immune response to two major bacterial infections in cattle. These data suggest that MAIT cells are functionally highly conserved and that cattle are an excellent large animal model to study the role of MAIT cells in important zoonotic infections.

Characterization of the domestic goat γδ T cell receptor gene loci and gene usage

Goats and cattle diverged 30 million years ago but retain similarities in immune system genes. Here, the caprine T cell receptor (TCR) gene loci and transcription of its genes were examined and compared to cattle. We annotated the TCR loci using an improved genome assembly (ARS1) of a highly homozygous San Clemente goat. This assembly has already proven useful for describing other immune system genes including antibody and leucocyte receptors. Both the TCRγ (TRG) and TCRδ (TRD) loci were similarly organized in goats as in cattle and the gene sequences were highly conserved. However, the number of genes varied slightly as a result of duplications and differences occurred in mutations resulting in pseudogenes. WC1⁺ γδ T cells in cattle have been shown to use TCRγ genes from only one of the six available cassettes. The structure of that Cγ gene product is unique and may be necessary to interact with WC1 for signal transduction following antigen ligation. Using RT-PCR and PacBio sequencing, we observed the same restriction for goat WC1⁺ γδ T cells. In contrast, caprine WC1⁺ and WC1^- γδ T cell populations had a diverse TCRδ gene usage although the propensity for particular gene usage differed between the two cell populations. Noncanonical recombination signal sequences (RSS) largely correlated with restricted expression of TCRγ and δ genes. Finally, caprine γδ T cells were found to incorporate multiple TRD diversity gene sequences in a single transcript, an unusual feature among mammals but also previously observed in cattle.

IMGT® Biocuration and Comparative Analysis of Bos taurus and Ovis aries TRA/TRD Loci

The adaptive immune response provides the vertebrate immune system with the ability to recognize and remember specific pathogens to generate immunity, and mount stronger attacks each time the pathogen is encountered. T cell receptors are the antigen receptors of the adaptive immune response expressed by T cells, which specifically recognize processed antigens, presented as peptides by the highly polymorphic major histocompatibility (MH) proteins. T cell receptors (TR) are divided into two groups, αβ and γδ, which express distinct TR containing either α and β, or γ and δ chains, respectively. The TRα locus (TRA) and TRδ locus (TRD) of bovine (Bos taurus) and the sheep (Ovis aries) have recently been described and annotated by IMGT^® biocurators. The aim of the present study is to present the results of the biocuration and to compare the genes of the TRA/TRD loci among these ruminant species based on the Homo sapiens repertoire. The comparative analysis shows similarities but also differences, including the fact that these two species have a TRA/TRD locus about three times larger than that of humans and therefore have many more genes which may demonstrate duplications and/or deletions during evolution.

The expansion of the TRB and TRG genes in domestic goats (Capra hircus) is characteristic of the ruminant species

BACKGROUND

Goats (Capra hircus), one of the first domesticated species, are economically important for milk and meat production, and their broad geographical distribution reflects their successful adaptation to diverse environmental conditions. Despite the relevance of this species, the genetic research on the goat traits is limited compared to other domestic species. Thanks to the latest goat reference genomic sequence (ARS1), which is considered to be one of the most continuous assemblies in livestock, we deduced the genomic structure of the T cell receptor beta (TRB) and gamma (TRG) loci in this ruminant species.

RESULTS

Our analyses revealed that although the organization of the goat TRB locus is broadly similar to that of the other artiodactyl species, with three in-tandem D-J-C clusters located at the 3' end, a complex and extensive series of duplications have occurred in the V genes at the 5' end, leading to a marked expansion in the number of the TRBV genes. This phenomenon appears to be a feature of the ruminant lineage since similar gene expansions have also occurred in sheep and cattle. Likewise, the general organization of the goat TRG genes is typical of ruminant species studied so far, with two paralogous TRG loci, TRG1 and TRG2, located in two distinct and distant positions on the same chromosome as result of a split in the ancestral locus. Each TRG locus consists of reiterated V-J-J-C cassettes, with the goat TRG2 containing an additional cassette relative to the corresponding sheep and cattle loci.

CONCLUSIONS

Taken together, these findings demonstrate that strong evolutionary pressures in the ruminant lineage have selected for the development of enlarged sets of TRB and TRG genes that contribute to a diverse T cell receptor repertoire. However, differences observed among the goat, sheep and cattle TRB and TRG genes indicate that distinct evolutionary histories, with independent expansions and/or contractions, have also affected each ruminant species.

Chromosome-level assembly of wild Bactrian camel genome reveals organization of immune gene loci

Camelids are characterized by their unique adaptive immune system that exhibits the generation of homodimeric heavy-chain immunoglobulins, somatic hypermutation of T-cell receptors, and low genetic diversity of major histocompatibility complex (MHC) genes. However, short-read assemblies are typically highly fragmented in these gene loci owing to their repetitive and polymorphic nature. Here, we constructed a chromosome-level assembly of wild Bactrian camel genome based on high-coverage long-read sequencing and chromatin interaction mapping. The assembly with a contig N50 of 5.37 Mb and a scaffold N50 of 76.03 Mb, represents the most contiguous camelid genome to date. The genomic organization of immunoglobulin heavy-chain locus was similar between the wild Bactrian camel and alpaca, and genes encoding for conventional and heavy-chain antibodies were intermixed. The organizations of two immunoglobulin light-chain loci and four T cell receptor loci were also fully deciphered using the new assembly. Additionally, the complete classical MHC region was resolved into a single contig. The high-quality assembly presented here provides an essential reference for future investigations examining the camelid immune system.

The T cell receptor (TRA) locus in the rabbit (Oryctolagus cuniculus): Genomic features and consequences for invariant T cells

The rabbit has been widely used in immunology and infectiology. Rabbit immunoglobulins have been extensively studied, leading to the discovery of their idiotypes, allotypic diversity, and of the diversification of the primary repertoire by hyperconversion. Much less is known about rabbit T cell receptors (TR), especially TRA. This isotype is particularly important for innate-like T cells, which typically express invariant TRA (iTRA). The presence of such cells in the rabbit remains an enigma. Rabbit NKT cells seem to be very rare, and lagomorphs lack MAIT cells. TRAV1, the variable gene expressed in the iTRA of these cells across most mammals, and MR1, the MH1-like receptor that present riboflavin derivatives to MAIT cells, are missing in rabbit. An alternative iTRA has been identified, that may be expressed by new innate-like T cells. To facilitate TRA repertoire analyses in rabbit, we report here a full description of TRA and TRD loci and a subgroup definition based on IMGT® classification. Rabbit TRA rearrangements follow the same temporal pattern that is observed in mouse and human. Rare transcripts expressing TRDV/TRDD/TRDJ rearrangements spliced to TRAC were detected. TRA and TRD genes have been made available in IMGT and IMGT/HighV-QUEST, allowing easy analysis of TRA/TRD RepSeq.

Divergent Evolution of TRC Genes in Mammalian Niche Adaptation

Mammals inhabit a wide variety of ecological niches, which in turn can be affected by various ecological factors, especially in relation to immunity. The canonical TRC repertoire (TRAC, TRBC, TRGC, and TRDC) codes C regions of T cell receptor chains that form the primary antigen receptors involved in the activation of cellular immunity. At present, little is known about the correlation between the evolution of mammalian TRC genes and ecological factors. In this study, four types canonical of TRC genes were identified from 37 mammalian species. Phylogenetic comparative methods (phyANOVA and PGLS) and selective pressure analyses among different groups of ecological factors (habitat, diet, and sociality) were carried out. The results showed that habitat was the major ecological factor shaping mammalian TRC repertoires. Specifically, trade-off between TRGC numbers and positive selection of TRAC and the balanced evolutionary rates between TRAC and TRDC genes were speculated as two main mechanisms in adaption to habitat and sociality. Overall, our study suggested divergent mechanisms for the evolution of TRCs, prompting mammalian immunity adaptions within diverse niches.

Computational detection and experimental validation of segmental duplications and associated copy number variations in water buffalo ( Bubalus bubalis )

Duplicated sequences are an important source of gene evolution and structural variation within mammalian genomes. Using a read depth approach based on next-generation sequencing, we performed a genome-wide analysis of segmental duplications (SDs) and associated copy number variations (CNVs) in the water buffalo (Bubalus bubalis). By aligning short reads of Olimpia (the reference water buffalo) to the UMD3.1 cattle genome, we identified 1,038 segmental duplications comprising 44.6 Mb (equivalent to ~1.73% of the cattle genome) of the autosomal and X chromosomal sequence in the buffalo genome. We experimentally validated 70.3% (71/101) of these duplications using fluorescent in situ hybridization. We also detected a total of 1,344 CNV regions across 14 additional water buffaloes, amounting to 59.8 Mb of variable sequence or the equivalent of 2.2% of the cattle genome. The CNV regions overlap 1,245 genes that are significantly enriched for specific biological functions including immune response, oxygen transport, sensory system and signal transduction. Additionally, we performed array Comparative Genomic Hybridization (aCGH) experiments using the 14 water buffaloes as test samples and Olimpia as the reference. Using a linear regression model, a high Pearson correlation (r = 0.781) was observed between the log₂ ratios between copy number estimates and the log₂ ratios of aCGH probes. We further designed Quantitative PCR assays to confirm CNV regions within or near annotated genes and found 74.2% agreement with our CNV predictions. These results confirm sub-chromosome-scale structural rearrangements present in the cattle and water buffalo. The information on genome variation that will be of value for evolutionary and phenotypic studies, and may be useful for selective breeding of both species.

Bovine T cell receptors and γδ WC1 co-receptor transcriptome analysis during the first month of life

Here we evaluated neonatal transcription of α, β, γ and δ TCR and the γδ T cell co-receptor family WC1 in peripheral blood mononuclear cells. A previous report showed a rapid and global shift in transcription of immunoglobulin genes in neonatal calves during the first month after birth but this was not found here for the T cell genes. Transcription frequency of genes within TRAV subgroups correlated with the number of members, indicating a stochastic choice. In contrast, of the approximately 60 TRDV genes those in two of eleven TRDV1 clades and TRDVb3 were transcribed significantly more than the others while those in only one TRBV subgroup were. Transcription of genes in the TRGV5-containing cassette predominated among TRGV genes as a result of their exclusive usage by the WC1⁺ γδ T cells with a preference for transcription of two of four TRGV genes in that cassette. Finally, we report no large differences in transcription frequencies among the 13 WC1 genes.

The occurrence of three D-J-C clusters within the dromedary TRB locus highlights a shared evolution in Tylopoda, Ruminantia and Suina

The αβ T cells are important components of the adaptive immune system and can recognize a vast array of peptides presented by MHC molecules. The ability of these T cells to recognize the complex depends on the diversity of the αβ TR, which is generated by a recombination of specific Variable, Diversity and Joining genes for the β chain, and Variable and Joining genes for the α chain. In this study, we analysed the genomic structure and the gene content of the TRB locus in Camelus dromedarius, which is a species belonging to the Tylopoda suborder. The most noteworthy result is the presence of three in tandem TRBD-J-C clusters in the dromedary TRB locus, which is similar to clusters found in sheep, cattle and pigs and suggests a common duplication event occurred prior to the Tylopoda/Ruminantia/Suina divergence. Conversely, a significant contraction of the dromedary TRBV genes, which was previously found in the TRG and TRD loci, was observed with respect to the other artiodactyl species.

Genomic and expression analyses of Tursiops truncatus T cell receptor gamma (TRG) and alpha/delta (TRA/TRD) loci reveal a similar basic public γδ repertoire in dolphin and human

BACKGROUND

The bottlenose dolphin (Tursiops truncatus) is a mammal that belongs to the Cetartiodactyla and have lived in marine ecosystems for nearly 60 millions years. Despite its popularity, our knowledge about its adaptive immunity and evolution is very limited. Furthermore, nothing is known about the genomics and evolution of dolphin antigen receptor immunity.

RESULTS

Here we report a evolutionary and expression study of Tursiops truncatus T cell receptor gamma (TRG) and alpha/delta (TRA/TRD) genes. We have identified in silico the TRG and TRA/TRD genes and analyzed the relevant mature transcripts in blood and in skin from four subjects. The dolphin TRG locus is the smallest and simplest of all mammalian loci as yet studied. It shows a genomic organization comprising two variable (V1 and V2), three joining (J1, J2 and J3) and a single constant (C), genes. Despite the fragmented nature of the genome assemblies, we deduced the TRA/TRD locus organization, with the recent TRDV1 subgroup genes duplications, as it is expected in artiodactyls. Expression analysis from blood of a subject allowed us to assign unambiguously eight TRAV genes to those annotated in the genomic sequence and to twelve new genes, belonging to five different subgroups. All transcripts were productive and no relevant biases towards TRAV-J rearrangements are observed. Blood and skin from four unrelated subjects expression data provide evidence for an unusual ratio of productive/unproductive transcripts which arise from the TRG V-J gene rearrangement and for a "public" gamma delta TR repertoire. The productive cDNA sequences, shared both in the same and in different individuals, include biases of the TRGV1 and TRGJ2 genes. The high frequency of TRGV1-J2/TRDV1- D1-J4 productive rearrangements in dolphins may represent an interesting oligo-clonal population comparable to that found in human with the TRGV9- JP/TRDV2-D-J T cells and in primates.

CONCLUSIONS

Although the features of the TRG and TRA/TRD loci organization reflect those of the so far examined artiodactyls, genomic results highlight in dolphin an unusually simple TRG locus. The cDNA analysis reveal productive TRA/TRD transcripts and unusual ratios of productive/unproductive TRG transcripts. Comparing multiple different individuals, evidence is found for a "public" gamma delta TCR repertoire thus suggesting that in dolphins as in human the gamma delta TCR repertoire is accompanied by selection for public gamma chain.

Genomic organization of the zebrafish (Danio rerio) T cell receptor alpha/delta locus and analysis of expressed products

In testing the hypothesis that all jawed vertebrate classes employ immunoglobulin heavy chain V (IgHV) gene segments in their T cell receptor (TCR)δ encoding loci, we found that some basic characterization was required of zebrafish TCRδ. We began by annotating and characterizing the TCRα/δ locus of Danio rerio based on the most recent genome assembly, GRCz10. We identified a total of 141 theoretically functional V segments which we grouped into 41 families based upon 70 % nucleotide identity. This number represents the second greatest count of apparently functional V genes thus far described in an antigen receptor locus with the exception of cattle TCRα/δ. Cloning, relative quantitative PCR, and deep sequencing results corroborate that zebrafish do express TCRδ, but these data suggest only at extremely low levels and in limited diversity in the spleens of the adult fish. While we found no evidence for IgH-TCRδ rearrangements in this fish, by determining the locus organization we were able to suggest how the evolution of the teleost α/δ locus could have lost IgHVs that exist in sharks and frogs. We also found evidence of surprisingly low TCRδ expression and repertoire diversity in this species.

An overview of the lagomorph immune system and its genetic diversity

Our knowledge of the lagomorph immune system remains largely based upon studies of the European rabbit (Oryctolagus cuniculus), a major model for studies of immunology. Two important and devastating viral diseases, rabbit hemorrhagic disease and myxomatosis, are affecting European rabbit populations. In this context, we discuss the genetic diversity of the European rabbit immune system and extend to available information about other lagomorphs. Regarding innate immunity, we review the most recent advances in identifying interleukins, chemokines and chemokine receptors, Toll-like receptors, antiviral proteins (RIG-I and Trim5), and the genes encoding fucosyltransferases that are utilized by rabbit hemorrhagic disease virus as a portal for invading host respiratory and gut epithelial cells. Evolutionary studies showed that several genes of innate immunity are evolving by strong natural selection. Studies of the leporid CCR5 gene revealed a very dramatic change unique in mammals at the second extracellular loop of CCR5 resulting from a gene conversion event with the paralogous CCR2. For the adaptive immune system, we review genetic diversity at the loci encoding antibody variable and constant regions, the major histocompatibility complex (RLA) and T cells. Studies of IGHV and IGKC genes expressed in leporids are two of the few examples of trans-species polymorphism observed outside of the major histocompatibility complex. In addition, we review some endogenous viruses of lagomorph genomes, the importance of the European rabbit as a model for human disease studies, and the anticipated role of next-generation sequencing in extending knowledge of lagomorph immune systems and their evolution.

Sheep (Ovis aries) T cell receptor alpha (TRA) and delta (TRD) genes and genomic organization of the TRA/TRD locus

Background

In mammals, T cells develop along two discrete pathways characterized by expression of either the αβ or the γδ T cell receptors. Human and mouse display a low peripheral blood γδ T cell percentage ("γδ low species") while sheep, bovine and pig accounts for a high proportion of γδ T lymphocytes ("γδ high species"). While the T cell receptor alpha (TRA) and delta (TRD) genes and the genomic organization of the TRA/TRD locus has been determined in human and mouse, this information is still poorly known in artiodactyl species, such as sheep.

Results

The analysis of the current Ovis aries whole genome assembly, Oar_v3.1, revealed that, as in the other mammalian species, the sheep TRD locus is nested within the TRA locus. In the most 5’ part the TRA/TRD locus contains TRAV genes which are intermingled with TRDV genes, then TRD genes which include seven TRDD, four TRDJ genes, one TRDC and a single TRDV gene with an inverted transcriptional orientation, and finally in the most 3’ part, the TRA locus is completed by 61 TRAJ genes and one TRAC gene.

Comparative sequence and analysis and annotation led to the identification of 66 TRAV genes assigned to 34 TRAV subgroups and 25 TRDV genes belonging to the TRDV1 subgroup, while one gene was found for each TRDV2, TRDV3 and TRDV4 subgroups. Multiple duplication events within several TRAV subgroups have generated the sheep TRAV germline repertoire, which is substantially larger than the human one. A significant proportion of these TRAV gene duplications seems to have occurred simultaneously with the amplification of the TRDV1 subgroup genes. This dynamic of expansion has also generated novel multigene subgroups, which are species-specific. Ovis aries TRA and TRD genes identified in this study were assigned IMGT definitive or temporary names and were approved by the IMGT/WHO-IUIS nomenclature committee.

The completeness of the genome assembly in the 3' part of the locus has allowed us to interpret rearranged CDR3 of cDNA from both TRA and TRD chain repertoires. The involvement of one up to four TRDD genes into a single transcript makes the potential sheep TRD chain much larger than any known TR chain repertoire.

Conclusions

The sheep genome, as the bovine genome, contains a large and diverse repertoire of TRA and TRD genes when compared to the “γδ T cell low” species genomes. The composition and length of the rearranged CDR3 in TRD V-delta domains influence the three-dimensional configuration of the antigen-combining site thus suggesting that in ruminants, γδ T cells play a more important and specific role in immune recognition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1790-z) contains supplementary material, which is available to authorized users.