Extensive analysis of D-J-C arrangements allows the identification of different mechanisms enhancing the diversity in sheep T cell receptor beta-chain repertoire

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut

Alternative splicing (AS), an important post-transcriptional regulation mechanism in eukaryotes, can significantly increase transcript diversity and contribute to gene expression regulation and many other complicated developmental processes. While plant gene AS events are well described, few studies have investigated the comprehensive regulation machinery of plant AS. Here, we use multi-omics to analyse peanut AS events. Using long-read isoform sequencing, 146 464 full-length non-chimeric transcripts were obtained, resulting in annotation corrections for 1782 genes and the identification of 4653 new loci. Using Iso-Seq RNA sequences, 271 776 unique splice junctions were identified, 82.49% of which were supported by transcriptome data. We characterized 50 977 polyadenylation sites for 23 262 genes, 12 369 of which had alternative polyadenylation sites. AS allows differential regulation of the same gene by miRNAs at the isoform level coupled with polyadenylation. In addition, we identified many long non-coding RNAs and fusion transcripts. There is a suppressed effect of 6mA on AS and gene expression. By analysis of chromatin structures, the genes located in the boundaries of topologically associated domains, proximal chromosomal telomere regions, inter- or intra-chromosomal loops were found to have more unique splice isoforms, higher expression, lower 6mA and more transposable elements (TEs) in their gene bodies than the other genes, indicating that chromatin interaction, 6mA and TEs play important roles in AS and gene expression. These results greatly refine the peanut genome annotation and contribute to the study of gene expression and regulation in peanuts. This work also showed AS is associated with multiple strategies for gene regulation.

The T Cell Receptor (TRB) Locus in Tursiops truncatus: From Sequence to Structure of the Alpha/Beta Heterodimer in the Human/Dolphin Comparison

The bottlenose dolphin (Tursiops truncatus) belongs to the Cetartiodactyla and, similarly to other cetaceans, represents the most successful mammalian colonization of the aquatic environment. Here we report a genomic, evolutionary, and expression study of T. truncatus T cell receptor beta (TRB) genes. Although the organization of the dolphin TRB locus is similar to that of the other artiodactyl species, with three in tandem D-J-C clusters located at its 3' end, its uniqueness is given by the reduction of the total length due essentially to the absence of duplications and to the deletions that have drastically reduced the number of the germline TRBV genes. We have analyzed the relevant mature transcripts from two subjects. The simultaneous availability of rearranged T cell receptor α (TRA) and TRB cDNA from the peripheral blood of one of the two specimens, and the human/dolphin amino acids multi-sequence alignments, allowed us to calculate the most likely interactions at the protein interface between the alpha/beta heterodimer in complex with major histocompatibility class I (MH1) protein. Interacting amino acids located in the complementarity-determining region according to IMGT numbering (CDR-IMGT) of the dolphin variable V-alpha and beta domains were identified. According to comparative modelization, the atom pair contact sites analysis between the human MH1 grove (G) domains and the T cell receptor (TR) V domains confirms conservation of the structure of the dolphin TR/pMH.

The Genomic Organisation of the TRA/TRD Locus Validates the Peculiar Characteristics of Dromedary δ-Chain Expression

The role of γδ T cells in vertebrate immunity is still an unsolved puzzle. Species such as humans and mice display a low percentage of these T lymphocytes (i.e., “γδ low species”) with a restricted diversity of γδ T cell receptors (TR). Conversely, artiodactyl species (i.e., “γδ high species”) account for a high proportion of γδ T cells with large γ and δ chain repertoires. The genomic organisation of the TR γ (TRG) and δ (TRD) loci has been determined in sheep and cattle, noting that a wide number of germline genes that encode for γ and δ chains characterise their genomes. Taking advantage of the current improved version of the genome assembly, we have investigated the genomic structure and gene content of the dromedary TRD locus, which, as in the other mammalian species, is nested within the TR α (TRA) genes. The most remarkable finding was the identification of a very limited number of variable germline genes (TRDV) compared to sheep and cattle, which supports our previous expression analyses for which the somatic hypermutation mechanism is able to enlarge and diversify the primary repertoire of dromedary δ chains. Furthermore, the comparison between genomic and expressed sequences reveals that D genes, up to four incorporated in a transcript, greatly contribute to the increased diversity of the dromedary δ chain antigen binding-site.

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.

IMGT® Biocuration and Comparative Study of the T Cell Receptor Beta Locus of Veterinary Species Based on Homo sapiens TRB

IMGT®, the international ImMunoGeneTics information system® is the global reference in immunogenetics and immunoinformatics. By its creation in 1989 by Marie-Paule Lefranc (Université de Montpellier and CNRS), IMGT® marked the advent of immunoinformatics, which emerged at the interface between immunogenetics and bioinformatics. IMGT® is specialized in the immunoglobulins (IG) or antibodies, T cell receptors (TR), major histocompatibility (MH), and proteins of the IgSF and MhSF superfamilies. T cell receptors are divided into two groups, αβ and γδ TR, which express distinct TR containing either α and β, or γ and δ chains, respectively. The TRβ locus (TRB) was recently described and annotated by IMGT® biocurators for several veterinary species, i.e., cat (Felis catus), dog (Canis lupus familiaris), ferret (Mustela putorius furo), pig (Sus scrofa), rabbit (Oryctolagus cuniculus), rhesus monkey (Macaca mulatta), and sheep (Ovis aries). The aim of the present study is to compare the genes of the TRB locus among these different veterinary species based on Homo sapiens. The results reveal that there are similarities but also differences including the number of genes by subgroup which may demonstrate duplications and/or deletions during evolution.

The Camel Adaptive Immune Receptors Repertoire as a Singular Example of Structural and Functional Genomics

The adaptive immune receptors repertoire is highly plastic, with its ability to produce antigen-binding molecules and select those with high affinity for their antigen. Species have developed diverse genetic and structural strategies to create their respective repertoires required for their survival in the different environments. Camelids, until now, considered as a case of evolutionary innovation because of their only heavy-chain antibodies, represent a new mammalian model particularly useful for understanding the role of diversity in the immune system function. Here, we review the structural and functional characteristics and the current status of the genomic organization of camel immunoglobulins (IG) or antibodies, α/ß and γ/δ T cell receptors (TR), and major histocompatibility complex (MHC). In camelid humoral response, in addition to the conventional antibodies, there are IG with “only-heavy-chain” (no light chain, and two identical heavy gamma chains lacking CH1 and with a VH domain designated as VHH). The unique features of these VHH offer advantages in biotechnology and for clinical applications. The TRG and TRD rearranged variable domains of Camelus dromedarius (Arabian camel) display somatic hypermutation (SHM), increasing the intrinsic structural stability in the γ/δ heterodimer and influencing the affinity maturation to a given antigen similar to immunoglobulin genes. The SHM increases the dromedary γ/δ repertoire diversity. In Camelus genus, the general structural organization of the TRB locus is similar to that of the other artiodactyl species, with a pool of TRBV genes positioned at the 5’ end of three in tandem D-J-C clusters, followed by a single TRBV gene with an inverted transcriptional orientation located at the 3’ end. At the difference of TRG and TRD, the diversity of the TRB variable domains is not shaped by SHM and depends from the classical combinatorial and junctional diversity. The MHC locus is located on chromosome 20 in Camelus dromedarius. Cytogenetic and comparative whole genome analyses revealed the order of the three major regions “Centromere-ClassII-ClassIII-ClassI”. Unexpectedly low extent of polymorphisms and haplotypes was observed in all Old World camels despite different geographic origins.

Overview of the Germline and Expressed Repertoires of the TRB Genes in Sus scrofa

The α/β T cell receptor (TR) is a complex heterodimer that recognizes antigenic peptides and binds to major histocompatibility complex (MH) molecules. Both α and β chains are encoded by different genes localized on two distinct chromosomal loci: TRA and TRB. The present study employed the recent release of the swine genome assembly to define the genomic organization of the TRB locus. According to the sequencing data, the pig TRB locus spans approximately 400 kb of genomic DNA and consists of 38 TRBV genes belonging to 24 subgroups located upstream of three in tandem TRBD-J-C clusters, which are followed by a TRBV gene in an inverted transcriptional orientation. Comparative analysis confirms that the general organization of the TRB locus is similar among mammalian species, but the number of germline TRBV genes varies greatly even between species belonging to the same order, determining the diversity and specificity of the immune response. However, sequence analysis of the TRB locus also suggests the presence of blocks of conserved homology in the genomic region across mammals. Furthermore, by analysing a public cDNA collection, we identified the usage pattern of the TRBV, TRBD, and TRBJ genes in the adult pig TRB repertoire, and we noted that the expressed TRBV repertoire seems to be broader and more diverse than the germline repertoire, in line with the presence of a high level of TRBV gene polymorphisms. Because the nucleotide differences seems to be principally concentrated in the CDR2 region, it is reasonable to presume that most T cell β-chain diversity can be related to polymorphisms in pig MH molecules. Domestic pigs represent a valuable animal model as they are even more anatomically, genetically and physiologically similar to humans than are mice. Therefore, present knowledge on the genomic organization of the pig TRB locus allows the collection of increased information on the basic aspects of the porcine immune system and contributes to filling the gaps left by rodent models.

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.

Data characterizing the genomic structure of the T cell receptor (TRB) locus in Camelus dromedarius

These data are presented in support of structural and evolutionary analysis of the published article entitled “The occurrence of three D-J-C clusters within the dromedary TRB locus highlights a shared evolution in Tylopoda, Ruminantia and Suina” (Antonacci et al., 2017) [1]. Here we describe the genomic structure and the gene content of the T cell receptor beta chain (TRB) locus in Camelus dromedarius. As in the other species of mammals, the general genomic organization of the dromedary TRB locus consists of a pool of TRBV genes located upstream of in tandem TRBD-J-C clusters, followed by a TRBV gene with an inverted transcriptional orientation. A peculiarity of the dromedary TRB locus structure is the presence of three TRBD-J-C clusters, which is a common feature of sheep, cattle and pig sequences.

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.