Genomic V exons from whole genome shotgun data in reptiles

Adaptive Immunity in Reptiles: Conventional Components but Unconventional Strategies

Recent studies have established that the innate immune system of reptiles is broad and robust, but the question remains: What role does the reptilian adaptive immune system play? Conventionally, adaptive immunity is described as involving T and B lymphocytes that display variable receptors, is highly specific, improves over the course of the response, and produces a memory response. While reptiles do have B and T lymphocytes that utilize variable receptors, their adaptive response is relatively non-specific, generates a prolonged antibody response, and does not produce a typical memory response. This alternative adaptive strategy may allow reptiles to produce a broad adaptive response that complements a strong innate system. Further studies into reptile adaptive immunity cannot only clarify outstanding questions on the reptilian immune system but can shed light on a number of important immunological concepts, including the evolution of the immune system and adaptive immune responses that take place outside of germinal centers.

Comparison of Reptilian Genomes Reveals Deletions Associated with the Natural Loss of γδ T Cells in Squamates

T lymphocytes or T cells are key components of the vertebrate response to pathogens and cancer. There are two T cell classes based on their TCRs, αβ T cells and γδ T cells, and each plays a critical role in immune responses. The squamate reptiles may be unique among the vertebrate lineages by lacking an entire class of T cells, the γδ T cells. In this study, we investigated the basis of the loss of the γδ T cells in squamates. The genome and transcriptome of a sleepy lizard, the skink Tiliqua rugosa, were compared with those of tuatara, Sphenodon punctatus, the last living member of the Rhynchocephalian reptiles. We demonstrate that the lack of TCRγ and TCRδ transcripts in the skink are due to large deletions in the T. rugosa genome. We also show that tuataras are on a growing list of species, including sharks, frogs, birds, alligators, and platypus, that can use an atypical TCRδ that appears to be a chimera of a TCR chain with an Ab-like Ag-binding domain. Tuatara represents the nearest living relative to squamates that retain γδ T cells. The loss of γδTCR in the skink is due to genomic deletions that appear to be conserved in other squamates. The genes encoding the αβTCR chains in the skink do not appear to have increased in complexity to compensate for the loss of γδ T cells.

Gouania willdenowi is a teleost fish without immunoglobulin genes

Immunoglobulin (Ig) genes encode antibodies in jawed vertebrates. They are essential elements of the adaptive immune response. Ig exists in soluble form or as part of the B cell membrane antigen receptor (BCR). Studies of Ig genes in fish genomes reveal the absence of Ig genes in Gouania willdenowi by deletion of the entire Ig locus from the canonical chromosomal region. The genes coding for integral BCR proteins, CD79a and CD79b, are also absent. Genes exist for T α/β lymphocyte receptors but not for the T γ/δ receptors. The results of the genomic analysis are independently corroborated with RNA-Seq transcriptomes from other Gobiesocidae species. From the transcriptome studies, Ig is also absent from these other Gobiesocidae species, Acyrtus sp. and Tomicodon sp. Present evidence suggests that Ig is missing from all species of the Gobiesocidae family.

Insights into the evolution of IG genes in Amphibians and reptiles

Immunoglobulins are essential proteins of the immune system to neutralize pathogens. Gene encoding B cell receptors and antibodies (Ig genes) first appeared with the emergence of early vertebrates having a jaw, and are now present in all extant jawed vertebrates, or Gnathostomata. The genes have undergone evolutionary changes. In particular, genomic structural changes corresponding to genes of the adaptive immune system were coincident or in parallel with the adaptation of vertebrates from the sea to land. In cartilaginous fish exist IgM, IgD/W, and IgNAR and in bony fish IgM, IgT, IgD. Amphibians and reptiles witnessed significant modifications both in the structure and orientation of IG genes. In particular, for these amphibians and Amniota that adapted to land, IgM and IgD genes were retained, but other isotypes arose, including genes for IgA(X)1, IgA(X)2, and IgY. Recent progress in high throughput genome sequencing is helping to uncover the IG gene structure of all jawed vertebrates. In this work, we review the work and present knowledge of immunoglobulin genes in genomes of amphibians and reptiles.

The reptilian perspective on vertebrate immunity: 10 years of progress

Ten years ago, ‘Understanding the vertebrate immune system: insights from the reptilian perspective’ was published. At the time, our understanding of the reptilian immune system lagged behind that of birds, mammals, fish and amphibians. Since then, great progress has been made in elucidating the mechanisms of reptilian immunity. Here, I review recent discoveries associated with the recognition of pathogens, effector mechanisms and memory responses in reptiles. Moreover, I put forward key questions to drive the next 10 years of research, including how reptiles are able to balance robust innate mechanisms with avoiding self-damage, how B cells and antibodies are used in immune defense and whether innate mechanisms can display the hallmarks of memory. Finally, I briefly discuss the links between our mechanistic understanding of the reptilian immune system and the field of eco-immunology. Overall, the field of reptile immunology is poised to contribute greatly to our understanding of vertebrate immunity in the next 10 years.

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.

Inter- and intraspecies comparison of phylogenetic fingerprints and sequence diversity of immunoglobulin variable genes

Protection and neutralization of a vast array of pathogens is accomplished by the tremendous diversity of the B cell receptor (BCR) repertoire. For jawed vertebrates, this diversity is initiated via the somatic recombination of immunoglobulin (Ig) germline elements. While it is clear that the number of these germline segments differs from species to species, the extent of cross-species sequence diversity remains largely uncharacterized. Here we use extensive computational and statistical methods to investigate the sequence diversity and evolutionary relationship between Ig variable (V), diversity (D), and joining (J) germline segments across nine commonly studied species ranging from zebrafish to human. Metrics such as guanine-cytosine (GC) content showed low redundancy across Ig germline genes within a given species. Other comparisons, including amino acid motifs, evolutionary selection, and sequence diversity, revealed species-specific properties. Additionally, we showed that the germline-encoded diversity differs across antibody (recombined V-D-J) repertoires of various B cell subsets. To facilitate future comparative immunogenomics analysis, we created VDJgermlines, an R package that contains the germline sequences from multiple species. Our study informs strategies for the humanization and engineering of therapeutic antibodies.

Iterative Variable Gene Discovery from Whole Genome Sequencing with a Bootstrapped Multiresolution Algorithm

In jawed vertebrates, variable (V) genes code for antigen-binding regions of B and T lymphocyte receptors, which generate a specific response to foreign pathogens. Obtaining the detailed repertoire of these genes across the jawed vertebrate kingdom would help to understand their evolution and function. However, annotations of V-genes are known for only a few model species since their extraction is not amenable to standard gene finding algorithms. Also, the more distant evolution of a taxon is from such model species, and there is less homology between their V-gene sequences. Here, we present an iterative supervised machine learning algorithm that begins by training a small set of known and verified V-gene sequences. The algorithm successively discovers homologous unaligned V-exons from a larger set of whole genome shotgun (WGS) datasets from many taxa. Upon each iteration, newly uncovered V-genes are added to the training set for the next predictions. This iterative learning/discovery process terminates when the number of new sequences discovered is negligible. This process is akin to “online” or reinforcement learning and is proven to be useful for discovering homologous V-genes from successively more distant taxa from the original set. Results are demonstrated for 14 primate WGS datasets and validated against Ensembl annotations. This algorithm is implemented in the Python programming language and is freely available at http://vgenerepertoire.org.

Tracing Antibody Repertoire Evolution by Systems Phylogeny

Antibody evolution studies have been traditionally limited to either tracing a single clonal lineage (B cells derived from a single V-(D)-J recombination) over time or examining bulk functionality changes (e.g., tracing serum polyclonal antibody proteins). Studying a single B cell disregards the majority of the humoral immune response, whereas bulk functional studies lack the necessary resolution to analyze the co-existing clonal diversity. Recent advances in high-throughput sequencing (HTS) technologies and bioinformatics have made it possible to examine multiple co-evolving antibody monoclonal lineages within the context of a single repertoire. A plethora of accompanying methods and tools have been introduced in hopes of better understanding how pathogen presence dictates the global evolution of the antibody repertoire. Here, we provide a comprehensive summary of the tremendous progress of this newly emerging field of systems phylogeny of antibody responses. We present an overview encompassing the historical developments of repertoire phylogenetics, state-of-the-art tools, and an outlook on the future directions of this fast-advancing and promising field.

A comprehensive analysis of the germline and expressed TCR repertoire in White Peking duck

Recently, many immune-related genes have been extensively studied in ducks, but relatively little is known about their TCR genes. Here, we determined the germline and expressed repertoire of TCR genes in White Peking duck. The genomic organization of the duck TCRα/δ, TCRγ and unconventional TCRδ2 loci are highly conserved with their counterparts in mammals or chickens. By contrast, the duck TCRβ locus is organized in an unusual pattern, (Vβ)_n-Dβ-(Jβ)₂-Cβ1-(Jβ)₄-Cβ2, which differs from the tandem-aligned clusters in mammals or the translocon organization in some teleosts. Excluding the first exon encoding the immunoglobulin domain, the subsequent exons of the two Cβ show significant diversity in nucleotide sequence and exon structure. Based on the nucleotide sequence identity, 49 Vα, 30 Vδ, 13 Vβ and 15 Vγ unique gene segments are classified into 3 Vα, 5 Vδ, 4 Vβ and 6 Vγ subgroups, respectively. Phylogenetic analyses revealed that most duck V subgroups, excluding Vβ1, Vγ5 and Vγ6, have closely related orthologues in chicken. The coding joints of all cDNA clones demonstrate conserved mechanisms that are used to increase junctional diversity. Collectively, these data provide insight into the evolution of TCRs in vertebrates and improve our understanding of the avian immune system.

A Comprehensive Analysis of the Phylogeny, Genomic Organization and Expression of Immunoglobulin Light Chain Genes in Alligator sinensis, an Endangered Reptile Species

Crocodilians are evolutionarily distinct reptiles that are distantly related to lizards and are thought to be the closest relatives of birds. Compared with birds and mammals, few studies have investigated the Ig light chain of crocodilians. Here, employing an Alligator sinensis genomic bacterial artificial chromosome (BAC) library and available genome data, we characterized the genomic organization of the Alligator sinensis IgL gene loci. The Alligator sinensis has two IgL isotypes, λ and κ, the same as Anolis carolinensis. The Igλ locus contains 6 Cλ genes, each preceded by a Jλ gene, and 86 potentially functional Vλ genes upstream of (Jλ-Cλ)n. The Igκ locus contains a single Cκ gene, 6 Jκs and 62 functional Vκs. All VL genes are classified into a total of 31 families: 19 Vλ families and 12 Vκ families. Based on an analysis of the chromosomal location of the light chain genes among mammals, birds, lizards and frogs, the data further confirm that there are two IgL isotypes in the Alligator sinensis: Igλ and Igκ. By analyzing the cloned Igλ/κ cDNA, we identified a biased usage pattern of V families in the expressed Vλ and Vκ. An analysis of the junctions of the recombined VJ revealed the presence of N and P nucleotides in both expressed λ and κ sequences. Phylogenetic analysis of the V genes revealed V families shared by mammals, birds, reptiles and Xenopus, suggesting that these conserved V families are orthologous and have been retained during the evolution of IgL. Our data suggest that the Alligator sinensis IgL gene repertoire is highly diverse and complex and provide insight into immunoglobulin gene evolution in vertebrates.

Coevolution of T-cell receptors with MHC and non-MHC ligands

The structure and amino acid diversity of the T-cell receptor (TCR), similar in nature to that of Fab portions of antibodies, would suggest that these proteins have a nearly infinite capacity to recognize antigen. Yet all currently defined native T cells expressing an α and β chain in their TCR can only sense antigen when presented in the context of a major histocompatibility complex (MHC) molecule. This MHC molecule can be one of many that exist in vertebrates, presenting small peptide fragments, lipid molecules, or small molecule metabolites. Here we review the pattern of TCR recognition of MHC molecules throughout a broad sampling of species and T-cell lineages and also touch upon T cells that do not appear to require MHC presentation for their surveillance function. We review the diversity of MHC molecules and information on the corresponding T-cell lineages identified in divergent species. We also discuss TCRs with structural domains unlike that of conventional TCRs of mouse and human. By presenting this broad view of TCR sequence, structure, domain organization, and function, we seek to explore how this receptor has evolved across time and been selected for alternative antigen-recognition capabilities in divergent lineages.

Evolution of V genes from the TRV loci of mammals

Information concerning the evolution of T lymphocyte receptors (TCR) can be deciphered from that part of the molecule that recognizes antigen presented by major histocompatibility complex (MHC), namely the variable (V) regions. The genes that code for these variable regions are found within the TCR loci. Here, we describe a study of the evolutionary origin of V genes that code for the α and β chains of the TCR loci of mammals. In particular, we demonstrate that most of the 35 TRAV and 25 TRBV conserved genes found in Primates are also found in other Eutheria, while in Marsupials, Monotremes, and Reptiles, these genes diversified in a different manner. We also show that in mammals, all TRAV genes are derived from five ancestral genes, while all TRBV genes originate from four such genes. In Reptiles, the five TRAV and three out of the four TRBV ancestral genes exist, as well as other V genes not found in mammals. We also studied the TRGV and TRDV loci from all mammals, and we show a relationship of the TRDV to the TRAV locus throughout evolutionary time.

V genes in primates from whole genome sequencing data

The adaptive immune system uses V genes for antigen recognition. However, the evolutionary diversification and selection processes within and across species and orders remain poorly understood. Here, we studied the amino acid (AA) sequences obtained from the translated in-frame V exons of immunoglobulins (IG) and T cell receptors (TR) from 16 primate species whose genomes have been sequenced. Multi-species comparative analysis supports the hypothesis that V genes in the IG loci undergo birth/death processes, thereby permitting rapid adaptability over evolutionary time. We also show that multiple cladistic groupings exist in the TRA (35 clades) and TRB (25 clades) V gene loci and that each primate species typically contributes at least one V gene to each of these clades. The results demonstrate that IG V genes and TR V genes have quite different evolutionary pathways; multiple duplications can explain the IG loci results, while coevolutionary pressures can explain the phylogenetic results of the TR V gene loci. Our results suggest that there exist evolutionary relationships between V gene clades in the TRA and TRB loci. Due to the long-standing preservation of these clades, such genes may have specific and necessary roles for the viability of a species.