Evolution and selection of primate T cell antigen receptor BV8 gene subfamily

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.

A novel locus on the X chromosome regulates post-maturity bone density changes in mice

Two mouse strains, AKR/J and SAMP6, were assessed longitudinally for bone mineral density of the spine. They displayed very different time courses of bone accrual, with the SAMP6 strain reaching a plateau for vertebral BMD at 3 months, whereas AKR/J mice continued to increase spine BMD for at least 8 months. Among 253 F(2) progeny of an AKR/JxSAMP6 cross, at 4 months of age, the BMD variance was 5-6% of the mean, vs. 15% for weight. Variance increased with age for every parameter measured, and was generally higher among males. The ratio of 6-month/4-month spine BMDs, termed DeltasBMD, had a normal distribution with 5.7% variance, and was largely independent of spine BMD (R=-0.23) or body weight (R=-0.12) at maturity. Heritability of the DeltasBMD trait was calculated at 0.59. Genetic mapping identified two significant loci, both distinct from those observed for BMD at maturity--implying that different genes regulate skeletal growth vs. remodeling. A locus on the X chromosome, replicated in two mouse F(2) populations (P<10(-4) for combined discovery and confirmation), affects age-dependent BMD change for both spine and the full skeleton. Its position agrees with a very narrow region identified by association mapping for effects on lumbar bone density in postmenopausal women [Parsons CA, Mroczkowski HJ, McGuigan FE, Albagha OM, Manolagas S, Reid DM, et al. Interspecies synteny mapping identifies a quantitative trait locus for bone mineral density on human chromosome Xp22. Hum Mol Genet 2005;14:3141-8]. A second locus, on chromosome 7, was observed in only one cross. Single-nucleotide polymorphisms (SNPs) are highly clustered near these loci, distinguishing the parental strains over only limited spans.

Sequence variation in the human T-cell receptor loci

Identifying common sequence variations known as single nucleotide polymorphisms (SNPs) in human populations is one of the current objectives of the human genome project. Nearly 3 million SNPs have been identified. Analysis of the relative allele frequency of these markers in human populations and the genetic associations between these markers, known as linkage disequilibrium, is now underway to generate a high-density genetic map. Because of the central role T cells play in immune reactivity, the T-cell receptor (TCR) loci have long been considered important candidates for common disease susceptibility within the immune system (e.g., asthma, atopy and autoimmunity). Over the past two decades, hundreds of SNPs in the TCR loci have been identified. Most studies have focused on defining SNPs in the variable gene segments which are involved in antigenic recognition. On average, the coding sequence of each TCR variable gene segment contains two SNPs, with many more found in the 5', 3' and intronic sequences of these segments. Therefore, a potentially large repertoire of functional variants exists in these loci. Association between SNPs (linkage disequilibrium) extends approximately 30 kb in the TCR loci, although a few larger regions of disequilibrium have been identified. Therefore, the SNPs found in one variable gene segment may or may not be associated with SNPs in other surrounding variable gene segments. This suggests that meaningful association studies in the TCR loci will require the analysis and typing of large marker sets to fully evaluate the role of TCR loci in common disease susceptibility in human populations.

Comparative genomics of the human and mouse T cell receptor loci

The availability of the complete genomic sequences of the human and mouse T cell receptor loci opens up new opportunities for understanding T cell receptors (TCRs) and their genes. The full complement of TCR gene segments is finally known and should prove a valuable resource for supporting functional studies. A rational nomenclature system has been implemented and is widely available through IMGT and other public databases. Systematic comparisons of the genomic sequences within each locus, between loci, and across species enable precise analyses of the various diversification mechanisms and some regulatory signals. The genomic landscape of the TCR loci provides fundamental insights into TCR evolution as highly localized and tightly regulated gene families.

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.

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

The diversity of T-cell receptors is generated primarily by the variable-region gene families, each of which is composed of a large number of member genes. The entire genomic sequence of the variable region (VB) of the T- cell receptor beta chain from humans and mice has become available. To understand the evolutionary dynamics of the VB gene family, we conducted a phylogenetic analysis of all VB genes from humans and mice, as well as a detailed analysis of internal DNA duplications in the human genomic VB region. The phylogenetic tree obtained shows that human and mouse VB genes intermingle extensively rather than forming two separate clusters and that many gene duplications occurred both before and after the divergence between primates and rodents. Analyzing the genomic maps of transposable elements (e.g., LINEs and SINEs) and relic VB genes in the VB gene region, we present evidence that a 20-kb VB region duplicated tandemly four times in the human lineage during the last 32 Myr, and 6 out of the 15 VB genes in this region have become nonfunctional during this period. Our results show that the VB gene family is subject to evolution by a birth-and-death process rather than to concerted evolution.