Evolutionary dynamics of the T-cell receptor VB gene family as inferred from the human and mouse genomic sequences

Evolution of T cells in the cancer-resistant naked mole-rat

Naked mole-rats (NMRs) are best known for their extreme longevity and cancer resistance, suggesting that their immune system might have evolved to facilitate these phenotypes. Natural killer (NK) and T cells have evolved to detect and destroy cells infected with pathogens and to provide an early response to malignancies. While it is known that NMRs lack NK cells, likely lost during evolution, little is known about their T-cell subsets in terms of the evolution of the genes that regulate their function, their clonotypic diversity, and the thymus where they mature. Here we find, using single-cell transcriptomics, that NMRs have a large circulating population of γδT cells, which in mice and humans mostly reside in peripheral tissues and induce anti-cancer cytotoxicity. Using single-cell-T-cell-receptor sequencing, we find that a cytotoxic γδT-cell subset of NMRs harbors a dominant clonotype, and that their conventional CD8 αβT cells exhibit modest clonotypic diversity. Consistently, perinatal NMR thymuses are considerably smaller than those of mice yet follow similar involution progression. Our findings suggest that NMRs have evolved under a relaxed intracellular pathogenic selective pressure that may have allowed cancer resistance and longevity to become stronger targets of selection to which the immune system has responded by utilizing γδT cells.

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 of Rare Thymic NKT Cell Precursors by Multiparameter Flow Cytometry

Mouse invariant natural killer T (NKT) cells are a subset of T lymphocytes which have been shown to play a significant role in innate and adaptive immune responses. Features of innate responses are attributed to these cells because they can be stimulated simultaneously with the same ligand to produce quickly and in large amount cytokines without prior immunization. Because these characteristics could be exploited for clinical applications, NKT cells have attracted considerable interest. Many studies have investigated the molecular mechanisms through which they are selected and differentiate. These studies are based on developmental models that serve as a scaffold to understand the specific roles played by various factors and to identify checkpoints during cellular development. Analysis of NKT cell precursors at the HSA^high stage, stage 0, can reveal potential selection defects, whereas analysis of NKT cells at the HSA^low stage can shed light on defects in the maturation/differentiation of the different NKT cell subsets (NKT1, 2, and 17). Unlike HSA^low NKT cell subsets, HSA^high NKT cell precursors are not accurately identified by flow cytometry because of their extreme rarity. Here, we describe an NKT cell enrichment strategy to identify unambiguously NKT cell precursors at the HSA^high stage that can be used to assess their distribution and characteristics by multicolor flow cytometry.

The Florida manatee (Trichechus manatus latirostris) T cell receptor loci exhibit V subgroup synteny and chain-specific evolution

The Florida manatee (Trichechus manatus latirostris) has limited diversity in the immunoglobulin heavy chain. We therefore investigated the antigen receptor loci of the other arm of the adaptive immune system: the T cell receptor. Manatees are the first species from Afrotheria, a basal eutherian superorder, to have an in-depth characterization of all T cell receptor loci. By annotating the genome and expressed transcripts, we found that each chain has distinct features that correlates to their individual functions. The genomic organization also plays a role in modulating sequence conservation between species. There were extensive V subgroup synteny blocks in the TRA and TRB loci between T. m. latirostris and human. Increased genomic locus complexity correlated to increased locus synteny. We also identified evidence for a VHD pseudogene for the first time in a eutherian mammal. These findings emphasize the value of including species within this basal eutherian radiation in comparative studies.

100 million years of multigene family evolution: origin and evolution of the avian MHC class IIB

BACKGROUND

Gene duplication has led to a most remarkable adaptation involved in vertebrates' host-pathogen arms-race, the major histocompatibility complex (MHC). However, MHC duplication history is as yet poorly understood in non-mammalian vertebrates, including birds.

RESULTS

Here, we provide evidence for the evolution of two ancient avian MHC class IIB (MHCIIB) lineages by a duplication event prior to the radiation of all extant birds >100 million years ago, and document the role of concerted evolution in eroding the footprints of the avian MHCIIB duplication history.

CONCLUSIONS

Our results suggest that eroded footprints of gene duplication histories may mimic birth-death evolution and that in the avian MHC the presence of the two lineages may have been masked by elevated rates of concerted evolution in several taxa. Through the presence of a range of intermediate evolutionary stages along the homogenizing process of concerted evolution, the avian MHCIIB provides a remarkable illustration of the erosion of multigene family duplication history.

High Throughput Sequencing of T Cell Antigen Receptors Reveals a Conserved TCR Repertoire

The T-cell receptor (TCR) repertoire is a mirror of the human immune system that reflects processes caused by infections, cancer, autoimmunity, and aging. Next-generation sequencing has become a powerful tool for deep TCR profiling. Herein, we used this technology to study the repertoire features of TCR beta chain in the blood of healthy individuals.Peripheral blood samples were collected from 10 healthy donors. T cells were isolated with anti-human CD3 magnetic beads according to the manufacturer's protocol. We then combined multiplex-PCR, Illumina sequencing, and IMGT/High V-QUEST to analyze the characteristics and polymorphisms of the TCR.Most of the individual T cell clones were present at very low frequencies, suggesting that they had not undergone clonal expansion. The usage frequencies of the TCR beta variable, beta joining, and beta diversity gene segments were similar among T cells from different individuals. Notably, the usage frequency of individual nucleotides and amino acids within complementarity-determining region (CDR3) intervals was remarkably consistent between individuals. Moreover, our data show that terminal deoxynucleotidyl transferase activity was biased toward the insertion of G (31.92%) and C (27.14%) over A (21.82%) and T (19.12%) nucleotides.Some conserved features could be observed in the composition of CDR3, which may inform future studies of human TCR gene recombination.

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.

Independent evolutionary origin of fem paralogous genes and complementary sex determination in hymenopteran insects

The primary signal of sex determination in the honeybee, the complementary sex determiner (csd) gene, evolved from a gene duplication event from an ancestral copy of the fem gene. Recently, other paralogs of the fem gene have been identified in several ant and bumblebee genomes. This discovery and the close phylogenetic relationship of the paralogous gene sequences led to the hypothesis of a single ancestry of the csd genetic system of complementary sex determination in the Hymenopteran insects, in which the fem and csd gene copies evolved as a unit in concert with the mutual transfers of sequences (concerted evolution). Here, we show that the paralogous gene copies evolved repeatedly through independent gene duplication events in the honeybee, bumblebee, and ant lineage. We detected no sequence tracts that would indicate a DNA transfer between the fem and the fem1/csd genes between different ant and bee species. Instead, we found tracts of duplication events in other genomic locations, suggesting that gene duplication was a frequent event in the evolution of these genes. These and other evidences suggest that the fem1/csd gene originated repeatedly through gene duplications in the bumblebee, honeybee, and ant lineages in the last 100 million years. Signatures of concerted evolution were not detectable, implicating that the gene tree based on neutral synonymous sites represents the phylogenetic relationships and origins of the fem and fem1/csd genes. Our results further imply that the fem1 and csd gene in bumblebees, honeybees, and ants are not orthologs, because they originated independently from the fem gene. Hence, the widely shared and conserved complementary sex determination mechanism in Hymenopteran insects is controlled by different genes and molecular processes. These findings highlight the limits of comparative genomics and emphasize the requirement to study gene functions in different species and major hymenopteran lineages.

The mammalian secreted RNases: Mechanisms of action in host defence

The mammalian ribonucleaseA family comprises a large group of structurally similar proteins which are secreted by a range of tissues and immune cells. Their physiological role is unclear. It has been suggested that some of these RNases contribute to host defence, notably eosinophil-derived neurotoxin, eosinophil cationic protein, eosinophil-associated RNases, RNase4, angiogenin (RNase5), RNase7, RNase8 and bovine seminal RNase. This review summarises data supporting the involvement of these proteins in host defence, focusing on their antimicrobial, cytotoxic and immunomodulatory activities. The extent to which the data support possible mechanisms of action for these proteins is discussed. This compilation of findings and current hypotheses on the physiological role of these RNases will provide a stimulus for further research and development of ideas on the contribution of the RNases to host defence.

Antigen-driven patterns of TCR bias are shared across diverse outcomes of human hepatitis C virus infection

Hepatitis C virus (HCV) infection causes significant morbidity and mortality worldwide. T cells play a central role in HCV clearance; however, there is currently little understanding of whether the disease outcome in HCV infection is influenced by the choice of TCR repertoire. TCR repertoires used against two immunodominant HCV determinants--the highly polymorphic, HLA-B*0801 restricted (1395)HSKKKCDEL(1403) (HSK) and the comparatively conserved, HLA-A*0101-restricted, (1435)ATDALMTGY(1443) (ATD)--were analyzed in clearly defined cohorts of HLA-matched, HCV-infected individuals with persistent infection and HCV clearance. In comparison with ATD, TCR repertoire selected against HSK was more narrowly focused, supporting reports of mutational escape in this epitope, in persistent HCV infection. Notwithstanding the Ag-driven divergence, T cell repertoire selection against either Ag was comparable in subjects with diverse disease outcomes. Biased T cell repertoires were observed early in infection and were evident not only in persistently infected individuals but also in subjects with HCV clearance, suggesting that these are not exclusively characteristic of viral persistence. Comprehensive clonal analysis of Ag-specific T cells revealed widespread use of public TCRs displaying a high degree of predictability in TRBV/TRBJ gene usage, CDR3 length, and amino acid composition. These public TCRs were observed against both ATD and HSK and were shared across diverse disease outcomes. Collectively, these observations indicate that repertoire diversity rather than particular Vβ segments are better associated with HCV persistence/clearance in humans. Notably, many of the anti-HCV TCRs switched TRBV and TRBJ genes around a conserved, N nucleotide-encoded CDR3 core, revealing TCR sequence mosaicism as a potential host mechanism to combat this highly variant virus.

Comparative genomics and evolution of the alpha-defensin multigene family in primates

Defensin genes encode small cationic antimicrobial peptides that form an important part of the innate immune system. They are divided into three families, alpha (α), beta (β), and theta (), according to arrangement of the disulfide bonding pattern between cysteine residues. Considering the functional importance of defensins, investigators have studied the evolution and the genomic organization of defensin genes. However, these studies have been restricted mainly to β-defensins. To understand the evolutionary dynamics of α-defensin genes among primates, we identified the α-defensin repertoires in human, chimpanzee, orangutan, macaque, and marmoset. The α-defensin genes in primates can be classified into three phylogenetic classes (class I, II, and III). The presence of all three classes in the marmoset indicates that their divergence occurred before the separation of New World and Old World monkeys. Comparative analysis of the α-defensin genomic clusters suggests that the makeup of the α-defensin gene repertoires between primates is quite different, as their genes have undergone dramatic birth-and-death evolution. Analysis of the encoded peptides of the α-defensin genes indicates that despite the overall high level of sequence divergence, certain amino acid residues or motifs are conserved within and between the three phylogenetic classes. The evolution of α-defensins in primates, therefore, appears to be governed by two opposing evolutionary forces. One force stabilizes specific amino acid residues and motifs to preserve the functional and structural integrity of the molecules and the other diversifies the sequences generating molecules with a wide range of activities against a large number of pathogens.

Genomic analysis reveals extensive gene duplication within the bovine TRB locus

Background

Diverse TR and IG repertoires are generated by V(D)J somatic recombination. Genomic studies have been pivotal in cataloguing the V, D, J and C genes present in the various TR/IG loci and describing how duplication events have expanded the number of these genes. Such studies have also provided insights into the evolution of these loci and the complex mechanisms that regulate TR/IG expression. In this study we analyze the sequence of the third bovine genome assembly to characterize the germline repertoire of bovine TRB genes and compare the organization, evolution and regulatory structure of the bovine TRB locus with that of humans and mice.

Results

The TRB locus in the third bovine genome assembly is distributed over 5 scaffolds, extending to ~730 Kb. The available sequence contains 134 TRBV genes, assigned to 24 subgroups, and 3 clusters of DJC genes, each comprising a single TRBD gene, 5–7 TRBJ genes and a single TRBC gene. Seventy-nine of the TRBV genes are predicted to be functional. Comparison with the human and murine TRB loci shows that the gene order, as well as the sequences of non-coding elements that regulate TRB expression, are highly conserved in the bovine. Dot-plot analyses demonstrate that expansion of the genomic TRBV repertoire has occurred via a complex and extensive series of duplications, predominantly involving DNA blocks containing multiple genes. These duplication events have resulted in massive expansion of several TRBV subgroups, most notably TRBV6, 9 and 21 which contain 40, 35 and 16 members respectively. Similarly, duplication has lead to the generation of a third DJC cluster. Analyses of cDNA data confirms the diversity of the TRBV genes and, in addition, identifies a substantial number of TRBV genes, predominantly from the larger subgroups, which are still absent from the genome assembly. The observed gene duplication within the bovine TRB locus has created a repertoire of phylogenetically diverse functional TRBV genes, which is substantially larger than that described for humans and mice.

Conclusion

The analyses completed in this study reveal that, although the gene content and organization of the bovine TRB locus are broadly similar to that of humans and mice, multiple duplication events have led to a marked expansion in the number of TRB genes. Similar expansions in other ruminant TR loci suggest strong evolutionary pressures in this lineage have selected for the development of enlarged sets of TR genes that can contribute to diverse TR repertoires.

Gene Duplication: A Drive for Phenotypic Diversity and Cause of Human Disease

Gene duplication is one of the key factors driving genetic innovation, i.e., producing novel genetic variants. Although the contribution of whole-genome and segmental duplications to phenotypic diversity across species is widely appreciated, the phenotypic spectrum and potential pathogenicity of small-scale duplications in individual genomes are less well explored. This review discusses the nature of small-scale duplications and the phenotypes produced by such duplications. Phenotypic variation and disease phenotypes induced by duplications are more diverse and widespread than previously anticipated, and duplications are a major class of disease-related genomic variation. Pathogenic duplications particularly involve dosage-sensitive genes with both similar and dissimilar over- and underexpression phenotypes, and genes encoding proteins with a propensity to aggregate. Phenotypes related to human-specific copy number variation in genes regulating environmental responses and immunity are increasingly recognized. Small genomic duplications containing defense-related genes also contribute to complex common phenotypes.

Adaptive evolution of multiple-variable exons and structural diversity of drug-metabolizing enzymes

BACKGROUND

The human genome contains a large number of gene clusters with multiple-variable-first exons, including the drug-metabolizing UDP glucuronosyltransferase (UGT1) and I-branching beta-1,6-N-acetylglucosaminyltransferase (GCNT2, also known as IGNT) clusters, organized in a tandem array, similar to that of the protocadherin (PCDH), immunoglobulin (IG), and T-cell receptor (TCR) clusters. To gain insight into the evolutionary processes that may have shaped their diversity, we performed comprehensive comparative analyses for vertebrate multiple-variable-first-exon clusters.

RESULTS

We found that there are species-specific variable-exon duplications and mutations in the vertebrate Ugt1, Gcnt2, and Ugt2a clusters and that their variable and constant genomic organizations are conserved and vertebrate-specific. In addition, analyzing the complete repertoires of closely-related Ugt2 clusters in humans, mice, and rats revealed extensive lineage-specific duplications. In contrast to the Pcdh gene clusters, gene conversion does not play a predominant role in the evolution of the vertebrate Ugt1, Gcnt2 and Ugt2 gene clusters. Thus, their tremendous diversity is achieved through "birth-and-death" evolution. Comparative analyses and homologous modeling demonstrated that vertebrate UGT proteins have similar three-dimensional structures each with N-terminal and C-terminal Rossmann-fold domains binding acceptor and donor substrates, respectively. Molecular docking experiments identified key residues in donor and acceptor recognition and provided insight into the catalytic mechanism of UGT glucuronidation, suggesting the human UGT1A1 residue histidine 39 (H39) as a general base and the residue aspartic acid 151 (D151) as an important electron-transfer helper. In addition, we identified four hypervariable regions in the N-terminal Rossmann domain that form an acceptor-binding pocket. Finally, analyzing patterns of nonsynonymous and synonymous nucleotide substitutions identified codon sites that are subject to positive Darwinian selection at the molecular level. These diversified residues likely play an important role in recognition of myriad xenobiotics and endobiotics.

CONCLUSION

Our results suggest that enormous diversity of vertebrate multiple variable first exons is achieved through birth-and-death evolution and that adaptive evolution of specific codon sites enhances vertebrate UGT diversity for defense against environmental agents. Our results also have interesting implications regarding the staggering molecular diversity required for chemical detoxification and drug clearance.

Concerted and birth-and-death evolution of multigene families

Until around 1990, most multigene families were thought to be subject to concerted evolution, in which all member genes of a family evolve as a unit in concert. However, phylogenetic analysis of MHC and other immune system genes showed a quite different evolutionary pattern, and a new model called birth-and-death evolution was proposed. In this model, new genes are created by gene duplication and some duplicate genes stay in the genome for a long time, whereas others are inactivated or deleted from the genome. Later investigations have shown that most non-rRNA genes including highly conserved histone or ubiquitin genes are subject to this type of evolution. However, the controversy over the two models is still continuing because the distinction between the two models becomes difficult when sequence differences are small. Unlike concerted evolution, the model of birth-and-death evolution can give some insights into the origins of new genetic systems or new phenotypic characters.

Selectionism and neutralism in molecular evolution

Charles Darwin proposed that evolution occurs primarily by natural selection, but this view has been controversial from the beginning. Two of the major opposing views have been mutationism and neutralism. Early molecular studies suggested that most amino acid substitutions in proteins are neutral or nearly neutral and the functional change of proteins occurs by a few key amino acid substitutions. This suggestion generated an intense controversy over selectionism and neutralism. This controversy is partially caused by Kimura's definition of neutrality, which was too strict (|2Ns|< or =1). If we define neutral mutations as the mutations that do not change the function of gene products appreciably, many controversies disappear because slightly deleterious and slightly advantageous mutations are engulfed by neutral mutations. The ratio of the rate of nonsynonymous nucleotide substitution to that of synonymous substitution is a useful quantity to study positive Darwinian selection operating at highly variable genetic loci, but it does not necessarily detect adaptively important codons. Previously, multigene families were thought to evolve following the model of concerted evolution, but new evidence indicates that most of them evolve by a birth-and-death process of duplicate genes. It is now clear that most phenotypic characters or genetic systems such as the adaptive immune system in vertebrates are controlled by the interaction of a number of multigene families, which are often evolutionarily related and are subject to birth-and-death evolution. Therefore, it is important to study the mechanisms of gene family interaction for understanding phenotypic evolution. Because gene duplication occurs more or less at random, phenotypic evolution contains some fortuitous elements, though the environmental factors also play an important role. The randomness of phenotypic evolution is qualitatively different from allele frequency changes by random genetic drift. However, there is some similarity between phenotypic and molecular evolution with respect to functional or environmental constraints and evolutionary rate. It appears that mutation (including gene duplication and other DNA changes) is the driving force of evolution at both the genic and the phenotypic levels.

CTL recognition of a bulged viral peptide involves biased TCR selection

MHC class I molecules generally present peptides of 8-10 aa long, forming an extended coil in the HLA cleft. Although longer peptides can also bind to class I molecules, they tend to bulge from the cleft and it is not known whether the TCR repertoire has sufficient plasticity to recognize these determinants during the antiviral CTL response. In this study, we show that unrelated individuals infected with EBV generate a significant CTL response directed toward an HLA-B*3501-restricted, 11-mer epitope from the BZLF1 Ag. The 11-mer determinant adopts a highly bulged conformation with seven of the peptide side chains being solvent-exposed and available for TCR interaction. Such a complex potentially creates a structural challenge for TCR corecognition of both HLA-B*3501 and the peptide Ag. Surprisingly, unrelated B*3501 donors recognizing the 11-mer use identical or closely related alphabeta TCR sequences that share particular CDR3 motifs. Within the small number of dominant CTL clonotypes observed, each has discrete fine specificity for the exposed side chain residues of the peptide. The data show that bulged viral peptides are indeed immunogenic but suggest that the highly constrained TCR repertoire reflects a limit to TCR diversity when responding to some unusual MHC peptide ligands.

The ribonuclease A superfamily of mammals and birds: identifying new members and tracing evolutionary histories

The RNase A superfamily has been important in biochemical, structural, and evolutionary studies and is believed to be the sole vertebrate-specific enzyme family. To understand the origin and diversification of the superfamily, we here determine its entire repertoire in the sequenced genomes of human, mouse, rat, and chicken. We report a previously unnoticed gene cluster in mouse chromosome 10 and a number of new genes, including mammalian RNases 11-13, which are close relatives of the recently identified RNases 9 and 10. Gene expression data imply male-reproductive functions for RNases 9-13, although their sequences suggest the lack of ribonucleolytic activities. In contrast to the presence of 13-20 functional genes in mammals, chicken has only 3 RNase genes, which are evolutionarily close to mammalian RNase 5, like other nonmammalian RNases. This and other evidence suggests that the RNase A superfamily originated from an RNase 5-like gene and expanded in mammals. Together with the fact that multiple lineages of the superfamily, including RNases 2, 3, 5, and 7, have antipathogenic activities, we suggest that the superfamily started off as a host-defense mechanism in vertebrates. Consistent with this hypothesis, all members of the superfamily exhibit high rates of amino acid substitution as is commonly observed in immunity genes.

Characterization of the T-cell receptor gamma locus and analysis of the variable gene segment expression in rabbit

The genomic organization and expression of genes of the T-cell receptor gamma (TRG) locus are described for mice and humans, but not for species such as rabbits (Oryctolagus cuniculus), in which gammadelta T cells compose a sizeable proportion of T cells in the periphery. We cloned 200 kb of the rabbit TRG locus and determined the TRGV gene usage in adult and newborn rabbits by RT-PCR. We identified two TRGJ genes, one TRGC gene, and 22 TRGV genes, all of which encoded functional variable regions. One TRGV gene is the unique member of the TRGV2 subgroup, whereas the other genes belong to the TRGV1 subgroup. Evolutionary analyses of TRGV1 genes identified three distinct groups that can be explained by separate duplication events in the rabbit genome. Evidence of gene conversion between TRGV1.1 and TRGV1.6 was observed. Both TRGV1 and TRGV2 subgroup genes were expressed in the spleen, intestine, and appendix of adult rabbits, and the repertoire of TRGV genes expressed in these tissues was similar. In these tissues from newborns, and in skin from adults, only the genes from the TRGV1 subgroup were expressed. Greater TRGV-J junctional diversity was found in tissues from adult compared to newborn rabbits. Our analyses indicate rabbits have a larger germ line encoded TRG repertoire compared with that of mice and humans. In addition, we found TRGV gene usage is alike in most tissues of rabbits similar to that found in humans but in contrast to that found in mice.

Comparative genomics and diversifying selection of the clustered vertebrate protocadherin genes

To explain the mechanism for specifying diverse neuronal connections in the brain, Sperry proposed that individual cells carry chemoaffinity tags on their surfaces. The enormous complexity of these connections requires a tremendous diversity of cell-surface proteins. A large number of neural transmembrane protocadherin (Pcdh) proteins is encoded by three closely linked human and mouse gene clusters (alpha, beta, and gamma). To gain insight into Pcdh evolution, I performed comprehensive comparative cDNA and genomic DNA analyses for the three clusters in the chimpanzee, rat, and zebrafish genomes. I found that there are species-specific duplications in vertebrate Pcdh genes and that additional diversity is generated through alternative splicing within the zebrafish "variable" and "constant" regions. Moreover, different codons (sites) in the mammalian Pcdh ectodomains (ECs) are under diversifying selection, with some under diversity-enhancing positive Darwinian selection and others, including calcium-binding sites, under strong purifying selection. Interestingly, almost all positively selected codon positions are located on the surface of ECs 2 and 3. These diversified residues likely play an important role in combinatorial interactions of Pcdh proteins, which could provide the staggering diversity required for neuronal connections in the brain. These results also suggest that adaptive selection is an additional evolutionary factor for increasing Pcdh diversity.

The T-cell receptor in primates: identifying and sequencing new owl monkey TRBV gene sub-groups

The New World primate Aotus nancymaae (owl monkey) has been shown to be an excellent experimental model when studying malarial parasites. Characterising the T-cell receptor (TR) alphabeta repertoire by means of the different variable beta (TRBV) genes displayed contributes to a better understanding of these lymphocytes' role in the response against several malarial antigens. This study describes identifying and characterising eleven new TRBV gene sub-groups in cDNA from Aotus nancymaae's peripheral blood lymphocytes; these 11 gene sequences displayed homology to the previously reported human TRBV3, TRBV10, TRBV11, TRBV14, TRBV18, TRBV19, TRBV20, TRBV25, TRBV27, TRBV29 and TRBV30 sub-groups, resulting in 83% overall homology at the amino acid level. An additional Aotus sequence was found having similarity with the human TRBJ-2-7*01 gene. Evolutionary relationships amongst these sequences and the homologous genes from both New and Old World primates have shown that the TRBV repertoire has been maintained in the species being studied, displaying varying association patterns and substitution rates, depending on the sub-group being studied. The degree of identity observed when comparing human and Aotus genes suggests that these species might have a similar TRBV repertoire.

Birth-and-Death Evolution with Strong Purifying Selection in the Histone H1 Multigene Family and the Origin of orphon H1 Genes

Histones are small basic nuclear proteins with critical structural and functional roles in eukaryotic genomes. The H1 multigene family constitutes a very interesting histone class gathering the greatest number of isoforms, with many different arrangements in the genome, including clustered and solitary genes, and showing replication-dependent (RD) or replication-independent (RI) expression patterns. The evolution of H1 histones has been classically explained by concerted evolution through a rapid process of interlocus recombination or gene conversion. Given such intriguing features, we have analyzed the long-term evolutionary pattern of the H1 multigene family through the evaluation of the relative importance of gene conversion, point mutation, and selection in generating and maintaining the different H1 subtypes. We have found the presence of an extensive silent nucleotide divergence, both within and between species, which is always significantly greater than the nonsilent variation, indicating that purifying selection is the major factor maintaining H1 protein homogeneity. The results obtained from phylogenetic analysis reveal that different H1 subtypes are no more closely related within than between species, as they cluster by type in the topologies, and that both RD and RI H1 variants follow the same evolutionary pattern. These findings suggest that H1 histones have not been subject to any significant effect of interlocus recombination or concerted evolution. However, the diversification of the H1 isoforms seems to be enhanced primarily by mutation and selection, where genes are subject to birth-and-death evolution with strong purifying selection at the protein level. This model is able to explain not only the generation and diversification of RD H1 isoforms but also the origin and long-term persistence of orphon RI H1 subtypes in the genome, something that is still unclear, assuming concerted evolution.

Characterizing T-cell receptor gamma-variable gene in Aotus nancymaae owl monkey peripheral blood

Gammadelta T lymphocytes have a heterodimeric complex formed by the association of gamma and delta chains as receptor. Proliferation of this lymphocyte population has been observed, when infection by several pathogens such as Mycobacterium tuberculosis and Plasmodium spp. occurs. The New World Monkey Aotus nancymaae has become a very good experimental model for the immunological and physiopathological study of these infectious agents. The A. nancymaae gamma-variable region was characterized from peripheral blood samples by using cDNA and genomic DNA polymerase chain reaction amplification, DNA sequencing, and dot-blot hybridization techniques. Seventeen different T-cell receptor gamma-variable (TCRGV) sequences were obtained. These sequences were distributed among TCRGV subsets 1, 2, or 3, according to human subset classification. Although no subset 4 amplification was obtained, this subset was detected by dot-blot hybridization. The presence of these 4 subsets resembles the behavior displayed by 'gammadelta-low species' (humans and mice), where high diversity among these lymphocytes can be observed. Homologies greater than 70% were found with respect to humans. Sequence convergence between human and A. nancymaae subsets 1 and 3 highlights Aotus as a promising model for studying these lymphocyte functions.

IMGT/PhyloGene: an on-line tool for comparative analysis of immunoglobulin and T cell receptor genes

IMGT/PhyloGene is an on-line software package for comparative analysis of immunoglobulin (IG) and T cell receptor (TR) variable genes of all vertebrate species, newly implemented in IMGT, the international ImMunoGeneTics information system ((R)). IMGT/PhyloGene is strongly associated with the IMGT gene and allele nomenclature and with the IMGT unique numbering for V-REGION, which directly creates standardized alignments from IMGT reference sequences. IMGT/PhyloGene is the first tool to use the IMGT expertized and standardized data for automated comparative analyses, and the first on-line software package for phylogenetic reconstruction to be integrated to a sequence database. Starting from a standardized alignment of selected sequences, IMGT/PhyloGene computes a matrix of evolutionary distances, builds a tree using the Neighbor-Joining (NJ) algorithm, and outputs various graphical tree representations. The resulting IMGT/PhyloGene tree is then used as a support for studying the evolution of particular subregions, such as the CDR-IMGT (Complementarity Determining Regions) or the V-RS (Variable gene Recombination Signals). IMGT/PhyloGene is freely available at http://imgt.cines.fr.

Analysis of the TCR beta variable gene repertoire in chimpanzees: identification of functional homologs to human pseudogenes

Chimpanzees are used for a variety of disease models such as hepatitis C virus (HCV) infection, where Ag-specific T cells are thought to be critical for resolution of infection. The variable segments of the TCR alphabeta genes are polymorphic and contain putative binding sites for MHC class I and II molecules. In this study, we performed a comprehensive analysis of genes that comprise the TCR beta variable gene (TCRBV) repertoire of the common chimpanzee Pan troglodytes. We identified 42 P. troglodytes TCRBV sequences representative of 25 known human TCRBV families. BV5, BV6, and BV7 are multigene TCRBV families in humans and homologs of most family members were found in the chimpanzee TCRBV repertoire. Some of the chimpanzee TCRBV sequences were identical with their human counterparts at the amino acid level. Notably four successfully rearranged TCRBV sequences in the chimpanzees corresponded to human pseudogenes. One of these TCR sequences was used by a cell line directed against a viral CTL epitope in an HCV-infected animal indicating the functionality of this V region in the context of immune defense against pathogens. These data indicate that some TCRBV genes maintained in the chimpanzee have been lost in humans within a brief evolutionary time frame despite remarkable conservation of the chimpanzee and human TCRBV repertoires. Our results predict that the diversity of TCR clonotypes responding to pathogens like HCV will be very similar in both species and will facilitate a molecular dissection of the immune response in chimpanzee models of human diseases.

Purifying selection and birth-and-death evolution in the histone H4 gene family

Histones are small basic proteins encoded by a multigene family and are responsible for the nucleosomal organization of chromatin in eukaryotes. Because of the high degree of protein sequence conservation, it is generally believed that histone genes are subject to concerted evolution. However, purifying selection can also generate a high degree of sequence homogeneity. In this study, we examined the long-term evolution of histone H4 genes to determine whether concerted evolution or purifying selection was the major factor for maintaining sequence homogeneity. We analyzed the proportion (p(S)) of synonymous nucleotide differences between the H4 genes from 59 species of fungi, plants, animals, and protists and found that p(S) is generally very high and often close to the saturation level (p(S) ranging from 0.3 to 0.6) even though protein sequences are virtually identical for all H4 genes. A small proportion of genes showed a low level of p(S) values, but this appeared to be caused by recent gene duplication. Our findings suggest that the members of this gene family evolve according to the birth-and-death model of evolution under strong purifying selection. Using histone-like genes in archaebacteria as outgroups, we also showed that H1, H2A, H2B, H3, and H4 histone genes in eukaryotes form separate clusters and that these classes of genes diverged nearly at the same time, before the eukaryotic kingdoms diverged.

Sequence variation and linkage disequilibrium in the human T-cell receptor beta (TCRB) locus

The T-cell receptor (TCR) plays a central role in the immune system, and > 90% of human T cells present a receptor that consists of the alpha TCR subunit (TCRA) and the beta subunit (TCRB). Here we report an analysis of 63 variable genes (BV), spanning 553 kb of TCRB that yielded 279 single-nucleotide polymorphisms (SNPs). Samples were drawn from 10 individuals and represent four populations-African American, Chinese, Mexican, and Northern European. We found nine variants that produce nonfunctional BV segments, removing those genes from the TCRB genomic repertoire. There was significant heterogeneity among population samples in SNP frequency (including the BV-inactivating sites), indicating the need for multiple-population samples for adequate variant discovery. In addition, we observed considerable linkage disequilibrium (LD) (r(2) > 0.1) over distances of approximately 30 kb in TCRB, and, in general, the distribution of r(2) as a function of physical distance was in close agreement with neutral coalescent simulations. LD in TCRB showed considerable spatial variation across the locus, being concentrated in "blocks" of LD; however, coalescent simulations of the locus illustrated that the heterogeneity of LD we observed in TCRB did not differ markedly from that expected from neutral processes. Finally, examination of the extended genotypes for each subject demonstrated homozygous stretches of >100 kb in the locus of several individuals. These results provide the basis for optimization of locuswide SNP typing in TCRB for studies of genotype-phenotype association.